网络流最大流 Dinic算法模板
生活随笔
收集整理的這篇文章主要介紹了
网络流最大流 Dinic算法模板
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
Dinic算法跟EK算法的優化地方就在于,EK是每次bfs記錄前驅結點,然后從匯點找到一條回到源點的路。
而dinic是bfs記錄深度,不斷從源點dfs到匯點并且更新路徑流量,直到dfs不到匯點。
鄰接矩陣代碼:
#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<set>
#include<stack>
#include<vector>
#include<map>
#include<queue>
#define myself i,l,r
#define lson i<<1
#define rson i<<1|1
#define Lson i<<1,l,mid
#define Rson i<<1|1,mid+1,r
#define half (l+r)/2
#define inff 0x3f3f3f3f
#define lowbit(x) x&(-x)
#define PI 3.14159265358979323846
#define min4(a,b,c,d) min(min(a,b),min(c,d))
#define min3(x,y,z) min(min(x,y),min(y,z))
#define pii make_pair
#define pr pair<int,int>
const int dir[4][2]={0,1,0,-1,1,0,-1,0};
typedef long long ll;
const ll inFF=9223372036854775807;
typedef unsigned long long ull;
using namespace std;
const int maxn=405;
int dis[maxn],c[maxn][maxn];
int n,m;
int bfs(int st,int ed)
{memset(dis,-1,sizeof(dis));//深度dis[st]=0;queue<int>q;q.push(1);while(!q.empty()){int u=q.front();q.pop();for(int v=1;v<=n;v++){if(dis[v]<0&&c[u][v]>0)//沒有遍歷過的點{dis[v]=dis[u]+1;q.push(v);}}}if(dis[ed]>0) return true;//如果都沒有跑的匯點,意思沒有增廣路return false;
}
int find(int u,int low)//其實是一個dfs的過程,上邊bfs記錄深度也是因為這里通過dfs跑出多次增廣路
{int a=0;if(u==n) return low;for(int i=1;i<=n;i++){if(c[u][i]>0&&dis[i]==dis[u]+1&&(a=find(i,min(low,c[u][i]))))//3個條件//1.這條路有殘留量2.這兩點相連,通俗的講就是這次找增廣路的過程中,水就是沿著這條路的,但是下次就不一定.3就是遞歸下一個點看能不能到匯點//同時在過程中找到這條路最小流量,這種方法在稠密的圖中有很大的優化作用{c[u][i]-=a;//更新殘留量c[i][u]+=a;return a;}}return 0;
}
int main()
{int x,y,z;while(cin>>m>>n){memset(c,0,sizeof(c));for(int i=0;i<m;i++)scanf("%d %d %d",&x,&y,&z),c[x][y]+=z;int ans=0,temp;while(bfs(1,n)){while(temp=find(1,inff))ans+=temp;}printf("%d\n",ans);}return 0;
}
鄰接表板子:
#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<set>
#include<stack>
#include<vector>
#include<map>
#include<queue>
#define myself i,l,r
#define lson i<<1
#define rson i<<1|1
#define Lson i<<1,l,mid
#define Rson i<<1|1,mid+1,r
#define half (l+r)/2
#define inff 0x3f3f3f3f
#define lowbit(x) x&(-x)
#define PI 3.14159265358979323846
#define min4(a,b,c,d) min(min(a,b),min(c,d))
#define min3(x,y,z) min(min(x,y),min(y,z))
#define pii make_pair
#define pr pair<int,int>
const int dir[4][2]={0,1,0,-1,1,0,-1,0};
typedef long long ll;
const ll inFF=9223372036854775807;
typedef unsigned long long ull;
using namespace std;
const int maxn=205;
int d[maxn];
int head[maxn],sign;
struct node
{int to,p,val;
}edge[maxn<<1];
int n,m;
void add(int u,int v,int val)
{edge[sign]=node{v,head[u],val};head[u]=sign++;
}
void init()
{sign=0;for(int i=0;i<=n;i++)head[i]=-1;
}
int find(int u,int low)
{int a=0;if(u==n) return low;for(int i=head[u];~i;i=edge[i].p){int v=edge[i].to;if(d[v]==d[u]+1&&edge[i].val>0&&(a=find(v,min(low,edge[i].val)))){edge[i].val-=a;edge[i^1].val+=a;return a;}}return 0;
}
int bfs(int st,int ed)
{memset(d,0,sizeof(d));d[st]=1;queue<int> q;q.push(st);while(!q.empty()){int u=q.front();q.pop();for(int i=head[u];~i;i=edge[i].p){int v=edge[i].to;if(!d[v]&&edge[i].val>0){d[v]=d[u]+1;q.push(v);}}}if(d[ed]>0) return true;return false;
}
int main()
{int x,y,z;while(cin>>m>>n){init();for(int i=0;i<m;i++){scanf("%d %d %d",&x,&y,&z);add(x,y,z),add(y,x,0);}int ans=0,temp;while(bfs(1,n)){while(temp=find(1,inff))ans+=temp;}printf("%d\n",ans);}
}
?
總結
以上是生活随笔為你收集整理的网络流最大流 Dinic算法模板的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 网络流最大流EK算法板子
- 下一篇: 网络流Dinic cur当前弧优化