网络流之最大流 模板

FF算法:有了Dinic，谁还用FF…
EK算法:有了Dinic，谁还用EK…
Dinic:

#include<bits/stdc++.h>
using namespace std;
int tot,st,ed,i,j,n,m,head[1111111],ans,dep[1111111],q[1111111];
struct eat{
	int wi,to,ne;
}e[1111111];
void add(int x,int y,int z)
{
	e[tot].wi=z;
	e[tot].to=y;
	e[tot].ne=head[x];
	head[x]=tot++;
}
bool bfs()
{
	memset(dep,0,sizeof dep);
	memset(q,0,sizeof q);
	int l=0,r=1;
	l++;
	q[l]=st;
	dep[st]=1;
	while(l<=r)
	{
		int u=q[l];
		l++;
		for(int i=head[u];i!=-1;i=e[i].ne)
		{
			int w=e[i].wi,v=e[i].to;
			if(w!=0&&dep[v]==0)
			{
				r++;
				q[r]=v;
				dep[v]=dep[u]+1;
				if(v==ed)
				{
					return true;
				}
			}
		}
	}
	return false;
}
int dfs(int u,int flow)
{
	if(u==ed) return flow;
	int res=flow,dd;
	for(int i=head[u];i!=-1;i=e[i].ne)
	{
		if(res==0) break;
		int v=e[i].to,w=e[i].wi;
		if(w!=0&&dep[v]==(dep[u]+1))
		{
			dd=dfs(v,min(res,w));
			if(dd==0) dep[v]=0;
			e[i].wi-=dd;
			e[i^1].wi+=dd;
			res-=dd;
		}
	}
	return flow-res;
}
void dinic()
{
	int f;
	while(bfs())
	{
		while(f=dfs(st,99999999))
		{
			ans+=f;
		}
	}
}
int main()
{
	memset(head,-1,sizeof head);
	memset(e,-1,sizeof e);
	scanf("%d%d%d%d",&n,&m,&st,&ed);
	for(i=1;i<=m;i++)
	{
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		add(x,y,z);
		add(y,x,0);
	}
	dinic();
	printf("%d\n",ans);
}

放道例题，以供复习。

【GDKOI2006】破解密文

Input

Output
输出文件包含一个整数，为你破译得到的最大价值

Sample Input
7
0 2 2 7
1 0
1 0
-1 3 1 6 7
1 2 1 6
1 2 2 7
-1 0

Sample Output
3

Data Constraint

这道题要求最大权闭合子图，可以用最小割求解，如果要用最大流的话，要用权值为正的总贡献和减去总流量即anflow即权值为负的总贡献流，这里要注意，连边要从被锁定的点连向锁定它的点，还有就是要将都点权设为正的，不过原来是正的要与超级源点相连，原来是负的要与超级汇点相连，还有就是要将原来题目给你的点的流量设为INF，直接跑一边Dinic即可，或者其他如ISAP等神奇诡异的高级算法均行。

代码如下：

#include<bits/stdc++.h>
using namespace std;
int inf=99999999;
int st,ed,n,i,j,a[1111111],sum,anflow,k;
int dep[1111111],q[1111111],l,r;
int tot,head[1111111];
struct ert{
	int ne,to,wi;
}e[1111111];
void add(int xx,int yy,int ww)
{
	e[tot].to=yy;
	e[tot].wi=ww;
	e[tot].ne=head[xx];
	head[xx]=tot++;
}
void init()
{
	tot=0;
	memset(e,-1,sizeof e);
	memset(head,-1,sizeof head);
}

bool bfs()
{
	memset(dep,0,sizeof dep);
	memset(q,0,sizeof q);
	l=0;
	r=0;
	r++;
	q[r]=st;
	dep[st]=1;
	while(l<r)
	{
		l++;
		int u=q[l],v,w;
		for(int i=head[u];i!=-1;i=e[i].ne)
		{
			v=e[i].to,w=e[i].wi;
			if(w!=0&&dep[v]==0)
			{
				dep[v]=dep[u]+1;
				r++;
				q[r]=v;
				if(v==ed) return true;
			}
		}
	}
	return false;
}
int dfs(int u,int flow)
{
//	cout<<u<<" "<<flow<<endl;
	if(u==ed) return flow;
	int res=flow,dd;
	for(int i=head[u];i!=-1;i=e[i].ne)
	{
		if(res==0) break;
		int v,w;
		v=e[i].to,w=e[i].wi;
		if(w!=0&&dep[v]==(dep[u]+1))
		{
			dd=dfs(v,min(w,res));
			if(dd==0) dep[v]=0;
			e[i].wi-=dd;
			e[i^1].wi+=dd;
			res-=dd;
		}
	}
	return (flow-res);
}
void dinic()
{
	int f;
	while(bfs()==1)
	{
		while(f=dfs(st,inf))
		{
//			cout<<f<<endl;
			anflow+=f;
		}
	}
}

int main()
{
	init();
	scanf("%d",&n);
	st=0;
	ed=n+1;
	for(i=1;i<=n;i++)
	{
		scanf("%d",&a[i]);
		scanf("%d",&k);
		for(j=1;j<=k;j++)
		{
			int asd;
			scanf("%d",&asd);
			add(i,asd,inf);
			add(asd,i,0);
		}
	}
	sum=0;
	for(i=1;i<=n;i++)
	{
		if(a[i]==1)
		{
			sum+=1;
			add(st,i,a[i]);
			add(i,st,0);
		}
		else if(a[i]==0)
		{
			add(st,i,0);
			add(i,st,0);
		}
		else
		{
//			cout<<a[i]<<endl;要设成正点权！！！ 
			add(i,ed,-a[i]);
			add(ed,i,0);
		}
	}
	dinic();
//	cout<<anflow<<endl;
	printf("%d\n",(sum-anflow));
}