树的统计 zjoi2008 —树链剖分
描述:
一棵树上有n个节点,编号分别为1到n,每个节点都有一个权值w。
我们将以下面的形式来要求你对这棵树完成一些操作:
I. CHANGE u t : 把结点u的权值改为t
II. QMAX u v: 询问从点u到点v的路径上的节点的最大权值
III. QSUM u v: 询问从点u到点v的路径上的节点的权值和
注意:从点u到点v的路径上的节点包括u和v本身
输入:
输入文件的第一行为一个整数n,表示节点的个数。
接下来n – 1行,每行2个整数a和b,表示节点a和节点b之间有一条边相连。
接下来一行n个整数,第i个整数wi表示节点i的权值。
接下来1行,为一个整数q,表示操作的总数。
接下来q行,每行一个操作,以“CHANGE u t”或者“QMAX u v”或者“QSUM u v”的形式给出。
输出:
对于每个“QMAX”或者“QSUM”的操作,每行输出一个整数表示要求输出的结果。
样例输入:
4
1 2
2 3
4 1
4 2 1 3
12
QMAX 3 4
QMAX 3 3
QMAX 3 2
QMAX 2 3
QSUM 3 4
QSUM 2 1
CHANGE 1 5
QMAX 3 4
CHANGE 3 6
QMAX 3 4
QMAX 2 4
QSUM 3 4
样例输出:
4
1
2
2
10
6
5
6
5
16
数据范围:
对于100%的数据,保证1<=n<=30000,0<=q<=200000;中途操作中保证每个节点的权值w在-30000到30000之间。
解题思路:
我也是够了,一道树链剖分的模板题用了我整整一个晚上,智商不够啊!而且还是被一个SB错误卡了下去。。。言归正传,正宗的树链剖分模板题,线段树上维护即可,判断轻重子节点拉成一条链,用num数组很巧妙地在树链和线段树之间转换,非常玄学的操作,废话不多说,上代码。
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
#define INF 0x3f3f3f3f
#define LD (o<<1)
#define RD ((o<<1)|1)
#define mem(a) memset(a,0,sizeof(a))
#define N 30010
struct Tree{int l,r,maxv,sum;}T[4*N];//线段树
struct Edge{int v,next;}E[N];//边
struct Node{
int dep,siz,fa,w,hson,num,top;
Node(){dep=siz=fa=w=hson=num=top=0;}
}P[N];//点
/*
dep深度
size节点数
fa父亲
w节点权值
hson重孩子
num在线段树编号
top表示每条链第一个位置
*/
int first[N],pred[N];
/*
first 邻接表头指针
pred线段树的编号是原来什么编号
*/
int n,m,last,tot;//tot记录线段树总编号
char c[10];
void init(){
last=tot=0;
mem(first);mem(pred);
}
void add(int u,int v){
E[++last].v=v;
E[last].next=first[u];
first[u]=last;
}
void dfs1(int s){//深搜出每个节点的重子节点,顺便建树
P[s].siz=1;//当前为根节点的子树大小
P[s].hson=0;//默认重孩子为0
for(int i=first[s];i;i=E[i].next){//邻接表
int &son=E[i].v;
if(son==P[s].fa)continue;
P[son].dep=P[s].dep+1;
P[son].fa=s;
dfs1(son);
P[s].siz+=P[son].siz;
if(P[P[s].hson].siz<P[son].siz)//打擂台判断重子节点
P[s].hson=son;
}
}
void dfs2(int s,int tp){//深搜出重链
P[s].top=tp;//每一条链的起点
P[s].num=++tot;//重链节点编号
pred[tot]=s;
if(P[s].hson)dfs2(P[s].hson,tp);//有重孩子继续深搜
for(int i=first[s];i;i=E[i].next){
int &son=E[i].v;
if(son!=P[s].fa&&son!=P[s].hson)
dfs2(son,son);//重新搜索一链
}
}
void pushup(int o){
T[o].sum=T[LD].sum+T[RD].sum;
T[o].maxv=max(T[LD].maxv,T[RD].maxv);
}
void build(int o,int l,int r){
T[o].l=l;T[o].r=r;
if(l==r){T[o].sum=T[o].maxv=P[pred[l]].w;return;}
int mid=(l+r)>>1;
build(LD,l,mid);
build(RD,mid+1,r);
pushup(o);
}
void update(int o,int p,int val){
if(T[o].l==p&&T[o].r==p){T[o].sum=T[o].maxv=val;return;}
int mid=(T[o].l+T[o].r)>>1;
if(p<=mid)update(LD,p,val);
if(p>=mid+1)update(RD,p,val);
pushup(o);
}
int query_sum(int o,int l,int r){
if(l<=T[o].l&&T[o].r<=r)return T[o].sum;
int mid=(T[o].l+T[o].r)>>1;
if(r<=mid)return query_sum(LD,l,r);
if(l>mid)return query_sum(RD,l,r);
if(l<=mid&&r>mid)return query_sum(LD,l,mid)+query_sum(RD,mid+1,r);
}
int query_maxv(int o,int l,int r){
if(l<=T[o].l&&T[o].r<=r)return T[o].maxv;
int mid=(T[o].l+T[o].r)>>1;
if(r<=mid)return query_maxv(LD,l,r);
if(l>mid)return query_maxv(RD,l,r);
if(l<=mid&&r>mid)return max(query_maxv(LD,l,mid),query_maxv(RD,mid+1,r));
}
//求链 目的是想办法跳到一条重链上来
//由于一颗树只有一个完整重莲,如果2个点在一颗子树上,那么先统计轻边,自然会转到重边上来。
int get_sum(int x,int y){
int ans=0;
while(P[x].top!=P[y].top){//不是一条重链
if(P[P[x].top].dep<P[P[y].top].dep)swap(x,y);
//链的首点深度大的在前
ans+=query_sum(1,P[P[x].top].num,P[x].num);
//统计x这条链
x=P[P[x].top].fa;//x变为这条链顶点的父亲
}
if(P[x].num>P[y].num)swap(x,y);
ans+=query_sum(1,P[x].num,P[y].num);
return ans;
}
int get_maxv(int x,int y){
int ans=-INF;
while(P[x].top!=P[y].top){//不是一条重链
if(P[P[x].top].dep<P[P[y].top].dep)swap(x,y);
//链的首点深度大的在前
ans=max(ans,query_maxv(1,P[P[x].top].num,P[x].num));
//统计x这条链
x=P[P[x].top].fa;//x变为这条链顶点的父亲
}
if(P[x].num>P[y].num)swap(x,y);
ans=max(ans,query_maxv(1,P[x].num,P[y].num));
return ans;
}
int main(){
init();
scanf("%d",&n);
for(int i=1;i<n;i++){
int a,b;scanf("%d%d",&a,&b);
add(a,b);
add(b,a);//双向建边
}
for(int i=1;i<=n;i++)scanf("%d",&P[i].w);
dfs1(1);
dfs2(1,1);
build(1,1,tot);//构建线段树
scanf("%d",&m);
for(int i=1;i<=m;i++){
int x,y;scanf("%s%d%d",c,&x,&y);
if(c[0]=='C'){P[x].w=y;update(1,P[x].num,y);}
else{
if(c[1]=='M')printf("%d\n",get_maxv(x,y));
if(c[1]=='S')printf("%d\n",get_sum(x,y));
}
}
return 0;
}