bzoj 1512 [POI2006]Pro-Professor Szu tarjan dp

先tarjan缩点，然后拓扑dp。
似乎有一些细节。。。

#include <bits/stdc++.h>
using namespace std;
#define N 1000010
#define inf 36501
int n,m,tot,top,cnt,scc,ans;
int head[N],nex[N],to[N];
int deep[N],low[N],st[N],bel[N];
bool ins[N],vis[N];
int du[N],size[N],f[N];
vector<int>vec[N];
void add(int x,int y)
{
    tot++;
    nex[tot]=head[x];head[x]=tot;
    to[tot]=y;
}
void tarjan(int x)
{
    deep[x]=low[x]=++cnt;
    st[++top]=x;ins[x]=1;
    for(int i=head[x];i;i=nex[i])
    {
        if(!deep[to[i]])
            tarjan(to[i]),low[x]=min(low[x],low[to[i]]);
        else if(ins[to[i]])
            low[x]=min(low[x],deep[to[i]]);
    }
    if(deep[x]==low[x])
    {
        int tmp=st[top--];
        ins[tmp]=0;bel[tmp]=++scc;
        while(tmp!=x)
        {
            tmp=st[top--];
            ins[tmp]=0;bel[tmp]=scc;
        }
    }
}
int main()
{
    //freopen("tt.in","r",stdin);
    scanf("%d%d",&n,&m);n++;
    for(int i=1,x,y;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        add(y,x);
        if(x==y)vis[x]=1;
    }
    for(int i=1;i<=n;i++)
        if(!deep[i])tarjan(i);

    for(int i=1;i<=n;i++)size[bel[i]]++;
    for(int i=1;i<=n;i++)
        for(int j=head[i];j;j=nex[j])
            if(bel[i]!=bel[to[j]])
                vec[bel[i]].push_back(bel[to[j]]),du[bel[to[j]]]++;
    for(int i=1;i<=n;i++)
        if(vis[i])size[bel[i]]++;
    int h=1,r=0;
    for(int i=1;i<=n;i++)
        if(!du[i])st[++r]=i;
    f[bel[n]]=1;
    if(size[bel[n]]>1)f[bel[n]]=inf;
    while(r>=h)
    {
        int tmp=st[h++];
        for(int i=0,t;i<vec[tmp].size();i++)
        {
            f[t=vec[tmp][i]]+=f[tmp];
            if(f[t]>inf||(size[t]>1&&f[t]))f[t]=inf;
            if(--du[t]==0)st[++r]=t;
        }
    }
    for(int i=1;i<n;i++)
        ans=max(ans,f[bel[i]]);
    if(ans!=inf)printf("%d\n",ans);
    else puts("zawsze");
    for(int i=1;i<n;i++)
        if(f[bel[i]]==ans)
            st[++top]=i;
    printf("%d\n",top);
    for(int i=1;i<=top;i++)
        printf("%d ",st[i]);
    puts("");
    return 0;
}