bzoj 2879: [Noi2012]美食节（费用流+动态加边）

Time Limit:10 SecMemory Limit:512 MB
Submit:1453Solved:778

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

题解：费用流+动态加边

刚开始写的单纯的费用流TLE了。

于是只能用一些高大上的东西——动态加点。

其实这道题除了数据范围，与scoi 2007 修车就是一个题。

建图

源点->所有厨师拆成的点（一个厨师sigma pi ）容量为1，费用为0

每种菜->汇点 容量为1，费用为0

本来应该一个厨师拆成的每个点都与不同的菜连容量为1，费用为t[i][j]×k 的边（k表示倒数第几个做，因为倒数第一做只有这个人需要等,做这个菜的时间不会其他菜） ，但是我们发现其实这些边最后一定有很多是不会走到的。

于是我们刚开始的时候只需要把倒数第一人，即t[i][j]连好就可以了。

每次增广的时候，其实只要一个是做菜的其余都是退菜，所以我们需要找出源点的后一个点，然后计算出他属于哪个厨师做的第几道菜，比如是第一个厨师的第2道菜，那么我们就把所有菜到第一个厨师的第3道菜的边联通，再继续增广就可以了。

 

 

#include  
#include  
#include  
#include  
#include    
#define N 100003  
#define M 1200000  
#define inf 1000000000  
using namespace std;  
int n,m,tot,mincost;  
int next[M],v[M],remain[M],c[M],sum,q[N];  
int point[N],can[N],dis[N],last[N],p[N],T[50][120];  
void add(int x,int y,int z,int k)  
{  
    tot++; next[tot]=point[x]; point[x]=tot; v[tot]=y; remain[tot]=z; c[tot]=k;  
    tot++; next[tot]=point[y]; point[y]=tot; v[tot]=x; remain[tot]=0; c[tot]=-k;  
}  
int addflow(int s,int t)  
{  
    int now=t; int ans=inf;  
    int k,a,b;  
    while (now!=s)  
    {  
        ans=min(ans,remain[last[now]]);  
        if (v[last[now]^1]==s)  
         {  
            k=v[last[now]];  
            a=(k-1)/sum+1; b=(k)%sum+1;  
         }  
        now=v[last[now]^1];  
    }  
    now=t;  
    while (now!=s)  
    {  
        remain[last[now]]-=ans;  
        remain[last[now]^1]+=ans;  
        now=v[last[now]^1]; 
    }  
    mincost+=ans*dis[t];
    for (int i=1;i<=m;i++)  
     add((a-1)*sum+b,n*sum+i,1,b*T[i][a]);  
}  
bool spfa(int s,int t)  
{  
    for (int i=s;i<=t;i++)  
     dis[i]=inf,can[i]=0;
    dis[s]=0; can[s]=1;  
    int head=0,tail=0;
    q[++tail]=s;
    while (headdis[now]+c[i]&&remain[i])  
         {  
            dis[v[i]]=dis[now]+c[i];  
            last[v[i]]=i;  
            if (!can[v[i]])  
             {  
                can[v[i]]=1;  
                q[++tail]=v[i];
             }  
         }  
        can[now]=0;  
    }  
    if (dis[t]==inf) return false;  
    return true;  
}  
int main()  
{   
    scanf("%d%d",&m,&n);  
    tot=-1;  sum=0;  
    memset(next,-1,sizeof(next)); memset(point,-1,sizeof(point));  
    for (int i=1;i<=m;i++) scanf("%d",&p[i]),sum+=p[i];  
    int num=100001;  
    for (int i=1;i<=m;i++)  
     for (int j=1;j<=n;j++) scanf("%d",&T[i][j]);  
	for(int i=1;i<=n*sum;i++)
        add(0,i,1,0);
    for(int i=1;i<=m;i++)
        add(n*sum+i,num,p[i],0);
    for(int i=1;i<=n;i++)
       for(int k=1;k<=m;k++)
          add((i-1)*sum+1,n*sum+k,1,T[k][i]);
    while (spfa(0,num))  
	  addflow(0,num);  
    printf("%d\n",mincost);  
    return 0;  
}