loj2322 「清华集训 2017」Hello world!
先吐槽一下,sb数据毁我青春败我前程。
首先,一个数开根开不了多少次。
当我们把它开到1的时候,我们以后就不需要开他了,我们可以利用并查集跳过他,这是套路。
但是这个每次走$k$步,让人很头痛。
于是乎……分块
首先,对于$k$比较大的情况,我们可以暴力跳。否则,
我们对于每个$k$建一棵树,第$i$棵树上$x$的父亲是原树上$x$的第$i$个祖先
然后树剖,用线段树或者树状数组维护一下就好了。
然后!!打开榜,rk1,啥子,并查集+暴力跳?!
诶,那就暴力跳吧。
//Serene
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
#define db double
#define For(i,a,b) for(int i=(a);i<=(b);++i)
#define Rep(i,a,b) for(int i=(a);i>=(b);--i)
const int maxn=1e5+7,maxs=200+7,maxt=23,W=19;
ll n,m,a[maxn],s; char cc;ll ff;
template<typename T>void read(T& aa) {
aa=0;cc=getchar();ff=1;
while((cc<'0'||cc>'9')&&cc!='-') cc=getchar();
if(cc=='-') ff=-1,cc=getchar();
while(cc>='0'&&cc<='9') aa=aa*10+cc-'0',cc=getchar();
aa*=ff;
} int fir[maxn],nxt[maxn],to[maxn],e=0;
void add(int x,int y) {
to[++e]=y;nxt[e]=fir[x];fir[x]=e;
to[++e]=x;nxt[e]=fir[y];fir[y]=e;
} int f[maxn];
int find(int x) {return x==f[x]? x:f[x]=find(f[x]);} int fa[maxn][maxt],jp[maxn][maxs],dep[maxn];
void dfs(int pos,int f) {
int y,z; dep[pos]=dep[f]+1;
fa[pos][0]=jp[pos][1]=f; jp[pos][0]=pos;
For(i,1,W) fa[pos][i]=fa[fa[pos][i-1]][i-1];
For(i,2,s) jp[pos][i]=jp[f][i-1];
for(y=fir[pos];y;y=nxt[y]) {
if((z=to[y])==f) continue;
dfs(z,pos);
}
} inline int get_fa(int x,int k) {
if(k<=s) return jp[x][k];
Rep(i,W,0) if(k>=(1<<i)) {x=fa[x][i]; k-=(1<<i);}
return x;
} int get_lca(int x,int y) {
if(dep[x]!=dep[y]) {
if(dep[x]<dep[y]) swap(x,y);
Rep(i,W,0) if(dep[fa[x][i]]>=dep[y]) x=fa[x][i];
}
if(x==y) return x;
Rep(i,W,0) if(fa[x][i]!=fa[y][i]) {x=fa[x][i]; y=fa[y][i];}
return fa[x][0];
} void chge(int x) {
if(a[x]==1) return;
a[x]=floor(sqrt(a[x]));
if(a[x]==1) f[x]=find(fa[x][0]);
} int q(int x,int y,int lca,int k) {
if(dep[y]-dep[lca]>=k) return get_fa(y,k);
k-=dep[y]-dep[lca]; y=lca;
return get_fa(x,dep[x]-dep[y]-k);
} inline int lst(int x,int k) {
if(k>s) return get_fa(x,k);
int y=find(fa[x][0]);
int p=(dep[x]-dep[y])%k;
if(p) p=k-p;
return jp[y][p];
} void get_chge(int x,int y,int k) {
int lca=get_lca(x,y),len=dep[x]+dep[y]-2*dep[lca];
if(len%k) chge(y),y=q(x,y,lca,len%k),lca=get_lca(x,y);
while(dep[x]>=dep[lca]) {chge(x); x=lst(x,k);}
while(dep[y]>dep[lca]) {chge(y); y=lst(y,k);}
} ll Yth(int x,int y,int k) {
int lca=get_lca(x,y),len=dep[x]+dep[y]-2*dep[lca]; ll rs=0;
if(len%k) rs+=a[y],y=q(x,y,lca,len%k),lca=get_lca(x,y);
rs+=(dep[x]+dep[y]-2*dep[lca])/k+1;
while(dep[x]>=dep[lca]) {rs+=a[x]-1; x=lst(x,k);}
while(dep[y]>dep[lca]) {rs+=a[y]-1; y=lst(y,k);}
return rs;
} int main() {
read(n); int op,x,y,k;
s=min(200,(int)sqrt(n)+1);
For(i,1,n) read(a[i]),f[i]=i;
For(i,1,n-1) {
read(x); read(y);
add(x,y);
}
dfs(1,0); read(m);
For(i,1,n) if(a[i]==1) f[i]=fa[i][0];
For(i,1,m) {
read(op); read(x); read(y); read(k);
if(op==0) get_chge(x,y,k);
else printf("%lld\n",Yth(x,y,k));
}
return 0;
}
loj2322 「清华集训 2017」Hello world!的更多相关文章
- Loj #2331. 「清华集训 2017」某位歌姬的故事
Loj #2331. 「清华集训 2017」某位歌姬的故事 IA 是一名会唱歌的女孩子. IOI2018 就要来了,IA 决定给参赛选手们写一首歌,以表达美好的祝愿.这首歌一共有 \(n\) 个音符, ...
- Loj #2324. 「清华集训 2017」小 Y 和二叉树
Loj #2324. 「清华集训 2017」小 Y 和二叉树 小Y是一个心灵手巧的OIer,她有许多二叉树模型. 小Y的二叉树模型中,每个结点都具有一个编号,小Y把她最喜欢的一个二叉树模型挂在了墙上, ...
- Loj #2321. 「清华集训 2017」无限之环
Loj #2321. 「清华集训 2017」无限之环 曾经有一款流行的游戏,叫做 *Infinity Loop***,先来简单的介绍一下这个游戏: 游戏在一个 \(n \times m\) 的网格状棋 ...
- Loj 2320.「清华集训 2017」生成树计数
Loj 2320.「清华集训 2017」生成树计数 题目描述 在一个 \(s\) 个点的图中,存在 \(s-n\) 条边,使图中形成了 \(n\) 个连通块,第 \(i\) 个连通块中有 \(a_i\ ...
- loj #2325. 「清华集训 2017」小Y和恐怖的奴隶主
#2325. 「清华集训 2017」小Y和恐怖的奴隶主 内存限制:256 MiB时间限制:2000 ms标准输入输出 题目类型:传统评测方式:文本比较 题目描述 "A fight? Co ...
- [LOJ#2330]「清华集训 2017」榕树之心
[LOJ#2330]「清华集训 2017」榕树之心 试题描述 深秋.冷风吹散了最后一丝夏日的暑气,也吹落了榕树脚下灌木丛的叶子.相识数年的Evan和Lyra再次回到了小时候见面的茂盛榕树之下.小溪依旧 ...
- [LOJ#2329]「清华集训 2017」我的生命已如风中残烛
[LOJ#2329]「清华集训 2017」我的生命已如风中残烛 试题描述 九条可怜是一个贪玩的女孩子. 这天她在一堵墙钉了 \(n\) 个钉子,第 \(i\) 个钉子的坐标是 \((x_i,y_i)\ ...
- [LOJ#2328]「清华集训 2017」避难所
[LOJ#2328]「清华集训 2017」避难所 试题描述 "B君啊,你当年的伙伴都不在北京了,为什么你还在北京呢?" "大概是因为出了一些事故吧,否则这道题就不叫避难所 ...
- [LOJ#2327]「清华集训 2017」福若格斯
[LOJ#2327]「清华集训 2017」福若格斯 试题描述 小d是4xx9小游戏高手. 有一天,小d发现了一个很经典的小游戏:跳青蛙. 游戏在一个 \(5\) 个格子的棋盘上进行.在游戏的一开始,最 ...
随机推荐
- 2018-8-10-win10-uwp-如何判断一个对象被移除
title author date CreateTime categories win10 uwp 如何判断一个对象被移除 lindexi 2018-08-10 19:16:50 +0800 2018 ...
- 宽域POST提交数据
小数据宽域可以使用jsonp,但是大数据跨域必须post那么有以下2种方式 1,传统方式 动态生成form var url = ''var $iframe = $("<iframe s ...
- LUOGU P1654 OSU! (概率期望)
传送门 解题思路 首先考虑对于一个点来说,如果这个点是1的话,那么对于答案来说 $(ans+1)^3=ans^3+3*ans^2+3*ans+1$,这对于上一个答案来说其实贡献了 $3*ans^2+3 ...
- Servlet和模本办法
抽象类HttpServlet,提供7个子方法的默认实现,当我们手写Servlet只需要继承HttpServlet并替换相应的do方法来完成自己的业务逻辑即可. 我们一般的实现方法: public My ...
- if __name__=='__main__'使用场景,彻底明白
本博中有一篇文章写了 if __name__=='__main__'的作用与原理http://www.cnblogs.com/fennudexiaoniao/p/7458324.html,但是好像似懂 ...
- 04_Hibernate检索方式
一.Hibernate检索方式概述 OID检索方式:按照对象的OID来检索对象(get/load) HQL检索方式:使用面向对象的HQL查询语言 QBC检索方式:使用QBC(Query By Crit ...
- selenium 三种断言以及异常类型
selenium 提供了三种模式的断言:assert .verify.waitfor Assert 失败时,该测试将终止. Verify 失败时,该测试将继续执行,并将错误记入日显示屏 .也就是说允许 ...
- Activiti流程实例管理
1.启动流程 在完成了流程定义部署后,就要启动流程实例了. /** * 1 启动流程 * 当流程到达一个节点时,会在act_ru_execution表中产生1条数据 * 如果当前节点是用户任务节点,这 ...
- 2019-2-18-VisualStudio-给项目添加特殊的-Nuget-的链接
title author date CreateTime categories VisualStudio 给项目添加特殊的 Nuget 的链接 lindexi 2019-02-18 15:56:48 ...
- HZOI20190823 C magic
数论板子合集... 我们要求: $N^{\sum\limits_{i=1}^{N}[gcd(i,N)==1]C_{n}^{i}}mod p$ 其中p为54184622,是个合数 指数是组合数,不能用快 ...