POJ-1273-最大流裸題。
阿新 • • 發佈:2018-12-30
用EK演算法水過
#include <cstdio> #include <cmath> #include <cstring> #include <string> #include <algorithm> #include <iostream> #include <queue> #include <set> #include <vector> #include<stack> using namespace std; #include <fstream> const int N=205; const int inf=2147483647; int min(int a,int b) {return a<b?a:b;} int n,m; int max_flow;//max_flow是最大流 int map[N][N],flow[N][N]; // map[i][j]是每條邊的容量,flow[i][j]是每條邊的流量 int res[N],pre[N];//res[]是每個點的剩餘流量,pre[]是每個點的父親 int EK(int st,int ed) { max_flow=0; queue<int> q; while(!q.empty()) q.pop(); memset(flow,0,sizeof(flow));//最開始每條邊的流量都是0 while(1) { memset(res,0,sizeof(res));//殘餘流量得變0,一開始所有點都沒流入 res[st]=inf; //源點,剩餘流量無限 q.push(st); //從源點開始進行BFS找增廣路 int u,v; while(!q.empty()) { u=q.front(); q.pop(); for (v=1;v<=m;v++) //遍歷所有點,找可行邊 { if (!res[v]&& map[u][v]>flow[u][v]) //!res[v]表示沒訪問過,且 容量大於流量 { pre[v]=u; q.push(v); res[v]=min(res[u],map[u][v]-flow[u][v]);儘量填滿 } } } if (res[ed]==0) //如果當前已經是最大流,匯點沒有殘餘流量 return max_flow; for (u=ed;u!=st;u=pre[u])//如果還能增廣,那麼回溯,從匯點往回更新每條走過的邊的流量 { flow[pre[u]][u]+=res[ed];//更新正向流量 flow[u][pre[u]]-=res[ed]; //更新反向流量 } max_flow+=res[ed]; //更新最大流 } } int main() { while(scanf("%d%d",&n,&m)!=EOF) { memset(pre,0,sizeof(pre)); memset(map,0,sizeof(map)); int u,v,w; for (int i=1;i<=n;i++) { scanf("%d%d%d",&u,&v,&w); map[u][v]+=w; //重邊看作一條邊 } printf("%d\n",EK(1,m)); } return 0; }