[xsy2123]毛毛虫
题意:有一棵带点权的树,链修改是把$(x,y)$这条链和与其相邻的节点都加上一个数,查询是问$(x,y)$这条链和与其相邻的节点的权值和
学到了一个新姿势?
考虑树链剖分,在剖重链时每次给当前节点的儿子顺次标号,那么一条重链上所有节点的儿子的标号是连续的,也就对应着题目的限制,再补上漏掉的$lca_{x,y}$和$fa_{lca_{x,y}}$即可
#include<stdio.h>
typedef long long ll;
void swap(int&a,int&b){a^=b^=a^=b;}
int h[100010],nex[200010],to[200010],v[100010],M;
void add(int a,int b){
M++;
to[M]=b;
nex[M]=h[a];
h[a]=M;
}
int siz[100010],dep[100010],fa[100010],son[100010],n;
void dfs(int x){
int i,k=0;
siz[x]=1;
dep[x]=dep[fa[x]]+1;
for(i=h[x];i;i=nex[i]){
if(to[i]!=fa[x]){
fa[to[i]]=x;
dfs(to[i]);
siz[x]+=siz[to[i]];
if(siz[to[i]]>siz[k])k=to[i];
}
}
son[x]=k;
}
int bl[100010],p[100010],pos[100010],L[100010],R[100010];
void dfs(int x,int chain){
int i;
bl[x]=chain;
L[x]=M+1;
for(i=h[x];i;i=nex[i]){
if(to[i]!=fa[x]){
p[++M]=to[i];
pos[to[i]]=M;
}
}
R[x]=M;
if(son[x])dfs(son[x],chain);
for(i=h[x];i;i=nex[i]){
if(to[i]!=fa[x]&&to[i]!=son[x])dfs(to[i],to[i]);
}
}
ll s[400010],d[400010];
void pushup(int x){s[x]=s[x<<1]+s[x<<1|1];}
void build(int l,int r,int x){
if(l==r){
s[x]=v[p[l]];
return;
}
int mid=(l+r)>>1;
build(l,mid,x<<1);
build(mid+1,r,x<<1|1);
pushup(x);
}
void ad(int x,int len,ll v){
d[x]+=v;
s[x]+=len*v;
}
void pushdown(int x,int ln,int rn){
if(d[x]){
ad(x<<1,ln,d[x]);
ad(x<<1|1,rn,d[x]);
d[x]=0;
}
}
void modify(int L,int R,int v,int l,int r,int x){
if(L>R)return;
if(L<=l&&r<=R)return ad(x,r-l+1,v);
int mid=(l+r)>>1;
pushdown(x,mid-l+1,r-mid);
if(L<=mid)modify(L,R,v,l,mid,x<<1);
if(mid<R)modify(L,R,v,mid+1,r,x<<1|1);
pushup(x);
}
ll query(int L,int R,int l,int r,int x){
if(L>R)return 0;
if(L<=l&&r<=R)return s[x];
int mid=(l+r)>>1;
ll s=0;
pushdown(x,mid-l+1,r-mid);
if(L<=mid)s+=query(L,R,l,mid,x<<1);
if(mid<R)s+=query(L,R,mid+1,r,x<<1|1);
return s;
}
int lca(int x,int y){
while(bl[x]!=bl[y]){
if(dep[bl[x]]<dep[bl[y]])swap(x,y);
x=fa[bl[x]];
}
return dep[x]<dep[y]?x:y;
}
void modify(int x,int y,int v){
int u=lca(x,y);
while(bl[x]!=bl[u]){
modify(L[bl[x]],R[x],v,1,n,1);
x=fa[bl[x]];
}
while(bl[y]!=bl[u]){
modify(L[bl[y]],R[y],v,1,n,1);
y=fa[bl[y]];
}
if(dep[x]>dep[y])swap(x,y);
modify(L[x],R[y],v,1,n,1);
modify(pos[x],pos[x],v,1,n,1);
if(fa[x])modify(pos[fa[x]],pos[fa[x]],v,1,n,1);
}
ll query(int x,int y){
int u=lca(x,y);
ll s=0;
while(bl[x]!=bl[u]){
s+=query(L[bl[x]],R[x],1,n,1);
x=fa[bl[x]];
}
while(bl[y]!=bl[u]){
s+=query(L[bl[y]],R[y],1,n,1);
y=fa[bl[y]];
}
if(dep[x]>dep[y])swap(x,y);
s+=query(L[x],R[y],1,n,1)+query(pos[x],pos[x],1,n,1);
if(fa[x])s+=query(pos[fa[x]],pos[fa[x]],1,n,1);
return s;
}
int main(){
int m,i,x,y,z;
scanf("%d",&n);
for(i=1;i<n;i++){
scanf("%d%d",&x,&y);
add(x,y);
add(y,x);
}
for(i=1;i<=n;i++)scanf("%d",v+i);
dfs(1);
M=1;
p[1]=1;
pos[1]=1;
dfs(1,1);
build(1,n,1);
scanf("%d",&m);
while(m--){
scanf("%d%d%d",&i,&x,&y);
if(i==1){
scanf("%d",&z);
modify(x,y,z);
}else
printf("%lld\n",query(x,y));
}
}
[xsy2123]毛毛虫的更多相关文章
- 【HAOI2009】【P1307】毛毛虫
感觉相比其他树归题简单多了,不过有点绕(也许是我的思路很奇怪一.一)(这是省选题啊,就算作为T1这题也太水了,HA好弱……) 原题: 对于一棵树,我们可以将某条链和与该链相连的边抽出来,看上去就象成一 ...
- [haoi2009]毛毛虫 树形dp
这道题细节处理不少,但要AC不难: 设以i节点为根节点的子树能形成的最大的毛毛虫长度为f[i],则f[i]=max(f[j])+i节点的孩子数: 答案需要f最大和次大的两个子树合并,而且若合并的位置不 ...
- [HAOI2009]毛毛虫
题目描述 对于一棵树,我们可以将某条链和与该链相连的边抽出来,看上去就象成一个毛毛虫,点数越多,毛毛虫就越大.例如下图左边的树(图 1 )抽出一部分就变成了右边的一个毛毛虫了(图 2 ). 输入输出格 ...
- 【HAOI2009】毛毛虫
题面 题目描述 对于一棵树,我们可以将某条链和与该链相连的边抽出来,看上去就象成一个毛毛虫,点数越多,毛毛虫就越大.例如下图左边的树(图 1 )抽出一部分就变成了右边的一个毛毛虫了(图 2 ). 输入 ...
- 【luogu3174】【HAOI2009】毛毛虫
Description 对于一棵树,我们可以将某条链和与该链相连的边抽出来,看上去就象成一个毛毛虫,点数越多,毛毛虫就越大. Input 在文本文件 worm.in 中第一行两个整数 N , M ,分 ...
- P3174 [HAOI2009]毛毛虫
题目描述 对于一棵树,我们可以将某条链和与该链相连的边抽出来,看上去就象成一个毛毛虫,点数越多,毛毛虫就越大.例如下图左边的树(图 1 )抽出一部分就变成了右边的一个毛毛虫了(图 2 ). 输入输出格 ...
- loj #6570. 毛毛虫计数
$ \color{#0066ff}{ 题目描述 }$ hsezoi 巨佬 olinr 喜欢 van 毛毛虫,他定义毛毛虫是一棵树,满足树上存在一条树链,使得树上所有点到这条树链的距离最多为 1. 给定 ...
- loj6570 毛毛虫计数(生成函数FFT)
link 巨佬olinr的题解 <-- olinr很强 考虑生成函数 考虑直径上点数>=4的毛毛虫的直径,考虑直径中间那些节点以及他上面挂的那些点的EGF \(A(x)=\sum_{i\g ...
- [洛谷P3174][HAOI2009]毛毛虫
题目大意:给一棵树,求其中最大的“毛毛虫”,毛毛虫的定义是一条链上分出几条边 题解:把每个点的权值定义为它的度数减一,跑带权直径即可,最后答案加二 卡点:无 C++ Code: #include &l ...
随机推荐
- html+js+node实现五子棋线上对战,五子棋最简易算法
首先附上我的github地址,https://github.com/jiangzhenfei/five,线上实例:http://47.93.103.19:5900/client/ 线上实例,你可以随意 ...
- 【HNOI】 期望面积
[题目描述]给定n个点,求这n个点组成凸包的期望面积.保证任意三点不共线. [数据范围]n<=100. 首先我们知道凸包面积的计算为所有在凸包上相邻的点的叉积和,那么我们可以枚举两个点,然后求出 ...
- php中的parse_ini_file函数
作用:parse_ini_file() 函数解析一个配置文件,并以数组的形式返回其中的设置 格式:parse_ini_file(file,true)// (第二个参数为可选参数.如果设置为 true, ...
- FindQQByProcess
看网上有许多通过进程寻找QQ号的例子,看了一下,里面涉及的知识点还是比较多,但网上的兼容性不太好,而且没有给出匹配字符的来源,所以自己动手写了一下,顺便给出一些我调试的结果. #include &qu ...
- Python3 xml模块的增删改查
xml数据示例 ? 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 <data> < ...
- wait与waitpid
1. 函数原型: #include <sys/wait.h> pid_t wait(int *statloc); pid_t waitpid(pid_t pid, int *statloc ...
- python基础===一道小学奥数题的解法
今早在博客园和大家分享了一道昨晚微博中看到的小学奥数题,后来有朋友给出了答案.然后我尝试用python解答它. 原题是这样的: 数学题:好事好 + 要做好 = 要做好事,求 “好.事.做.要”的值分别 ...
- 异步网络模块之aiohttp的使用(一)
异步网络模块之aiohttp的使用(一) 平时我们也许用的更多的是requests模块,或者是requests_hml模块,但是他们都属于阻塞类型的不支持异步,速度很难提高,于是后来出现了异步的gre ...
- 2014ACM/ICPC亚洲区广州站题解
这一场各种计算几何,统统没有做. HDU 5129 Yong Zheng's Death HDU 5136 Yue Fei's Battle
- 关于ueditor在Java中文件上传问题,404问题
问题困扰了两天,部署要求导入到webcontent下,我导入到了整个项目目录下,自己粗心犯错,导致页面访问不到404. 解决了上面的问题,试着进行文件上传,却一直找不到图片: 调出浏览器控制台: 刚开 ...