标签:log out sed cst csdn 记录 sdn print size
【算法】网络流-最大流(dinic)
【题解】网络流:http://m.blog.csdn.net/article/details?id=9401909
当前弧优化是因为DFS过程中访问x点时一旦流入量=流出量就退出,所以可以记录下此时正在考虑的弧,下次从此处继续考虑即可。
当前弧之前的弧,不能使流入量-流出量=0,那么一定该弧以及该弧之后的弧中有断裂,那么下次再考虑就没有意义了。
当前弧本身,使流入量-流出量=0,也就是使该点前面的弧中最小的一条断裂了,当前弧以及当前弧连出去之后的弧只是有可能断裂,那么下次就应该从这条当前弧开始考虑。
我的写法中直接用flow表示剩余流量,思路和程序结构参考紫书。
#include<cstdio> #include<algorithm> #include<cstring> using namespace std; const int maxn=210,inf=0x3f3f3f3f; struct edge{int from,v,flow;}e[maxn*3]; int n,m,tot=2,S,T,d[maxn],q[510],first[maxn],cur[maxn]; void insert(int u,int v,int flow) {tot++;e[tot].v=v;e[tot].flow=flow;e[tot].from=first[u];first[u]=tot;} bool bfs() { memset(d,-1,sizeof(d)); int head=0,tail=1;q[head]=S; d[S]=0; while(head!=tail) { int x=q[head++];if(head>=501)head=0; for(int i=first[x];i;i=e[i].from) if(d[e[i].v]==-1&&e[i].flow>0) { q[tail++]=e[i].v;d[e[i].v]=d[x]+1; if(tail>=501)tail=0; } } if(d[T]==-1)return 0;else return 1; } int dfs(int x,int a) { if(x==T||a==0)return a; int flow=0,f; for(int& i=cur[x];i;i=e[i].from)//当前弧优化 if(d[e[i].v]==d[x]+1&&(f=dfs(e[i].v,min(a,e[i].flow)))>0) { e[i].flow-=f; e[i^1].flow+=f; flow+=f; a-=f; if(a==0)break; } return flow; } int main() { scanf("%d%d",&m,&n); S=1;T=n; for(int i=1;i<=m;i++) { int u,v,w; scanf("%d%d%d",&u,&v,&w); insert(u,v,w); insert(v,u,0); } int ans=0; while(bfs()) { for(int i=1;i<=n;i++)cur[i]=first[i]; ans+=dfs(S,inf); } printf("%d",ans); return 0; }
标签:log out sed cst csdn 记录 sdn print size
原文地址:http://www.cnblogs.com/onioncyc/p/6421652.html