题外话

(其实模板题没必要在这里水题解的……主要是想说这个qwq)

小常数的快乐.jpg

我也不知道我为啥常数特别小跑得飞快……不加快读就能在 luogu 的最优解上跑到 rank5 (

说不定深夜提交个 fread 版本能上rk3

题目链接

题意

给定一棵带点权的树,边权为1,在线处理m次操作

  1. 0 x k 询问所有与 x 距离不超过 k 的点权和
  2. 1 x y 修改 x 的点权为 y

思路

首先建点分树。对于每个点维护两个树状数组,以距离为下标,权值为内容。第一个维护子树中距离该点为k的权值和,第二个维护父亲的内容(同第一个)。

对于修改操作,暴力爬树高,\(log^2\) 复杂度(树高加上树状数组)。

对于查询操作,一样爬树高,进行容斥(把当前子树 \(k\) 的先加起来,往祖先上爬,如果距离小于 \(k\) ,假设为 \(d\),到祖先上去求一个 \(k-d\),再容斥掉原来这棵子树里被计算过的)。

顺便,把最开始的建树当成 modify 能省很多事qwq

代码

#include <bits/stdc++.h>
#define lowbit(x) ((x)&(-x))
using namespace std;
const int N=1e5+10,inf=1e9+7;
struct edge
{
int to,nxt;
}e[N<<1];
int n,m,head[N],val[N],fa[N][20],dis[N][20],sz[N],f[N],dep[N],rt,siz,tot=0;
bool vis[N];
vector<int> sbit[N],fbit[N]; void add( int u,int v )
{
e[++tot]=(edge){v,head[u]}; head[u]=tot;
e[++tot]=(edge){u,head[v]}; head[v]=tot;
} void get_rt( int x,int fat )
{
sz[x]=1; f[x]=0;
for ( int i=head[x]; i; i=e[i].nxt )
{
int y=e[i].to;
if ( vis[y] || fat==y ) continue;
get_rt( y,x ); sz[x]+=sz[y]; f[x]=max( f[x],sz[y] );
}
f[x]=max( f[x],siz-sz[x] );
if ( f[x]<f[rt] ) rt=x;
} void pre_build( int x,int ancestor,int fat,int d )
{
for ( int i=head[x]; i; i=e[i].nxt )
{
int y=e[i].to;
if ( vis[y] || y==fat ) continue;
fa[y][++dep[y]]=ancestor; dis[y][dep[y]]=d;
pre_build( y,ancestor,x,d+1 );
}
} void build( int x )
{
vis[x]=1; pre_build( x,x,0,1 );
int sav=siz; sbit[x].resize( siz+1 ); fbit[x].resize( siz+1 );
for ( int i=head[x]; i; i=e[i].nxt )
{
int y=e[i].to; if ( vis[y] ) continue;
siz=sz[y]; if ( siz>sz[x] ) siz=sav-sz[x];
rt=0; get_rt( y,x ); build( rt );
}
} int qsum_s( int x,int k )
{
int res=val[x],lim=sbit[x].size()-1; k=min( k,lim );
for ( int i=k; i; i-=lowbit(i) ) res+=sbit[x][i];
return res;
} int qsum_f( int x,int k )
{
int res=0,lim=fbit[x].size()-1; k=min( k,lim );
for ( int i=k; i; i-=lowbit(i) ) res+=fbit[x][i];
return res;
} void modify( int x,int v )
{
int d=dis[x][dep[x]],lim=sbit[x].size()-1;
for ( int j=d; j<=lim && j; j+=lowbit(j) ) fbit[x][j]+=v;
for ( int i=dep[x]; i; i-- )
{
d=dis[x][i]; lim=sbit[fa[x][i]].size()-1;
for ( int j=d; j<=lim; j+=lowbit(j) ) sbit[fa[x][i]][j]+=v;
d=dis[x][i-1];
for ( int j=d; j<=lim && j; j+=lowbit(j) ) fbit[fa[x][i]][j]+=v;
}
} int query( int x,int k )
{
int res=qsum_s( x,k );
for ( int i=dep[x]; i; i-- )
if ( dis[x][i]<=k ) res+=qsum_s( fa[x][i],k-dis[x][i] )-qsum_f( fa[x][i+1],k-dis[x][i] );
return res;
} int main()
{
scanf( "%d%d",&n,&m );
for ( int i=1; i<=n; i++ )
scanf( "%d",&val[i] );
for ( int i=1,u,v; i<n; i++ )
scanf( "%d%d",&u,&v ),add( u,v ); f[0]=inf; siz=n; get_rt( 1,0 ); build( rt );
for ( int i=1; i<=n; i++ )
fa[i][dep[i]+1]=i;
for ( int i=1; i<=n; i++ )
modify( i,val[i] );
int las=0;
while ( m-- )
{
int opt,x,y; scanf( "%d%d%d",&opt,&x,&y );
x^=las; y^=las;
if ( opt==0 ) las=query( x,y ),printf( "%d\n",las );
else modify( x,y-val[x] ),val[x]=y;
}
}

【题解】P6329 【模板】点分树 | 震波的更多相关文章

  1. P6329-[模板]点分树 | 震波

    正题 题目链接:https://www.luogu.com.cn/problem/P6329 解题思路 给出\(n\)个点的一棵树,每个点有权值,有\(m\)次操作 修改一个点\(x\)的权值为\(y ...

  2. 【bzoj3435】[Wc2014]紫荆花之恋 替罪点分树套SBT

    题目描述 强强和萌萌是一对好朋友.有一天他们在外面闲逛,突然看到前方有一棵紫荆树.这已经是紫荆花飞舞的季节了,无数的花瓣以肉眼可见的速度从紫荆树上长了出来.仔细看看的话,这个大树实际上是一个带权树.每 ...

  3. 线段树学习笔记(基础&进阶)(一) | P3372 【模板】线段树 1 题解

    什么是线段树 线段树是一棵二叉树,每个结点存储需维护的信息,一般用于处理区间最值.区间和等问题. 线段树的用处 对编号连续的一些点进行修改或者统计操作,修改和统计的复杂度都是 O(log n). 基础 ...

  4. 题解 P3372 【【模板】线段树1 】(zkw)

    看了一下题解里的zkw线段树,感觉讲的不是很清楚啊(可能有清楚的但是我没翻到,望大佬勿怪). 决定自己写一篇...希望大家能看明白... zkw线段树是一种优秀的非递归线段树,速度比普通线段树快两道三 ...

  5. 题解 P3372 【【模板】线段树1 】

    看了一下题解里的zkw线段树,感觉讲的不是很清楚啊(可能有清楚的但是我没翻到,望大佬勿怪). 决定自己写一篇...希望大家能看明白... zkw线段树是一种优秀的非递归线段树,速度比普通线段树快两道三 ...

  6. 洛谷 P3373 【模板】线段树 2 题解

    P3373 [模板]线段树 2 题目描述 如题,已知一个数列,你需要进行下面三种操作: 1.将某区间每一个数乘上x 2.将某区间每一个数加上x 3.求出某区间每一个数的和 输入格式 第一行包含三个整数 ...

  7. 洛谷P3373 【模板】线段树 2 && P2023 [AHOI2009]维护序列——题解

    题目传送: P3373 [模板]线段树 2  P2023 [AHOI2009]维护序列 该题较传统线段树模板相比多了一个区间乘的操作.一提到线段树的区间维护问题,就自然想到了“懒标记”:为了降低时间复 ...

  8. 一篇自己都看不懂的点分治&点分树学习笔记

    淀粉质点分治可真是个好东西 Part A.点分治 众所周知,树上分治算法有$3$种:点分治.边分治.链分治(最后一个似乎就是树链剖分),它们名字的不同是由于分治方式的不同的.点分治,顾名思义,每一次选 ...

  9. 【洛谷3345_BZOJ3924】[ZJOI2015]幻想乡战略游戏(点分树)

    大概有整整一个月没更博客了 -- 4 月为省选爆肝了一个月,最后压线进 B 队,也算给 NOIP2018 翻车到 316 分压线省一这个折磨了五个月的 debuff 画上了一个不算太差的句号.结果省选 ...

随机推荐

  1. uboot——do_bootm

    do_bootm |----------根据参数得到 image的起始地址 |----------比较header的 magic_num 是否为 zImage | |是 | | zImage路线 | ...

  2. flink1.10版local模式提交job流程分析

    1.WordCount程序实例 2.本地监听9000端口后测试结果 3.job提交流程 4.local模式执行StreamGraph任务 5.流程分析 flink job提交流程个人理解可以大致分为定 ...

  3. 对图片进行Base64转码和解码

    Base64代码 base64.c #include <stdlib.h> #include <stdio.h> #include <string.h> #incl ...

  4. Collectors工具类

    Collector是专门用来作为Stream的collect方法的参数的:而Collectors是作为生产具体Collector的工具类. Collectors是一个工具类,是JDK预实现Collec ...

  5. JavaScript复习大纲

    1. HTML.CSS和JavaScript各自在网页设计中的作用. 1.HTML生成结构. 2.CSS样式美化. 3.JavaScript的作用: (1) 操作HTML及CSS,让网页具有动态行为. ...

  6. BurpSuite多重代理的情形

    有时候为了隐藏自己的真实IP,访问目标网站需要走代理.比如说,通过代理IP访问目标网站并且代理IP可以随时切换,这样可以避免IP被封堵后无法访问目标网站的尴尬. 一. 首先把代理切换为全局模式 操作步 ...

  7. JVM字节码执行引擎

    一.概述 在不同的虚拟机实现里面,执行引擎在执行Java代码的时候可能会有解释执行(通过解释器执行)和编译器执行(通过即时编译器产生本地代码执行)两种选择,所有的Java虚拟机的执行引擎都是一致的:输 ...

  8. 总是说spring难学?来看完这些spring的注解及其解释,真香!

    前言 用过spring的人都知道,spring简单的通过注解就可以完成很多事情,但这些东西是如何实现的呢以及如何应用到我们自己的代码中?接下来,让我们一起开启注解的旅程. 1. @Controller ...

  9. Vegas让人物回眸更有韵味的方法分享

    "回眸一笑百媚生,六宫粉黛无颜色",是白居易在<长恨歌>中描述杨贵妃美貌的名句,这一句运用夸张的手法,反映了杨贵妃回眸时的娇媚横生,百般娇媚. 接下来,小编就教你用视频 ...

  10. jQuery 第六章 实例方法 动画

    .show() .hide() .toggle() .fadeIn() .fadeout() .fadeToggle() .fadeTo() .slideDown() .slideUp() .slid ...