【NOI2012】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

输出仅一行包含一个整数，为总等待时间的最小值。

每个厨师建$p$个点，每个点代表做一个菜，然后和菜形成二分图。倒数第$k$个菜的费用是$k*t_{i,j}$。但是这样总共有$O(mp)$个点，过不了。
注意到总共只做$p$道菜，可以在某个厨师做完倒数第$i$道菜后再加入倒数第$i+1$道菜的点。这样点数减少到了$O(m+p)$。
这样显得有问题，因为成了贪心每次优先选择费用最小的，而本来应该先做费用大的再做费用小的。但是实际上反向边退流进行了修正。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int s=50005,t=50006,oo=0x3f3f3f3f,mod=50007;
int fir[50010],ne[500010],to[500010],w[500010],c[500010],
ti[45][110],
bel[50010],id[50010],
que[50010],in[50010],minw[50010],len[50010],fa[50010],
n,m,num,tot,ans;
int node(int u,int x)
{
    tot++;
    bel[tot]=u;
    id[tot]=x;
    return tot;
}
void add(int u,int v,int x,int y)
{
    num++;
    ne[num*2]=fir[u];
    fir[u]=num*2;
    to[num*2]=v;
    w[num*2]=x;
    c[num*2]=y;
    ne[num*2+1]=fir[v];
    fir[v]=num*2+1;
    to[num*2+1]=u;
    w[num*2+1]=0;
    c[num*2+1]=-y;
}
void check()
{
    for (int i=1;i<=t;i++)
        for (int j=fir[i];j;j=ne[j])
            if (w[j]) printf("%d->%d:%d,%d\n",i,to[j],w[j],c[j]);
    printf("------------------------------\n");
}
bool solve()
{
    //check();
    int u,v,hd=0,tl=1,x,y;
    que[0]=s;
    in[s]=1;
    for (int i=1;i<=tot;i++)
    {
        minw[i]=0;
        len[i]=oo;
    }
    minw[t]=0;
    len[t]=oo;
    minw[s]=oo;
    len[s]=0;
    while (hd!=tl)
    {
        u=que[hd++];
        hd%=mod;
        for (int i=fir[u];i;i=ne[i])
            if (w[i]&&len[v=to[i]]>len[u]+c[i])
            {
                len[v=to[i]]=len[u]+c[i];
                fa[v]=i;
                minw[v]=min(minw[u],w[i]);
                if (!in[v])
                {
                    in[v]=1;
                    que[tl++]=v;
                    tl%=mod;
                }
            }
        in[u]=0;
    }
    if (!minw[t]) return 0;
    ans+=minw[t]*len[t];
    for (int i=fa[t];i;i=fa[to[i^1]])
    {
        w[i]-=minw[t];
        w[i^1]+=minw[t];
        if (bel[to[i^1]]) x=to[i^1];
    }
    //x=to[fa[to[fa[t]^1]]^1];
    y=node(bel[x],id[x]+1);
    //printf("node:%d chef:%d num:%d food:%d cost:%d\n",x,bel[x],id[x],to[fa[t]^1],len[t]);
    add(s,y,1,0);
    for (int i=1;i<=n;i++)
        add(y,i,1,ti[i][bel[y]]*id[y]);
    return 1;
}
int main()
{
    int x;
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
    {
        scanf("%d",&x);
        add(i,t,x,0);
    }
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            scanf("%d",&ti[i][j]);
    tot=n;
    for (int i=1;i<=m;i++)
    {
        x=node(i,1);
        add(s,x,1,0);
        for (int j=1;j<=n;j++)
            add(x,j,1,ti[j][i]);
    }
    while (solve());
    printf("%d\n",ans);
}