网络流入门

网络流入门，网络一共有三种算法：分别是EK算法，dinic算法。但是一般用的是dinic算法 EK算法是用bfs暴力找增广路，找到增广路后回溯一下。dinic算法是先将图分个层，然后进行dfs，这样会快很多。

多源多汇问题

问题描述：源有多个，汇也有多个，流可以从任意一个源中流出，最终可以流向任意一个汇。总流量等于所有源流出的总流量，也等于所有汇流入的总流量。

建立一个超级源点s’和超级汇点t’。从s‘向每一个源点引一条有向弧，容量为无穷大；从每一个汇点向t’引一条有向弧，容量为无穷大。

结点容量

问题描述：每个结点都有一个允许通过的最大流量，称为结点容量。

把每个结点u分裂成两个结点u1和u2。从u1向u2引一条有向弧，容量为结点u的结点容量。对于原先指向u的有向边，改成指向u1；对于原先从u指出的有向边，改成从u2指出；容量都不变。

无源无汇有容量下界网络的可行流

问题描述：每条弧除了有容量上界c之外，还有一个容量下界b。 
原始的最大流问题的容量下界全等于0。

建立附加源s和汇点t，首先添加弧（t，s，INF），然后对弧进行改造，改造成下界为0的弧。对于上界为c、下界为b的有向弧(u,v,c),拆成三条有向弧（u，v，c-b）、（u，t，b）、（s，v，b），然后合并。 
求改造后的网络的s-t最大流即可。当且仅当所有附加弧满载时原网络有可行流。

           我觉得网络流精髓所在是每两个点之间设立一对权值为0，x的边。这样的作用是可以为之前的行为一个反悔的操作。如果有条边占用了d，那么边的权值为x-d d。如果要执行反向边k，那么就是说这条边应该退回K的大小。EK算法每次都要初始化pre，还有flow数组，flow的值是选取这条路之前和自己这条边的最小值。然后EK就是找到一条边就算了，找的每一次并不是最优解。然后走多次bfs知道没有增广路为止。

#include <iostream>
#include <queue>
#include<string.h>
#include <cstdio>
using namespace std;
#define arraysize 201
int maxData = 0x7fffffff;
int capacity[arraysize][arraysize]; //记录残留网络的容量
int flow[arraysize];                //标记从源点到当前节点实际还剩多少流量可用
int pre[arraysize];                 //标记在这条路径上当前节点的前驱,同时标记该节点是否在队列中
int n,m;
queue<int> myqueue;
int BFS(int src,int des)
{
    int i,j;
    while(!myqueue.empty())       //队列清空
        myqueue.pop();
    for(i=1;i<m+1;++i)
    {
        pre[i]=-1;
    }
    pre[src]=0;
    flow[src]= maxData;
    myqueue.push(src);
    while(!myqueue.empty())
    {
        int index = myqueue.front();
        myqueue.pop();
        if(index == des)            //找到了增广路径
            break;
        for(i=1;i<m+1;++i)
        {
            if(i!=src && capacity[index][i]>0 && pre[i]==-1)
            {
                 printf("%d %d\n",index,i);
                 pre[i] = index; //记录前驱
                 flow[i] = min(capacity[index][i],flow[index]);   //关键：迭代的找到增量
                 myqueue.push(i);
            }
        }
    }
    if(pre[des]==-1)      //残留图中不再存在增广路径
        return -1;
    else
        return flow[des];
}
int maxFlow(int src,int des)
{
    int increasement= 0;
    int sumflow = 0;
    while((increasement=BFS(src,des))!=-1)
    {
         int k = des;          //利用前驱寻找路径
         while(k!=src)
         {
              int last = pre[k];
              capacity[last][k] -= increasement; //改变正向边的容量
              capacity[k][last] += increasement; //改变反向边的容量
              k = last;
              printf("\n");
         }
         sumflow += increasement;
    }
    return sumflow;
}

然后讲的是dinic算法。dinic是使用bfs dfs交错使用。dfs是将这个图分个层次。那么dfs时候遍历的点就变少了，dinic中flow是指的是star点到现在该点增广路上最小值所剩余的值。然后int dfs(int u,int flow)，中的flow指的是增广路径上的最小值。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
int const INF =0x3f3f3f3f;
int const maxe=1e5+1000;
int const maxv =1650;
using namespace std;
int N,F,D;
int S,T;//源点和汇点
struct edge
{
    int to, w, next;
} e[maxe<<1];
int head[maxv<<1],depth[maxv],cnt;
void init ()
{
    memset(head,-1,sizeof(head));
    cnt=-1;
}
void add_edge(int u,int to,int w)
{
    e[++cnt].w=w;
    e[cnt].to=to;
    e[cnt].next=head[u];
    head[u]=cnt;
}
void _add(int u,int v,int w)
{
    add_edge(u,v,w);
    add_edge(v,u,0);
}
bool bfs()
{
    queue<int>q;
    while(!q.empty()) q.pop();
    memset(depth,0,sizeof(depth));
    depth[S]=1;
    q.push(S);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(int i=head[u]; i!=-1; i=e[i].next)
        {
            if(!depth[e[i].to]&&e[i].w>0)
            {
                depth[e[i].to]=depth[u]+1;
                q.push(e[i].to);
            }
        }
    }
    if(!depth[T]) return false;
    return true;
}
int dfs(int u,int flow)
{
    if(u==T) return flow;
    int ret =0;
    for(int i=head[u]; i!=-1&& flow; i=e[i].next)
    {
        if(depth[u]+1==depth[e[i].to]&& e[i].w!=0)
        {
            int temp=dfs(e[i].to,min(e[i].w,flow));
            if(temp>0)
            {
                flow=flow-temp;
                ret=ret+temp;
                e[i].w=e[i].w-temp;
                e[i^1].w=e[i^1].w+temp;
            }
        }
    }
    return ret;
}
int Dinic()
{
    int ans=0;
    while(bfs())
        ans=ans+dfs(S,INF);
    printf("%d\n",ans);
}
int main()
{
    init();
    S=0,T=600;
    scanf("%d%d%d",&N,&F,&D);
    for(int i=1; i<=F; i++) _add(S,i,1); //源点到食物
    for(int i=1; i<=D; i++) _add(100+i,T,1); //饮料到汇点
    for(int i=1; i<=N; i++) //i+200和i+300表示牛的编号
    {
        int f,d;
        scanf("%d%d",&f,&d);
        _add(i+200,i+300,1);//牛到牛
        while(f--)//食物
        {
            int u;
            scanf("%d",&u);
            _add(u,i+200,1);//食物到牛
        }
        while(d--)//饮料
        {
            int u;
            scanf("%d",&u);
            _add(i+300,u+100,1);
        }
    }
    Dinic();
    return 0;
}