可达性统计

题目描述
给定一张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=30010;
int n,m;
int p[N];//通向y的节点的个数
vector<int> e[N];//存储图
vector<int> tuopo;//存储拓扑排序的结果
bitset<N> f[N];
void tuoposort(){
    tuopo.clear();
    queue<int> q;
    for(int i=1;i<=n;++i){
        if(p[i]==0)//如果前面没有节点就直接进入拓扑序列
        q.push(i);
        }
    while(q.size()){
        int w=q.front();
        q.pop();
        tuopo.push_back(w);
        for(int i=0;i<e[w].size();++i)//遍历队列，看看有没有只有一个节点的
        if(--p[e[w][i]]==0)
        q.push(e[w][i]);
    }
}
void init()
{
    cin>>n>>m;
    memset(p,0,sizeof(p));
    for(int i=1;i<=m;++i)
    {
        int x,y;
        cin>>x>>y;
        p[y]++;//记录入度的个数
        e[x].push_back(y);
    }
}
void find(){
    tuoposort();
    for(int i=tuopo.size()-1;i>=0;i--)
    {
        int x=tuopo[i];
        f[x].reset(),f[x][x]=1;
        for(int k=0;k<e[x].size();++k)
        f[x]|=f[e[x][k]];
    }
    for(int i=1;i<=n;++i)
    cout<<f[i].count()<<endl;
} int main()
{
     ios::sync_with_stdio(false);
    init();

    find();

    return 0;
}