[BZOJ2879] [Noi2012] 美食节 (费用流 & 动态加边)

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

Solution

　　和$BZOJ1070$其实是一道题，只不过这道题强制动态加边，并不知道为什么第$2$个点会$TLE$，总之$BZOJ$可以$AC$

#include <bits/stdc++.h>
using namespace std;
struct edge
{
    int u, v, w, c, nxt;
}e[80005];
int n, fst[1005], t[45][105], etot = 1, level[1005];
int q[1005], belong[1005], fin[105], pth[1005], use;
bool inq[1005];
  
void addedge(int u, int v, int w, int c)
{
    e[++etot] = (edge){u, v, w, c, fst[u]}, fst[u] = etot;
    e[++etot] = (edge){v, u, 0, -c, fst[v]}, fst[v] = etot;
}
  
bool SPFA()
{
    int front = 0, back, u;
    memset(level, 63, sizeof(level));
    level[n + 1] = 0;
    q[back = 1] = n + 1, inq[n + 1] = true;
    while(front != back)
    {
        u = q[(++front % 1005)];
        front %= 1005, inq[u] = false;
        for(int i = fst[u]; i; i = e[i].nxt)
            if(e[i].w && level[e[i].v] > level[u] + e[i].c)
            {
                level[e[i].v] = level[u] + e[i].c;
                pth[e[i].v] = i;
                if(!inq[e[i].v])
                {
                    q[(++back % 1005)] = e[i].v;
                    back %= 1005, inq[e[i].v] = true;
                }
            }
    }
    return level[n + 2] < 1000000000;
}
  
int Edmond_Karp()
{
    for(int i = pth[n + 2]; i; i = pth[e[i].u])
    {
        --e[i].w, ++e[i ^ 1].w;
        if(e[i].u == n + 1) use = e[i].v;
    }
    return level[n + 2];
}
  
int main()
{
    int m, tmp, ans = 0, ptot;
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= n; ++i)
    {
        scanf("%d", &tmp);
        addedge(i, n + 2, tmp, 0);
    }
    for(int i = 1; i <= n; ++i)
        for(int j = 1; j <= m; ++j)
            scanf("%d", &t[i][j]);
    ptot = n + 2;
    for(int i = 1; i <= m; ++i)
    {
        belong[++ptot] = i, ++fin[i];
        for(int j = 1; j <= n; ++j)
            addedge(ptot, j, 1, t[j][i]);
        addedge(n + 1, ptot, 1, 0);
    }
    while(SPFA())
    {
        ans += Edmond_Karp();
        tmp = belong[++ptot] = belong[use], ++fin[tmp];
        for(int i = 1; i <= n; ++i)
            addedge(ptot, i, 1, t[i][tmp] * fin[tmp]);
        addedge(n + 1, ptot, 1, 0);
    }
    printf("%d
", ans);
    return 0;
}