题目描述
一棵树上有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”的操作,每行输出一个整数表示要求输出的结果。
输入输出样例
输入 #1 复制
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
输出 #1 复制
4
1
2
2
10
6
5
6
5
16
说明/提示
对于100%的数据,保证1<=n<=30000,0<=q<=200000;中途操作中保证每个节点的权值w在-30000到30000之间。
解题思路:
直接使用树链剖分+线段树就可以,这是板题
#include <iostream>
#include<cstdio>
#include<fstream>
#include<cstring>
#include<queue>
#include<cmath>
#include<vector>
#include<algorithm>
using namespace std;
const int maxn = 30005;
const int inf = 0x3f3f3f3f;
struct Edge{
int to,nxt;
}edges[maxn<<1];
int head[maxn<<1];
int tot;
void addEdge(int u, int v)
{
edges[tot] = Edge{v,head[u]};
head[u] = tot++;
}
void init()
{
memset(head,-1,sizeof(head));
tot = 1;
}
int n;
int val[maxn];
/*-----------开始树链剖分--------------*/
int fa[maxn],depth[maxn],sz[maxn],son[maxn];
void dfs1(int u, int pre) //当前结点,当前结点的父亲
{
fa[u] = pre;
depth[u] = depth[pre]+1;
sz[u] = 1;
for(int i = head[u]; ~i; i = edges[i].nxt)
{
int v = edges[i].to;
if(v == pre) continue;
dfs1(v,u);
sz[u]+=sz[v];
if(sz[v]>sz[son[u]])
son[u] = v; //找重儿子
}
}
int top[maxn], id[maxn], rk[maxn];
int num; //这个需要初始化
void dfs2(int u, int t)
{
top[u] = t;
id[u] = ++num;
rk[num] = u;
if(son[u]) dfs2(son[u],t);
for(int i = head[u]; ~i; i = edges[i].nxt)
{
int v = edges[i].to;
if(v == fa[u]||v == son[u]) continue;
dfs2(v,v);
}
}
void chain_init()
{
memset(son,0,sizeof(son));
num = 0;
}
/*-----------树链剖分结束--------------*/
/*-----------线段树开始--------------*/
int Sum[maxn<<2], Max[maxn<<2];
void PushUp(int rt)
{
Sum[rt] = Sum[rt<<1]+Sum[rt<<1|1];
Max[rt] = max(Max[rt<<1],Max[rt<<1|1]);
}
void Build(int l, int r, int rt)
{
if(l==r)
{
Sum[rt] = val[rk[l]];
Max[rt] = val[rk[l]];
return ;
}
int mid = (l+r)>>1;
Build(l,mid,rt<<1);
Build(mid+1,r,rt<<1|1);
PushUp(rt);
}
void update(int L, int C, int l, int r, int rt)
{
if(l==r)
{
Sum[rt] = C;
Max[rt] = C;
return ;
}
int mid = (l+r)>>1;
if(L<=mid) update(L,C,l,mid,rt<<1);
else update(L,C,mid+1,r,rt<<1|1);
PushUp(rt);
}
int Query_Sum(int L, int R, int l, int r, int rt)
{
if(l>=L&&r<=R)
return Sum[rt];
int tmp = 0;
int mid = (l+r)>>1;
if(mid>=L)
tmp+=Query_Sum(L,R,l,mid,rt<<1);
if(mid<R)
tmp+=Query_Sum(L,R,mid+1,r,rt<<1|1);
return tmp;
}
int Query_Max(int L, int R, int l, int r, int rt)
{
if(l>=L&&r<=R)
return Max[rt];
int tmp = -inf;
int mid = (l+r)>>1;
if(mid>=L)
tmp = max(Query_Max(L,R,l,mid,rt<<1),tmp);
if(mid<R)
tmp = max(Query_Max(L,R,mid+1,r,rt<<1|1),tmp);
return tmp;
}
/*-----------线段树结束--------------*/
/*-----------链上查询开始--------------*/
int qsum(int x, int y)
{
int res = 0;
while(top[x]!=top[y])
{
if(depth[top[x]] < depth[top[y]]) swap(x,y);
res+=Query_Sum(id[top[x]],id[x],1,n,1);
x = fa[top[x]];
}
if(depth[x]>depth[y])swap(x,y);
res+=Query_Sum(id[x],id[y],1,n,1);
return res;
}
int qmax(int x, int y)
{
int res = -inf;
while(top[x]!=top[y])
{
if(depth[top[x]] < depth[top[y]]) swap(x,y);
res = max(Query_Max(id[top[x]],id[x],1,n,1),res);
x = fa[top[x]];
}
if(depth[x]>depth[y])swap(x,y);
res = max(Query_Max(id[x],id[y],1,n,1),res);
return res;
}
/*-----------链上查询结束--------------*/
int main()
{
scanf("%d", &n);
init();
int u,v;
for(int i = 1; i < n; i++)
{
scanf("%d%d", &u, &v);
addEdge(u,v);
addEdge(v,u);
}
chain_init();
dfs1(1,0);
dfs2(1,1);
for(int i = 1; i <= n; i++)
scanf("%d", &val[i]);
Build(1,n,1);
int q;
scanf("%d", &q);
char op[10];
while(q--)
{
scanf("%s%d%d",op,&u, &v);
if(op[1] == 'M')
printf("%d\n",qmax(u,v));
else if(op[1] == 'S')
printf("%d\n",qsum(u,v));
else
update(id[u],v,1,n,1);
}
return 0;
}