网络流

最大流
1.FF方法
2.EK算法
时间复杂度：O(n*m^2)
实际运用中远远达不到这个上界，一般能处理10^3 - 10^4规模的网络

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;
const int inf=1<<29,N=2010,M=20010;
int head[N],ver[M],edge[M],ne[M],v[N],incf[N],pre[N];
int n,m,s,t,tot,maxflow;
void add(int x,int y,int z)
{
	ver[++tot]=y;edge[tot]=z;ne[tot]=head[x];head[x]=tot;
	ver[++tot]=x;edge[tot]=0;ne[tot]=head[y];head[y]=tot;
}
bool bfs()
{
	memset(v,0,sizeof v);
	queue<int> q;
	q.push(s);
	v[s]=1;
	incf[s]=inf;
	while(q.size())
	{
		int x=q.front();q.pop();
		for(int i=head[x];i;i=ne[i])
		{
			if(edge[i])
			{
				int y=ver[i];
				if(v[y]) continue;
				incf[y]=min(incf[x],edge[i]);
				pre[y]=i;
				q.push(y);v[y]=1;
				if(y==t) return 1;
			}
		}
		
	}
	return 0; 
}
void update()
{
	int x=t;
	while(x!=s)
	{
		int i=pre[x];
		edge[i]-=incf[t];
		edge[i^1]+=incf[t];
		x=ver[i^1];
	}
	maxflow+=incf[t];
}
int main()
{
	cin>>n>>m>>s>>t;
	memset(head,0,sizeof head);
	tot=1;maxflow=0;
	for(int i=1;i<=m;i++)
	{
		int x,y,c;
		scanf("%d%d%d",&x,&y,&c);
		add(x,y,c);
	}
	
	while(bfs()) update();
	cout<<maxflow<<endl;
}

3.Sap算法
4.iSap算法
5.Dinic算法
时间复杂度：O(n^2m)
实际运用中远远达不到这个上界，可以说是比较容易实现的效率最高的网络流算法之一，一般能处理10^4 ~ 10^5规模的网络。
Dinic算法求解二分图最大匹配的时间复杂度为O(msqrt(n)),实际表现则更快

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;
const int N=1e6+100;
const int M=N*2;
const int inf=1e8;
int n,m,s,t;
int h[N],edge[M],ver[M],ne[M],tot;
int cur[N],d[N];
void add(int a,int b,int c)
{
	ver[tot]=b;edge[tot]=c;ne[tot]=h[a];h[a]=tot++;
	ver[tot]=a;edge[tot]=0;ne[tot]=h[b];h[b]=tot++;
}
int dfs()
{
	memset(d,-1,sizeof d);
	queue<int> q;
	q.push(s);d[s]=0;cur[s]=h[s];
	while(q.size())
	{
		int x=q.front();q.pop();
		for(int i=h[x];i!=-1;i=ne[i])
		{
			int y=ver[i];
			if(d[y]==-1&&edge[i])
			{
				d[y]=d[x]+1;
				cur[y]=h[y];
				q.push(y);
				if(y==t) return 1;
			}
		}
	}
	return 0;

}
int find(int x,int limit)
{
	if(x==t) return limit;
	int flow=0;
	for(int i=cur[x];i!=-1&&flow<limit;i=ne[i])
	{
		cur[x]=i;//当前弧优化
		int y=ver[i];
		if(d[y]==d[x]+1&&edge[i])
		{
			int k=find(y,min(limit-flow,edge[i]));
			if(!k) cur[y]=-1;
			edge[i]-=k;
			edge[i^1]+=k;
			flow+=k;
		}
	}
	return flow;
}
int dinic()
{
	int flow;
	int maxflow=0;
	while(dfs())
	{
		while(flow=find(s,inf))
		{
			maxflow+=flow;
		}
	}
	return maxflow;
}
int main()
{
	memset(h,-1,sizeof h);
	scanf("%d%d%d%d",&n,&m,&s,&t);
	for(int i=1;i<=m;i++)
	{
		int a,b,c;
		scanf("%d%d%d",&a,&b,&c);
		add(a,b,c);
	}
	printf("%d\n",dinic());	
} 

费用流

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=5010,M=100010,inf=1e8;
int n,m,S,T;
int h[N],e[M],f[M],w[M],ne[M],idx;
int q[N],d[N],pre[N],incf[N];
bool st[N];
void add(int a,int b,int c,int d)
{
	e[idx]=b;f[idx]=c;w[idx]=d;ne[idx]=h[a];h[a]=idx++;
	e[idx]=a;f[idx]=0;w[idx]=-d;ne[idx]=h[b];h[b]=idx++;
}
bool spfa()
{
	int hh=0,tt=1;
	memset(d,0x3f,sizeof d);
	memset(incf,0,sizeof incf);
	q[0]=S;d[S]=0,incf[S]=inf;
	while(hh!=tt)
	{
		int t=q[hh++];
		if(hh==N) hh=0;
		st[t]=0;
		for(int i=h[t];i!=-1;i=ne[i])
		{
			int ver=e[i];
			if(f[i]&&d[ver]>d[t]+w[i])
			{
				d[ver]=d[t]+w[i];
				pre[ver]=i;
				incf[ver]=min(f[i],incf[t]);
				if(!st[ver])
				{
					q[tt++]=ver;
					if(tt==N) tt=0;
					st[ver]=1;
				}
			}
		}
	}
	return incf[T]>0;
}
void EK(int &flow,int &cost)
{
	flow=cost=0;
	while(spfa())
	{
		int t=incf[T];
		flow+=t,cost+=t*d[T];
		for(int i=T;i!=S;i=e[pre[i]^1])
		{
			f[pre[i]]-=t;
			f[pre[i]^1]+=t;
		}
	}
}
int main()
{
	scanf("%d%d%d%d",&n,&m,&S,&T);
	memset(h,-1,sizeof h);
	while(m--)
	{
		int a,b,c,d;
		scanf("%d%d%d%d",&a,&b,&c,&d);
		add(a,b,c,d);
	}
	int flow,cost;
	EK(flow,cost);
	printf("%d %d\n",flow,cost);
	return 0;
}