#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<queue>
#define INF 0x3f3f3f3f
using namespace std;
int x[510][510],n,m;//存容量
struct zp//标号
{
int head;//上一个节点
int value;//最多传来的流量
} w[1010];
int maxriver()
{
int vis[1010],sum=0;
memset(vis,0,sizeof(vis));//标记该点是否找过
w[1].head=0;//1为原点n为汇点
w[1].value=INF;
queue<int> q;
q.push(1);
w[n].value=-1;
vis[1]=1;
int k=1;
while(!q.empty())
{
for(int i=1; i<=n; i++)//从源点开始广搜
if(!vis[i]&&x[k][i]>0)//通过有流量的边到可到达的点
{
w[i].head=k;//记录上个节点
w[i].value=min(x[k][i],w[k].value);//更新传过来的流量
q.push(i);
vis[i]=1;//标记已访问过
}
q.pop();
k=q.front();
if(w[n].value!=-1)//汇点有流量更新各边的流量
{
sum+=w[n].value;
int p=w[n].value;
k=n;
while(k!=1)//根据各节点回溯回去更新值
{
int flag=w[k].head;
x[flag][k]-=p;
x[k][flag]+=p;
k=flag;
}
memset(vis,0,sizeof(vis));
vis[1]=1;
while(!q.empty())
q.pop();
q.push(1);
w[n].value=-1;
k=1;
}
}
printf("%d\n",sum);//输出最大流
}
int main()
{
// int ncase;
//scanf("%d",&ncase);
// while(ncase--)
// {
while(~scanf("%d%d",&n,&m))
{
memset(x,0,sizeof(x));
int a,b,c;
for(int i=0; i<m; i++)//读入边及其容量
{
scanf("%d%d%d",&a,&b,&c);
x[a][b]+=c;//重边
}
maxriver();
}
//}
}
网络流EK算法
最新推荐文章于 2022-08-27 15:06:08 发布