Wikioi 2822爱在心中（强连通缩点+dfs）

2822 爱在心中（强连通缩点+遍历）
2822 爱在心中
时间限制: 1 s
空间限制: 128000 KB
题目等级 : 钻石 Diamond
题解
题目描述 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

program df;
var i,j,n,m,x,y,z,s,max,min,t:longint;
a,e:array[0..1000,0..1000] of longint;
low,dis,tm,x1,y1,c,f:array[0..1000000] of longint;
zhan,re:array[0..10000000] of boolean;

function find(x:longint):longint;
begin
if f[x]<>x then f[x]:=find(f[x]);
find:=f[x];
end;

procedure dfs(x:longint); //求强连通分量
var i,y:longint;
begin
inc(t);
low[x]:=t;
zhan[x]:=true;
re[x]:=true;
for i:=1 to a[x,0] do
begin
y:=a[x,i];
if not re[y] then dfs(y);
y:=find(y);
if (low[x]>low[y]) and (zhan[y]) then
begin
f[x]:=y;
low[x]:=low[y];
end;
end;
zhan[x]:=false;
end;

procedure spfa(x:longint);
var i,y:longint;
begin
for i:=1 to a[x,0] do
begin
y:=a[x,i];
if zhan[y] then
begin
zhan[y]:=false;
inc(t);
spfa(y);
end;
end;
if t=n then exit;
end;

procedure sea(l,r:longint);
var i,j,n,m:longint;
begin
i:=l; j:=r;
m:=e[max,(l+r) div 2];
repeat
while e[max,i]小于m do inc(i);
while e[max,j]>m do dec(j);
if i<=j then
begin
n:=e[max,i]; e[max,i]:=e[max,j]; e[max,j]:=n;
inc(i); dec(j);
end;

until i>j;
if r>i then sea(i,r);
if j>l then sea(l,j);
end;

begin
readln(n,m);
for i:=1 to n do
f[i]:=i;
for i:=1 to m do
begin
readln(x1[i],y1[i]);
x:=x1[i];
y:=y1[i];
inc(a[x,0]);
a[x,a[x,0]]:=y;
end;
for i:=1 to n do
if not re[i] then
begin
t:=0;
fillchar(low,sizeof(low),0);
dfs(i);
end;

for i:=1 to n do
if find(i)<>i then
begin
x:=find(i);
if re[x] then
begin
inc(s); //记录多少个环（即多少个爱心天使）
c[s]:=x; //每个环的根节点
re[x]:=false;
inc(e[x,0]);
e[x,e[x,0]]:=x; //每个环的代表(同样属于环，最后用于排序输出)
end;
inc(e[x,0]);
e[x,e[x,0]]:=i; //每个环的各个元素（用于最后排序输出）
end;
writeln(s);
if s=0 then
begin
writeln(-1);
halt;
end;
fillchar(a,sizeof(a),0);
for i:=1 to m do
begin
inc(a[y1[i],0]);
a[y1[i],a[y1[i],0]]:=x1[i]; //反向建边（作用在下面）
end;

for i:=1 to s do
begin
t:=1;
fillchar(zhan,sizeof(zhan),true);
zhan[c[i]]:=false;
spfa(c[i]); //每一个爱心天使都遍历一遍（反向建边时才能实现）
if t=n then //若恰好能够能遍历完，则排序输出（注：不可能出现两个，否则图将无法连通）
begin
max:=c[i];
sea(1,e[max,0]);
for j:=1 to e[max,0] do
write(e[max,j],’ ‘);
halt;
end;
end;
writeln(-1); //不存在则输出-1
end.