可达性统计

题目描述

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

输入描述

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

输出描述

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

#include<iostream>
#include<cstdio>
#include<queue>
#include<bitset>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=105;
const int M=N*N;
const int INF=0x3f3f3f3f;
queue<int>Q;
bitset<N>c[N];
int head[N],tot=0;
struct node{
	int ver;
	int nxt;
}tu[M];
int n,m,cnt;
int deg[N];
int f[N];
void init(){
	tot=0;
	memset(head,-1,sizeof(head));
}
void add(int x,int y){
	tu[++tot].ver=y;
	tu[tot].nxt=head[x];
	head[x]=tot;
	deg[y]++;
}
void topsort(){
	for(int i=1;i<=n;i++){
		if(deg[i]==0){
			Q.push(i);
		}
	}
	while(Q.size()){
		int x=Q.front();
		Q.pop();
		f[++cnt]=x;
		for(int i=head[x];~i;i=tu[i].nxt){
			int y=tu[i].ver;
			if(--deg[y]==0){
				Q.push(y);
			}
		}
	}
}
void solve(){
	for(int i=cnt;i>=1;i--){
		int x=f[i];
		c[x][x]=1;
		for(int j=head[x];~j;j=tu[j].nxt){
			int y=tu[j].ver;
			c[x]|=c[y];
		}
		
	}
}
signed main(){
	init();
	scanf("%d %d",&n,&m);
	for(int i=1;i<=m;i++){
		int x,y;
		scanf("%d %d",&x,&y);
		add(x,y);
	}
	topsort();
	solve();
	for(int i=1;i<=n;i++){
		printf("%d\n",c[i].count());
	}
	return 0;
}