poj2516（最小费用最大流模版）

首先就这题而言，各k互不影响，可以分开来算，累加即可

对于构图来说，以后可以源点s=0，汇点t=n+m，其中一组点是1~n，另一组n+1~n+m，注意可以这样构图

对于最小费用最大流，要用spfa来做最短路径，只需记录路径即可

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#define N 130
using namespace std;
const int inf=0x3f3f3f3f;
int n,m,k,offer[N][N],need[N][N],map[N][N],cost[N][N],sumn[N],sumo[N],pre[N],s,t;
bool spfa()
{
	bool vis[N]={false};
	int dis[N];
	memset(dis,inf,sizeof(dis));//对了，不能int dis【1000】=｛inf｝这样付初值，这样绝对不可以！！！！！！！！！！！！！
	queue<int> q;
	vis[s]=true;
	dis[s]=0;
	q.push(s);
	while (!q.empty())
	{
		int u=q.front();
		q.pop();
		for (int i=0;i<=t;i++)
		if (map[u][i]>0&&dis[i]>dis[u]+cost[u][i])
		{
			dis[i]=dis[u]+cost[u][i];
			pre[i]=u;//注意记录路径！！
			if (!vis[i])
			{
				q.push(i);
				vis[i]=true;
			}
		}
		vis[u]=false;
	}/spfa，
	if (dis[t]==inf) return false;else return true;//判断是否有增广路，然后返回即可。
}
int work()
{
	int mi=inf,ans=0;
	while (spfa())
	{
		for (int i=t;i;i=pre[i])
		mi=min(mi,map[pre[i]][i]);//这种方式找最小
		for (int i=t;i;i=pre[i])
		ans+=mi*cost[pre[i]][i],map[pre[i]][i]-=mi,map[i][pre[i]]+=mi;//费用流的关键，注意反向边，别忘了减去流量
	}
	return ans;//返回最小费用！
}
int main()
{
	while (scanf("%d%d%d",&n,&m,&k))
	{
		memset(sumn,0,sizeof(sumn));
		memset(sumo,0,sizeof(sumo));
		memset(need,0,sizeof(need));
		memset(offer,0,sizeof(offer));
		if (n==0) break;
		int sign=0;
		for (int i=1;i<=n;i++)
			for (int j=1;j<=k;j++) scanf("%d",&need[i][j]),sumn[j]+=need[i][j];
		for (int i=1;i<=m;i++)
			for (int j=1;j<=k;j++) scanf("%d",&offer[i][j]),sumo[j]+=offer[i][j];
		for (int j=1;j<=k;j++)
		if (sumo[j]<sumn[j])
		{
			sign=1;
			break;
		}	
		s=0;
		t=n+m+1;
		int ans=0;
		for (int kk=1;kk<=k;kk++)
		{
			memset(map,0,sizeof(map));
			memset(cost,0,sizeof(cost));
			
			for (int i=m+1;i<=n+m;i++)
			for (int j=1;j<=m;j++) 
			{
				scanf("%d",&cost[j][i]);
				cost[i][j]=-cost[j][i];
			}
			if (sign) continue;
			for (int i=1;i<=m;i++) map[s][i]=offer[i][kk];
			for (int i=m+1;i<=n+m;i++) map[i][t]=need[i-m][kk];
			for (int i=1;i<=m;i++)
			for (int j=m+1;j<=n+m;j++) map[i][j]=offer[i][kk];
			ans+=work();
		}	
		
		if (sign)printf("-1\n");else printf("%d\n",ans);
	}
	return 0;
}