【noi 2012】美食节

题目描述
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）。他想知道的是最小的总等待时间是多少。
输入格式：
输入文件的第1行包含两个正整数n和m，表示菜品的种数和厨师的数量。 第2行包含n个正整数，其中第i个数为pi，表示点第i种菜品的人数。 接下来有n行，每行包含m个非负整数，这n行中的第i行的第j个数为ti,j，表示第j个厨师制作第i种菜品所需的时间。 输入文件中每行相邻的两个数之间均由一个空格隔开，行末均没有多余空格。
输出格式：
输出仅一行包含一个整数，为总等待时间的最小值。

这道题跟SCOI2007修车很像，就是数据范围很大，80000个点，3200000条边，直接跑费用流会T飞。
考虑一下，其实大部分边是用不上的。由于我们跑一次spfa只能找出一次增广路，所以我们可以暂时不连不需要的边。一开始把所有厨师做倒数第１道菜与所有菜连好，然后找一条增广路，这条增广路上一定经过了一个点，表示第j个厨师做倒数第１道菜，于是我们添加点（第j个厨师做倒数第２道菜），与汇点和所有菜连边，以此类推。
可以发现每次增广需要的边都已经加好了。
代码（人丑常数大）：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
#define S 80100
#define T 80101
using namespace std;
struct edge{
    int pre,next,to,flow,cost;
}ed[20000010];
pair <int,int> yuan[81010];
int n,m;
int ind[110][810],t[500][1100],vis[81100],dis[81100],from[81100],p[500],size=-1,head[81100];
void add(int from,int to,int flow,int cost)
{
    size++;
    ed[size].to=to;
    ed[size].next=head[from];
    ed[size].flow=flow;
    ed[size].cost=cost;
    ed[size].pre=from;
    head[from]=size;
}
void graph()
{
    int tmp=ed[from[T]].pre;
    int x=yuan[tmp].first;
    int y=yuan[tmp].second;
    y++;
    for(int i=1;i<=n;i++)
        add(i,ind[x][y],1,t[i][x]*y),add(ind[x][y],i,0,-t[i][x]*y);
}
bool spfa()
{
    queue<int> q;
    memset(dis,0x7f,sizeof(dis));
    int inf=dis[0];
    q.push(S);
    vis[S]=1;
    dis[S]=0;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        vis[u]=0;
        for(int i=head[u];~i;i=ed[i].next)
        {
            int v=ed[i].to;
            if(dis[v]>dis[u]+ed[i].cost&&ed[i].flow)
            {
                dis[v]=dis[u]+ed[i].cost;
                from[v]=i;
                if(!vis[v])
                {
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
    if(dis[T]<inf) return true;
    else return false;
}
int mcf()
{
    int ans=0;
    while(spfa())
    {
        int x=1e9;
        for(int i=from[T];~i;i=from[ed[i].pre])
            x=min(x,ed[i].flow);
        for(int i=from[T];~i;i=from[ed[i].pre])
        {
            ed[i].flow-=x;
            ed[i^1].flow+=x;
            ans+=x*ed[i].cost;
        }
        graph();
    }
    return ans;
}
int main()
{
    memset(from,-1,sizeof(from));
    memset(head,-1,sizeof(head));
    scanf("%d%d",&n,&m);
    int cnt=n+1,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",&t[i][j]);
    for(int i=1;i<=m;i++)
        for(int j=1;j<=sum;j++)
            ind[i][j]=++cnt,yuan[cnt]=make_pair(i,j);
    for(int i=1;i<=n;i++) add(S,i,p[i],0),add(i,S,0,0);
    for(int i=1;i<=m;i++)
        for(int j=1;j<=sum;j++)
            add(ind[i][j],T,1,0),add(T,ind[i][j],0,0);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            add(i,ind[j][1],1,t[i][j]),add(ind[j][1],i,0,-t[i][j]);
    printf("%d",mcf());
    return 0;
}