Codevs 2822 爱在心中

http://codevs.cn/problem/2822/

思路
Tarjin算法找出图中所有的节点数>1的环即爱心天使，将每个爱心天使缩成一个点，重新建图，在新图中找到一个入度为0的点（至多一个），只有此点可能被所有人爱。
那么我们在建图的时候反向建边，看此点能否走到其他所有点即可。

代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
int n,m,ru,rv,tot,cnt,num,sum,ans,s;
int first[100010],nxt[100010],dfn[100010],low[100010],stak[100010],siz[100010],sccno[100010],Ed[100010],As[100010];
bool flag;
bool b[100010],arv[100010];
struct edge
{
    int u,v;
}l[100010];
inline void build(int f,int t)
{
    l[++tot]=(edge){f,t};
    nxt[tot]=first[f];
    first[f]=tot;
}
void dfs(int k)
{
    dfn[k]=low[k]=++tot;
    stak[++num]=k;
    for(int i=first[k];i!=-1;i=nxt[i])
    {
        int x=l[i].v;
        if(!dfn[x])
        {
            dfs(x);
            low[k]=min(low[k],low[x]);
        }
        else if(!sccno[x])
        low[k]=min(low[k],dfn[x]);
    }
    if(low[k]==dfn[k])
    {
        cnt++;
        while(1)
        {
            sccno[stak[num]]=cnt;
            siz[cnt]++;
            num--;
            if(stak[num+1]==k)
            break;
        }
    }
}
void bfs(int k)
{
    for(int i=first[k];i!=-1;i=nxt[i])
    {
        int x=l[i].v;
        arv[x]=1;
        bfs(x);
    }
}
int main()
{
    memset(first,-1,sizeof(first));
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&ru,&rv);
        build(ru,rv);
    }
    sum=tot;
    tot=0;
    for(int i=1;i<=n;i++)
    {
        if(!dfn[i])
        dfs(i);
    }
    ans=cnt;
    for(int i=1;i<=cnt;i++)
    if(siz[i]<=1)
    ans--;
    printf("%d\n",ans);
    num=0;
    for(int i=1;i<=sum;i++)
    {
        if(sccno[l[i].u]!=sccno[l[i].v])
        {
            Ed[++num]=sccno[l[i].u];
            As[num]=sccno[l[i].v];
            b[sccno[l[i].u]]=1;
        }
    }
    for(int i=1;i<=cnt;i++)
    {
        if(b[i]==0)
        {
            s=i;
            break;
        }
    }
    if(s!=0&&ans!=0)
    {
        memset(first,-1,sizeof(first));
        memset(nxt,0,sizeof(nxt));
        memset(l,0,sizeof(l));
        tot=0;
        for(int i=1;i<=num;i++)
        build(As[i],Ed[i]);
        bfs(s);
        arv[s]=1;
        for(int i=1;i<=cnt;i++)
        {
            if(arv[i]==0)
            {
                flag=1;
                break;
            }
        }
        tot=0;
        if(!flag)
        {
            for(int i=1;i<=n;i++)
            {
                if(sccno[i]==s)
                Ed[++tot]=i;
            }
            sort(Ed+1,Ed+tot+1);
            for(int i=1;i<=tot;i++)
            printf("%d ",Ed[i]);
        }
        else printf("-1");
    }
    else printf("-1");
    return 0;
}