P2590 树的统计(树链剖分+线段树)

题目描述

一棵树上有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;
}