[bzoj2879][网络流-费用流]美食节

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

HINT

【样例说明】 厨师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

题解

这题不是DP想DP的可以看路牌了。。
不考虑数据范围的情况下，这题其实和scoi2007的修车一模一样。设总共需要做的菜为sum个，每个厨师拆成sum个点，第i个点表示这个厨师做倒数第(sum-i+1)个菜。
于是菜向拆点的厨师连边，设菜为i，拆点的厨师原来是第j个，这个菜是倒数第k个做
那么i->(j-1)*sum+(sum-k+1)，流量为1，费用为k*t[i][j]
为什么？让这个菜倒数第k个做的时候，前面对他的影响在前面已经算过了，这里算的是这个菜对后面的影响以及对他自己的影响
所以正确性显然
但是考虑数据范围，这题比那题的数据多了很多。所以要动态加边。每一次费用流跑完之后，由于只会找出一条路径，那么这条路径的倒数第k+1个就建边

#include
#include
#include
#include
#include
#include
using namespace std;
struct node
{
    int x,y,c,d,next,other;
}a[2110000];int len,last[210000];
void ins(int x,int y,int c,int d)
{
    int k1,k2;
    k1=++len;
    a[len].x=x;a[len].y=y;a[len].c=c;a[len].d=d;
    a[len].next=last[x];last[x]=len;
    k2=++len;
    a[len].x=y;a[len].y=x;a[len].c=0;a[len].d=-d;
    a[len].next=last[y];last[y]=len;
    a[k1].other=k2;
    a[k2].other=k1;
}
int d[210000],st,ed;
int v[210000],tt;
queue<int> q;
int pre[210000],tmp[210000];
bool spfa()
{
    memset(d,63,sizeof(d));
    d[st]=0;
    tt++;q.push(st);
    while(!q.empty())
    {
        int x=q.front();
        for(int k=last[x];k;k=a[k].next)
        {
            int y=a[k].y;
            if(a[k].c>0 && d[y]>d[x]+a[k].d)
            {
                d[y]=d[x]+a[k].d;
                pre[y]=k;tmp[y]=x;
                if(v[y]!=tt)
                {
                    v[y]=tt;
                    q.push(y);  
                }
            }
        }
        q.pop();v[x]=tt-1;
    }
    if(d[ed]<2147483647)return true;
    return false;
}
int p[51],ct[51][110];
int n,m,sum;
void mcf()
{
    int cnt=0,ans=0;
    while(spfa())
    {
        int t=ed,minn=999999999,aa,bb;
        while(t!=st)
        {
            minn=min(minn,a[pre[t]].c);
            if(t==ed){aa=(tmp[t]-n-1)/sum+1;bb=(tmp[t]-n-1)%sum+1;}
            t=tmp[t];
        }
        t=ed;
        ans+=d[t]*minn;
        while(t!=st)
        {
            a[pre[t]].c-=minn;a[a[pre[t]].other].c+=minn;
            t=tmp[t];
        }
        cnt+=minn;
        if(cnt==sum)break;
        for(int i=1;i<=n;i++)ins(i,n+sum*(aa-1)+bb-1,1,ct[i][aa]*(sum-(bb-1)+1));
        ins(n+sum*(aa-1)+bb-1,ed,1,0);
    }
    printf("%d\n",ans);
}
int main()
{
    scanf("%d%d",&n,&m);//n种菜 m种厨师 
    sum=0;
    for(int i=1;i<=n;i++)scanf("%d",&p[i]),sum+=p[i];
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)scanf("%d",&ct[i][j]);
    st=n+sum*m+1;ed=st+1;
    for(int i=1;i<=n;i++)ins(st,i,p[i],0);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            ins(i,n+(j-1)*sum+sum,1,ct[i][j]);
    for(int i=1;i<=m;i++)ins(n+(i-1)*sum+sum,ed,1,0);
    mcf();
    return 0;
}