BZOJ 2879 NOI2012美食节

链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2879

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）。他想知道的是最小的总等待时间是多少。

思路:动态建点，每次如果能有增广路，就新建点。

 1 #include<algorithm>
 2 #include<cstdio>
 3 #include<cmath>
 4 #include<cstring>
 5 #include<iostream>
 6 #define N 200005
 7 int tot,go[3000005],first[200005],next[3000005],flow[3000005],cost[3000005];
 8 int op[3000005],from[N],edge[N],dis[N],vis[N],c[2*N];
 9 int a[505][505],n,m,p[N],ans,S,T,all;
10 int read(){
11     int t=0,f=1;char ch=getchar();
12     while (ch<‘0‘||ch>‘9‘){if (ch==‘-‘) f=-1;ch=getchar();}
13     while (‘0‘<=ch&&ch<=‘9‘){t=t*10+ch-‘0‘;ch=getchar();}
14     return t*f;
15 }
16 void insert(int x,int y,int z,int l){
17     tot++;go[tot]=y;next[tot]=first[x];first[x]=tot;flow[tot]=z;cost[tot]=l;
18 }
19 void add(int x,int y,int z,int l){
20     insert(x,y,z,l);op[tot]=tot+1;insert(y,x,0,-l);op[tot]=tot-1;
21 }
22 bool spfa(){
23     for (int i=0;i<=T;i++) dis[i]=0x3f3f3f3f,vis[i]=0;
24     int h=1,t=1;dis[S]=0;c[1]=S;
25     while (h<=t){
26         int now=c[h++];
27         for (int i=first[now];i;i=next[i]){
28             int pur=go[i];
29             if (flow[i]&&dis[pur]>dis[now]+cost[i]){
30                 dis[pur]=dis[now]+cost[i];
31                 from[pur]=now;
32                 edge[pur]=i;
33                 if (vis[pur]) continue;
34                 vis[pur]=1;
35                 c[++t]=pur;
36             }
37         }
38         vis[now]=0;
39     }
40     return dis[T]!=0x3f3f3f3f;
41 }
42 void updata(){
43     int mn=0x7fffffff,x,y;
44     for (int i=T;i!=S;i=from[i]){
45         mn=std::min(mn,flow[edge[i]]);
46         if (from[i]==S){
47             x=(i-1)/all+1;
48             y=(i)%all+1;
49         }
50     }
51     for (int i=T;i!=S;i=from[i]){
52         ans+=mn*cost[edge[i]];
53         flow[edge[i]]-=mn;
54         flow[op[edge[i]]]+=mn;
55     }
56     for (int i=1;i<=m;i++)
57      add((x-1)*all+y,n*all+i,1,y*a[i][x]);
58 }
59 int main(){
60     m=read();n=read();
61     S=0;
62     for (int i=1;i<=m;i++){
63         p[i]=read();
64         all+=p[i];
65     }
66     T=m+n*all+1;
67     for (int i=1;i<=n*all;i++)
68      add(S,i,1,0);
69     for (int i=1;i<=m;i++)
70      add(i+n*all,T,p[i],0);
71     for (int i=1;i<=m;i++)
72      for (int j=1;j<=n;j++)
73       a[i][j]=read();
74     for (int i=1;i<=n;i++)
75      for (int j=1;j<=m;j++)
76       add((i-1)*all+1,n*all+j,1,a[j][i]);
77     while (spfa()) updata();
78     printf("%d\n",ans);
79 }