洛谷 P3376 【模板】网络最大流 ISAP

题目链接

模板题。

ISAP：Dinic的优化，Dinic的每次bfs重新求分层图较为浪费时间，而ISAP就对此进行了优化，对于每次增广实际每个点的深度并没有太多改变，所以ISAP定义了每个点到汇点的最短距离，对于每个点的最短距离，每次增广之后是不会变小(若变小则一开始就存在更短的增广路)。

GAP优化：当与汇点距离为当前距离的结点数为0时，说明已经求得最大流，可以直接跳出dfs。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<string>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define PI acos(-1)
#define INF 0x3f3f3f3f
#define NUM 100100
#define debug true
#define lowbit(x) ((-x)&x)
#define ffor(i,d,u) for(int i=d;i<=u;++i)
#define _ffor(i,u,d) for(int i=u;i>=d;--i)
#define mst(array,Num) memset(array,Num,sizeof(array))
const int p = 1e9+7;
int n,m,maxflow=0,s,t;
int head[NUM/10],d[NUM/10]={},num[NUM/10]={},ednum=-1;
struct edge
{
    int next,to,w;
}e[NUM<<1];
void read(int &x){
    char ch = getchar();x = 0;
    for (; ch < '0' || ch > '9'; ch = getchar());
    for (; ch >='0' && ch <= '9'; ch = getchar()) x = x * 10 + ch - '0';
}
template <typename T>
void write(T x)
{
    int len=0;char c[21];
    if(x<0)putchar('-'),x*=(-1);
    do{++len;c[len]=(x%10)+'0';}while(x/=10);
    _ffor(i,len,1)putchar(c[i]);
}
inline void bfs()
{
    int x,y;
    ffor(i,1,n)d[i]=INF;
    queue < int > q;
    q.push(t),++num[d[t]=1];
    while(!q.empty())
    {
        x=q.front();q.pop();
        for(int i=head[x];i!=-1;i=e[i].next)
        {
            y=e[i].to;
            if(d[y]!=INF||e[i].w!=0)continue;
            ++num[d[y]=d[x]+1];
            q.push(y);
        }
    }
}
int dfs(int vertex,int minx)
{
    if(!minx||vertex==t)return minx;
    int f,x,flow=0;
    for(int i=head[vertex];i!=-1;i=e[i].next)
    {
        x=e[i].to;
        if(d[vertex]==d[x]+1&&(f=dfs(x,min(minx,e[i].w))))
        {
            minx-=f,flow+=f,e[i].w-=f,e[i^1].w+=f;
            if(!minx)return flow;
        }
    }
    if(!(--num[d[vertex]]))d[s]=n+1;
    ++num[++d[vertex]];
    return flow;
}
inline void ISAP()
{
    bfs();
    while(d[s]<=n)maxflow+=dfs(s,INF);
    write(maxflow);
}
inline void AC()
{
    read(n),read(m),read(s),read(t);
    int x,y,wei;
    mst(head,-1);
    ffor(i,1,m)
    {
        read(x),read(y),read(wei);
        e[++ednum].to=y,e[ednum].w=wei,e[ednum].next=head[x],head[x]=ednum;
        e[++ednum].to=x,e[ednum].w=0,e[ednum].next=head[y],head[y]=ednum;
    }
    ISAP();
}
int main()
{
    AC();
    return 0;
}