Acwing164. 可达性统计

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

输入格式

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

输出格式

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

数据范围

1≤N,M≤30000
1≤x,y≤N

输入样例：

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

难度：中等

时/空限制：1s / 256MB

总通过数：7123

总尝试数：12259

来源：《算法竞赛进阶指南》

算法标签

#include<bits/stdc++.h>
using namespace std;
const int N = 3e4 + 10;
int h[N], e[N], ne[N], d[N];
int n, m, idx,q[N];
bitset<N> f[N];
/*
* 将每个点可以到达的点用二进制数表示
* 二进制数位为1即可以到达
*/
void add(int a, int b)
{
	e[idx] = b;
	ne[idx] = h[a];
	h[a] = idx++;
}
void topsort()
{
	int hh = 0, tt = -1;
	for (int i = 1; i <= n; i++)
	{
		if (d[i] == 0)
			q[++tt] = i;
	}
	while (hh <= tt)
	{
		int t = q[hh++];
		for (int i = h[t]; i != -1; i = ne[i])
		{
			int j = e[i];
			if (--d[j] == 0)
				q[++tt] = j;
		}
	}
}
int main()
{
	memset(h, -1, sizeof h);
	cin >> n >> m;
	while (m--)
	{
		int a, b;
		cin >> a >> b;
		add(a, b);
		d[b]++;
	}
	topsort();
	for (int i = n - 1; i >= 0; i--)
	{
		int j = q[i];//获取拓扑序列中的元素
		f[j][j] = 1;
		for (int k = h[j]; k != -1; k = ne[k])
			f[j] |= f[e[k]];//遍历子节点中可以到达的点
	}
	for (int i = 1; i <= n; i++)
	{
		cout << f[i].count() << endl;
	}
	return 0;
}