POJ

题意：

给你n个点，和m条单向边，问你有多少点满足(G)={v∈V|∀w∈V:(v→w)⇒(w→v)}关系，并把这些点输出（要注意的是这个关系中是蕴含关系而不是且（&&）关系）

题解：

单独一个强连通分量中的所有点是满足题目要求的
但如果它连出去到了其他点那里，要么成为新的强连通分量，要么失去原有的符合题目要求的性质
所以只需tarjan缩点求出所有强连通分量，再O(E)枚举所有边，是否会成为连接一个分量与另一个分量的边——即一条出度——即可
如果一个分量没有出度，那么他中间的所有点都是符合题目要求的点

代码（多组输入记得每次初始化）：

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<queue>
#include<map>
#include<vector>
#include<math.h>
#define mem(a,x) memset(a,x,sizeof(a))
using namespace std;
const int maxn=50005;
const int mod=26;
const int INF=0x3f3f3f3f;
struct edge
{
    int u,v,next;
    bool sign;
} e[maxn*2];
int head[maxn],cnt;
void add_edge(int x,int y)
{
    e[cnt].u=x;
    e[cnt].v=y;
    e[cnt].next=head[x];
    head[x]=cnt++;
}
int dfn[maxn],low[maxn],stacks[maxn],top,tot;
int taj;
int belong[maxn],visit[maxn],result[maxn];
vector<int>w[maxn];
void tarjan(int x,int fx)
{
    dfn[x]=low[x]=++tot;
    stacks[top++]=x;
    visit[x]=1;
    for(int i=head[x]; i!=-1; i=e[i].next)
    {
        int v=e[i].v;
        if(!dfn[v])
        {
            tarjan(v,x);
            low[x]=min(low[x],low[v]);
            if(dfn[x]<low[v])
            {
                e[i].sign=1;
            }
        }
        else if(visit[v])
        {
            low[x]=min(low[x],dfn[v]);
        }
    }
    if(low[x]==dfn[x])
    {
        int now;
        taj++;
        w[taj].clear();
        do
        {
            now=stacks[--top];
            visit[now]=0;
            belong[now]=taj;
            w[taj].push_back(now);
        }
        while(now!=x);
    }
}
void tarjan_init(int n)
{
    memset(visit,0,sizeof(visit));
    memset(low,0,sizeof (low));
    memset(dfn,0,sizeof (dfn));
    memset(stacks,0,sizeof (stacks));
    memset(belong,0,sizeof belong);
    top=tot=taj=0;
    for(int i=1; i<=n; ++i)
    {
        if(!dfn[i]) tarjan(i,i);
    }
}
vector<int>g[maxn];
int cu[maxn],ru[maxn];
void suodian()
{
    memset(cu,0,sizeof(cu));
    memset(ru,0,sizeof(ru));
    for(int i=1; i<=taj; ++i)
        g[i].clear();
    for(int i=0; i<cnt; ++i)
    {
        int u=belong[e[i].u];
        int v=belong[e[i].v];
        if(u!=v) g[u].push_back(v),cu[u]++,ru[v]++;//printf("%d %d
",u,v);
    }
}
int vis[maxn];
void init()
{
    memset(head,-1,sizeof(head));
    cnt=0;
}
int main()
{
    int n,m;
    while(~scanf("%d",&n))
    {
        if(n==0) return 0;
        scanf("%d",&m);
        init();
        for(int i=1; i<=m; ++i)
        {
            int a,b;
            scanf("%d%d",&a,&b);
            add_edge(a,b);
        }
        tarjan_init(n);
        suodian();
        int index=0;
        for(int i=1;i<=n;++i)
        {
            if(cu[belong[i]]==0)
                result[index++]=i;
        }
        if(index==0)
        {
            printf("
");
            continue;
        }
        sort(result,result+index);
        for(int i=0;i<index;++i)
            if(i==index-1)
                printf("%d
",result[i]);
            else printf("%d ",result[i]);
    }
    return 0;
}