首页 > 代码库 > BZOJ2879: [Noi2012]美食节
BZOJ2879: [Noi2012]美食节
题解：
果然是神题，orz！！！
搬运题解：
做过修车就好办了。这个题仅仅是数据范围变大了一坨。建图就不说了，主要是动态加边。倒过来做，因为一个厨师最后一个菜做的时间是不会影响到其他菜的时间 的。
而且每一个厨师确定了最后一个菜才能去想倒数第二个菜是什么。所以每跑一次SPFA，就回来看看是哪个厨师做的菜，然后在多加一个点限制一下流量，将 这个点连向所有的菜。
实现很有技巧！！！我参考了lyd的代码
代码：
1 #include<cstdio> 2 3 #include<cstdlib> 4 5 #include<cmath> 6 7 #include<cstring> 8 9 #include<algorithm> 10 11 #include<iostream> 12 13 #include<vector> 14 15 #include<map> 16 17 #include<set> 18 19 #include<queue> 20 21 #include<string> 22 23 #define inf 1000000000 24 25 #define maxn 200000+5 26 27 #define maxm 2000000+5 28 29 #define eps 1e-10 30 31 #define ll long long 32 33 #define pa pair<int,int> 34 35 #define for0(i,n) for(int i=0;i<=(n);i++) 36 37 #define for1(i,n) for(int i=1;i<=(n);i++) 38 39 #define for2(i,x,y) for(int i=(x);i<=(y);i++) 40 41 #define for3(i,x,y) for(int i=(x);i>=(y);i--) 42 43 #define for4(i,x) for(int i=head[x],y;i;i=e[i].next) 44 45 #define for5(n,m) for(int i=1;i<=n;i++)for(int j=1;j<=m;j++) 46 47 #define mod 1000000007 48 49 using namespace std; 50 51 inline int read() 52 53 { 54 55 int x=0,f=1;char ch=getchar(); 56 57 while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();} 58 59 while(ch>=‘0‘&&ch<=‘9‘){x=10*x+ch-‘0‘;ch=getchar();} 60 61 return x*f; 62 63 } 64 int n,m,k,cnt[maxn],now[maxn],cntcnt,mincost,maxflow,tot=1,s,t,a[1000],b[1000][1000],sum,c[maxn],head[maxn],d[maxn],from[2*maxm]; 65 66 bool v[maxn]; 67 68 queue<int>q; 69 70 struct edge{int from,go,next,v,c;}e[2*maxm]; 71 72 void add(int x,int y,int v,int c) 73 74 { 75 76 e[++tot]=(edge){x,y,head[x],v,c};head[x]=tot; 77 78 e[++tot]=(edge){y,x,head[y],0,-c};head[y]=tot; 79 80 } 81 82 bool spfa() 83 84 { 85 86 for (int i=0;i<=cntcnt;i++){v[i]=0;d[i]=inf;} 87 88 q.push(s);d[s]=0;v[s]=1; 89 90 while(!q.empty()) 91 92 { 93 94 int x=q.front();q.pop();v[x]=0; 95 96 for (int i=head[x],y;i;i=e[i].next) 97 98 if(e[i].v&&d[x]+e[i].c<d[y=e[i].go]) 99 100 {101 102 d[y]=d[x]+e[i].c;from[y]=i;103 104 if(!v[y]){v[y]=1;q.push(y);}105 106 }107 108 }109 110 return d[t]!=inf;111 112 }113 114 void update()115 116 {117 118 for(int i=from[t];i!=s;i=from[e[i].from])119 {120 e[i].v--;e[i^1].v++;121 }122 mincost+=d[t];sum--;123 124 }125 126 int main()127 128 {129 130 freopen("input.txt","r",stdin);131 132 freopen("output.txt","w",stdout);133 134 n=read();m=read();135 for1(i,n)a[i]=read(),sum+=a[i];136 for1(i,n)for1(j,m)b[i][j]=read();137 s=0;t=n+1;cntcnt=n+m+1;138 for1(i,n)add(s,i,a[i],0);139 for1(i,m)add(n+1+i,t,1,0),cnt[i]=1,now[i]=tot;140 for1(i,n)for1(j,m)add(i,n+1+j,1,b[i][j]);141 while(sum)142 {143 if(spfa())update();144 for1(i,m)145 if(!e[now[i]^1].v)146 {147 add(++cntcnt,t,1,0),now[i]=tot,cnt[i]++;148 for1(j,n)add(j,cntcnt,1,cnt[i]*b[j][i]);149 break;150 }151 }152 cout<<mincost<<endl;153 154 return 0;155 156 }
2879: [Noi2012]美食节
Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 733 Solved: 411
[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
BZOJ2879: [Noi2012]美食节