可达性统计------------------------------思维(拓扑排序+dp)

解析:
由题意可知是一个有向无环图 DAG，所以我们做一遍拓扑排序。
然后根据dp的思想
f(u)= f(u1) ∪ f(u2) ∪ f(u3)
对于想要求解出u的情况。必须把u1,u2,u3求解出来。所以我们倒着遍历拓扑排序后的数组。然后借助bitset数组统计1的情况即可

#include<bits/stdc++.h>
using namespace std;
const int N=3e4+10;
bitset<N> f[N];
vector<int> G[N],v;
int in[N];
int n,m,a,b;
void topsort()
{
    queue<int> q;
    for(int i=1;i<=n;i++)
    {
        if(!in[i]) q.push(i);
    }
    while(q.size())
    {
        int t=q.front();q.pop();
        v.push_back(t);
        for(auto x:G[t])
        {
            if(--in[x]==0) q.push(x);
        }
    }
}
int main()
{
    cin>>n>>m;
    for(int i=1;i<=m;i++)
    {
        cin>>a>>b;
        G[a].push_back(b);
        in[b]++;
    }
    topsort();
    for(int i=v.size()-1;i>=0;i--)
    {
        int x=v[i];
        f[x][x]=1;
        for(auto it:G[x])
        {
            f[x]|=f[it];
        }
    }
    for(int i=1;i<=n;i++)cout<<f[i].count()<<endl;
}