关于构图:
首先,这是一道最大流记录路径的题目,构图不难,记录路径才是重点。
那么,我们就借鉴二分图的方法用网络流做(类似于试题库问题)。
构图:
1.有一个超级源点和汇点(st=0,ed=n+m+1)。
2.每一个单位连向源点,流量为该单位的人数。
3.每一个单位连向每一个餐桌,流量为1,表示这个单位的人可以坐任何一个餐桌,但为了使相同单位的人做不到一个餐桌上,所以流量为1。
4.每一个餐桌连向汇点,流量为该餐桌的人数,当全部都留满时,表示有符合要求的情况。
以上不难理解,很裸。
重点来了:
for(int i=1;i<=n;i++)
{
for(int j=last[i];j>=0;j=a[j].next)
{
y=a[j].y;
if(y!=st&&!a[j].z) printf("%d ",y-n); /*路径的记录,将每一个单位过一遍,找出它所流向的餐桌(当流量为0时,表示它有流过去),且由于反向边的缘故,要判断没有流向源点。注意:在构图时,我们将餐桌放在右边,加上了n,此时就得减去n。*/
}
printf("\n");
}
关于代码:
#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
#define INF 2147483647
using namespace std;
queue<int> f;
struct node{int x,y,z,next;} a[100000];
int last[10000];
int n,m,len=-1,st,ed;
void ins(int x,int y,int z)
{
a[++len].x=x;a[len].y=y;a[len].z=z;a[len].next=last[x];last[x]=len;
}
int h[100000];
bool bfs()
{
memset(h,0,sizeof(h));
h[st]=1;
f.push(st);
while(!f.empty())
{
int x=f.front();
for(int i=last[x];i>=0;i=a[i].next)
{
int y=a[i].y;
if(a[i].z>0&&h[y]==0)
{
h[y]=h[x]+1;
f.push(y);
}
}
f.pop();
}
if(h[ed]) return true; else return false;
}
int dfs(int x,int f)
{
int s=0,t;
if(x==ed) return f;
for(int i=last[x];i>=0;i=a[i].next)
{
int y=a[i].y;
if(a[i].z>0&&h[y]==h[x]+1&&f>s)
{
s+=(t=(dfs(y,min(f-s,a[i].z))));
a[i].z-=t;
a[i^1].z+=t;
}
}
if(!s) h[x]=0;
return s;
}
int dinic()
{
int sum=0;
while(bfs())
sum+=dfs(st,INF);
return sum;
}
int main()
{
int x,y,sum=0;
scanf("%d %d",&n,&m);
st=0,ed=n+m+1;
memset(last,-1,sizeof(last));
for(int i=1;i<=n;i++)
{
scanf("%d",&x);
sum+=x;
ins(st,i,x),ins(i,st,0);
for(int j=1;j<=m;j++)
ins(i,j+n,1),ins(j+n,i,0);
}
for(int i=1;i<=m;i++)
{
scanf("%d",&x);
ins(i+n,ed,x),ins(ed,i+n,0);
}
if(dinic()<sum)
{
printf("0");
return 0;
}
printf("1\n");
for(int i=1;i<=n;i++)
{
for(int j=last[i];j>=0;j=a[j].next)
{
y=a[j].y;
if(y!=st&&!a[j].z) printf("%d ",y-n);
}
printf("\n");
}
}