【模板】网络最大流

题目描述:点击进入

思路

最大流模板（小蓝书）（未完待续~~~）

代码

#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstring>
using namespace std;
typedef long long ll;
const int maxn=1e4+10;
struct node
{
    int to;
	int next;
	int val;
}edge[maxn];
int head[maxn],tot;
void add(int u,int v,int val)//建边 
{
    edge[tot].to=v;
	edge[tot].val=val;
    edge[tot].next=head[u];
	head[u]=tot++;
}
int deep[maxn];//记录深度 
int cur[maxn];//当前弧优化 
bool vis[maxn];
int n,m;
ll res;
int s;//　源点：只有流出去的点
int t;//　汇点：只有流进来的点
int dfs(int pos/*当前点*/,int flow)//dfs找增广路
{
    if(pos==t) return flow;
    int rest=flow,k,i;
    //i=cur[pos]当前弧优化 
    for(i=cur[pos];i!=-1&&rest;i=edge[i].next)
    {
    	int v=edge[i].to;
        int val=edge[i].val;
        if(deep[v]==deep[pos]+1&&val)
        {
            k=dfs(v,min(val,rest));//递归求最小值
            if(!k) deep[v]=0;//剪枝，去掉增广完毕的点 
                //回溯的时候修改反向边权值和正向边权值
            edge[i].val-=k;//修改反正向边权值
            edge[i^1].val+=k;//异或修改反向边
            rest-=k; 
        }
	}
	cur[pos]=i;
    return flow-rest;
}
bool bfs()//bfs处理深度 
{
	memset(deep,0,sizeof(deep));
    deep[s]=1;
    queue<int>q; 
    q.push(s);
    cur[s]=head[s];
    while(q.size())
    {
        int now=q.front();
		q.pop();
		for(int i=head[now];i!=-1;i=edge[i].next)
		{
			int v=edge[i].to;
            if(edge[i].val&&!deep[v])
            {
            	q.push(v);
            	cur[v]=head[v];
                deep[v]=deep[now]+1;
				if(v==t) return 1;
            }
		}
    }
    return 0;
}
void dinic()
{
	while(bfs())
        res+=dfs(s,0x3f3f3f3f);
}
int main()
{
	memset(head,-1,sizeof(head));
    scanf("%d%d%d%d",&n,&m,&s,&t);
    for(int i=1;i<=m;i++)//建图
    {
        int u,v,w;
        scanf("%d%d%d",&u,&v,&w);
        add(u,v,w);//正边 
		add(v,u,0);//反边 
		//注意一定要连续塞进去
    }
    dinic();
    printf("%lld\n",res);//输出答案
	return 0;
}