最大流模板

#include<bits/stdc++.h>
#define mem(a,b) memset(a,b,sizeof(a))
#define For(a,b,c) for(int a = b;a <= c;a++)
using namespace std;
typedef long long ll;
const int maxn = 1000005;
const int INF = 0x3f3f3f3f;
const int e_maxn = 1000 * 4;
const int v_maxn = 500;
struct ppp
{
	int v,nex,cap,flow,c;//分别为下一个点，下一条边，当前边容量和流量，还有费用
}e[e_maxn];
int head[v_maxn],dis[v_maxn],cur[v_maxn];//链表头节点，层次图中每个点所属层次，dfs所用的临时链表头
int tole,N,M,s,t;//总加入边数，总输入点数和总输入边数，源点和汇点
void make_edge(int u,int v,int cap,int c)//求最大流时费用不必考虑，可以省去，加上也没关系
{
	e[tole].v = v;e[tole].flow = 0;e[tole].cap = cap;e[tole].c = c;e[tole].nex = head[u];head[u] = tole++;
}
void add_edge(int u,int v,int cap,int c)
{
	make_edge(u,v,cap,c);//正向边容量为cap，费用为c
	make_edge(v,u,0,-c);//反向边容量为0，费用为-c
}
int bfs()//其中s为源点，t为汇点
{
	queue<int> que;
	que.push(s);
	mem(dis,-1);
	dis[s] = 0;
	int temp,v;
	while(!que.empty())
	{
		temp = que.front();
		que.pop();
		for(int i = head[temp];~i;i = e[i].nex)
		{
			v = e[i].v;
			if(dis[v] == -1 && e[i].cap > e[i].flow)
			{
				dis[v] = dis[temp] + 1;
				que.push(v);
			}
		}
	}
	return dis[t] != -1;
}

int dfs(int x,int a)
{
	if(x == t || !a) return a;
	int v,f,ret = 0;
	for(int &i = cur[x];~i;i = e[i].nex)
	{
		v = e[i].v;
		if((dis[v] == dis[x] + 1) && (f = dfs(v,min(a,e[i].cap - e[i].flow))) > 0)
		{
			e[i].flow += f;
			e[i ^ 1].flow -= f;
			a -= f;
			ret += f;
			if(!a)break;
		}
	}
	return ret;
}

void Dinic(int &ans)//ans为最大流
{
	while(bfs())
	{
		for(int i = 0;i <= t;i++)cur[i] = head[i];
		ans += dfs(s,INF);
	}
}