AcWing164. 可达性统计（拓扑排序&&位运算）

输入样例：

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;
typedef long long ll;
const int N=3e4+5,M=3e4+5;
int n,m,d[N],p[N];
int h[N],e[M],ne[M],idx;
void add(int a,int b){
	e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
bitset<N>f[N];
int main(){
	scanf("%d%d",&n,&m);
	memset(h,-1,sizeof h);
	for(int i=1;i<=m;i++){
		int x,y;
		scanf("%d%d",&x,&y);
		d[y]++;
		add(x,y);
		
	}
	vector<int>v;
	queue<int>q;
	for(int i=1;i<=n;i++) if(!d[i]) q.push(i);	
	int cnt=0;
	while(q.size()){
		int t=q.front();
		q.pop();
		v.push_back(t);
		for(int i=h[t];~i;i=ne[i]){
			int j=e[i];
			p[j]=max(p[j],p[t]+1);
			if(--d[j]==0) q.push(j);
		}
	}
	for(int i=v.size()-1;i>=0;i--){
		int j=v[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++) printf("%d\n",f[i].count());
	return 0;
}