[BZOJ2879][NOI2012]美食节（时间费用流+动态加边）

题目：

我是超链接

题解：

这题目一眼费用流了，以前好像有一道基本一样的题目
简单说也就是每个菜建一个点和源点连容量1，厨师*菜数个点表示第i个厨师倒数第j个菜做了啥
但是这个东西数据范围有点大，光点就是80000个，边多的上天了，根本不行。
考虑动态加边，首先对于每一道菜给每个厨师的倒数第一个菜连边，因为倒数第k个菜肯定比倒数第k+1个菜便宜。
每次增广路都只找一条最短路，所以对于这条找的最短路的厨师，我们再给所有的菜连他的倒数第二个点，同时连汇点

代码：

#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
#define INF 1e9
const int N=90000;
int tot,point[N],dis[N],pre[N],v[N*2],c[N*2],remind[N*2],nxt[N*2],a[45] [105],p[45];bool vis[N];
int n,m,sum;
void addline(int x,int y,int z,int vv)
{
    ++tot; nxt[tot]=point[x]; point[x]=tot; v[tot]=y; remind[tot]=z; c[tot]=vv;
    ++tot; nxt[tot]=point[y]; point[y]=tot; v[tot]=x; remind[tot]=0; c[tot]=-vv;
}
bool spfa(int s,int t)
{
    memset(dis,0x7f,sizeof(dis));
    memset(vis,0,sizeof(vis));
    dis[s]=0;queue<int>q;
    q.push(s);
    while (!q.empty())
    {
        int now=q.front(); q.pop(); vis[now]=0;
        for (int i=point[now];i!=-1;i=nxt[i])
          if (dis[v[i]]>dis[now]+c[i] && remind[i])
          {
            dis[v[i]]=dis[now]+c[i]; pre[v[i]]=i; 
            if (!vis[v[i]]) vis[v[i]]=1,q.push(v[i]);
          }
    }
    return dis[t]<INF;
}
int work(int s,int t)
{
    int ans=0;
    while (spfa(s,t))
    {
        int x=t,minn=INF;
        while (x!=s)
        {
            int i=pre[x]; minn=min(minn,remind[i]);
            x=v[i^1];
        }
        x=t;
        while (x!=s)
        {
            int i=pre[x]; 
            remind[i]-=minn; remind[i^1]+=minn;
            x=v[i^1];
        } 
        ans+=dis[t]*minn;
        x=v[pre[t]^1]; addline(x+1,t,1,0);
        for (int i=1;i<=n;i++) addline(i+sum*m,x+1,1,a[i][x/sum+1]*(x%sum+1));
    }
    return ans;
}
int main()
{
    tot=-1;memset(point,-1,sizeof(point));
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++){scanf("%d",&p[i]);sum+=p[i];}
    int s=0,t=m*sum+n+1;
    for (int i=1;i<=n;i++)
    {
        addline(s,i+sum*m,p[i],0);
        for (int j=1;j<=m;j++)
        {
            scanf("%d",&a[i][j]);
            addline(i+sum*m,(j-1)*sum+1,1,a[i][j]);
        }
    }  
    for (int i=1;i<=m;i++) addline((i-1)*sum+1,t,1,0);
    printf("%d",work(s,t));
}