「BZOJ2879」[Noi2012] 美食节（费用流+动态加点）

题目

传送门

题解

这个和[SCOI2007] 修车实际上是一道题，用到了动态加点；
所谓的动态加点，刚开始每个厨师都连倒数第一个做的菜。跑完最大流后，对于找到增广路的厨师再给安排倒数第二个菜……这可以放在MCMF函数中实现。

代码

#include<iostream>  
#include<cstdio>  
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
#define inf 0x7fffffff
#define N 100005
#define M 7000000
using namespace std;

int n,m,p,ans,cnt,last[N],vis[N],pre[N],dis[N],a[50][120],sum[50],s,t;
struct edge{int to,from,op,c,w,next;}e[M];
queue <int> q;

void insert(int u,int v,int c,int w)
{  
    e[++cnt].to=v;e[cnt].from=u;e[cnt].c=c;e[cnt].w=w;e[cnt].op=cnt+1;e[cnt].next=last[u];last[u]=cnt;  
    e[++cnt].to=u;e[cnt].from=v;e[cnt].c=0;e[cnt].w=-w;e[cnt].op=cnt-1;e[cnt].next=last[v];last[v]=cnt;  
}
bool spfa()  
{  
    for (int i=s;i<=t;i++)  
    {  
        dis[i]=inf/3;  
        vis[i]=pre[i]=0;  
    }  
    vis[s]=1;  
    dis[s]=0;  
    while (!q.empty()) q.pop();  
    q.push(s);  
    while (!q.empty())  
    {  
        int u=q.front();  
        q.pop();  
        for (int i=last[u];i;i=e[i].next)  
            if (e[i].c&&dis[u]+e[i].w<dis[e[i].to])  
            {  
                pre[e[i].to]=i;  
                dis[e[i].to]=e[i].w+dis[u];  
                if (!vis[e[i].to])  
                {  
                    q.push(e[i].to);  
                    vis[e[i].to]=1;  
                }  
            }  
        vis[u]=0;  
    }  
    if (dis[t]==inf/3) return 0;  
    return 1;  
}  

void mcf()  
{  
    int mnf=inf,i=t;  
    while (pre[i])  
    {  
        mnf=min(mnf,e[pre[i]].c);  
        i=e[pre[i]].from;  
    }  
    ans+=mnf*dis[t];  
    i=t;  
    while (pre[i])  
    {  
        e[pre[i]].c-=mnf;  
        e[e[pre[i]].op].c+=mnf;  
        if (e[pre[i]].to==t&&(e[pre[i]].from-n)%p>0)  
        {  
            insert(e[pre[i]].from+1,t,1,0);  
            int kind=(e[pre[i]].from-n)/p+1,num=(e[pre[i]].from-n)%p+1;  
            for (int j=1;j<=n;j++)  
                insert(j,e[pre[i]].from+1,1,a[j][kind]*num);
        }
        i=e[pre[i]].from;
    }
}

int main()
{  
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
    {
        scanf("%d",&sum[i]);
        p+=sum[i];
    }
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            scanf("%d",&a[i][j]);
    s=0;t=n+m*p+1;
    for (int i=1;i<=n;i++)
        insert(s,i,sum[i],0);
    for (int i=1;i<=m;i++)
    {
        for (int j=1;j<=n;j++)
            insert(j,n+(i-1)*p+1,1,a[j][i]);
        insert(n+(i-1)*p+1,t,1,0);
    }
    while (spfa()) mcf();
    printf("%d",ans);
    return 0;
}

总结

说起来容易做起来难，码力太差