Acwing 164. 可达性统计（拓扑排序+bitset）

由于给定图是DAG所以可以对图进行拓扑排序，这样排序结束后，从入度最大的点开始开始计算可达点，然后一层一层向前统计即可。为了防止同一节点被反复累加，我们采用bitset进行或运算，来保证一个节点之后被计入一次。

// #pragma GCC optimize (2)
// #pragma G++ optimize (2)
#include <iostream>
#include <cstring>
#include <algorithm>
#include <queue>
#include <vector>
#include <unordered_map>
#include <map>
#include <bitset>
#include <set>
#define endl '\n'
#define int long long
#define lowbit(x) x &(-x)
#define rep(i, a, n) for (int i = a; i <= n; i++)
using namespace std;
void TLE(){ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);}
typedef long long ll;
const int N = 3e4 + 10;
int n, m;
int h[N], e[N], ne[N], idx;
int d[N], seq[N], cnt;
bitset<N> f[N];
void topsort()//拓扑排序
{
    queue<int>q;
    for(int i=1;i<=n;i++)
    {
        if(d[i]==0)q.push(i);
    }
    while(q.size())
    {
        auto t=q.front();
        q.pop();
        seq[cnt++]=t;
        for(int i=h[t];~i;i=ne[i])
        {
            int j=e[i];
            if(--d[j]==0)q.push(j);
        }
    }
}
void add(int a, int b)
{
    e[idx] = b, ne[idx] = h[a], h[a] = idx++;
}
void solve()
{
    memset(h, -1, sizeof h);
    cin >> n >> m;
    while (m--)
    {
        int a, b;
        cin >> a >> b;
        add(a, b);
        d[b]++;
    }

    topsort();

    for(int i=n-1;i>=0;i--)
    {
        int j=seq[i];
        f[j][j]=1;
        for(int k=h[j];~k;k=ne[k])
        {
            f[j]|=f[e[k]];
        }
    }
    for(int i=1;i<=n;i++)cout<<f[i].count()<<endl;
}
signed main()
{
    TLE();
    int T;T = 1;
    // int T;cin>>T;
    while (T--)
        solve();
    return 0;
}