Dinic
//建立反向边是因为贪心不最优,给自己反悔的机会
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<queue> #include<vector> #include<math.h> using namespace std; #define M 9499999 int c[10005][10005]; int dep[10001]; int n,m,s,t; int dfs(int s,int t,int f) { if(s==t) return f; int d; for(int v=1;v<=n;v++) if(c[s][v]>0&&dep[v]==dep[s]+1&&(d=dfs(v,t,min(f,c[s][v])))) { c[s][v]-=d; c[v][s]+=d; return d; } return 0; } int bfs(int s,int t) { queue<int>q;q.push(s); memset(dep,-1,sizeof(dep)); dep[s]=0; while(!q.empty()){ int u=q.front();q.pop(); for(int v=1;v<=m;v++) if(c[u][v]>0&&dep[v]==-1) { dep[v]=dep[u]+1; q.push(v); } } return dep[t]!=-1; } int dinic() { int flow=0,f; while(bfs(s,t)) while(1){ f=dfs(s,t,M); if(!f) break; flow+=f; } return flow; } int main() { scanf("%d%d%d%d",&n,&m,&s,&t); for(int i=1,u,v,w;i<=m;i++) { scanf("%d%d%d",&u,&v,&w); c[u][v]+=w; } cout<<dinic(); return 0; }
链表的(较快)
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<queue> #include<vector> #include<math.h> using namespace std; #define M 94999999 int h[10002],nex[200002],na[200002],w[200002]; int dep[10001]; int n,m,s,t; int dfs(int s,int t,int f)//当前点,目标,可用流量 { if(s==t) return f; int d; for(int i=h[s];i;i=nex[i]) { int v=na[i]; if(w[i]>0&&dep[v]==dep[s]+1&&(d=dfs(v,t,min(f,w[i])))) { w[i]-=d; w[i^1]+=d; return d; } } } int bfs(int s,int t)//分层 { queue<int>q;q.push(s); memset(dep,-1,sizeof(dep)); dep[s]=0; while(!q.empty()){ int u=q.front();q.pop(); for(int i=h[u];i;i=nex[i]) { int v=na[i]; if(w[i]>0&&dep[v]==-1) { dep[v]=dep[u]+1; q.push(v); } } } return dep[t]!=-1; } int dinic() { int flow=0,f; while(bfs(s,t)) while(1){ f=dfs(s,t,M); if(!f) break; flow+=f; } return flow; } int main() { scanf("%d%d%d%d",&n,&m,&s,&t); int tot=1; for(int i=1,u,v,ww;i<=m;i++) { scanf("%d%d%d",&u,&v,&ww); na[++tot]=v;nex[tot]=h[u];h[u]=tot;w[tot]=ww; na[++tot]=u;nex[tot]=h[v];h[v]=tot; } cout<<dinic(); return 0; }
Dinic的核心是对搜索进行诱导(横向),降低复杂度。