bzoj2879 [Noi2012]美食节

这个题一眼看出是修车，但数据范围有点卡，我傻傻地以为写个zkw就万事大吉了，结果发现zkw还没有spfa跑的快......

这个题的关键是动态加边，就是每次找到已经增广到最新的点的厨师给他加一个新点，这样一共需要增广P次，就可以跑出答案了。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#define maxn 100000
#define maxm 100000
#define inf 2147483647
using namespace std;
struct et
{
    int s,t,val,cost,next;
}e[maxm];
int fir[maxn],dis[maxn],pre[maxn],q[maxm],t[200][200],now[maxn],cnt[maxn];
bool inque[maxn];
int n,m,tot,st,ed,sum,fare,k,ans;

bool find()
{
    int head=0,tail=1;
    for (int i=st;i<=ed;i++) dis[i]=inf;
    q[1]=st; dis[st]=0; inque[st]=1;
    while (head<tail)
    {
        int now=q[++head];
        for (int j=fir[now];j;j=e[j].next)
        {
            int k=e[j].t;
            if (e[j].val&&dis[now]+e[j].cost<dis[k])
            {
                dis[k]=dis[now]+e[j].cost;
                pre[k]=j;
                if (!inque[k]) q[++tail]=k,inque[k]=1;
            }
        }
        inque[now]=0;
    }
}

void add(int x,int y,int z,int w)
{
    e[++tot].s=x; e[tot].t=y; e[tot].val=z; e[tot].cost=w; e[tot].next=fir[x]; fir[x]=tot;
    e[++tot].s=y; e[tot].t=x; e[tot].val=0; e[tot].cost=-w; e[tot].next=fir[y]; fir[y]=tot;
}

void adjust()
{
    for (int j=1;j<=m;j++)
        if (e[now[j]].val==0) {
            cnt[j]++;
            now[j]=tot+1;
            add(n+sum*(j-1)+cnt[j],ed,1,0);
            for (int i=1;i<=n;i++) add(i,n+sum*(j-1)+cnt[j],1,t[i][j]*cnt[j]);
            return ;
        }
}

void fare_flow()
{
    find();
    fare+=dis[ed];
    for (int j=pre[ed];j;j=pre[e[j].s]) e[j].val--,e[j^1].val++;
    adjust();
}

int main()
{
    freopen("delicacy.in","r",stdin);
    freopen("delicacy.out","w",stdout);
    scanf("%d%d",&n,&m);
    tot=1;
    int x;
    st=0; 
    for (int i=1;i<=n;i++) {
        scanf("%d",&x);
        add(st,i,x,0);
        sum+=x;
    }
    ed=m*sum+n+1; 
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
        {
            scanf("%d",&t[i][j]);
            add(i,n+(j-1)*sum+1,1,t[i][j]);
        }
    for (int i=1;i<=m;i++) 
    {
        cnt[i]=1;
        now[i]=tot+1;
        add(n+(i-1)*sum+1,ed,1,0);
    }
    for (int i=sum;i>0;i--) fare_flow();
    printf("%d\n",fare);
    return 0;
}

 

常数大......