bzoj2879: [Noi2012]美食节 最大流

题目大意：
有n种菜，需求为p_i。有m名厨师，第j名厨师做第i种菜的时间为t_ij。各厨师独立做菜，最小化上菜等待时间之和。
简单题解：
费用流（跪翱犇）。
一名厨师做的最后一个菜的时间被算了一次，第二个两次，依次类推。使用这种方式来统计答案。
初始时，菜向厨师连边，表示最后一次做菜。

对当前增广的厨师，新建一个点表示其前一次做菜，此时费用应为两倍；再次新建即为三倍，依次类推。

注意last数组可以省很多的时间，即记录第一个没有满流的地点，形成一个凸包。

#include <cstdio>
#include <queue>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define maxm 3000005
#define maxn 100010
#define INF 0x3f3f3f3f
struct node
{
    int u,v,f,c,next;
}e[maxm];
int k,head[maxn],pre[maxn],dist[maxn],vis[maxn],last[maxn];
int next[maxn],cnt[maxn],to[maxn],cost[maxn],flow[maxn];
int en=1,s,t,maxflow,mincost,m,n,num,sum;
int save[500][500],need[500];
void add(int u,int v,int f,int c)
{
    en++;
    e[en].u=u;
    e[en].v=v;
    e[en].c=c;
    e[en].f=f;
    e[en].next=head[u];
    head[u]=en;

    en++;
    e[en].u=v;
    e[en].v=u;
    e[en].c=-c;
    e[en].f=0;
    e[en].next=head[v];
    head[v]=en;
}
int q[maxm];
int spfa()
{
    int i,u,v;
    for(i=0;i<=t;i++)
        dist[i]=INF;
    dist[s]=0;
    vis[s]=1;
    int H=0,T=0;
    T++;q[T]=s;
    while(H<T)
    {
        H++;
        u=q[H];
        for(i=head[u];i;i=e[i].next)
        {
            v=e[i].v;
            if(e[i].f>0&&dist[u]+e[i].c<dist[v])
            {
                dist[v]=dist[u]+e[i].c;
                pre[v]=i;
                if(!vis[v])
                {
                    vis[v]=1;
                    T++;
                    q[T]=v;
                }
            }
        }
        vis[u]=0;
    }
}
void add2()
{
    int v,x,y,a,b;
    int maxf=INF;
    for(v=pre[t];v;v=pre[e[v].u])
    {
        e[v].f--;
        e[v^1].f++;
    }
    mincost+=dist[t];
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&need[i]);
        sum+=need[i];
    }
    for(int i=1;i<=n;i++)
        add(s,i,need[i],0);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&save[i][j]);
            add(i,n+j,1,save[i][j]);
        }
    }
        t=n+m+sum+1;
        for(int i=1;i<=m;i++)
        {
                cnt[i]=1;
                add(n+i,t,1,0);
                last[i]=en-1;
        }
        for(int j;sum;sum--)
        {
                spfa();
                add2();
                for(j=1;j<=m && e[last[j]].f;j++);
                
                ++cnt[j];
                for(int i=1;i<=n;++i)
                        add(i,n+m+sum,1,save[i][j]*cnt[j]);
                add(n+m+sum,t,1,0);
                last[j]=en-1;
        }
        printf("%d\n",mincost);
        return 0;
}
/*
3 10 2 3 3 2
5
6
7
*/