bzoj3887 [Usaco2015 Jan]Grass Cownoisseur tarjan+拓补排序

题意不说了。
%%%popoqqq
既然没有要求一条边不能被经过两次，那么在强连通分量内的所有点都有贡献，所以先缩点，这个挺显然的。
问题是这个逆边要怎么选，，有一个明显结论
将某条边反向后 缩点之后的图形成了一个包含1号节点所在强连通分量的环 这样才能使答案增加
把这个环从反向的边和1号节点所在的强连通分量断开，一条路是从1到n，另外一条是从n到1.
那么我们缩点以后拓补排序一下，把正图和反图中1号节点所在的强连通分量的dis算出来然后枚举一下看看那条边反过来就好了。

#include
#include
#include
#include
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=2e5+5;
int n,m,ans;
int bel[N],go[N],next[N],head[N];
int tot,low[N],bl,dfn[N],cnt,stack[N],top,size[N];
bool vis[N],ins[N];

struct toplogysort
{
    int f[N],du[N],q[N],t,w;
    int head[N],go[N],next[N],tot;
    inline void add(int x,int y) 
    {
        du[y]++;
        go[++tot]=y;
        next[tot]=head[x];
        head[x]=tot;
    }
    inline void tsort()
    {
        //using namespace Origin_graph;
        memset(f,0xef,sizeof(f));
        //printf("%d\n",f[0]);
        f[bel[1]]=size[bel[1]];
        fo(i,1,bl)if (!du[i])q[++w]=i;
        while (tint x=q[++t];
            for(int i=head[x];i;i=next[i])
            {
                int v=go[i];
                f[v]=max(f[v],f[x]+size[v]);
                if (!--du[v])q[++w]=v;
            }
        }
    }
}pos,neg;
inline void add(int x,int y)
{
    go[++tot]=y;
    next[tot]=head[x];
    head[x]=tot;
}
inline void tarjan(int x)
{
    vis[x]=1;
    dfn[x]=low[x]=++cnt;
    ins[x]=1;
    stack[++top]=x;
    for(int i=head[x];i;i=next[i])
    {
        int v=go[i];
        if (!vis[v])
        {
            tarjan(v);
            low[x]=min(low[x],low[v]);
        }
        else if (ins[v])low[x]=min(low[x],dfn[v]);
    }
    if (dfn[x]==low[x])
    {
        ++bl;
        int t=-1;
        while (t!=x)
        {
            t=stack[top--];
            ins[t]=0;
            bel[t]=bl;
            size[bl]++;
        }
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    fo(i,1,m)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        add(x,y);

    }
    fo(i,1,n)if (!vis[i])tarjan(i);
    fo(x,1,n)
    for(int i=head[x];i;i=next[i])
    {
        int v=go[i];
        if (bel[x]!=bel[v])
        {
            pos.add(bel[x],bel[v]);
            neg.add(bel[v],bel[x]);
        } 
    }
    pos.tsort();
    neg.tsort();
    ans=2*size[bel[1]];
    fo(x,1,n)
    for(int i=head[x];i;i=next[i])
    {
        int v=go[i];
        ans=max(ans,pos.f[bel[v]]+neg.f[bel[x]]);
    }
    printf("%d\n",ans-size[bel[1]]);
}