题意:给你n个点,m条边,让你找一些点满足要求:从一个点A出发能够到达的点集S,也必须能够从S出发到达点A。
分析:强联通缩点,然后找出度为0的点。
代码:
#pragma comment(linker,"/STACK:102400000,102400000")
#include <iostream>
#include <string.h>
#include <stdio.h>
#include <algorithm>
#include <vector>
#include <string>
#include <math.h>
#include <queue>
#include <stack>
#include <map>
#include <set>
using namespace std;
typedef long long ll; //ŒÇµÃ±ØÒªµÄʱºòžÄ³ÉÎÞ·ûºÅ
const int maxn=5105; //µãÊý
const int maxm=1000005; //±ßÊý
const int INF=1000000000;
struct EdgeNode
{
int from;
int to;
int next;
}edge[maxm];
int head[maxn],cnt;
void add(int x,int y)
{
edge[cnt].from=x;edge[cnt].to=y;edge[cnt].next=head[x];head[x]=cnt++;
}
void init()
{
cnt=0;
memset(head,-1,sizeof(head));
}
int dfn[maxn],low[maxn],ins[maxn],cixu,scc_count,n,ru[maxn],chu[maxn],sd[maxn];
stack<int>S;
void dfs(int u)
{
int v,x;
dfn[u]=low[u]=++cixu;
S.push(u);
ins[u]=1;
for(int i=head[u];i!=-1;i=edge[i].next)
{
v=edge[i].to;
if(!dfn[v]){
dfs(v);
low[u]=min(low[u],low[v]);
}
else if(ins[v]) //Èç¹ûÒÑŸËѹý£¬²¢ÇÒÔÚÕ»ÀïÃ棬Èç¹û²»ÔÚÕ»ÀïÃ棬˵Ã÷žÃµãÊôÓÚÆäËûµÄÇ¿Á¬Íš·ÖÁ¿ÀïÃæ
low[u]=min(low[u],dfn[v]); //ŽËŽŠlowŽúÌædfnÒ²ÊÇ¿ÉÒԵģ¬µ«ÊÇÔÚtarjanÇóÎÞÏòÍŒµÄʱºò²»ÄÜ
}
if(dfn[u]==low[u]) //±íÊŸuΪŽËÇ¿Á¬Íš·ÖÁ¿µÚÒ»žö±»·¢Ïֵĵ㣬ÊôÓÚžÃÇ¿Á¬Íš·ÖÁ¿µÄËùÓеÄÆäËûµã¶ŒÔÚÕ»ÖУ¬ËûµÄÉÏÃæ
{
scc_count++;
do
{
x=S.top();
S.pop();
sd[x]=scc_count;
ins[x]=0;
}while(x!=u);
}
}
void tarjan()
{
for(int i=1;i<=n;i++)sd[i]=i; //Ëõµã±êŒÇÊý×é
memset(ins,0,sizeof(ins));
memset(dfn,0,sizeof(dfn));
memset(low,0,sizeof(low));
cixu=scc_count=0;
for(int i=1;i<=n;i++) //¿ÉÄÜËùžøµÄ͌Ϊ²»Á¬ÍšµÄÍŒ£¬·Ö¿é£¬ÒòŽËÐèÒª±éÀúËùÓеĵã
{
if(!dfn[i])
dfs(i);
}
}
void SD()
{
int x,y,len=cnt;
init();
memset(ru,0,sizeof(ru));
memset(chu,0,sizeof(chu));
for(int i=0;i<len;i++){
x=sd[edge[i].from]; y=sd[edge[i].to];
if(x==y)continue;
add(x,y);
chu[x]++; ru[y]++;
}
}
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
int m,i,x,y,c[5005];
while(~scanf("%d",&n)&&n)
{
scanf("%d",&m);
init();
for(i=1;i<=m;i++){
scanf("%d%d",&x,&y);
add(x,y);
}
tarjan();
SD();
int ans=0;
for(i=1;i<=n;i++){
if(chu[sd[i]]==0)c[++ans]=i;
}
if(ans==0)printf("\n");
else{
for(i=1;i<=ans;i++){
if(i==1)printf("%d",c[i]);
else printf(" %d",c[i]);
}
printf("\n");
}
}
return 0;
}