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

题解：

本题由于数据范围原因，不能像SCOI修车那道题那样直接一下子全部把边建好，所以此题的关键是要动态加边。

下面说一下本题的建边方法。首先从源点向所有的菜建流量为p，费用为0的边。记sigma pi为sum，然后把每个厨师拆成sum个点（显然一个厨师最多做sum道菜，也就是所有菜都是他一个人做），对于每个厨师，拆出来的第i个点表示这个厨师倒数第i次做菜（与修车类似）。

每个点表示的是倒数第几次做菜，因为我们发现，对于一个厨师做的倒数第i道菜，会有i个人需要等待这道菜的时间，也就是要这个厨师做倒数第i到倒数第1道菜的这i个人需要花费等待这道菜的时间。所以第i个点与每道菜的费用应是i*pi，但是我们先不急着连，因为本题是要动态加边的嘛。

一开始并不把所有拆出来的表示厨师的点与汇点连边，我们一开始只把每个厨师倒数第一次做菜的点与汇点连流量为1费用为0的边（每次只能做一道菜），并且把每道菜与每个厨师连流量为1，费用为pi*1的边。

这样之后，我们变开始增广，每增广完一次，我们就找到是哪个厨师做了菜，当前这个厨师做完了倒数第i道菜，我们就把表示他倒数第i+1次做菜的点与汇点连流量为1，费用为0的边，并把每一道菜与这个点连流量为1，费用为（i+1）*pi的边。一直增广下去就行。

下面是代码：

#include <bits/stdc++.h>
using namespace std;

int n,m,st,ed,hed[100001],cnt,ti[501][501],sum,ji[10001]; 
struct node
{
    int next,from,to,c,cost;
}a[5000001];
int v[81001],w[81001],f[100001],inq[100001],q[100001],h,t,book,ans;
inline void add(int from,int to,int c,int cost)
{
    a[++cnt].from=from;
    a[cnt].to=to;
    a[cnt].c=c;
    a[cnt].cost=cost;
    a[cnt].next=hed[from];
    hed[from]=cnt;
    a[++cnt].from=to;
    a[cnt].to=from;
    a[cnt].c=0;
    a[cnt].cost=-cost;
    a[cnt].next=hed[to];
    hed[to]=cnt;
}
inline void bfs()
{
    memset(v,0x3f,sizeof(v));
    memset(w,0,sizeof(w));
    q[1]=st;
    w[st]=2e9;
    v[st]=0;
    h=1;t=2;
    while(h!=t)
    {
        int x=q[h];
        inq[x]=0;
        for(int i=hed[x];i;i=a[i].next)
        {
            int y=a[i].to;
            if(a[i].c&&v[y]>v[x]+a[i].cost)
            {
                w[y]=min(a[i].c,w[x]);
                v[y]=v[x]+a[i].cost;
                f[y]=i;
                if(!inq[y])
                {
                    q[t++]=y;
                    inq[y]=1;
                }
            }
        }
        h++;
    }
    book=a[f[ed]].from;
    for(int i=f[ed];i;i=f[a[i].from])
    {
        a[i].c-=w[ed];
        a[i^1].c+=w[ed];
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    cnt=1;
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&ji[i]);
        sum+=ji[i];
    }
    st=n+sum*m+1;
    ed=st+1;
    for(int i=1;i<=n;i++)
    add(st,i,ji[i],0);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&ti[i][j]);
            add(i,(j-1)*sum+n+1,1,ti[i][j]);
        }
    }
    for(int i=1;i<=m;i++)
    {
        add((i-1)*sum+n+1,ed,1,0);
    }
    while(1)
    {	
        bfs();
        if(w[ed]>0)
        {
            ans+=w[ed]*v[ed];
            if(book+1!=st)//注意这个地方，我被这个卡了很久。。。 
            add(book+1,ed,1,0);
            for(int i=1;i<=n;i++)
            {
                add(i,book+1,1,ti[i][(book+1-n-1)/sum+1]*((book+1-n-1)%sum+1));	
            }
        }
        else
        break;
    }
    printf("%d\n",ans);
    return 0;
}