1. 程式人生 > >POJ-1273-最大流裸題。

POJ-1273-最大流裸題。

用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;

}