CH 2101 可达性统计(进阶指南, 位运算,拓扑排序)

算法竞赛进阶指南, 99 页

本题要点:
1、有向无环图,拓扑排序之后,得到拓扑序列 top_order[MaxN]。
top_order[i] 可以到达的点 假设为 f(top_order[i]),等于 top_order[i]所有相邻的点 (假设其中任意一点为j)
可以到达的点f(j)的并集。因此,从后往前扫描拓扑序列,依次求并集。

2、 用c++的 bitset 模板来处理位运算。参考 https://www.cnblogs.com/magisk/p/8809922.html
构造函数 : bit.set(i), 表示第i位置1;
统计1的个数:bit.count()

#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#include <bitset>
using namespace std;
const int MaxN = 30010;
const int MaxM = 30010;
bitset<MaxN> bit[MaxN];
int top_order[MaxN];
int top_cnt = 0;
int n, m, tot;
int ver[MaxN], head[MaxN], ne[MaxN], deg[MaxN];

void add(int x, int y)
{
	ver[++tot] = y, ne[tot] = head[x], head[x] = tot;
	deg[y]++;
}

void topsort()
{
	queue<int> q;
	for(int i = 1; i <= n; ++i)	//所有点的号码 1 ~ n
	{
		if(0 == deg[i])
		{
			q.push(i);
		}
	}
	while(q.size())
	{
		int x = q.front();
		q.pop();
		top_order[top_cnt++] = x;
		for(int i = head[x]; i; i = ne[i])
		{
			int y = ver[i];
			if(--deg[y] == 0)
			{
				q.push(y);
			}
		}
	}
}

void solve()
{
	for(int i = 1; i <= n; ++i)
	{
		bit[i].set(i);	
	}
	for(int i = top_cnt - 1; i >= 0; --i)
	{
		int x = top_order[i];
		for(int j = head[x]; j; j = ne[j])	
		{
			int y = ver[j];	
			bit[x] |= bit[y];
		}
	}
	for(int i = 1; i <= n; ++i)
	{
		printf("%d\n", (int)(bit[i].count()));
	}
}

int main()
{
	int a, b;
	scanf("%d%d", &n, &m);
	for(int i = 1; i <= m; ++i)
	{
		scanf("%d%d", &a, &b);
		add(a, b);
	}
	topsort();
	solve();
	return 0;
}

/*
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
*/


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值