ISAP复习

随手瞎写，只为自己印象深刻，所以写得不好请多多包涵。。
ISAP，也就是improved SAP，是SAP的升级版本，在应对一些特殊情况时的运行速度非常优秀，比SAP快，但是在平均统计下速度应该和SAP差不多，一般来说出题人也不会故意去卡SAP（我觉得这卡不了吧）。

我们知道无论是dinic，SAP还是ISAP都是建立在寻找增广路的基础上。不同之处只不过是在对于图的处理上。

前置结构：
d[i]:表示第i个点到汇点t的距离。
cur[i]:当前弧
remain[i]:剩余流量

一句话概括最大流，就是不断寻找最短增广路进行增广直到不能增广为止。
那么，对于ISAP来说，我们先bfs倒着预处理出d,然后，明显可以发现，如果d[i]=d[j]+1，走j->i这样的路径一定是最短路，我们把符合这样的边叫做允许弧。

那么明显，我们只要在残量网络（注意不是原图）中不断走允许弧，就可以达到最短增广路的效果。
当我们找到沿着残量网络找到一条增广路，增广后，残量网络肯定会变化（至少少了一条边），因此决定允许弧的d数组要进行相应的更新。然而，Dinic 的做法是每次增广都重新计算d数组，这样虽然容易理解，但是速度慢了很多。ISAP改进的地方之一就是，其实没有必要马上更新d数组。这是因为，去掉一条边只可能令路径变得更长，而如果增广之前的残量网络存在另一条最短路，并且在增广后的残量网络中仍存在，那么这条路径毫无疑问是最短的。所以，ISAP 的做法是继续增广，直到遇到死路，才执行重标号操作。
那么，其实重标号就是对d数组的更新，我们如何重标号呢？
我们假设在x点没有找到增广路，那么就可以把d[x]更新为x的所有相邻点中最小残余流量+1。这是为什么呢？因为在这里，其实x和t之间已经不能通过允许弧相通，那么实际上在残余网络中x和t之间的距离是多少呢（不走允许弧），其实不就是相邻点的答案加上1吗？
一种特殊情况就是，如果x没有相邻的点，就把d[x]设为inf就好，表示这个点已经与网络脱节。
最后，实现的时候注意加上gap与当前弧优化。
模板（bzoj3396）

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#define fo(i,a,b) for(ll i=a;i<=b;i++)
#define fd(i,a,b) for(ll i=a;i>=b;i--)
using namespace std; 
typedef long long ll;
const int N=1e5+5;
const int inf=1e9+5;
int n,maxn,u,t,maxflow;
int tot;
int d[N],cur[N],last[N],num[N];
int head[N],next[N],v[N],remain[N];
queue<int> q;
inline void add(int x,int y,int z)
{
    v[++tot]=y;remain[tot]=z;next[tot]=head[x];head[x]=tot;
    v[++tot]=x;remain[tot]=0;next[tot]=head[y];head[y]=tot;
}
inline void bfs(int t)
{
    fo(i,1,maxn)d[i]=maxn;
    d[t]=0;
    fo(i,1,maxn)cur[i]=head[i];
    while (!q.empty()) q.pop();
    q.push(t);  
    while (!q.empty())
    {
        int now=q.front();
        q.pop();
        for (int i=head[now];i!=-1;i=next[i])
        {
            if (d[v[i]]==maxn&&remain[i^1])
            {
                d[v[i]]=d[now]+1;
                q.push(v[i]);
            }

        }
    }
}
inline int augment(int s,int t)
{
    int now=t,ans=inf;
    while (now!=s)
    {
        ans=min(ans,remain[last[now]]);
        now=v[last[now]^1];
    }
    now=t;
    while (now!=s)
    {
        remain[last[now]]-=ans;
        remain[last[now]^1]+=ans;
        now=v[last[now]^1];
    }
    return ans;
}
inline void isap(int s,int t)
{
    bfs(t);
    fo(i,1,maxn)++num[d[i]];
    int now=s;
    while (d[s]<maxn)
    {
        if (now==t)
        {
            maxflow+=augment(s,t);
            now=s;
        }
        bool flag=false;
        int i=cur[now];
        while (i!=-1)
        {
            if (d[v[i]]+1==d[now]&&remain[i])
            {
                flag=true;
                cur[now]=i;
                last[v[i]]=i;
                now=v[i];
                break;
            }
            i=next[i];
        }
        if (!flag)
        {
            int minn=maxn-1;
            int i=head[now];
            while(i!=-1)
            {
                if (remain[i])
                minn=min(minn,d[v[i]]);
                i=next[i];
            }
            if (!(--num[d[now]]))break;
            num[d[now]=minn+1]++;
            cur[now]=head[now];
            if (now!=s)now=v[last[now]^1];
        }
    }
}
int main()
{
    tot=-1;scanf("%d",&n);
    memset(head,-1,sizeof(head));
    memset(next,-1,sizeof(next));
    maxn=60;
    fo(i,1,n)
    {
        char x,y;
        int z;
        cin>>x>>y>>z;
        //scanf("%d%d%d",&x,&y,&z);
        u=x-'A'+2,t=y-'A'+2;
        add(u,t,z); 
    }
    add(1,2,inf);
    add(27,maxn,inf);
    isap(1,maxn);
    printf("%d\n",maxflow);
    return 0;
}