hdu 3667 Transportation

最小费用流的变形，这题边上的权值为a*f^2，a为已知系数，f为流值

一开始的想法是，用spfa进行最小权值增广时，判断条件为 dis[j]>dis[i]+a*(2*f+1),wa了（2*f+1=(f+1)^2-f^2)

后来注意到每条边上的容量c<=5，想到可以把一条边拆成c条容量为1的边，边权一次为a,3a,5a...(2*c-1)a

 

今天又想了想，终于发现昨天那个判断条件哪里错了，错在如果产生回流的情况，a为负值，f也为负值，dis值的增量应为 -a*(2*f+1)，修改后AC了，速度还快了一倍

 

自己出的数据不够强，没有出现回流的情况

 

总结：设计网络流方面的算法时，一定要注意回流的情况!还有就是重边！

另外，最小费用流解决问题时，如果因为流值不同会导致单位流量的费用不同，可以将边进行拆分，将单位流量相同的流量拆成一条边

POJ3422也运用到了拆边的思想，可惜当时构图方法是看的别人的，没有从本质上理解

 

代码:(注释部分为拆边的代码)

#include<iostream> #include<memory.h> #include<string> #include<cstdio> #include<algorithm> #include<math.h> #include<queue> using namespace std; const int MAX=105; const int inf=0x7fffffff; struct node { int u,v,b,c,f,next; }g[MAX*5000]; int adj[MAX],pre[MAX],vis[MAX]; __int64 dis[MAX]; int n,m,k,e,s,t,flag; void add(int u,int v,int b,int c) { g[e].u=u; g[e].v=v; g[e].b=b; g[e].c=c; g[e].f=0; g[e].next=adj[u]; adj[u]=e++; g[e].u=v; g[e].v=u; g[e].b=-b; g[e].c=0; g[e].f=0; g[e].next=adj[v]; adj[v]=e++; } void spfa() { int i,j,f,cur; memset(vis,0,sizeof(vis)); queue<int>q; q.push(s); vis[s]=1; while(!q.empty()) { i=q.front(); q.pop(); vis[i]=0; for(cur=adj[i];cur!=-1;cur=g[cur].next) { j=g[cur].v; f=g[cur].f; int tmp; if(f<0) tmp=2*(-f)-1; else tmp=2*f+1; //tmp=1; if(g[cur].c>g[cur].f&&dis[j]>dis[i]+g[cur].b*tmp) { dis[j]=dis[i]+g[cur].b*tmp; pre[j]=cur; if(!vis[j]) { vis[j]=true; q.push(j); } } } } } __int64 mcmf() { int i,minf; while(true) { memset(pre,-1,sizeof(pre)); for(i=0;i<=t;i++) dis[i]=inf; dis[s]=0; spfa(); if(pre[t]==-1) break; for(i=pre[t];i!=-1;i=pre[g[i].u]) { g[i].f+=1; g[i^1].f-=1; } } int j=0; for(i=adj[t];i!=-1;i=g[i].next) { j+=-g[i].f; } if(j<k) { return -1; } else { __int64 ans=0; for(i=0;i<e;i+=2) { ans+=g[i].f*g[i].f*g[i].b; } return ans; } } int main() { int i,j,w,c,x; while(scanf("%d%d%d",&n,&m,&k)!=EOF) { memset(adj,-1,sizeof(adj)); e=0; s=0; t=n+1; /*for(i=0;i<k;i++) { add(s,1,0,1); add(n,t,0,1); }*/ add(s,1,0,k); add(n,t,0,k); while(m--) { scanf("%d%d%d%d",&i,&j,&w,&c); /*for(x=0;x<c;x++) { add(i,j,w*(2*x+1),1); }*/ add(i,j,w,c); } printf("%I64d/n",mcmf()); } return 0; }