1、什么是最小费用最大流问题
上篇文章我们讲解了最大流问题，那什么是最小费用最大流呢？听名字就可以看出，我们要在满足最大流的同时找到达成最大流的最小费用。
对于一个网络流，最大流是一定的，但是组成最大流的费用是可以不同的，这里就有了在最大流网络上产生的费用流网络，就有了最小花费问题。
简单来说，就是满足最大流的路径可能有多条，我们要从这多条路径中找到一条花费代价最小的路径。所以最大流是解决这类问题的前提
2、EK算法 + SPFA 最短路
我们用每条边单位流量的花费作为边权，假如一条合法路径上每条边的花费分别为 c1,c2,.......ck ， 并且这条边上的最小流量为flow，
那么这条路径上的花费为 ： c1 * flow + c2*flow + ..... + ck*flow = （c1+ c2 + c3 + .... + ck）*flow = dis [ci] * flow
这里的 dis[ci] 就是我们要求的最短路！
3、算法思想
采用贪心的思想，每次找到一条从源点到达汇点的路径，增加流量，且该条路径满足使得增加的流量的花费最小，直到无法找到一条从源点到达汇点的路径，算法结束。
由于最大流量有限，每执行一次循环流量都会增加，因此该算法肯定会结束，且同时流量也必定会达到网络的最大流量；同时由于每次都是增加的最小的花 费，即当前的最小花费是所有到达当前流量flow时的花费最小值，因此最后的总花费最小。
4、求解步骤
1、找到一条从源点到达汇点的花费最小的路径，该花费 = 使用该路径上的边的单位费用之和。2、然后找出这条路径上的边的容量的最小值f，则当前最大流 max_flow 扩充f（求最大流的过程），同时当前最小费用 min_cost 扩充 f*min_dist(s,t)。
3、更新路径流量 flow[][] ，将这条路径上的每条正向边的容量都减少f，每条反向边的容量都增加f。
4、重复（1~3）直到无法找到从源点到达汇点的路径。
5、需要注意几点
1、注意超级源点和超级终点的建立。这里的超级源点/汇点的建立是把所有边统一起来，构成一张网络图，因为有些题目不会直接给你网络图，
需要你自己去建立，例如hdu 1533
2、初始化时，正向边的单位流量费用为cost[u][v]，那么反向边的单位流量费用就为 -cost[u][v]。正向边和反向边的费用互为相反数，因为回流费用减少。
3、费用cost数组和容量cap数组每次都要初始化为0。
4、求解从源点到汇点的“最短”路径时，由于网络中存在负权边，因此使用SPFA来实现。
6、代码
int cap[500][500],cost[500][500],flow[500][500];//cap是容量，cost是花费，flow是流量
int pre[500],dis[500],vis[500],cnt[500];//pre是记录增广路的前驱节点，dis[i]是记录源点到节点i的最小距离
//vis[i]标记点是否在队列中，cnt[i]记录点i进入队列的次数
int n,m;
int st,endd;
int mn_cost,mx_flow;
int spfa()
{
for(int i=0;i<n;i++)
dis[i]=mx;
memset(vis,0,sizeof(vis));
memset(cnt,0,sizeof(cnt));
queue<int>p;
p.push(st);//st是起点
vis[st]=1;
cnt[st]=1;
dis[st]=0;
while(!p.empty())
{
int now=p.front();
p.pop();
vis[now]=0;
for(int i=0;i<n;i++)
{
if(cap[now][i]>flow[now][i]&&dis[i]>dis[now]+cost[now][i])//这里将费用看成是长度，求源点到汇点的最小距离
{
dis[i]=dis[now]+cost[now][i];
pre[i]=now;
if(vis[i]==0)
{
p.push(i);
vis[i]=1;
cnt[i]++;
if(cnt[i]>n)
return 0;
}
}
}
}
if(dis[endd]>=mx)
return 0;
return 1;
}
void bfs(int n)//顶点数
{
memset(flow,0,sizeof(flow));
mx_flow=0;
mn_cost=0;
while(spfa())
{
int f=mx;
for(int i=endd;i!=st;i=pre[i])
f=min(f,cap[pre[i]][i]-flow[pre[i]][i]);
for(int i=endd;i!=st;i=pre[i])
{
flow[pre[i]][i]+=f;
flow[i][pre[i]]-=f;
}
mn_cost+=dis[endd]*f; mx_flow+=f;
}
}