codeforces1307F 思维

探讨了在无向图中,存在r个关键点,每次能走k步的情况。利用并查集与多源BFS策略,解决两点u和v之间是否可达的问题。通过拆点优化,将边转化为u->n+i->v形式,再进行BFS。对于询问,若两点距离超过2k,则判断它们各自向LCA走k步后是否在同一个集合内。
摘要由CSDN通过智能技术生成

无向图上有r个关键点,一个人每次只能连续走k步,多次询问u和v两点能否到达.

做法

如果两点间距离<=k,就用并查集合并,所以考虑对r个关键点进行多源bfs,但是注意到这样会使得任意两个距离<=2k的点都可达,考虑拆点,将每条边拆成u->n+i->v,然后多源bfs k步即可.
对于每个询问,如果u和v距离<=2
k显然可以,否则两个点各自向对方走k步(u->lca->v),此时两个点肯定不会相遇(因为距离>2*k),最后看两个点是否在一个集合中即可.

#include <bits/stdc++.h>
#define mp make_pair
#define pb push_back
#define sz(x) (int)x.size()
#define all(x) begin(x), end(x)
#define fi first
#define se second
#define debug(x) cerr << #x << " " << x << '\n'
using namespace std;
using ll = long long;
using pii = pair<int,int>;
using pli = pair<ll,int>;
const int INF = 0x3f3f3f3f, N = 4e5 + 5;
const ll LINF = 1e18 + 5;
constexpr int mod = 1e9 + 7;
constexpr int LOG = 19;
int n, k, r, sp[N], step[N], f[N];
vector <int> G[N];
int fa[N][22], dep[N]; 
void dfs(int u)
{	
	for(int i=1; i<=LOG; i++) 
		fa[u][i] = fa[fa[u][i-1]][i-1];
	for(int v : G[u])
	{
		if(v==fa[u][0]) continue;
		dep[v] = dep[u] + 1;
		fa[v][0] = u; 
		dfs(v);
	}
}
int LCA(int x,int y) 
{
	if(dep[x]>dep[y]) swap(x,y);
	for(int i=LOG; i>=0; i--)
		if(dep[fa[y][i]]>=dep[x]) y = fa[y][i];
	if(x==y) return x;
	for(int i=LOG; i>=0; i--)
		if(fa[x][i]!=fa[y][i]) x = fa[x][i], y = fa[y][i]; 
	return fa[x][0]; 
}
int up(int x,int d)
{
	for(int i=LOG; i>=0; i--)
		if((d>>i)&1) x = fa[x][i];
	return x;
}
int find(int x)
{
	if(x==f[x]) return x;
	return f[x] = find(f[x]);
}
void merge(int x, int y)
{
	x = find(x), y = find(y);
	if(x!=y) f[x] = y;
}
void bfs()
{
	queue <int> q;
	for(int i=1; i<=2*n; i++) f[i] = i, step[i] = -1;
	for(int i=1; i<=r; i++) q.push(sp[i]), step[sp[i]] = 0;
	while(sz(q))
	{
		int u = q.front(); q.pop();
		if(step[u]>=k) break;
		for(int v : G[u])
		{
			merge(u, v);
			if(step[v]==-1)
			{
				step[v] = step[u] + 1;
				q.push(v);
			}
		}
	}
}
void walk(int &u, int v, int w, int k)
{
	if(dep[u]-dep[w]>=k) u = up(u, k);
	else u = up(v, dep[v]-k+dep[u]-2*dep[w]);
}
bool ok(int u, int v)
{
	int w = LCA(u, v);
	if(dep[u]+dep[v]-2*dep[w]<=2*k) return 1;
	walk(u, v, w, k); walk(v, u, w, k);
	if(find(u)==find(v)) return 1;
	return 0;
}
int main()
{
	scanf("%d%d%d", &n, &k, &r);
	for(int i=1; i<n; i++)
	{
		int u, v; scanf("%d%d", &u, &v);
		G[u].pb(n+i); G[v].pb(n+i);
		G[n+i].pb(u); G[n+i].pb(v);
	}
	for(int i=1; i<=r; i++) scanf("%d", &sp[i]);
	dep[1] = 1; dfs(1);
	bfs();
	int q; scanf("%d", &q);
	while(q--)
	{
		int u, v; scanf("%d%d", &u, &v);
		if(ok(u, v)) puts("YES");
		else puts("NO");
	}
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值