【边权树链剖分】BZOJ1984 月下“毛景树”

题面在这里

题目很水,但是细节很多

首先是用树链剖分维护边权的姿势:每个点代表连向父亲的边

这样就要注意两点的LCA不能算入答案中

然后是线段树两个标记,只能同时存在一个

具体写法自己看代码

示例程序:

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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
#include<cstdio>
#include<algorithm>
using namespace std;

const int maxn=100005,maxe=200005,maxs=4*maxn;
int n,N,w[maxn],id[maxn];
int son[maxe],nxt[maxe],lnk[maxn],tot;
inline void add(int x,int y){
son[++tot]=y;nxt[tot]=lnk[x];lnk[x]=tot;
}
int mx[maxs],ad[maxs],cv[maxs];
inline void pushup(int x) {mx[x]=max(mx[x<<1],mx[x<<1|1]);}
inline void addad(int x,int w) {mx[x]+=w,ad[x]+=w;if (cv[x]) cv[x]+=w;}
inline void addcv(int x,int w) {mx[x]=w;cv[x]=w;ad[x]=0;}
inline void pushdown(int x){
if (ad[x]) {addad(x<<1,ad[x]);addad(x<<1|1,ad[x]);ad[x]=0;}
if (cv[x]) {addcv(x<<1,cv[x]);addcv(x<<1|1,cv[x]);cv[x]=0;}
}
void ist(int x,int l,int r,int ql,int qr,int w){ //+w:add -w:cover
if (ql<=l&&r<=qr) {if (w>0) addad(x,w);else addcv(x,-w);return;}
if (qr<l||r<ql) return;
int mid=l+r>>1; pushdown(x);
ist(x<<1,l,mid,ql,qr,w); ist(x<<1|1,mid+1,r,ql,qr,w);
pushup(x);
}
int qry(int x,int l,int r,int ql,int qr){
if (ql<=l&&r<=qr) return mx[x];
if (qr<l||r<ql) return 0;
int mid=l+r>>1; pushdown(x);
return max(qry(x<<1,l,mid,ql,qr),qry(x<<1|1,mid+1,r,ql,qr));
}
int H_son[maxn],siz[maxn],fa[maxn],dep[maxn],in[maxn],out[maxn],top[maxn],times;
void getH(int x){
siz[x]=1;dep[x]=dep[fa[x]]+1;H_son[x]=0;
for (int j=lnk[x];j;j=nxt[j])
if (son[j]!=fa[x]){
id[j+1>>1]=son[j];fa[son[j]]=x; getH(son[j]); siz[x]+=siz[son[j]];
if (!H_son[x]||siz[H_son[x]]<siz[son[j]]) H_son[x]=son[j];
}
}
void dfs(int x,int lst){
top[x]=lst;in[x]=++times;
if (H_son[x]) dfs(H_son[x],lst);
for (int j=lnk[x];j;j=nxt[j])
if (son[j]!=fa[x]&&son[j]!=H_son[x]) dfs(son[j],son[j]);
out[x]=times;
}
void cao(int x,int y,int w){
while (top[x]!=top[y]){
if (dep[top[x]]<dep[top[y]]) swap(x,y);
ist(1,1,N,in[top[x]],in[x],w); x=fa[top[x]];
}
if (in[x]>in[y]) swap(x,y);
ist(1,1,N,in[x]+1,in[y],w);
}
int ask(int x,int y){
int res=0;
while (top[x]!=top[y]){
if (dep[top[x]]<dep[top[y]]) swap(x,y);
res=max(res,qry(1,1,N,in[top[x]],in[x])); x=fa[top[x]];
}
if (in[x]>in[y]) swap(x,y);
return max(res,qry(1,1,N,in[x]+1,in[y]));
}
char s[20];
int main(){
scanf("%d",&n); N=n;
for (int i=1,x,y;i<n;i++) scanf("%d%d%d",&x,&y,&w[i]),add(y,x),add(x,y);
fa[1]=1;getH(1); dfs(1,1);
for (int i=1;i<n;i++) ist(1,1,N,in[id[i]],in[id[i]],-w[i]);
for (scanf("%s",s);s[0]!='S';scanf("%s",s))
if (s[1]=='h'){
int x,w;scanf("%d%d",&x,&w);
ist(1,1,N,in[id[x]],in[id[x]],-w);
}else
if (s[1]=='o'){
int x,y,w;scanf("%d%d%d",&x,&y,&w);
cao(x,y,-w);
}else
if (s[1]=='d'){
int x,y,w;scanf("%d%d%d",&x,&y,&w);
cao(x,y,w);
}else{
int x,y;scanf("%d%d",&x,&y);
printf("%d\n",ask(x,y));
}
return 0;
}