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

正解：费用流。

这题除了数据范围跟修车完全一样啊。。

然而我在考场上知道要动态加边，但是没考虑顺序，每次都加了m个厨师，所以跑不出来，然后乖乖打了暴力。。

还有费用流。。这题如果spfa不加标记是会T飞的。。我以前从来没有打过标记。。

思路跟修车的一样，然后我们只要加上一个优化，那就是每次增广以后记录一下增广的是哪个厨师，然后对于那个厨师就新建一个时间加1的点，连到图里，再跑一遍增光路就好了。如果流量已经到了p了，就直接输出，不要再跑了。

//It is made by wfj_2048~
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <vector>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <set>
#define inf (1<<30)
#define il inline
#define RG register
#define ll long long
#define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)

using namespace std;

struct edge{ int nt,to,flow,cap,cost; }g[5000010];

int head[1000010],dis[1000010],fa[1000010],f[1000010],vis[1000010],vi[1000010],p[1000010],q[10000010],tim[1010],pp[50],v[1010][1010],c[1010][1010],n,m,S,T,sz,tot,flow,cost,last,num=1;

il int gi(){
    RG int x=0,q=1; RG char ch=getchar(); while ((ch<'0' || ch>'9') && ch!='-') ch=getchar();
    if (ch=='-') q=-1,ch=getchar(); while (ch>='0' && ch<='9') x=x*10+ch-48,ch=getchar(); return q*x;
}

il void insert(RG int from,RG int to,RG int cap,RG int cost){ g[++num]=(edge){head[from],to,0,cap,cost},head[from]=num; }

il int spfa(RG int S,RG int T){
    for (RG int i=1;i<=tot;++i) dis[i]=inf; dis[S]=0,dis[T]=inf;
    RG int h=0,t=1; q[t]=S,f[S]=inf,vis[S]=1;
    while (h<t){
    RG int x=q[++h];
    for (RG int i=head[x];i;i=g[i].nt){
        RG int v=g[i].to;
        if (dis[v]>dis[x]+g[i].cost && g[i].cap>g[i].flow){
        dis[v]=dis[x]+g[i].cost,fa[v]=x,pp[v]=i;
        f[v]=min(f[x],g[i].cap-g[i].flow);
        if (!vis[v]) vis[v]=1,q[++t]=v;
        }
    }
    vis[x]=0;
    }
    if (dis[T]==inf) return 0; flow+=f[T],cost+=f[T]*dis[T],last=0;
    for (RG int x=T;x!=S;x=fa[x]){
    g[pp[x]].flow+=f[T],g[pp[x]^1].flow-=f[T];
    if (!last) last=vi[x];
    }
    return 1;
}

il void work(){
    n=gi(),m=gi(),tot=n,S=1000001,T=1000002;
    for (RG int i=1;i<=n;++i) p[i]=gi(),sz+=p[i],insert(S,i,p[i],0),insert(i,S,0,0);
    for (RG int i=1;i<=n;++i) for (RG int j=1;j<=m;++j) v[i][j]=gi();
    for (RG int i=1;i<=m;++i) c[i][1]=++tot,vi[tot]=i,insert(c[i][1],T,1,0),insert(T,c[i][1],0,0);
    for (RG int i=1;i<=n;++i) for (RG int j=1;j<=m;++j) insert(i,c[j][1],1,v[i][j]),insert(c[j][1],i,0,-v[i][j]);
    spfa(S,T); for (RG int i=1;i<=m;++i) tim[i]=1;
    while (flow<sz){
    c[last][++tim[last]]=++tot,vi[tot]=last,insert(c[last][tim[last]],T,1,0),insert(T,c[last][tim[last]],0,0);
    for (RG int i=1;i<=n;++i)
        insert(i,c[last][tim[last]],1,tim[last]*v[i][last]),insert(c[last][tim[last]],i,0,-tim[last]*v[i][last]);
    spfa(S,T);
    }
    printf("%d
",cost); return;
    return;
}

int main(){
    File("delicacy");
    work();
    return 0;
}