hdu 1532 网络流-最大流 DINIC

最新推荐文章于 2020-07-25 00:50:30 发布

原创最新推荐文章于 2020-07-25 00:50:30 发布 · 557 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#hdu 1532 #网络流-最大流 #DINIC

hdu 同时被 2 个专栏收录

40 篇文章

订阅专栏

网络流

2 篇文章

订阅专栏

本文介绍如何使用DINIC模板解决网络流问题，通过实例代码演示了从理解概念到实现过程的逐步学习路径。从网络流的基础概念出发，逐步深入探讨DINIC算法的原理和应用，旨在帮助初学者掌握网络流问题的求解方法。

//	hdu 1532 网络流-最大流 DINIC
//
//	网络流,刚开始接触感觉很好玩,但一直不得要领,现在呢
//	还是不得要领,学习的过程就是慢慢熟悉吧,学习着刘老师
//	的DINIC模板,好好加油

#include <cstdio>
#include <cmath>
#include <map>
#include <queue>
#include <iostream>
#include <vector>
#include <cstring>
#include <algorithm>
#define cls(x,a) memset(x,a,sizeof(x))
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int MAX_N = 2000 + 8;
const int MOD = 1e9 + 7;
const double PI = acos(-1);
const int INF = 1e9 + 7;
int n,m;


struct Edge{
	int from;
	int to;
	int cap;
	int flow;
	Edge(){

	}
	Edge(int u,int v,int c,int f):from(u),to(v),
		cap(c),flow(f){

		}
};

struct Dinic{
	int n,m,s,t;

	vector<Edge> edge;
	vector<int> G[MAX_N];
	bool vis[MAX_N];
	int d[MAX_N];
	int cur[MAX_N];

	void init(){
		edge.clear();
		for (int i = 0;i < MAX_N;i ++)
			G[i].clear();
		memset(vis,0,sizeof(vis));
		memset(cur,0,sizeof(cur));

	}

	void add_edge(int u,int v,int c){
		edge.push_back(Edge(u,v,c,0));
		edge.push_back(Edge(v,u,0,0));
		int m = edge.size();
		G[u].push_back(m-2);
		G[v].push_back(m-1);
	}
	bool BFS(){
		memset(vis,0,sizeof(vis));
		queue<int> Q;
		d[s] = 0;
		vis[s] = 1;
		Q.push(s);
		int u;
		while(!Q.empty()){
			u = Q.front();
			Q.pop();
			for (int i = 0;i < G[u].size();i ++){
				Edge& e = edge[G[u][i]];

				if (!vis[e.to] && e.cap > e.flow){
				 	vis[e.to] = 1;
					d[e.to] = d[u] + 1;
					Q.push(e.to);
				}
			}
		}
		return vis[t];
	}

	int DFS(int x,int a){
		if (x == t || a == 0)
			return a;
		int flow = 0,f;
		for (int& i = cur[x];i < G[x].size();i ++){
			Edge& e = edge[G[x][i]];

			if (d[x] + 1 == d[e.to] && (f = DFS(e.to,min(a,e.cap - e.flow))) > 0){
				e.flow += f;
				edge[G[x][i]^1].flow -= f;
				flow += f;
				a -= f;
				if (a == 0)
					break;
			}
		}
		return flow;
	}

	int Maxflow(int s,int t){
		this->s = s;
		this->t = t;
		int flow = 0;
		while(BFS()){
			memset(cur,0,sizeof(cur));
			flow += DFS(s,INF);
		}
		return flow;
	}

}di;

void input(){
	int u,v,c;
	di.init();
	for (int i = 1;i <= m;i ++){
		scanf("%d%d%d",&u,&v,&c);
		di.add_edge(u,v,c);
	}
}

void solve(){
	printf("%d\n",di.Maxflow(1,n));
}

int main(){
//	freopen("te.txt","r",stdin);
	while(scanf("%d%d",&m,&n)!=EOF){
		input();
		solve();	
	}
	return 0;
}