最大流dinic模板

以hdu1532为例

 

#include<cstring>  
#include<cstdio>  
#include<cmath>
#include<iostream>
#include<algorithm>
#include<queue>
#define N 405 //估算点
#define maxn 1<<26  
using namespace std;  

struct edge{
	int u,v,w,next;
}e[20*N];

int n,m;
int edn;  //边数 
int head[N];
int d[N];
int sp,tp; //源点、汇点

void add(int u,int v,int w){
	e[edn].u=u;e[edn].v=v;e[edn].w=w;
	e[edn].next=head[u];head[u]=edn++;
	e[edn].u=v;e[edn].v=u;e[edn].w=0;
	e[edn].next=head[v];head[v]=edn++;
} 
int bfs(){
	queue<int>q;
	int i,j;
	memset(d,-1,sizeof(d));
	d[sp]=0;
	q.push(sp);
	while(!q.empty()){
		int cur=q.front();
		q.pop();
		for(i=head[cur];~i;i=e[i].next){
			int v=e[i].v;
			if(d[v]==-1&&e[i].w>0){
				d[v]=d[cur]+1;
				q.push(v);
			}
		}
	}
	return d[tp]!=-1;//返回是否能够到达汇点 
}
int dfs(int a,int b){
	int r=0;
	int i;
	if(a==tp)return b;
	for(i=head[a];~i&&r<b;i=e[i].next){
		int v=e[i].v;
		if(e[i].w>0&&d[v]==d[a]+1){
			int x=min(e[i].w,b-r);
			x=dfs(v,x);
			r+=x;
			e[i].w-=x;
			e[i^1].w+=x;
		}
	}
	if(!r)d[a]=-2;
	return r;
}
int dinic(int sp,int tp){
	int total=0;
	int t;
	while(bfs()){
		while(t=dfs(sp,maxn)){
			total+=t;
		}
	}
	return total;
}
int main(){
	while(scanf("%d%d",&m,&n)!=EOF){
		int i;
		edn=0;
		memset(head,-1,sizeof(head));
		for(i=1;i<=m;i++){
			int u,v,w;
			scanf("%d%d%d",&u,&v,&w);
			add(u,v,w);
		}
		sp=1;	
		tp=n;
		printf("%d\n",dinic(sp,tp));
	} 
	return 0;
} 

 

 

 

*/

5 4 1 2 40 1 4 20 2 4 20 2 3 30 3 4 10

50*/

 

 

 

 

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值