ACwing-164. 可达性统计 拓扑排序

给定一张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

以某点x出发能够到达的点集合f(x):
f(x) = {x} ∪ ( ∪ f(y): 存在有向边(x,y) )
拓扑排序可以保证任意一条边(x,y) x排在y前
按倒序拓扑序的顺序计算即可
可以用bitset存储集合状态
求并集就是按位或

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<vector>
#include<queue>
#include<bitset>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=3e4+3;
struct Edge{
	int to,next;
	Edge(){}
	Edge(int to,int next):to(to),next(next){}
}E[N];
int head[N],tot;
inline void addedge(int u,int v){
	E[tot]=Edge(v,head[u]),head[u]=tot++;
}
int in[N],a[N];
int n,m,cnt,u,v;
void topo(){
	queue<int>q;
	for(int i=1;i<=n;++i)if(!in[i])q.push(i);
	while(!q.empty()){
		u=q.front(),q.pop();
		a[++cnt]=u;
		for(int i=head[u];v=E[i].to,~i;i=E[i].next)if(!--in[v])q.push(v);
	}
}
bitset<N>f[N];
int main(){
	scanf("%d%d",&n,&m);
	memset(head,-1,sizeof(head));
	for(int i=1;i<=m;++i)scanf("%d%d",&u,&v),addedge(u,v),++in[v];
	topo();
	for(int i=n;u=a[i],f[u][u]=1,i;--i)
		for(int j=head[u];v=E[j].to,~j;j=E[j].next)f[u]|=f[v];
	for(int i=1;i<=n;++i)printf("%d\n",f[i].count());
	return 0;
}