BZOJ2879 [Noi2012]美食节
本文版权归ljh2000和博客园共有，欢迎转载，但须保留此声明，并给出原文链接，谢谢合作。
本文作者：ljh2000
作者博客：http://www.cnblogs.com/ljh2000-jump/
转载请注明出处，侵权必究，保留最终解释权！
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
//It is made by ljh2000
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <ctime>
#include <vector>
#include <queue>
#include <map>
#include <set>
#include <string>
#include <complex>
using namespace std;
typedef long long LL;
const int MAXN = 45;
const int MAXM = 120;
const int inf = (1<<30);
const int maxn = 100011;
const int maxm = 3200011;
int n,m,S,T,cook_cnt,tot,t[45][150],head,tail;//n表示菜品,m表示厨师数量
int num[MAXN],first[maxn],ecnt,ans,dis[maxn],pre[maxn],ppre[maxn],dui[maxn*30];
bool in[maxn];
//感觉很像我以前做的一道省选题，叫修车来着...正着不好做，考虑算每次的贡献，倒数第k次对答案的总贡献就是k*tim
//但是这题的数据范围很大，我们需要优化。考虑因为假设一个厨师做的倒数第k道菜尚未决策，那么没有必要考虑倒数第k+1次。所以我们动态加边，当倒数第k道菜做完才加入倒数第k+1道菜的位置，每次修改增广路时很容易得到。
//这样的话，复杂度是可以保证的。因为每次图的点数都是n+m，边数是n*m...
struct edge{
int to,f,w,next;
}e[maxm];
inline void link(int x,int y,int f,int z){
e[++ecnt].next=first[x]; first[x]=ecnt; e[ecnt].to=y; e[ecnt].f=f; e[ecnt].w=z;
e[++ecnt].next=first[y]; first[y]=ecnt; e[ecnt].to=x; e[ecnt].f=0; e[ecnt].w=-z;
}
inline int getint(){
int w=0,q=0; char c=getchar(); while((c<'0'||c>'9') && c!='-') c=getchar();
if(c=='-') q=1,c=getchar(); while (c>='0'&&c<='9') w=w*10+c-'0',c=getchar(); return q?-w:w;
}
inline bool SPFA(){
for(int i=0;i<=T;i++) dis[i]=inf,in[i]=false,pre[i]=ppre[i]=-1;
head=tail=0; dui[++tail]=0; in[0]=true; dis[0]=0; int u;
while(head<tail) {
u=dui[++head]; in[u]=false;
for(int i=first[u];i;i=e[i].next) {
if(e[i].f==0) continue; int v=e[i].to;
if(dis[v]>dis[u]+e[i].w) {
dis[v]=dis[u]+e[i].w;
pre[v]=u; ppre[v]=i;
if(!in[v]) { in[v]=true; dui[++tail]=v; }
}
}
}
if(dis[T]==inf) return false;
return true;
}
inline void solve(){
int u=T,ff=inf;
for(;u!=S;u=pre[u]) ff=min(ff,e[ppre[u]].f);
for(u=T;;u=pre[u]) { ans+=ff*e[ppre[u]].w,e[ppre[u]].f-=ff,e[ppre[u]^1].f+=ff; if(pre[u]==S) break; }
int cook,cc; cook=(u-1)/tot+1; cc=(u-1)%tot+1; if(cc==tot) return ;
cc++; u++;
link(S,u,1,0);
for(int i=1;i<=n;i++) link(u,cook_cnt+i,1,t[i][cook]*cc);
}
inline void work(){
n=getint(); m=getint(); for(int i=1;i<=n;i++) num[i]=getint(),tot+=num[i];//tot表示需要做的菜的总数
for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) t[i][j]=getint();
cook_cnt=m*tot;//厨师拆成点的总数
S=0; T=cook_cnt+n+1; ecnt=1;
for(int i=1;i<=n;i++) link(cook_cnt+i,T,num[i],0);
for(int i=1;i<=m;i++) link(S,(i-1)*tot+1,1,0);
for(int i=1;i<=m;i++)//从每个厨师的第一个点出发连边
for(int j=1;j<=n;j++)
link((i-1)*tot+1,cook_cnt+j,1,t[j][i]);
while(SPFA())
solve();
printf("%d",ans);
}
int main()
{
work();
return 0;
}
