AcWing 164. 可達性統計 (拓撲序列+bitset)

整理的算法模板：ACM算法模板總結（分類詳細版）

給定一張N個點M條邊的有向無環圖，分別統計從每個點出發能夠到達的點的數量。

輸入格式

第一行兩個整數N,M，接下來M行每行兩個整數x,y，表示從x到y的一條有向邊。

輸出格式

輸出共N行，表示每個點能夠到達的點的數量。

數據範圍

1≤N,M≤300001≤N,M≤30000

輸入樣例：

10 10
3 8
2 3
2 5
5 9
5 9
2 3
3 9
4 8
2 10
4 9

輸出樣例：

1
6
3
3
2
1
1
1
1
1

首先有向無環圖是一個拓撲圖，所以只需要按照拓撲圖來做就行；

要求i的狀態，那麼後面所有集合的狀態都要事先求出來；所以，需要逆拓撲序來處理；

然後藉助STL裏面的bitset容器實現集合的合併；

#include <bits/stdc++.h>
using namespace std;
const int N=30007,M=N;
int h[N],e[M],ne[M],idx;
int n,m;
int din[N];
bitset<N> f[N];
vector<int> v;
void add(int a,int b)
{
    e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
int main()
{
    memset(h,-1,sizeof h);
    cin >>n>>m;
    for(int i=0;i<m;i++)
    {
        int a,b;
        cin >>a>>b;
        add(a,b);
        din[b]++;
    }
    queue<int> q;
    for(int i=1;i<=n;i++)
    {
        if(!din[i]) q.push(i);
    }
    while(q.size())
    {
        int x=q.front();
        v.push_back(x);
        q.pop();
        for(int i=h[x];~i;i=ne[i])
        {
            int j=e[i];
            if(--din[j]==0) q.push(j);
        }
    }
    reverse(v.begin(),v.end());
    for(auto i : v)
    {
        f[i][i]=1;
        for(int j=h[i];~j;j=ne[j])
        {
            int k=e[j];
            f[i]|=f[k];
        }
    }
    for(int i=1;i<=n;i++) cout <<f[i].count()<<endl;
}