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

给定一张N个点M条边的有向无环图，分别统计从每个点出发能够到达的点的数量。

输入格式
第一行两个整数N,M，接下来M行每行两个整数x,y，表示从x到y的一条有向边。

输出格式
输出共N行，表示每个点能够到达的点的数量。

数据范围
1≤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

#include <bits/stdc++.h>
using namespace std;
const int N = 30005;
bitset<N> f[N];
struct node {
    int v, ne;
} s[N];
int h[N], vis[N], a[N], b[N];
int cnt = 0, k = 0, u, v, n, m;
void add(int u, int v) {
    s[++cnt].v = v;
    s[cnt].ne = h[u];
    h[u] = cnt;
}
void tuopu() {
    queue<int> q;
    for (int i = 1; i <= n; i++)
        if (!b[i]) q.push(i);
    while (!q.empty()) {
        int u = q.front();
        a[++k] = u;
        q.pop();
        for (int i = h[u]; i; i = s[i].ne) {
            int v = s[i].v;
            b[v]--;
            if (!b[v]) q.push(v);
        }
    }
}
void solve() {
    for (int i = k; i; i--) {
        int u = a[i];
        f[u][u] = 1;
        for (int j = h[u]; j; j = s[j].ne) {
            v = s[j].v;
            f[u] |= f[v];
        }
    }
}
int main() {
    ios::sync_with_stdio(false);
    std::cin.tie(0);
    cin >> n >> m;
    while (m--) {
        cin >> u >> v;
        add(u, v);
        b[v]++;
    }
    tuopu();
    solve();
    for (int i = 1; i <= n; i++) 
        cout << f[i].count() << '\n';
}