[题解]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

Solution

这题是bzoj1070(SCOI2007)修车（修车的题解）的加强版，数据超级变态，方法基本相同，但是要动态加边。

我们不是把每个厨师拆成sigama(p)个点嘛，为了提高时间效率我们先只能把每个厨师的第一个点连上边，然后每一次增广完之后找到做菜的厨师，激活他的下一个节点（向他的下一个节点连边）。这样可以减少很多不需要的边的计算。即使如此，这道题也有TLE的风险。注意常数！而且一开始不能把所有厨师节点到汇点的边连好，必须在动态加边时再连，否则在除了bzoj以外我试过的题库都要T 2到3个点。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

inline int read(){
    int xx=0,f=1;char ch=getchar();
    for(;ch<'0'||ch>'9';ch=getchar())if(ch=='-')f=-1;
    for(;ch>='0'&&ch<='9';ch=getchar())xx=xx*10+ch-'0';
    return xx*f;
}

const int maxn=100010,oo=0x7fffffff;
struct edge{
    int to,next,size,cost;
}e[maxn*30];
int n,m,num=1,sum=0,top=0,pre[maxn],head[maxn],que[maxn],dis[maxn],p[maxn],T;
int f[45][105];
bool exist[maxn];

void adde(int u,int v,int w,int f){
    e[++num].to=v;e[num].size=w;
    e[num].cost=f;e[num].next=head[u];
    head[u]=num;
}
void add(int u,int v,int w,int f){
    adde(u,v,w,f);adde(v,u,0,-f);
}
bool spfa(int &cost){
    memset(exist,false,sizeof exist);
    memset(dis,127,sizeof dis);
    int maxx=dis[0],h=0,t=1;
    exist[que[t]=0]=true;
    dis[0]=0;
    while(h!=t){
        if((++h)>maxn-5)h=1;
        int x=que[h];
        exist[x]=false;
        for(int i=head[x];i;i=e[i].next){
            if(e[i].size&&dis[e[i].to]>dis[x]+e[i].cost){
                dis[e[i].to]=dis[x]+e[i].cost;
                pre[e[i].to]=i;
                if(!exist[e[i].to]){
                    if((++t)>maxn-5)t=1;
                    que[t]=e[i].to;
                    exist[e[i].to]=true;
                }
            }
        }
    }
    if(dis[T]==maxx)return false;
    int fl=oo;
    for(int i=T;i;i=e[pre[i]^1].to){
        fl=min(fl,e[pre[i]].size);
    }
    for(int i=T;i;i=e[pre[i]^1].to){
        e[pre[i]].size-=fl;
        e[pre[i]^1].size+=fl;
    }
    cost+=dis[T]*fl;
    int x=e[pre[T]^1].to;
    int a=(x-n-1)/sum+1,b=(x-n)%sum+1;
    if(b<=sum){
        for(int i=1;i<=n;i++){
            add(i,n+(a-1)*sum+b,1,b*f[i][a]);
        }
    }
    add(n+(a-1)*sum+b,T,1,0);
    return true;
}

int main(){
    n=read();m=read();
    for(int i=1;i<=n;i++){
        sum+=(p[i]=read());
        add(0,i,p[i],0);
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            f[i][j]=read();
        }
    }
    top=n+sum*m;T=top+1;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            add(i,n+(j-1)*sum+1,1,f[i][j]);
        }
    }
    for(int i=1;i<=m;i++){
        add(n+(i-1)*sum+1,T,1,0);
    }
    int cost=0;
    while(spfa(cost));
    printf("%d\n",cost);
    return 0;
}