BZOJ 2879: [Noi2012]美食节

2879: [Noi2012]美食节

Time Limit: 10 Sec Memory Limit: 512 MB
Submit: 1834 Solved: 969
[Submit][Status][Discuss]

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

HINT

Source

[Submit][Status][Discuss]

建图基本同BZOJ修车一题，但是需要动态加点。每次找出新满流的边，加入这个厨师的下一个点。

  1 #include <cstdio>
  2 
  3 inline int nextChar(void) {
  4     static const int siz = 1024;
  5     static char buf[siz];
  6     static char *hd = buf + siz;
  7     static char *tl = buf + siz;
  8     if (hd == tl)
  9         fread(hd = buf, 1, siz, stdin);
 10     return *hd++;
 11 }
 12 
 13 inline int nextInt(void) {
 14     register int ret = 0;
 15     register int neg = false;
 16     register int bit = nextChar();
 17     for (; bit < 48; bit = nextChar())
 18         if (bit == '-')neg ^= true;
 19     for (; bit > 47; bit = nextChar())
 20         ret = ret * 10 + bit - 48;
 21     return neg ? -ret : ret;
 22 }
 23 
 24 template <class T>
 25 inline T min(T a, T b) {
 26     return a < b ? a : b;
 27 }
 28 
 29 template <class T>
 30 inline T max(T a, T b) {
 31     return a > b ? a : b;
 32 }
 33 
 34 const int inf = 1e9;
 35 const int siz = 1000005;
 36 
 37 int tot;
 38 int s, t;
 39 int edges;
 40 int hd[siz];
 41 int to[siz];
 42 int nt[siz];
 43 int fl[siz];
 44 int vl[siz];
 45 
 46 inline void add(int u, int v, int f, int w) { 
 47     nt[edges] = hd[u]; to[edges] = v; fl[edges] = f; vl[edges] = +w; hd[u] = edges++;
 48     nt[edges] = hd[v]; to[edges] = u; fl[edges] = 0; vl[edges] = -w; hd[v] = edges++;
 49 }
 50 
 51 int dis[siz];
 52 int pre[siz];
 53 
 54 inline bool spfa(void) {
 55     static int que[siz];
 56     static int inq[siz];
 57     static int head, tail;
 58     for (int i = 0; i <= tot; ++i)
 59         dis[i] = inf, inq[i] = 0;
 60     dis[s] = 0;
 61     pre[s] = -1;
 62     head = 0, tail = 0;
 63     inq[que[tail++] = s] = 1;
 64     while (head != tail) {
 65         int u = que[head++], v; inq[u] = 0; 
 66         for (int i = hd[u]; ~i; i = nt[i])
 67             if (fl[i] && dis[v = to[i]] > dis[u] + vl[i]) { 
 68                 dis[v] = dis[u] + vl[i], pre[v] = i ^ 1;
 69                 if (!inq[v])inq[que[tail++] = v] = 1;
 70             }
 71     }
 72     return dis[t] < inf;
 73 }
 74 
 75 inline int flow(void) {
 76     int flow = inf;
 77     for (int i = pre[t]; ~i; i = pre[to[i]])
 78         flow = min(flow, fl[i ^ 1]);
 79     for (int i = pre[t]; ~i; i = pre[to[i]])
 80         fl[i] += flow, fl[i ^ 1] -= flow;
 81     return dis[t] * flow;
 82 }
 83 
 84 int n, m, sum;
 85 int tim[50][150];
 86 
 87 int kth[siz];
 88 int lst[siz];
 89 int edg[siz];
 90 
 91 signed main(void) 
 92 {
 93     for (int i = 0; i < siz; ++i)hd[i] = -1;
 94     
 95     n = nextInt();
 96     m = nextInt();
 97     
 98     s = 0, t = tot = n + 1;
 99     
100     for (int i = 1; i <= n; ++i) {
101         int cnt = nextInt();
102         add(s, i, cnt, 0);
103         sum += cnt;
104     }
105     
106     for (int i = 1; i <= n; ++i)
107         for (int j = 1; j <= m; ++j)
108             tim[i][j] = nextInt();
109             
110     for (int i = 1; i <= m; ++i) {
111         int p = ++tot;
112         for (int j = 1; j <= n; ++j)
113             add(j, p, 1, tim[j][i]);
114         add(p, t, 1, 0);
115         kth[i] = 1;
116         lst[i] = p;
117         edg[p] = edges - 2;
118     }
119     
120     int answer = 0;
121             
122     while (sum--) {
123         if (spfa())
124             answer += flow();
125         for (int i = 1; i <= m; ++i)
126             if (!fl[edg[lst[i]]]) {
127                 int p = ++tot; ++kth[i];
128                 for (int j = 1; j <= n; ++j)
129                     add(j, p, 1, tim[j][i] * kth[i]);
130                 add(p, t, 1, 0);
131                 lst[i] = p;
132                 edg[p] = edges - 2;
133                 break;
134             }
135     }
136     
137     printf("%d
", answer);
138 }

@Author: YouSiki