[bzoj1984]月下“毛景树”
Description
毛毛虫经过及时的变形,最终逃过的一劫,离开了菜妈的菜园。毛毛虫经过千山万水,历尽千辛万苦,最后来到了小小的绍兴一中的校园里。爬啊爬~爬啊爬~~毛毛虫爬到了一颗小小的“毛景树”下面,发现树上长着他最爱吃的毛毛果~~~
“毛景树”上有$N$个节点和$N-1$条树枝,但节点上是没有毛毛果的,毛毛果都是长在树枝上的。但是这棵“毛景树”有着神奇的魔力,他能改变树枝上毛毛果的个数:
$Change\;k\;w$:将第$k$条树枝上毛毛果的个数改变为$w$个。
$Cover\;u\;v\;w$:将节点$u$与节点$v$之间的树枝上毛毛果的个数都改变为$w$个。
$Add\;u\;v\;w$:将节点$u$与节点$v$之间的树枝上毛毛果的个数都增加$w$个。
由于毛毛虫很贪,于是他会有如下询问:
$Max\;u\;v$:询问节点$u$与节点$v$之间树枝上毛毛果个数最多有多少个。
Input
第一行一个正整数$N$。
接下来$N-1$行,每行三个正整数$U_i,V_i$和$W_i$,第$i+1$行描述第$i$条树枝。表示第$i$条树枝连接节点$U_i$和节点$V_i$,树枝上有$W_i$个毛毛果。
接下来是操作和询问,以“$Stop$”结束。
Output
对于毛毛虫的每个询问操作,输出一个答案。
Sample Input
4
1 2 8
1 3 7
3 4 9
Max 2 4
Cover 2 4 5
Add 1 4 10
Change 1 16
Max 2 4
Stop
Sample Output
9
16
HINT
$1\;\leq\;N\;\leq\;10^5$,操作+询问数目$\;\leq\;10^5$。
保证在任意时刻,所有树枝上毛毛果的个数都不会超过$10^9$个。
Solution
把边权转化成边的两端点中深度较大的那个点的点权.
树剖+线段树(注意下传$lazy$的细节,以及没被覆盖时值不该是$0$,因为存在全覆盖成$0$的情况)
#include<cmath>
#include<ctime>
#include<queue>
#include<stack>
#include<cstdio>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define N 100005
#define M 300005
using namespace std;
struct linetree{
int l,r,m,c,a;
}lt[M];
struct graph{
int nxt,to,w,n;
}e[M];
char c[10];
int g[N],a[N],b[N],w[N],n,u,v,x,cnt;
int f[N],p[N],dep[N],siz[N],son[N],top[N];
inline void addedge(int x,int y,int w,int n){
e[++cnt].nxt=g[x];g[x]=cnt;
e[cnt].to=y;e[cnt].w=w;e[cnt].n=n;
}
inline void dfs1(int u){
int m=0;siz[u]=1;
for(int i=g[u];i;i=e[i].nxt)
if(!dep[e[i].to]){
f[e[i].to]=u;
a[e[i].to]=e[i].w;
b[e[i].n]=e[i].to;
dep[e[i].to]=dep[u]+1;
dfs1(e[i].to);
siz[u]+=siz[e[i].to];
if(siz[e[i].to]>m){
m=siz[e[i].to];
son[u]=e[i].to;
}
}
}
inline void dfs2(int u,int tp){
top[u]=tp;p[u]=++cnt;w[cnt]=a[u];
if(son[u]) dfs2(son[u],tp);
for(int i=g[u];i;i=e[i].nxt)
if(f[u]!=e[i].to&&son[u]!=e[i].to)
dfs2(e[i].to,e[i].to);
}
inline void build(int u,int l,int r){
lt[u].l=l;lt[u].r=r;lt[u].c=-1;
if(lt[u].l<lt[u].r){
int lef=u<<1,rig=u<<1|1;
int mid=(lt[u].l+lt[u].r)>>1;
build(lef,l,mid);build(rig,mid+1,r);
lt[u].m=max(lt[lef].m,lt[rig].m);
}
else lt[u].m=w[lt[u].l];
}
inline void push(int u){
int lef=u<<1,rig=u<<1|1;
if(lt[u].c>=0){
lt[lef].a=lt[rig].a=0;
lt[lef].c=lt[rig].c=lt[u].c;
lt[lef].m=lt[rig].m=lt[u].c;
lt[u].c=-1;
}
if(lt[u].a){
lt[lef].a+=lt[u].a;
lt[rig].a+=lt[u].a;
lt[lef].m+=lt[u].a;
lt[rig].m+=lt[u].a;
lt[u].a=0;
}
}
inline void cover(int u,int l,int r,int k){
if(lt[u].l>=l&<[u].r<=r){
lt[u].m=lt[u].c=k;lt[u].a=0;
}
else if(lt[u].l<lt[u].r){
int lef=u<<1,rig=u<<1|1;
int mid=(lt[u].l+lt[u].r)>>1;
push(u);
if(l<=mid) cover(lef,l,r,k);
if(r>mid) cover(rig,l,r,k);
lt[u].m=max(lt[lef].m,lt[rig].m);
}
}
inline void add(int u,int l,int r,int k){
if(lt[u].l>=l&<[u].r<=r){
lt[u].m+=k;lt[u].a+=k;
}
else if(lt[u].l<lt[u].r){
int lef=u<<1,rig=u<<1|1;
int mid=(lt[u].l+lt[u].r)>>1;
push(u);
if(l<=mid) add(lef,l,r,k);
if(r>mid) add(rig,l,r,k);
lt[u].m=max(lt[lef].m,lt[rig].m);
}
}
inline int ask(int u,int l,int r){
if(lt[u].l>=l&<[u].r<=r)
return lt[u].m;
if(lt[u].l<lt[u].r){
int lef=u<<1,rig=u<<1|1,ret=0;
int mid=(lt[u].l+lt[u].r)>>1;
push(u);
if(l<=mid) ret=max(ret,ask(lef,l,r));
if(r>mid) ret=max(ret,ask(rig,l,r));
return ret;
}
}
inline void cov(int x,int y,int k){
int lt;
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]]){
lt=x;x=y;y=lt;
}
cover(1,p[top[x]],p[x],k);
x=f[top[x]];
}
if(p[x]>p[y]){
lt=x;x=y;y=lt;
}
if(p[x]<p[y]) cover(1,p[x]+1,p[y],k);
}
inline void ad(int x,int y,int k){
int lt;
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]]){
lt=x;x=y;y=lt;
}
add(1,p[top[x]],p[x],k);
x=f[top[x]];
}
if(p[x]>p[y]){
lt=x;x=y;y=lt;
}
if(p[x]<p[y]) add(1,p[x]+1,p[y],k);
}
inline int q(int x,int y){
int ret=0,lt;
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]]){
lt=x;x=y;y=lt;
}
ret=max(ret,ask(1,p[top[x]],p[x]));
x=f[top[x]];
}
if(p[x]>p[y]){
lt=x;x=y;y=lt;
}
if(p[x]<p[y]) ret=max(ret,ask(1,p[x]+1,p[y]));
return ret;
}
inline void Aireen(){
scanf("%d",&n);
for(int i=1,j,k,l;i<n;++i){
scanf("%d%d%d",&j,&k,&l);
addedge(j,k,l,i);addedge(k,j,l,i);
}
dep[1]=1;dfs1(1);
cnt=0;dfs2(1,1);
build(1,1,n);
while(true){
scanf("%s",c);
if(c[1]=='t') return;
if(c[1]=='h'){
scanf("%d%d",&u,&x);
cover(1,p[b[u]],p[b[u]],x);
}
else if(c[1]=='o'){
scanf("%d%d%d",&u,&v,&x);
cov(u,v,x);
}
else if(c[1]=='d'){
scanf("%d%d%d",&u,&v,&x);
ad(u,v,x);
}
else if(c[1]=='a'){
scanf("%d%d",&u,&v);
printf("%d\n",q(u,v));
}
}
}
int main(){
freopen("moon.in","r",stdin);
freopen("moon.out","w",stdout);
Aireen();
fclose(stdin);
fclose(stdout);
return 0;
}
[bzoj1984]月下“毛景树”的更多相关文章
- BZOJ1984: 月下“毛景树”
1984: 月下“毛景树” Time Limit: 20 Sec Memory Limit: 64 MBSubmit: 713 Solved: 245[Submit][Status] Descri ...
- [BZOJ1984]月下“毛景树”解题报告|树链剖分
Description 毛毛虫经过及时的变形,最终逃过的一劫,离开了菜妈的菜园. 毛毛虫经过千山万水,历尽千辛万苦,最后来到了小小的绍兴一中的校园里.爬啊爬~爬啊爬~~毛毛虫爬到了一颗小小的“毛景树” ...
- 【树链剖分】【分块】【最近公共祖先】【块状树】bzoj1984 月下“毛景树”
裸题,但是因为权在边上,所以要先把边权放到这条边的子节点上,然后进行链更新/查询的时候不能更新/查询其lca. #include<cstdio> #include<cmath> ...
- 2018.10.27 bzoj1984: 月下“毛景树”(树链剖分)
传送门 唉蒟蒻又退化了,这道sb题居然做了20min,最后发现是updcovupdcovupdcov写成了updaddupdaddupdadd我还能说什么233233233 就是让你转边权为点权之后, ...
- 【BZOJ1984】月下“毛景树” 树链剖分+线段树
[BZOJ1984]月下"毛景树" Description 毛毛虫经过及时的变形,最终逃过的一劫,离开了菜妈的菜园. 毛毛虫经过千山万水,历尽千辛万苦,最后来到了小小的绍兴一中的校 ...
- 【BZOJ-1984】月下“毛景树” 树链剖分
1984: 月下“毛景树” Time Limit: 20 Sec Memory Limit: 64 MBSubmit: 1314 Solved: 416[Submit][Status][Discu ...
- 树剖+线段树||树链剖分||BZOJ1984||Luogu4315||月下“毛景树”
题面:月下“毛景树” 题解:是道很裸的树剖,但处理的细节有点多(其实是自己线段树没学好).用一个Dfs把边权下移到点权,用E数组记录哪些边被用到了:前三个更新的操作都可以合并起来,可以发现a到b节点间 ...
- BZOJ 1984: 月下“毛景树” [树链剖分 边权]
1984: 月下“毛景树” Time Limit: 20 Sec Memory Limit: 64 MBSubmit: 1728 Solved: 531[Submit][Status][Discu ...
- Bzoj 1984: 月下“毛景树” 树链剖分
1984: 月下“毛景树” Time Limit: 20 Sec Memory Limit: 64 MBSubmit: 1282 Solved: 410[Submit][Status][Discu ...
随机推荐
- js删除数据的几种方法
js 删除数组几种方法 var arr=['a','b','c']; 若要删除其中的'b',有两种方法: 1.delete方法:delete arr[1] 这种方式数组长度不变,此时arr[1]变为u ...
- Symbol ES6 新增的一种值类型数据,表示一种绝不重复的值
let s1 = Symbol(33); let s2 = Symbol(33); alert(typeof(s1)); //数据类型判断 // alert(s1.toString()); //可把一 ...
- Smoothing in fMRI analysis (FAQ)
Source: http://mindhive.mit.edu/node/112 1. What is smoothing? "Smoothing" is generally us ...
- Basic: Fisher's transform
来源:http://bbs.chinahrd.net/thread-709742-1-1.html,Kenneth的回答. z = 0.5 * ln [ (1+r)/(1-r) ]" C0 ...
- Spring TestContext测试框架搭建
同样是测试,JUnit和Spring TestContext相比,Spring TestContext优势如下: 1.Spring TestContext可以手动设置测试事务回滚,不破坏数据现场 2. ...
- Java7并发编程实战(一) 线程的管理
1:线程的创建 1:继承Thread类,并且覆盖run()方法 2:创建一个实现Runnable接口的类.使用带参数的Thread构造器来构造 2:example-->计算打印乘法表 首先 ...
- 实验三 敏捷开发与XP实践
实验内容 1. XP基础 2. XP核心实践 3. 相关工具 实验要求 1.没有Linux基础的同学建议先学习<Linux基础入门(新版)><Vim编辑器> 课程 2.完成实验 ...
- windows7下启动mysql服务出现服务名无效
出现提示: WIN 7 cmd命令行下,net start mysql,出现 服务名无效提示: 问题原因: mysql服务没有安装. 解决办法: 在 mysql bin目录下 以管理员的权限 执行 m ...
- mvc5+ef6+Bootstrap 项目心得--WebGrid
1.mvc5+ef6+Bootstrap 项目心得--创立之初 2.mvc5+ef6+Bootstrap 项目心得--身份验证和权限管理 3.mvc5+ef6+Bootstrap 项目心得--WebG ...
- js中字符串和数组相互转化的方法
p.p1 { margin: 0.0px 0.0px 0.0px 0.0px; font: 12.0px "Helvetica Neue"; color: #e4af0a } p. ...