codevs 2822 爱在心中 Tarjan 解题报告

题目描述 Description

“每个人都拥有一个梦，即使彼此不相同，能够与你分享，无论失败成功都会感动。爱因为在心中，平凡而不平庸，世界就像迷宫，却又让我们此刻相逢Our Home。”

在爱的国度里有N个人，在他们的心中都有着一个爱的名单，上面记载着他所爱的人（不会出现自爱的情况）。爱是具有传递性的，即如果A爱B，B爱C，则A也爱C。
如果有这样一部分人，他们彼此都相爱，则他们就超越了一切的限制，用集体的爱化身成为一个爱心天使。
现在，我们想知道在这个爱的国度里会出现多少爱心天使。而且，如果某个爱心天使被其他所有人或爱心天使所爱则请输出这个爱心天使是由哪些人构成的，否则输出-1。

输入描述 Input Description

第1行，两个数N、M，代表爱的国度里有N个人，爱的关系有M条。
第2到第M+1行，每行两个数A、B，代表A爱B。

输出描述 Output Description

第1行，一个数，代表爱的国度里有多少爱心天使。
第2行，如果某个爱心天使被其他所有人和爱心天使所爱则请输出这个爱心天使是由哪些人构成的（从小到大排序），否则输出-1。

样例输入 Sample Input

样例输入1：

6 7
1 2
2 3
3 2
4 2
4 5
5 6
6 4

样例输入2：

3 3
1 2
2 1
2 3

样例输出 Sample Output

样例输出1：

2
2 3

样例输出2：

1
-1

数据范围及提示 Data Size & Hint

各个测试点1s

思路

这道题就是找强连通分量然后缩点。输出每个belong数组。

代码

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=100000+5;
int s,n,m,x,y,tim,num=0,top,sum=0,ans,flag,vis[N];
int du[N],dfn[N],low[N],size[N],head[N],belong[N],stack[N];
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
struct edge
{
    int u,v,next;
};
edge ed[N];
void build(int u,int v)
{
    ed[++num].v=v;
    ed[num].next=head[u];
    head[u]=num;
}
void tarjan(int now)
{
    dfn[now]=low[now]=++tim;
    stack[++top]=now;
    vis[now]=1;
    for (int i=head[now];i!=-1;i=ed[i].next)
    {
        int v=ed[i].v;
        if (vis[v]) low[now]=min(low[now],dfn[v]);
        else if (!dfn[v]) {tarjan(v);low[now]=min(low[now],low[v]);}
    }
    if (low[now]==dfn[now])
    {
        sum++;
        belong[now]=sum;size[sum]++;
        for (;stack[top]!=now;top--)
        {
            belong[stack[top]]=sum;
            vis[stack[top]]=0;
            size[sum]++;
        }
        vis[now]=0;
        top--;
    }
}
void suopoint()
{
    for (int i=1;i<=n;i++)
    for (int j=head[i];j!=-1;j=ed[j].next) 
    if (belong[ed[j].v]!=belong[i]) du[belong[i]]++;
}
void init() 
{
    memset(dfn,0,sizeof(dfn));
    memset(du,0,sizeof(du));
    memset(head,-1,sizeof(head));
}
int main()
{
    n=read();m=read();
    init();
    for (int i=1;i<=m;i++)
    {
        x=read();y=read();
        build(x,y);
    }
    for (int i=1;i<=n;i++)
    if (!dfn[i]) tarjan(i);
    suopoint();
    for (int i=1;i<=sum;i++)
    {
        if (size[i]>1) ans++;
        if (!du[i]) {s++;x=1;}
    }
    printf("%d\n",ans);
    if (s==1&&size[x]!=1)
    {
        for (int i=1;i<=n;i++)
        if (belong[i]==x) printf("%d ",i);
    }
    else printf("-1\n");
    return 0;
}