网络流最大流EK算法板子
生活随笔
收集整理的這篇文章主要介紹了
网络流最大流EK算法板子
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
剛剛看了個板子,自己能夠理解的也不是太深,但是大體上懂什么意思了。
#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 flow[maxn][maxn],c[maxn][maxn];//flow表示當前路已經裝的流量,c表示容量
int res[maxn],pre[maxn];//res表示殘留量,pre記錄前驅點
int m,n;
int EK(int st,int ed)
{int max_flow=0;memset(flow,0,sizeof(flow));//一開始所有流量都為0while(1){memset(res,0,sizeof(res));//res賦值為0起到一個標記的作用,當res有值時表示殘留,即可以用的量res[st]=inff;//源點的可用的量是無盡的queue<int> q;q.push(st);while(!q.empty()){int u=q.front();q.pop();for(int v=1;v<=n;v++)//遍歷所有的邊,找到可行的增廣路{if(!res[v]&&flow[u][v]<c[u][v])//只要這個點沒有被訪問過且流量小于容量{pre[v]=u;//前一個點記錄res[v]=min(res[u],c[u][v]-flow[u][v]);//記錄殘留量,盡量裝的最多,q.push(v);}}}if(res[ed]==0) //如果沒有殘留量,即為最大流return max_flow;//更新正邊和反邊流量,因為某條反邊也可能是另一條路的反邊for(int i=ed;i!=st;i=pre[i])//從匯點往會跑更新每條路的流,{flow[pre[i]][i]+=res[ed];//更新正邊流量flow[i][pre[i]]-=res[ed];//更新反邊流量}max_flow+=res[ed];//每次加上增廣的那條邊的值}
}
int main()
{int x,y,z;while(cin>>m>>n){memset(c,0,sizeof(c));memset(pre,0,sizeof(pre));while(m--){scanf("%d %d %d",&x,&y,&z);c[x][y]+=z;//重邊就當一條}printf("%d\n",EK(1,n));}return 0;
}
更新:另一個版本不記錄flow[][]流量,而是直接將c[][]容量進行改變,因為前邊跑的增廣路已經是確定的了,所以直接用c[][]容量判斷條件更加簡單了。
#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 c[maxn][maxn];
int pre[maxn],res[maxn];
int n,m;
int EK(int st,int ed)
{memset(pre,0,sizeof(pre));int max_flow=0;while(1){queue<int> q;q.push(1);memset(res,0,sizeof(res));res[1]=inff;while(!q.empty()){int u=q.front();q.pop();for(int v=1;v<=n;v++){if(!res[v]&&c[u][v])//判斷條件直接改成c{pre[v]=u;res[v]=min(c[u][v],res[u]);q.push(v);}}}if(res[ed]==0) return max_flow;for(int u=ed;u!=st;u=pre[u]){c[pre[u]][u]-=res[ed];//同時更新的也是容量c[u][pre[u]]+=res[ed];}max_flow+=res[ed];}
}
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;}printf("%d\n",EK(1,n));}
}
鄰接表版本:
#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 res[maxn];
int head[maxn],sign;
int n,m;
struct node
{int to,p,w;
}edge[maxn<<1];
struct nod
{int v,id;
}pre[maxn];//pre數組不僅要找前驅點,還要標記鄰接表對應下標
void init()
{sign=0;for(int i=0;i<=n;i++)head[i]=-1;
}
void add(int u,int v,int val)
{edge[sign]=node{v,head[u],val};head[u]=sign++;
}
int EK(int st,int ed)
{int max_flow=0;memset(pre,0,sizeof(pre));while(1){queue<int> q;q.push(st);memset(res,0,sizeof(res));res[st]=inff;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(!res[v]&&edge[i].w){pre[v].v=u;//標記前驅pre[v].id=i;//標記鄰接表下標res[v]=min(edge[i].w,res[u]);q.push(v);}}}if(res[ed]==0) return max_flow;max_flow+=res[ed];for(int u=ed;u!=st;u=pre[u].v){edge[pre[u].id].w-=res[ed];//因為是連著兩個加入,剛好是^1就是對應的正反邊edge[pre[u].id^1].w+=res[ed];}}
}
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);}printf("%d\n",EK(1,n));}return 0;
}
?
總結
以上是生活随笔為你收集整理的网络流最大流EK算法板子的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: HDU 3306 Another kin
- 下一篇: 网络流最大流 Dinic算法模板