可达性（强连通+入度）

题目描述 
给出一个 0 ≤ N ≤ 105 点数、0 ≤ M ≤ 105 边数的有向图，
输出一个尽可能小的点集，使得从这些点出发能够到达任意一点，如果有多个这样的集合，输出这些集合升序排序后字典序最小的。
输入描述:
第一行为两个整数 1 ≤ n, m ≤ 105，
接下来 M 行，每行两个整数 1 ≤ u, v ≤ 105 表示从点 u 至点 v 有一条有向边。
数据保证没有重边、自环。
输出描述:
第一行输出一个整数 z，表示作为答案的点集的大小；
第二行输出 z 个整数，升序排序，表示作为答案的点集。
示例1
输入


7 10
4 5
5 1
2 5
6 5
7 2
4 2
1 2
5 3
3 5
3 6
输出


2

4 7

思路：

就是进行类似于缩点的操作，然后进行入度的判断。

具体详见代码：

#include<bits/stdc++.h>
#define MAXN 100000+10
#define MAXM 100000+10
#define INF 100000+10
using namespace std;
struct Edge
{
    int from, to, next;
}edge[MAXM];

int head[MAXN], edgenum;
int dfs_clock;//时间戳
int sccno[MAXN], scc_cnt;//sccno[i]表示i属于哪个SCC scc_cnt是SCC计数器
bool Instack[MAXN];
int low[MAXN], dfn[MAXN];
int minn_scc_num[MAXN];
int in[MAXN];
int n, m;
stack<int> S;


void init()
{
    edgenum = 0;
    memset(head, -1, sizeof(head));
}
void addEdge(int u, int v)
{
    Edge E ={u, v, head[u]};
    edge[edgenum] = E;
    head[u] = edgenum++;
}
void getMap()
{
    int a, b;
    while(m--)
    {
        scanf("%d%d", &a, &b);
        addEdge(a, b);
    }
}
void tarjan(int u,int fa)
{
    int v;
    low[u] = dfn[u] = ++dfs_clock;
    S.push(u);
    Instack[u] = true;
    for(int i = head[u]; i != -1; i = edge[i].next)
    {
        v = edge[i].to;
        if(!dfn[v])
        {
            tarjan(v, u);
            low[u] = min(low[u], low[v]);
        }
        else if(Instack[v])
        low[u] = min(low[u], dfn[v]);//反向边
    }
    if(low[u] == dfn[u])
    {
        scc_cnt++;//记录SCC数目
        minn_scc_num[scc_cnt]=MAXN;
        for(;;)
        {
            v = S.top();
            S.pop();
            sccno[v]=scc_cnt;
            Instack[v] = false;
            minn_scc_num[scc_cnt]=min(v,minn_scc_num[scc_cnt]);
            if(v == u)
                break;
        }
    }
}
void find_cut(int l, int r)
{
    memset(Instack, false, sizeof(Instack));
    memset(low, 0, sizeof(low));
    memset(dfn, 0, sizeof(dfn));
    memset(sccno, 0, sizeof(sccno));
    scc_cnt = dfs_clock = 0;
    for(int i = l; i <= r; i++)
    if(!dfn[i])
        tarjan(i, -1);
}
int main()
{
        cin>>n>>m;
        init();
        getMap();
        find_cut(1, n);
        for(int i=1;i<=n;i++)
        {
            in[i]=0;
        }
        for(int i=0;i<edgenum;i++)
        {
            if(sccno[edge[i].from]!=sccno[edge[i].to])
                in[sccno[edge[i].to]]++;
        }
        int chu[MAXN];
        int ans=0;
        for(int i=1;i<=scc_cnt;i++)
        {
            if(!in[i])
            {
                chu[ans++]=minn_scc_num[i];
            }
        }
        sort(chu,chu+ans);
        printf("%d\n",ans);
        for(int i=0;i<ans;i++)
        {
            if(i==0)
                printf("%d",chu[i]);
            else
                printf(" %d",chu[i]);
        }
        printf("\n");

    return 0;
}