网络流/费用流

　　跟 BZOJ 1070 修车几乎是一道题，只是这题“要修的车”（即菜）多了很多……几乎是从$n$变成了$n^2$，所以建图的时候就得动态加点……

　　也就是说，当一个厨师已经确定了他的后k道菜的时候，再加入表示倒数第k+1道菜的节点。

　　我个xx一开始写成了每次加一层……不 T 才怪……

　　加点过程看min_cost_flow（即mcf）这个函数吧，每次通过找流过来的节点，来确定要加的点。

 /**************************************************************

     Problem: 2879

     User: Tunix

     Language: C++

     Result: Accepted

     Time:7416 ms

     Memory:72872 kb

 ****************************************************************/

 //BZOJ 2879

 #include<cmath>

 #include<vector>

 #include<cstdio>

 #include<cstring>

 #include<cstdlib>

 #include<iostream>

 #include<algorithm>

 #define rep(i,n) for(int i=0;i<n;++i)

 #define F(i,j,n) for(int i=j;i<=n;++i)

 #define D(i,j,n) for(int i=j;i>=n;--i)

 #define pb push_back

 #define CC(a,b) memset(a,b,sizeof(a))

 using namespace std;

 int getint(){

     int v=,sign=; char ch=getchar();

     while(!isdigit(ch)) {if(ch=='-') sign=-; ch=getchar();}

     while(isdigit(ch))  {v=v*+ch-''; ch=getchar();}

     return v*sign;

 }

 const int N=,M=,INF=~0u>>;

 const double eps=1e-;

 /*******************template********************/

 int n,m,ans,tot,tim[][],p[];

 struct edge{int from,to,v,c;};

 struct Net{

     edge E[M];

     int head[N],next[M],cnt;

     void ins(int x,int y,int z,int c){

         E[++cnt]=(edge){x,y,z,c};

         next[cnt]=head[x]; head[x]=cnt;

     }

     void add(int x,int y,int z,int c){

         ins(x,y,z,c); ins(y,x,,-c);

     }

     int from[N],Q[M],d[N],S,T,ed;

     bool inq[N],sign;

     bool spfa(){

         int l=,r=-;

         F(i,,T) d[i]=INF;

         d[S]=; Q[++r]=S; inq[S]=;

         while(l<=r){

             int x=Q[l++];

             inq[x]=;

             for(int i=head[x];i;i=next[i])

                 if(E[i].v> && d[x]+E[i].c<d[E[i].to]){

                     d[E[i].to]=d[x]+E[i].c;

                     from[E[i].to]=i;

                     if (!inq[E[i].to]){

                         Q[++r]=E[i].to;

                         inq[E[i].to]=;

                     }

                 }

         }

         return d[T]!=INF;

     }

     void mcf(){

         int x=INF,a,b;

         for(int i=from[T];i;i=from[E[i].from]){

             x=min(x,E[i].v);

             if (E[i].to==T){

                 a=E[i].from%tot+;

                 b=(E[i].from-)/tot;

             }

         }

         for(int i=from[T];i;i=from[E[i].from]){

             E[i].v-=x;

             E[i^].v+=x;

         }

         ans+=x*d[T];

         if (a!=)

             F(i,,n)

                 add(i,b*tot+a,,a*tim[i][b]);

     }

     void init(){

         n=getint(); m=getint(); cnt=ed=; tot=;

         F(i,,n) {p[i]=getint(); tot+=p[i];}

         S=; T=tot+tot*m+;

         F(i,,n) F(j,,m) tim[i][j]=getint();

         F(i,,n) add(S,i,p[i],);

         F(i,tot+,tot+tot*m) add(i,T,,);

         F(i,,n)

             F(j,,m)

                 add(i,j*tot+,,tim[i][j]);

         ans=;

         while(spfa()) mcf();

         printf("%d\n",ans);

     }

 }G1;

 int main(){

 #ifndef ONLINE_JUDGE

     freopen("2879.in","r",stdin);

     freopen("2879.out","w",stdout);

 #endif

     G1.init();

     return ;

 }

确定

【BZOJ】【2879】【NOI2012】美食节

网络流/费用流