【bzoj2879】【网络流】【NOI2012】【美食节】

2879: [Noi2012]美食节

Time Limit: 10 Sec   Memory Limit: 512 MB
Submit: 794   Solved: 442
[ Submit][ Status]

Description

CZ市为了欢迎全国各地的同学，特地举办了一场盛大的美食节。作为一个喜欢尝鲜的美食客，小M自然不愿意错过这场盛宴。他很快就尝遍了美食节所有的美食。然而，尝鲜的欲望是难以满足的。尽管所有的菜品都很可口，厨师做菜的速度也很快，小M仍然觉得自己桌上没有已经摆在别人餐桌上的美食是一件无法忍受的事情。于是小M开始研究起了做菜顺序的问题，即安排一个做菜的顺序使得同学们的等待时间最短。小M发现，美食节共有n种不同的菜品。每次点餐，每个同学可以选择其中的一个菜品。总共有m个厨师来制作这些菜品。当所有的同学点餐结束后，菜品的制作任务就会分配给每个厨师。然后每个厨师就会同时开始做菜。厨师们会按照要求的顺序进行制作，并且每次只能制作一人份。此外，小M还发现了另一件有意思的事情: 虽然这m个厨师都会制作全部的n种菜品，但对于同一菜品，不同厨师的制作时间未必相同。他将菜品用1, 2, ..., n依次编号，厨师用1, 2, ..., m依次编号，将第j个厨师制作第i种菜品的时间记为 ti,j 。小M认为：每个同学的等待时间为所有厨师开始做菜起，到自己那份菜品完成为止的时间总长度。换句话说，如果一个同学点的菜是某个厨师做的第k道菜，则他的等待时间就是这个厨师制作前k道菜的时间之和。而总等待时间为所有同学的等待时间之和。现在，小M找到了所有同学的点菜信息: 有 pi 个同学点了第i种菜品（i=1, 2, ..., n）。他想知道的是最小的总等待时间是多少。

Input

 输入文件的第1行包含两个正整数n和m，表示菜品的种数和厨师的数量。 第2行包含n个正整数，其中第i个数为pi，表示点第i种菜品的人数。 接下来有n行，每行包含m个非负整数，这n行中的第i行的第j个数为ti,j，表示第j个厨师制作第i种菜品所需的时间。 输入文件中每行相邻的两个数之间均由一个空格隔开，行末均没有多余空格。

Output

 输出仅一行包含一个整数，为总等待时间的最小值。

Sample Input

3 2
3 1 1
5 7
3 6
8 9

Sample Output

47

【样例说明】
厨师1先制作1份菜品2，再制作2份菜品1。点这3道菜的3个同学的等待时间分别为3，3+5=8，3+5+5=13。
厨师2先制作1份菜品1，再制作1份菜品3。点这2道菜的2个同学的等待时间分别为7，7+9=16。
总等待时间为3+8+13+7+16=47。
虽然菜品1和菜品3由厨师1制作更快，如果这些菜品都由厨师1制作，总等待时间反而更长。如果按上述的做法，将1份菜品1和1份菜品3调整到厨师2制作，这样厨师2不会闲着，总等待时间更短。
可以证明，没有更优的点餐方案。


【数据规模及约定】
对于100%的数据，n <= 40, m <= 100, p <= 800, ti,j <= 1000（其中p = ∑pi，即点菜同学的总人数）。
每组数据的n、m和p值如下：
测试点编号 n m p
1 n = 5 m = 5 p = 10
2 n = 40 m = 1 p = 400
3 n = 40 m = 2 p = 300
4 n = 40 m = 40 p = 40
5 n = 5 m = 40 p = 100
6 n = 10 m = 50 p = 200
7 n = 20 m = 60 p = 400
8 n = 40 m = 80 p = 600
9 n = 40 m = 100 p = 800
10 n = 40 m = 100 p = 800

思路：这道题的思路跟修车很像，建立一个源点，汇点，和三排点，从源点向第一排点连容量为这种产品的需求量，费用为0的边；从第一排向第二排连容量为这种产品需求量，费用为0的边；从第二排向第三排连容量为1，费用为k*cost[i][j]的边；从第三排边向汇点连容量为1，费用为0的边。

这个题如果就单单这样去做，会T，只能得60分，因为数据范围太大了，所以要加一些优化：第三排的边我们要动态去加，每当我们求出一条增光路时，我们看看这条路上是经过了哪个厨师，这个厨师做了几个菜，假如是第j个厨师做了第k个菜，然后我们就再连上第j个厨师做k+1道菜的边，这样动态加边才会过这个题。（就是在动态加边的过程中，推公式推了好久，没恶心死我！！）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int point[100000],next[10000000],pre[100000],l[10000000],h,t,tot=1,dis[100000],n,m,cost[41][101];
bool f[100000];
char * ptr=(char *)malloc(2000000);
inline void in(int &x){
	while(*ptr<'0'||*ptr>'9')++ptr;
	x=0;
	while(*ptr>47&&*ptr<58)x=x*10+*ptr++-'0';
}
struct S{
	int st,en,va,co;
}aa[10000000];
inline void add(int i,int j,int x,int y){
	tot+=1;next[tot]=point[i];point[i]=tot;
	aa[tot].st=i;aa[tot].en=j;aa[tot].va=x;aa[tot].co=y;
	tot+=1;next[tot]=point[j];point[j]=tot;
	aa[tot].st=j;aa[tot].en=i;aa[tot].va=0;aa[tot].co=-y;
}
inline int SPFA(int x,int y){
	memset(f,1,sizeof(f));
	memset(dis,127,sizeof(dis));
	memset(pre,0,sizeof(pre));
	int i,j,u;
	dis[x]=0;h=t=1;l[h]=x,pre[x]=x;
	while(h<=t){
		u=l[h];f[u]=true;
		for(i=point[u];i;i=next[i]){
			if(aa[i].va>0&&dis[aa[i].en]>dis[u]+aa[i].co){
				dis[aa[i].en]=dis[u]+aa[i].co;
				pre[aa[i].en]=i;
				if(f[aa[i].en]){
					t=(t+1)%10000000;
					f[aa[i].en]=false;
					l[t]=aa[i].en;
				}
			}
		}
		h=(h+1)%10000000;
	}
	return dis[y]>2100000000?0:dis[y];
} 
inline int ISAP(int x,int y){
	int i,j,minn=y,num=1,ii;
	for(i=y;i!=x;i=aa[pre[i]].st){
		if(num==2) ii=i;
		num+=1;
		minn=min(minn,aa[pre[i]].va);
	}
	for(i=y;i!=x;i=aa[pre[i]].st){
		aa[pre[i]].va-=minn;
		aa[pre[i]^1].va+=minn;
	}
	return ii;
}
int main()
{
	freopen("noi12_delicacy.in","r",stdin);
	freopen("noi12_delicacy.out","w",stdout);
	fread(ptr,1,2000000,stdin);
	int i,j,num=0,x,k,minn=1,ans=0,del;
	in(n);in(m);
	for(i=1;i<=n;++i){
		in(x);
		num+=x;
		add(1,i+1,x,0);
		for(j=1;j<=m;++j){
			add(i+1,1+j*n+i,x,0);//add(i+1,1+n+(j-1)*n+i,x,0)
		}
	}	
	for(i=1;i<=n;++i){
		for(j=1;j<=m;++j){
			in(cost[i][j]);
			add(1+n*j+i,1+n*(m+1)+(j-1)*num+1,1,cost[i][j]);
			if(i==1) add(1+n*(m+1)+(j-1)*num+1,2+n*(m+1)+num*m,1,0);
		}
	}	
/*	int u;
	for(i=1;i<=2+n*(m+1)+num*m;++i){
		u=point[i]; cout<<i<<endl;
		while(u){
			cout<<aa[u].st<<' '<<aa[u].en<<' '<<aa[u].va<<' '<<aa[u].co<<endl;
			u=next[u];
		}
	}*/
	minn=1;
	while(minn){
		minn=SPFA(1,2+n*(m+1)+num*m);
		if(minn) ans+=minn;
		else break;
		del=ISAP(1,2+n*(m+1)+num*m);
		j=((del-1-n*(m+1))/num)+1;
		k=((del-1-n*(m+1))%num)+1;
		//if(!(del-1-n*(m+1))%num){
		    for(i=1;i<=n;++i){
			    add(1+n+(j-1)*n+i,1+n+n*m+(j-1)*num+k,1,k*cost[i][j]);
			    if(i==1) add(1+n+n*m+(j-1)*num+k,2+n*(m+1)+num*m,1,0);
		    }
	//	}
	//	cout<<ans<<endl;
	}
	printf("%d",ans);
}