2879: [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）。他想知道的是最小的总等待时间是多少。

---

首先回忆一下1070: [SCOI2007]修车

把每个厨师拆成个点\(p\)，第i个点表示做这个厨师的倒数第\(i\)道菜

那么这个点对于答案的贡献就是做这道菜所花时间的\(i\)倍

然后喜提TLE（...）

可以发现在倒数第一个还没有流过的情况下倒数第二个是不可能被流到的！

动态加点就行了

---

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#define MP make_pair
#define TS top().second
#define M 100000
#define N 100000
using namespace std;

priority_queue<pair<int,int> >q;
int cn,a[1000][1000],cc[1000],sum,t,n,m,k,ver[M],edge[M],head[N],nex[M],cnt=1,d[N],h[N],c[M],g[N],x,y,z,s,b[N],cur[N],ans,cost,w[M],bb[M];
void add(int x,int y,int co,int z)
{
    ver[++cnt]=y; nex[cnt]=head[x]; head[x]=cnt; edge[cnt]=z; c[cnt]=co;
    ver[++cnt]=x; nex[cnt]=head[y]; head[y]=cnt; edge[cnt]=0; c[cnt]=-co;
}

bool dji()
{
    while(q.size()) q.pop();
    memset(d,0,sizeof(d)); memset(g,0x3f,sizeof(g)); memset(b,0,sizeof(b));
    d[0]=1; g[0]=0; q.push(MP(0,0));
    while(q.size())
    {
        while(q.size() && b[q.TS]) q.pop();
        if(!q.size()) break;
        int x=q.TS; q.pop(); b[x]=1;
        for(int i=head[x];i;i=nex[i])
        if(edge[i] && g[ver[i]]>g[x]+c[i]+h[x]-h[ver[i]])
        {
            g[ver[i]]=g[x]+c[i]+h[x]-h[ver[i]];
            d[ver[i]]=d[x]+1;
            q.push(MP(-g[ver[i]],ver[i]));
        }
    }
    if(g[t]<0x3f3f3f3f) return 1;
    return 0;
}

int dinic(int x,int flow)
{
    if(x==t || !flow) return flow;
    int re=flow, k;
    for(int& i=cur[x];i && re;i=nex[i])
    if(edge[i] && d[ver[i]]==d[x]+1 && g[ver[i]]==g[x]+c[i]+h[x]-h[ver[i]])
    {
        k=dinic(ver[i],min(re,edge[i]));
        if(k && ver[i]==t) bb[++bb[0]]=w[x];
        re-=k; edge[i]-=k; edge[i^1]+=k;
        if(!k) d[ver[i]]=0;
    } 
    return flow-re;
}

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&k),add(0,i,0,k), sum+=k;
    t=sum*m+n+1; cn=n+m;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&a[i][j]);
            add(i,j+n,a[i][j],1);
        }
    for(int i=1;i<=m;i++) cc[i]=1,add(i+n,t,0,1),w[i+n]=i;
    while(dji())
    {
        bb[0]=0;
        memcpy(cur,head,sizeof(head)); z=ans;
        while(k=dinic(0,0x3f3f3f3f)) ans+=k;
        for(int i=1;i<=cn;i++) {h[i]=g[i]+h[i]; if(g[i]==0x3f3f3f3f) h[i]=0;}
        h[t]+=g[t]; 
        cost+=(ans-z)*h[t];
        for(int i=1;i<=bb[0];i++)
        {
            int x=bb[i]; w[++cn]=bb[i];
            for(int j=1;j<=n;j++) 
                add(j,cn,a[j][x]*(cc[x]+1),1);
            add(cn,t,0,1);
            cc[x]+=1;
        }
    }
    printf("%d",cost);
}