网络流

dinic
网络流

先利用bfs进行分层,然后再在上面跑dfs.

用的是dinic

具体还是看代码把

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
struct node
{
	int point;
	int value;
	int nxt;
};
node line[200100];
int head[50100],tail=-1;//利用的yi
int cur[50100],maxflow;
int dep[50100];
int inf=0x7fffffff;
int n,m,s,t;
queue<int>q;
void add(int x,int y,int dis,int base)
{
	line[++tail].point=y;
	if(base)	
		line[tail].value=dis;
	line[tail].nxt=head[x];
	head[x]=tail;
}
bool bfs(int start,int end)
{
	memset(dep,0,sizeof(dep));
	while(!q.empty())
		q.pop();
	for(int i=1;i<=n;i++)
		cur[i]=head[i];
	dep[start]=1;
	q.push(start);
	while(!q.empty())
	{
		int now=q.front();
		q.pop();
		for(int i=head[now];i!=-1;i=line[i].nxt)
		{
			if(!dep[line[i].point]&&line[i].value)
			{
				dep[line[i].point]=dep[now]+1;
				q.push(line[i].point);
			}
		}
	}
	if(dep[end])
		return true;
	else
		return false;
}
int dfs(int now,int aim,int limte)
{
	if(now==aim||!limte)
		return limte;
	int flow=0,f;
	for(int i=cur[now];i!=-1;i=line[i].nxt)
	{
		cur[now]=i;
		if(dep[line[i].point]==dep[now]+1&&(f=dfs(line[i].point,aim,min(limte,line[i].value))))
		{
			flow+=f;
			limte-=f;
			line[i].value-=f;
			line[i^1].value+=f;
			if(!limte)
				break;
		}
	}
	return flow;
	
}
int dinic(int s,int e)
{
	int ans=0;
	while(bfs(s,e))//利用bfs分层
		ans+=dfs(s,e,inf);//dfs寻找增广路
	return ans;
}
int main()
{
	scanf("%d%d%d%d",&n,&m,&s,&t);//读入数据
        int a,b,c;
	for(int i=1;i<=n;i++)
		head[i]=-1;//初始化链表
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&a,&b,&c);
		add(a,b,c,1);//最后一个参数是判断是否是正向边还是反向边
		add(b,a,c,0);
	}
	printf("%d",dinic(s,t));
}
posted @ 2018-04-27 15:57  Lance1ot  阅读(126)  评论(0编辑  收藏  举报