可达性统计

bitset topsort

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<bitset>
#include<queue>
using namespace std;
int deg[30010];
int path[30010],ans=0;
int h[30010],ver[30010],ne[30010],cnt=0;
bitset<30010> f[30010];int n,m;
void add(int a,int b)
{
	ver[cnt]=b;
	ne[cnt]=h[a];
	h[a]=cnt++;
}
void tuosort()
{
	queue<int> q;
	for(int i=1;i<=n;i++)
	{
		if(!deg[i]) q.push(i);
	}
	while(q.size())
	{
		int t=q.front();
		q.pop();
		path[ans++]=t;
		for(int i=h[t];i!=-1;i=ne[i])
		{
			int y=ver[i];
			if(--deg[y]==0) q.push(y);
		}
	}
}
int main()
{
	
	cin>>n>>m;
	memset(h,-1,sizeof h);
	while(m--)
	{
		int x,y;
		scanf("%d%d",&x,&y);
		add(x,y);
		deg[y]++;
	}
	tuosort();
	for(int i=n-1;i>=0;i--)
	{
		int j=path[i];
		f[j][j]=1;
		for(int k=h[j];k!=-1;k=ne[k])
		{
			f[j]|=f[ver[k]];
		}
	}
	for(int i=1;i<=n;i++) cout<<f[i].count()<<endl;
}