网络流总结

基本模型：
1.最大流
2.最小割
3.费用流
4.上下界网络流
无源汇的上下界可行流
无源汇的最小费用上下界可行流
有源汇的上下界最大流：t向s连一条inf，求一个可行流，ans+=t向s这条边的实际流量，然后删去这条边，再跑一个s到t的最大流，ans+=这个最大流。
有源汇的上下界最小流：t向s连一条inf，求一个可行流，ans+=t向s这条边的实际流量，然后删去这条边，再跑一个t到s的最大流，ans-=这个最大流。
5.无向图的点连通度：求最少去掉多少个点能使无向图不连通
做法：拆点求最小割
6.最小割的可行边和必须边
两者的大前提:u v满流
可行边:tarjan完了之后不在同一个强连通分量里
必须边：dfs
7.最大权闭合子图:V中所有点的出边都指向V里的点 则V是一个闭合子图
（对二元关系的理解嘛 只算额外代价）
8.限制距离模型
9.混合图的欧拉回路

模板：
最大流

int bfs()
{
    memset(dist,-1,sizeof(dist));
    dist[S]=0;
    q.push(S);
    while(!q.empty())
    {
        int x=q.front();
        q.pop();
         for(int i=head[x];i!=-1;i=next[i])
        {
           int y=key[i];
           if(w[i]&&dist[y]==-1)
           dist[y]=dist[x]+1;
        }
    }
    if(dist[T]==-1)return 0;
    else return 1;
}
int zeng(int x,int b)
{
    if(x==T)return b;
    int now=0;
    for(int i=head[x];~i&&now<b;i=next[i])
    {
        int y=key[i];
        int tmp=zeng(y,min(b-now,w[i]));
        now+=tmp;
        w[i]-=tmp;
        w[i^1]+=tmp;
    }
    if(now==0)dist[x]=-1;
    return now;
}

当前弧优化 在边多的情况下优化速度真的很有效

int bfs()
{
    memset(dist,-1,sizeof(dist));
    dist[S]=0;
    q.push(S);
    while(!q.empty())
    {
        int x=q.front();
        q.pop();
         for(int i=head[x];i!=-1;i=next[i])
        {
           int y=key[i];
           if(w[i]&&dist[y]==-1)
           dist[y]=dist[x]+1;
        }
    }
    if(dist[T]==-1)return 0;
    else return 1;
}
int zeng(int x,int b)
{
    if(x==T)return b;
    int now=0;
    for(int i=cur[x];~i&&now<b;i=next[i])
    {
        int y=key[i];
        int tmp=zeng(y,min(b-now,w[i]));
        now+=tmp;
        w[i]-=tmp;
        w[i^1]+=tmp;
        if(w[i])cur[x]=i;
    }
    if(now==0)dist[x]=-1;
    return now;
}
while(bfs())
{
    for(int i=S;i<=T;i++)
    cur[i]=head[i];
    ans+=zneg(S,inf);
}

费用流

int head[N],next[M],key[M],w[M],tot,cost[M];
long long dist[N],inq[N],mn[N],fr[N];
queue<int>q;
int S,T;

void add(int x,int y,int z,int v)
{
tot++;next[tot]=head[x];head[x]=tot;key[tot]=y;w[tot]=z;cost[tot]=v;
tot++;next[tot]=head[y];head[y]=tot;key[tot]=x;w[tot]=0;cost[tot]=-v;
}
int spfa()
{
   memset(mn,0x3f,sizeof(mn));
   memset(dist,0x3f,sizeof(dist));
   dist[S]=0;q.push(S);
   while(!q.empty())
   {
     int x=q.front();inq[x]=0;q.pop();
     for(int i=head[x];~i;i=next[i])
     {
        int y=key[i];
        if(dist[y]>dist[x]+cost[i]&&w[i])
        {
            dist[y]=dist[x]+cost[i];
            fr[y]=i;
            mn[y]=min(mn[x],w[i]);
            if(inq[y]==0)inq[y]=1,q.push(y);
        }
     }
   }
   if(dist[T]==inf)return 0;
   else return 1;
}
long long zeng()
{
   for(int i=T;i!=S;i=key[fr[i]^1])
   w[fr[i]]-=mn[T],w[fr[i]^1]+=mn[T];
   return (long long)dist[T]*mn[T];
}