最大流模板

<2015年1月30日更新>

#include <cstdio>
#include <fstream>
#include <iostream>

#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>

#include <queue>
#include <vector>
#include <map>
#include <set>
#include <stack>

typedef unsigned int uint;
typedef long long int ll;
typedef unsigned long long int ull;
typedef double db;

#define DBG printf("*")

using namespace std;

const int INF=(1<<28)-1;
const ll LINF=((ll)1<<52)-1;

struct edge
{
	int in;
	ll c;
	edge*nxt;
	edge*ptr;
}pool[100000];
edge*et=pool;
edge*eds[200];
void addedge(int i,int j,ll c)
{
	et->ptr=et+1;
	et->c=c; et->in=j; et->nxt=eds[i]; eds[i]=et++;
	et->ptr=et-1;
	et->c=0; et->in=i; et->nxt=eds[j]; eds[j]=et++; 
}
#define FOREACH_EDGE(i,j) for(edge*i=eds[j];i!=NULL;i=i->nxt)

int n;

int st,ed;
int dep[200];
ll DFS(int x,ll mi)
{
	if(x==ed) return mi;
	ll res=0,c;
	FOREACH_EDGE(i,x)
	if(i->c>0 && dep[x]+1==dep[i->in] && (c=DFS(i->in,min(mi,i->c))) )
	{
		res+=c;
		i->c-=c;
		i->ptr->c+=c;
		mi-=c;
		if(mi<=0) break;
	}
	if(res<=0) dep[x]=-1;
	return res;
}

int qh,qt;
int q[400];
ll DINIC()
{
	ll res=0;
	
	while(true)
	{
		memset(dep,0xFF,sizeof(int)*(n+1));
		
		qh=qt=0;
		q[qt++]=st;
		dep[st]=0;
		while(qh!=qt)
		{
			int&cur=q[qh];
			FOREACH_EDGE(i,cur)
			if(i->c>0 && dep[i->in]==-1)
			{
				dep[i->in]=dep[cur]+1;
				q[qt++]=i->in;
			}
			qh++;
		}
		
		if(dep[ed]==-1) break;
	
		res=max((ll)0,res);
		res+=DFS(st,LINF);
	}
	
	return res;
}


//blocks define
#define ST(i) (i)
#define ED(i) ((i)+ntot)

int ntot,m;

int main()
{
	int T;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d%d",&ntot,&m);
		et=pool;
		
		ntot++;
		n=ntot*2;
		st=ED(0);
		ed=ST(ntot-1);
		
		memset(eds,0,sizeof(int)*(n+1));
		
		for(int i=1;i<ntot-1;i++)
		{
			ll c;
			scanf("%I64d",&c);
			addedge(ST(i),ED(i),c);
		}
		
		addedge(ST(0),ED(0),LINF);
		addedge(ST(ntot-1),ED(ntot-1),LINF);
		
		for(int i=0;i<m;i++)
		{
			int a,b;
			scanf("%d%d",&a,&b);
			addedge(ED(b),ST(a),LINF);
			addedge(ED(a),ST(b),LINF);
		}

		ll res=DINIC();
		
		if(res==0) printf("Min!\n");
		else if(res>=LINF) printf("Max!\n");
		else printf("%I64d\n",res);
	}
	
	return 0;
}

AC VIJOS P1590

把整个程序copy过来了....

现在打一个DINIC已经没有什么问题了....

构图方面还是比较坑...

然后是类似于DINIC的最小费用流.多路增广.存图同上,这里省略了.

int n;

int st,ed;
int dist[205];
int cost_add;
bool used[205];
int DFS(int x,int mi)
{
	if(x==ed)return mi;
	
	used[x]=true;
	int res=0;
	int c;
	FOREACH_EDGE(i,x)
	if(!used[i->in] && i->c>0 && dist[i->in]==dist[x]+i->v && (c=DFS(i->in,min(i->c,mi))))
	{
		i->c-=c;
		i->ptr->c+=c;
		mi-=c;
		res+=c;
		cost_add+=c*i->v;
		if(mi<=0) break;
	}
	used[x]=false;
	if(res<=0) dist[x]=INF;
	return res;
}

int q[1000000];
int qt,qh;
int DINIC()
{
	int res=0;
	
	while(true)
	{
		for(int i=0;i<n;i++)
		dist[i]=INF;
		
		qh=qt=0;
		dist[st]=0;
		q[qt++]=st;
		while(qh!=qt)
		{
			int&cur=q[qh];
			FOREACH_EDGE(i,cur)
			{
				int&nxt=i->in;
				if(i->c>0 && i->v + dist[cur] < dist[nxt])
				{
					dist[nxt]=dist[cur]+i->v;
					q[qt++]=nxt;
				}
			}
			qh++;
		}
		
		if(dist[ed]==INF) break;
		
		cost_add=0;
		DFS(st,INF);
		res+=cost_add;
	}
	
	return res;
}

注意DFS中使用了used在深搜的时候避免搜到环(最短路图中可能存在0环).

以后写最大流/最小费用最大流就果断背模板了...

反正最大流的题目重不在算法,重在建模.....

而且DINIC又好理解又好写(相对于SAP而言)