bzoj 2879: [Noi2012]美食节 费用流+动态加边 详解

 

题意：有n种菜和m个厨师，第i种菜有sum[i]个，每种菜给每个厨师做有一个花费时间a[i,j]，每个厨师一次只能做一道菜。每道菜的等待时间为这道菜完成的时间。求最小的等待时间之和。p= sum[i]  之和.

PS：如果还没有做过：[SCOI2007]修车 可以先去做这个题目会比较好 附加题解报告「BZOJ1070」[SCOI2007] 修车 详解

 

试了一下直接修改 修车 的代码，建成2+n+p*m个点，n*m*p条边的有向图，对其求费用流，时间复杂度O((2+n+p*m)^2*(n*m*p))，实际测试中得到40分。

　　既然不能做这么庞大的图的网络流，我们考虑动态加点。

　　刚开始连边：

　现在，所有厨师都闲着，那我们把每位厨师做倒数第一道菜的边连向N个点。

　　

　　然后做费用流。

　　在增广的过程中，如果k厨师做倒数第i个菜连到汇点满流，则每个菜对k厨师做倒数第i+1个菜连一条边，流量为1，权值为(i+1) * t[][k]。这样做到动态加边，避免因为边数过多导致超时。

谢谢Africamonkey大神在其博客中详细地写到如何动态加点

 

代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
#define inf 0x7fffffff
#define N 100005
#define M 7000000
using namespace std;

int n,m,p,ans,cnt,last[N],vis[N],pre[N],dis[N],a[50][120],sum[50],s,t;
struct edge{int to,from,op,c,w,next;}e[M];
queue <int> q;

void add_(int u,int v,int c,int w)
{
	e[++cnt].to=v;e[cnt].from=u;e[cnt].c=c;e[cnt].w=w;e[cnt].op=cnt+1;e[cnt].next=last[u];last[u]=cnt;
	e[++cnt].to=u;e[cnt].from=v;e[cnt].c=0;e[cnt].w=-w;e[cnt].op=cnt-1;e[cnt].next=last[v];last[v]=cnt;
}

bool spfa()
{
	for (int i=s;i<=t;i++)
	{
		dis[i]=inf;
		vis[i]=pre[i]=0;//初始化
	}
	vis[s]=1;
	dis[s]=0;
	while (!q.empty()) q.pop();
	q.push(s);
	while (!q.empty())
	{
		int u=q.front();
		q.pop();vis[u]=0;
		for (int i=last[u];i;i=e[i].next)
			if (e[i].c&&dis[u]+e[i].w<dis[e[i].to])
			{
			    int y=e[i].to;
				dis[y]=e[i].w+dis[u];
				pre[y]=i;//几号边
				if (!vis[y])
				{
					q.push(y);
					vis[y]=1;
				}
			}
	}
	return dis[t]!=inf;
}

void mcf()
{
	int mnf=inf,i=t;
	while (pre[i])
	{
		mnf=min(mnf,e[pre[i]].c);//最小流
		i=e[pre[i]].from;//回溯
	}
	ans+=mnf*dis[t];//花费
	i=t;
	while (pre[i])
	{
	    int it=pre[i];
		e[pre[i]].c-=mnf;
		e[e[pre[i]].op].c+=mnf;//其反边的流量
		if (e[pre[i]].to==t&&(e[pre[i]].from-n)%p>0)//该点属于厨师  -n是减去左边源点与n个点连的边
		{
			add_(e[pre[i]].from+1,t,1,0);//  再在下一个点上建立新边 通道流量为1 花费为0；
			int kind=(e[pre[i]].from-n)/p+1,num=(e[pre[i]].from-n)%p+1;//属于那一种类的菜  第几道菜
			for (int j=1;j<=n;j++)//加边 该厨师 连接每一道菜  流量为1 花费为 a[j][kind]*num；
				add_(j,e[pre[i]].from+1,1,a[j][kind]*num);
		}
		i=e[pre[i]].from;
	}
}

int main()
{
	scanf("%d%d",&n,&m);
	for (int i=1;i<=n;i++)
	{
		scanf("%d",&sum[i]);//每一种有几个人
		p+=sum[i];
	}
	for (int i=1;i<=n;i++)
    for (int j=1;j<=m;j++)
        scanf("%d",&a[i][j]);
	s=0;t=n+m*p+1;
	for (int i=1;i<=n;i++)
     add_(s,i,sum[i],0);//源点连向 n 种菜，费用 0，流量为这种菜的 数量
	for (int i=1;i<=m;i++)
	{
		for (int j=1;j<=n;j++)//只将每种菜连向每个拆点的第一个点
        add_(j,n+(i-1)*p+1,1,a[j][i]);
		add_(n+(i-1)*p+1,t,1,0);//将这 m*p 个点连向汇点，费用 0，流量 1
	}
	while (spfa()) mcf();
	printf("%d",ans);
	return 0;
}