[NOI2012]美食节(费用流)

题目描述

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

题解

一开始按照修车那题的思路搞了个费用流，就是把厨师拆成tim个点，表示他是倒数第几个做的，然后边权就是对后面的人的代价和，然后T到飞起。

原因是直接把图建出来导致图的规模过大，然后spfa就死了。

所以我们可以动态加边，每增广一次就加一个点。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#define N 49
#define M 109
#define P 809
#define R register
using namespace std;
queue<int>q;
int p[N],head[M*P],tot=1,dis[M*P],pre[M*P],fl[M*P],ans,n,m,sum,a[N][M],tim[M*P],id[M],be[M*P];
bool vis[M*P];
inline int rd(){
    int x=0;char c=getchar();bool f=0;
    while(!isdigit(c)){if(c=='-')f=1;c=getchar();}
    while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
    return f?-x:x;
} 
struct edge{int n,to,l,f;}e[N*M*P*2];
inline void add(int u,int v,int l,int f){
    e[++tot].n=head[u];e[tot].to=v;head[u]=tot;e[tot].l=l;e[tot].f=f;
    e[++tot].n=head[v];e[tot].to=u;head[v]=tot;e[tot].l=0;e[tot].f=-f;
}
inline bool spfa(int s,int t){
    memset(dis,0x3f,sizeof(dis));
    q.push(s);dis[s]=0;fl[s]=2e9;
    while(!q.empty()){
        int u=q.front();q.pop();vis[u]=0;
    //    cout<<u<<" "<<dis[u]<<endl;
        for(R int i=head[u];i;i=e[i].n){
            int v=e[i].to;
            if(e[i].l&&dis[v]>dis[u]+e[i].f){
                dis[v]=dis[u]+e[i].f;pre[v]=i;fl[v]=min(fl[u],e[i].l);
                if(!vis[v]){vis[v]=1;q.push(v);}
            }
        } 
    }
    return dis[t]!=0x3f3f3f3f;
}
inline void calc(int s,int t){
    int x=t;
    while(x!=s){
        int i=pre[x];
        e[i].l-=fl[t];e[i^1].l+=fl[t];x=e[i^1].to;
    }
    ans+=dis[t]*fl[t];
}
int main(){
    n=rd();m=rd();
    for(R int i=1;i<=n;++i)p[i]=rd(),sum+=p[i];
    for(R int i=1;i<=n;++i)for(R int j=1;j<=m;++j)a[i][j]=rd();
    int now=n;
    for(R int i=1;i<=m;++i)for(R int j=1;j<=sum;++j){
      ++now;tim[now]=j==1?1:tim[now-1]+1;if(j==1)id[i]=now;be[now]=i;
    }
    for(R int i=1;i<=n;++i)for(R int j=1;j<=m;++j)add(i,id[j],1,a[i][j]);
    for(R int i=1;i<=m;++i)add(id[i],now+1,1,0);
    for(R int i=1;i<=n;++i)add(0,i,p[i],0);
    while(spfa(0,now+1)){
      calc(0,now+1);
      int x=e[pre[now+1]^1].to;x++;
      add(x,now+1,1,0);
      for(int i=1;i<=n;++i)add(i,x,1,tim[x]*a[i][be[x]]);
    }
    cout<<ans;
    return 0;
}