【dfs序+二分】ABC202 E Count Descendants

33 篇文章 0 订阅

E - Count Descendants (atcoder.jp)

题意:

思路:

一开始的想法是,树形DP预处理出在子树u中,离结点u距离为j的结点个数,然后直接输出答案

#include <bits/stdc++.h>

#define int long long

using namespace std;

const int mxn=2e5+10;
const int mxe=2e5+10;

struct ty{
	int to,next;
}edge[mxe<<2];

int N,Q,u,d;
int tot=0;
int p[mxn];
int head[mxn],dep[mxn],dp[mxn][100];

void G_init(){
	tot=0;
	for(int i=0;i<=N;i++){
		head[i]=-1;
	}
}
void add(int u,int v){
	edge[tot].to=v;
	edge[tot].next=head[u];
	head[u]=tot++;
}
void dfs1(int u,int fa){
	dep[u]=dep[fa]+1;
	dp[u][0]=1;
	for(int i=head[u];~i;i=edge[i].next){
		if(edge[i].to==fa) continue;
		dfs1(edge[i].to,u);
	}
}
void dfs2(int u,int fa){
	for(int i=head[u];~i;i=edge[i].next){
		if(edge[i].to==fa) continue;
		dfs2(edge[i].to,u);
		for(int j=1;j<=100;j++){
			dp[u][j]+=dp[edge[i].to][j-1];
		}
	}
}
void solve(){
	cin>>N;
	G_init();
	for(int i=2;i<=N;i++) cin>>p[i];
	for(int i=2;i<=N;i++){
		add(i,p[i]);
		add(p[i],i);
	}
	dfs1(1,0);
	for(int i=1;i<=N;i++) dep[i]--;
	dfs2(1,0);
	cin>>Q;
	for(int i=1;i<=Q;i++){
		cin>>u>>d;
		if(dep[u]>d) cout<<0<<'\n';
		else{
			int delta=d-dep[u];
			cout<<dp[u][delta]<<'\n';
		}
	}
}
signed main(){
	ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
	int __=1;//cin>>__;
	while(__--)solve();return 0;
}

显然,题目没有说它是一棵二叉树,因此树高不止100,不能用树形DP预处理

正解是dfs序

先预处理出每个结点的入序和出序,对于每一次询问,查询在深度d的所有结点中dfs序在[In[u],Out[u]]之间的结点个数即可,所以我们需要预处理出对于每个深度,这个深度的所有结点的dfs序,把它们放进vector就好了

当我们要在子树中找点时,考虑在dfs序中这个序列里面某个区间里面找点

Code:

#include <bits/stdc++.h>

#define int long long

using namespace std;

const int mxn=2e5+10;
const int mxe=2e5+10;

struct ty{
	int to,next;
}edge[mxe<<2];

vector<int> V[mxn];

int N,Q,u,d;
int tot=0,idx=0;
int head[mxn];
int p[mxn],In[mxn],Out[mxn],dep[mxn];

void G_init(){
	tot=0;
	for(int i=0;i<=N;i++){
		head[i]=-1;
	}
}
void add(int u,int v){
	edge[tot].to=v;
	edge[tot].next=head[u];
	head[u]=tot++;
}
void dfs(int u,int fa){
	dep[u]=dep[fa]+1;
	In[u]=++idx;
	V[dep[u]].push_back(In[u]);
	for(int i=head[u];~i;i=edge[i].next){
		if(edge[i].to==fa) continue;
		dfs(edge[i].to,u);
	}
	Out[u]=++idx;
}
void solve(){
	cin>>N;
	G_init();
	for(int i=2;i<=N;i++) cin>>p[i];
	for(int i=2;i<=N;i++){
		add(i,p[i]);
		add(p[i],i);
	}
	dfs(1,0);
	cin>>Q;
	for(int i=1;i<=Q;i++){
		cin>>u>>d;
		d++;
		cout<<lower_bound(V[d].begin(),V[d].end(),Out[u])-lower_bound(V[d].begin(),V[d].end(),In[u])<<'\n';
	}
}
signed main(){
	ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
	int __=1;//cin>>__;
	while(__--)solve();return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值