P2050 [NOI2012]美食节（费用流）

P2050 [NOI2012]美食节

P2053 [SCOI2007]修车的加强版

因为数据较大，一次性把所有边都加完会T

于是我们每次只连需要的边跑费用流

就是开始先连所有厨师做倒数第1道菜

跑费用流后会选一个厨师做菜，于是就把该厨师做第2道菜的情况加进去

每跑一次都这样动态加边

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
#define N 100005
#define M 4000005
int n,m,w[],e[][],v,tot;
int d[N],a[N],p[N],S,T,tC;
queue <int> h; bool inh[N];
int cnt=,hd[N],nxt[M],ed[N],poi[M],val[M],cst[M];
inline void adde(int x,int y,int v1,int v2){
    nxt[ed[x]]=++cnt, hd[x]=hd[x]?hd[x]:cnt,
    ed[x]=cnt, poi[cnt]=y, val[cnt]=v1, cst[cnt]=v2;
}
inline void link(int x,int y,int v1,int v2){adde(x,y,v1,v2),adde(y,x,,-v2);}
bool bfs(){
    memset(d,,sizeof(d)); int inf=d[];
    h.push(S); inh[S]=; d[S]=; a[S]=inf;
    while(!h.empty()){
        int x=h.front(); h.pop(); inh[x]=;
        for(int i=hd[x];i;i=nxt[i]){
            int to=poi[i];
            if(val[i]>&&d[to]>d[x]+cst[i]){
                d[to]=d[x]+cst[i]; p[to]=i;
                a[to]=min(a[x],val[i]);
                if(!inh[to]) h.push(to),inh[to]=;
            }
        }
    }if(d[T]==inf) return ;
    tC+=a[T]*d[T]; int id=poi[p[T]^]+m;//做菜的厨师一定是poi[p[T]^1]%m ，它一定是通过正向边到T
    for(int i=T;i!=S;i=poi[p[i]^]) val[p[i]]-=a[T],val[p[i]^]+=a[T];
    link(id,T,,);
    for(int i=;i<=n;++i) link(tot*m+i,id,,((id-)/m+)*e[i][(id-)%m+]);
    return ;
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=;i<=n;++i) scanf("%d",&w[i]),tot+=w[i];
    S=tot*m+n+; T=S+;
    for(int i=;i<=n;++i) link(S,tot*m+i,w[i],);
    for(int i=;i<=m;++i) link(i,T,,);
    for(int i=;i<=n;++i)
        for(int j=;j<=m;++j)
            scanf("%d",&e[i][j]),link(tot*m+i,j,,e[i][j]);
    while(bfs());
    printf("%d",tC);
    return ;
}