可達性-牛客（tarjan縮點）

題意： 給出一個 0 ≤ N ≤ 10⁵ 點數、0 ≤ M ≤ 10⁵ 邊數的有向圖，輸出一個儘可能小的點集，使得從這些點出發能夠到達任意一點，如果有多個這樣的集合，輸出這些集合升序排序後字典序最小的。

解：縮點，如果只有一個點直接輸出1，整張圖連通。

縮點過程，維護兩個陣列fa[]點集字典序最小的店和num[]儲存新點的點集,縮點所代表的點有哪些

直接tarjan演算法求即可。

樣例模擬：

回溯到3，回溯到5

沿5往下走5->1->2->5

回溯到2,更改low2=1；回溯到1，出棧；

出棧1 2 5 3 6之後，回溯到2，1…

程式碼：

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+10;
const int mod=142857;
const int inf=0x3f3f3f3f;
typedef long long ll;
typedef pair<int,int> pii;
const int N=5e5+10;
 
struct node
{
    int v,next;
}G[maxn];
int head[maxn],cnt;
 
void add(int x,int y)
{
    G[++cnt]=node{y,head[x]};
    head[x]
 
=cnt;
}
 
int dfn[maxn],low[maxn],tot;
stack<int> st;
int vis[maxn];
int fa[maxn],num[maxn],ncnt;
 
void tarjan(int u)
{
    dfn[u]=low[u]=++tot;
    st.push(u);
    vis[u]=1;
    for(int i=head[u]; i!=-1; i=G[i].next)
    {
        int v=G[i].v;
        if(!dfn[v])
        {
            tarjan(v);
            low[u]
 
=min(low[u],low[v]);
        }
        else if(vis[v])
            low[u]=min(low[u],dfn[v]);
    }
    if(dfn[u]==low[u])
    {
        ncnt++;
        int t;
        do{
            t=st.top();
            st.pop();
            num[t]=ncnt;
            fa[ncnt]=min(t,fa[cnt]);//初始置無限大;
            vis[t]=0;
        }while(t!=u);
    }
}
int n,m;
int ans[maxn],x;
int in[maxn];
int main()
{
    cin>>n>>m;
    memset(fa,inf,sizeof(fa));
    memset(head,-1,sizeof(head));
    for(int i=1; i<=m; i++)
    {
        int u,v;
        cin>>u>>v;
        add(u,v);
    }
    for(int i=1; i<=n; i++)
    {
        if(!dfn[i])
         tarjan(i);
    }
    if(ncnt==1)
     cout<<1<<endl;
    else
    {
        for(int i=1; i<=n; i++)
        {
            for(int j=head[i]; j!=-1; j=G[j].next)
            {
                int u=num[i];
                int v=num[G[j].v];
                if(u!=v)
                 in[v]++;
            }
        }
        for(int i=1; i<=ncnt; i++)
        {
            if(!in[i])
             ans[x++]=fa[i];
        }
        sort(ans,ans+x);
        cout<<x<<endl;
        for(int i=0; i<x; i++) cout<<ans[i]<<" ";
        cout<<"\n";
    }
    system("pause");
 
}