hdu 6430 TeaTree 线段树合并

博客围绕树上节点问题展开,已知树上每个节点权值,需计算每个节点以其为LCA的两个节点GCD的最大值即heard值。采用线段树合并方法,为每个节点建线段树,保存该点权值因数,将子节点线段树与父节点合并来维护答案。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题意:一棵树上每个节点权值为v[i],每个节点的heard值是:以它为LCA的两个节点的GCD的最大值,要求输出每个节点的heard值

分析:线段树合并,对每个节点建一棵线段树,第x个位置就表示v[i]有值为x的因子,只要2个节点有相同的因子就维护答案,每次都把所有子节点的线段树与父节点合并。写法是建立n个线段树,每个线段树保存的是该点权值的所有因数,线段树区间[l,r]表示区间l-r范围内最大的因数,首先将1e5个数的所有因数打表,然后再将n个点的权值建线段树。

代码:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
const int maxn=1e5+10;
int rt[maxn],ls[maxn*400],rs[maxn*400],sum[maxn*400];
vector<int>G[maxn],fac[maxn];
int cnt=0,sz=1e5,ans[maxn];
void init(){
	for(int i=1;i<maxn;i++)
        for(int j=i;j<maxn;j+=i)
            fac[j].push_back(i);//将每个数的所有因数保存起来
}

void update(int k,int l,int r,int &o){
	if(!o)o=++cnt;
	if(l==r){
		sum[o]=k;
		return;
	}
	int m=(l+r)/2;
	if(k<=m)update(k,l,m,ls[o]);
	else update(k,m+1,r,rs[o]);
	sum[o]=max(sum[ls[o]],sum[rs[o]]);
}

int merge(int o,int pre,int& v){
	if(!o||!pre)return o|pre;//将父线段树的因数和儿子线段树因数合并到父线段树,
	if(sum[o]==sum[pre])v=max(v,sum[o]);
	if(ls[o]|ls[pre])ls[o]=merge(ls[o],ls[pre],v);//父线段树和儿子线段树左节点都是空的,就不用遍历合并
	if(rs[o]|rs[pre])rs[o]=merge(rs[o],rs[pre],v);
	return o;
}

void dfs(int u){
	ans[u]=-1;
	for(int i=0;i<G[u].size();i++){
		int v=G[u][i];
		dfs(v);
		merge(rt[u],rt[v],ans[u]);
	}
}

int main(){
	int n,v;
	init();
	scanf("%d",&n);
	for(int i=2;i<=n;i++){
		scanf("%d",&v);
		G[v].push_back(i);
	}
	for(int i=1;i<=n;i++){
		scanf("%d",&v);
		for(int j=0;j<fac[v].size();j++)
		update(fac[v][j],1,sz,rt[i]);
	}
	dfs(1);
	for(int i=1;i<=n;i++)
        printf("%d\n",ans[i]);
    return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值