题意:一棵树,每个点都有自己val(1 <= val <= 1e5),而任意两个点u,v可以对lca(u,v) 产生gcd(valu,valv)的贡献,求每个点能接受到来自子树贡献的最大值。

分析:一个数w和其整数数倍的数gcd值还是w。记录每个值对应树的下标,枚举1- max{val}中的数,遍历其整数倍中两两结点的lca,检查其是否需要被更新。

其中遍历LCA的操作并不用二重循环枚举顶点,可以将顶点按dfs序排序后遍历其中相邻两点的LCA。

这里使用了RMQ的方式O(1) 查询LCA。

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=1e5+;
struct node{
int v,next;
}edges[N<<];
int head[N],e;
int id[N];
int RMQ[N*][];
int curID;
int F[N*],B[N*];
int n,m,Q,root;
int val[N];
int res[N];
vector<int> num[N];
int mx; void init(){
e = ; curID = ;
memset(head,-,sizeof(head));
for(int i=;i<N;++i) num[i].clear();
} void Add (int u,int v)
{
edges[e].v=v;
edges[e].next=head[u];
head[u]=e++;
} void DFS (int u,int p,int Dep)
{
int i,v;
curID++;
F[curID]=u;
B[curID]=Dep;
id[u]=curID;
for (i=head[u];i!=-;i=edges[i].next){
v=edges[i].v;
if (v==p) continue;
DFS(v,u,Dep+);
curID++;
F[curID]=u;
B[curID]=Dep;
}
} void initRMQ()
{
int i,j,x,y;
for (i=;i<=curID;i++)
RMQ[i][]=i;
for (j=;(<<j)<=curID;j++)
for (i=;i+(<<j)-<=curID;i++){
x=RMQ[i][j-];
y=RMQ[i+(<<(j-))][j-];
RMQ[i][j]=B[x]<B[y]?x:y;
}
} int getLCA (int a,int b)
{
int k,x,y;
a=id[a];b=id[b];
if (a>b) k=a,a=b,b=k;
k = - __builtin_clz(b - a + );
x=RMQ[a][k];
y=RMQ[b-(<<k)+][k];
return B[x]<B[y]?F[x]:F[y];
} void gao(int w)
{
vector<int> p;
for(int i=w;i<=mx;i+=w){
for(auto &v:num[i]){
p.push_back(v);
}
}
sort(p.begin(),p.end(),[&](int x,int y){return id[x]<id[y];});
for(int i=;i<p.size();++i){
int lca= getLCA(p[i],p[i-]);
res[lca] = max(res[lca],w);
}
} int main()
{
#ifndef ONLINE_JUDGE
freopen("in.txt","r",stdin);
freopen("out.txt","w",stdout);
#endif
int N;
while(~scanf("%d",&N)){
init();
memset(res,-,sizeof(res));
int u,v,tmp;
for(int i=;i<=N;++i){
scanf("%d",&u);
Add(u,i);
Add(i,u);
}
DFS(,-,);
initRMQ();
mx = -;
for(int i=;i<=N;++i){
scanf("%d",&val[i]);
mx = val[i]>mx? val[i] : mx;
num[val[i]].push_back(i);
}
for(int i=mx;i>=;--i)
gao(i);
for(int i=;i<=N;++i)
printf("%d\n",res[i]);
}
return ;
}

HDU - 6430 Problem E. TeaTree 2018 Multi-University Training Contest 10 (LCA+枚举因子)的更多相关文章

  1. HDU 6430 Problem E. TeaTree(虚树)

    Problem E. TeaTree Problem Description Recently, TeaTree acquire new knoledge gcd (Greatest Common D ...

  2. HDU - 6435 Problem J. CSGO 2018 Multi-University Training Contest 10 (二进制枚举+思维)

    题意:有N个主武器(MW)和M个副武器(SW),每个武器都有自己的S值,和K个附加属性xi.要选取一对主副武器搭配,搭配后获得的性能由该公式得出: 求获得最大的性能为多少. 分析:由于|xm - xs ...

  3. 2018 Nowcoder Multi-University Training Contest 10

    Practice Link J. Rikka with Nickname 题意: 给出\(n\)个字符串,要求依次合并两个串\(s, t\),满足将\(t\)合并到\(s\)中变成\(r\),使得\( ...

  4. hdu 5416 CRB and Tree(2015 Multi-University Training Contest 10)

    CRB and Tree                                                             Time Limit: 8000/4000 MS (J ...

  5. 2015 Multi-University Training Contest 10 hdu 5406 CRB and Apple

    CRB and Apple Time Limit: 12000/6000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)To ...

  6. 2015 Multi-University Training Contest 10 hdu 5412 CRB and Queries

    CRB and Queries Time Limit: 12000/6000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Other ...

  7. 2018 Multi-University Training Contest 10 Solution

    A - Problem A.Alkane 留坑. B - Problem B. Beads 留坑. C - Problem C. Calculate 留坑. D - Problem D. Permut ...

  8. HDU 6141 - I am your Father! | 2017 Multi-University Training Contest 8

    思路来自 FXXL 最小树形图模板用kuangbin的 /* HDU 6141 - I am your Father! [ 最小树形图 ] | 2017 Multi-University Traini ...

  9. HDU - 6437 Problem L.Videos 2018 Multi-University Training Contest 10 (最小费用最大流)

    题意:M个影片,其属性有开始时间S,结束时间T,类型op和权值val.有K个人,每个人可以看若干个时间不相交的影片,其获得的收益是这个影片的权值val,但如果观看的影片相邻为相同的属性,那么收益要减少 ...

随机推荐

  1. iframe宽高百分百显示

    <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/ ...

  2. Python笔记8:网络编程

    python内置封装了很多常见的网络协议的库,因此python成为了一个强大的网络编程工具,这里是对python的网络方面编程的一个简单描述. urllib 和 urllib2模块 urllib 和u ...

  3. js计算百分比

    function GetPercent(num, total) { /// <summary> /// 求百分比 /// </summary> /// <param na ...

  4. Linux文件的打包与压缩

    打包命令: tar tar 的选项与参数非常的多!我们只讲几个常用的选项,更多选项您可以自行 man tar 查询罗! [root@www ~]# tar [-j|-z] [cv] [-f 创建的档名 ...

  5. jvm之 国际酒店 一次报表 load数据死循环导致的FULLGC

    从目前这几天的运行情况看,每天晚上0点到3点,应该是你们那边有大数据任务处理,在这个时间段,SOA会频繁报ZK断开重联的告警,见下图.   类似的问题,在CRM系统上也有出现过,我们dump内存,分析 ...

  6. C# Expression 树转化为SQL语句(一)

    sql有有四中基本语句,分别是增删改查,在建立model后如何生成这四中sql语句,降低开发时间. 我们先模拟出一张学生表: public class Student { public int id ...

  7. poj 2923(状态压缩+背包)

    比较巧妙的一道题目,拿到题目就想用暴力直接搜索,仔细分析了下发现复杂度达到了2^n*n! ,明显不行,于是只好往背包上想. 于是又想二分找次数判断可行的方法,但是发现复杂度10^8还是很悬... 然后 ...

  8. springboot2.1.1 中集成websocket 单元测试异常

    单元测试在没有集成websocket之前是好好的,当集成websocket之后就出现了下面的异常(只贴出来关键信息): 2019-01-11 10:05:42 [ERROR] [org.springf ...

  9. git 四个基本对象、分支、三个存储区、reset-revert-变基、cherry-pick

    1:git四个基本对象     2:工作区.缓存去.历史区   3:Git 分支介绍 https://blog.csdn.net/wh_19910525/article/details/7470964 ...

  10. grunt 相关知识

    /** * Created by lee on 2014.07.02 002. */ module.exports = function (grunt) { // 项目配置 grunt.initCon ...