网络流24题之圆桌问题

二分图多重最大匹配：

在二分图最大匹配中，每个点（不管是X方点还是Y方点）最多只能和一条匹配边相关联，然而，我们经常遇到这种问题，即二分图匹配中一个点可以和多条匹配边相关联，但有上限，或者说，L_i表示点i最多可以和多少条匹配边相关联。

用网络流解决，建图方法：

在原图上建立源点S和汇点T，S向每个X方点连一条容量为该X方点L值的边，每个Y方点向T连一条容量为该Y方点L值的边，原来二分图中各边在新的网络中仍存在，容量为1（若该边可以使用多次则容量大于1），求该网络的最大流，就是该二分图多重最大匹配的值。

代码：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<string>
#include<climits>
#include<queue>
#include<stack>
#include<map>
#include<set>
#define N 1001
#define M N*N
#define inf 1<<30
using namespace std;
int n,m,tot,s,t;
int head[N],pos=-1,cur[N];
struct edge{int to,next,c;}e[M];
void add(int a,int b,int c)
{pos++;e[pos].to=b,e[pos].next=head[a],e[pos].c=c,head[a]=pos;}
queue<int>Q;bool vis[N];int dis[N];
bool bfs()
{
	for(int i=s;i<=t;i++)dis[i]=-1,vis[i]=0;
	dis[s]=0,vis[s]=1;Q.push(s);
	while(!Q.empty())
	{
		int u=Q.front();Q.pop();
		for(int i=head[u];i!=-1;i=e[i].next)
		{
			int v=e[i].to;
			if(vis[v]||e[i].c<=0)continue;
			dis[v]=dis[u]+1;
			vis[v]=1;Q.push(v);
		}
	}return vis[t];
}
int dfs(int u,int a)
{
	if(u==t||!a)return a;
	int f,flow=0;
	for(int &i=cur[u];i!=-1;i=e[i].next)
	{
		int v=e[i].to;
		if(dis[v]==dis[u]+1&&(f=dfs(v,min(a,e[i].c)))>0)
		{
			flow+=f,a-=f;
			e[i].c-=f,e[i^1].c+=f;
			if(!a)break;
		}
	}return flow;
}
int dinic()
{
	int ret=0;
	while(bfs())
	{
		for(int i=s;i<=t;i++)cur[i]=head[i];
		ret+=dfs(s,inf);
	}return ret;
}
void init(){memset(head,-1,sizeof(head));}
void print()
{
	for(int i=1;i<=n;i++)
	{
		for(int j=head[i];j!=-1;j=e[j].next)
		{
			int v=e[j].to;
			if(e[j].c>0)continue;
			if(v>n&&v<=n+m)
				printf("%d ",v-n);
		}printf("\n");
	}
}
int main()
{
//	freopen("table.in","r",stdin);
//	freopen("table.out","w",stdout);
	scanf("%d%d",&n,&m);
	s=0,t=n+m+1;init();
	for(int i=1;i<=n;i++)
	{
		int x;scanf("%d",&x);
		tot=tot+x;
		add(s,i,x);add(i,s,0);
	}
	for(int i=1;i<=m;i++)
	{
		int x;scanf("%d",&x);
		add(i+n,t,x);add(t,i+n,0);
	}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			add(i,j+n,1),add(j+n,i,0);
	if(dinic()==tot)
	{
		printf("1\n");
		print();
	}
	else printf("0\n");
}