树链剖分 更新点 树的统计Count

一棵树上有n个节点,编号分别为1到n,每个节点都有一个权值w。我们将以下面的形式来要求你对这棵树完成
一些操作: I. CHANGE u t : 把结点u的权值改为t II. QMAX u v: 询问从点u到点v的路径上的节点的最大权值 I
II. QSUM u v: 询问从点u到点v的路径上的节点的权值和 注意:从点u到点v的路径上的节点包括u和v本身

Input
  输入的第一行为一个整数n,表示节点的个数。接下来n – 1行,每行2个整数a和b,表示节点a和节点b之间有
一条边相连。接下来n行,每行一个整数,第i行的整数wi表示节点i的权值。接下来1行,为一个整数q,表示操作
的总数。接下来q行,每行一个操作,以“CHANGE u t”或者“QMAX u v”或者“QSUM u v”的形式给出。
对于100%的数据,保证1<=n<=30000,0<=q<=200000;中途操作中保证每个节点的权值w在-30000到30000之间。

Output
  对于每个“QMAX”或者“QSUM”的操作,每行输出一个整数表示要求输出的结果。

Sample Input
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
Sample Output
4
1
2
2
10
6
5
6
5
16

#include<bits/stdc++.h>
using namespace std;
const int maxn=35010;
struct edge
{
    int to,nxt;
}edge[maxn*2];
int head[maxn],tot;
int ee[maxn];
int top[maxn],fa[maxn],deep[maxn],num[maxn],p[maxn],fp[maxn],son[maxn],pos;
void init()
{
    tot=0;
    memset(head,-1,sizeof(head));
    pos=0;
    memset(son,-1,sizeof(son));
}
void add_edge(int u,int v)
{
    edge[tot].to=v;
    edge[tot].nxt=head[u];
    head[u]=tot++;
}
void dfs1(int u,int pre,int d) //???dfs??fa,deep,num,son
{
    deep[u] = d;
    fa[u] = pre;
    num[u] = 1;
    for(int i = head[u];i != -1;i = edge[i].nxt)
    {
        int v = edge[i].to;
        if(v != pre)
        {
            dfs1(v,u,d+1);
            num[u] += num[v];
            if(son[u] == -1 || num[v] > num[son[u]])
                son[u] = v;
        }
    }
}


void getpos(int u,int sp)
{
    top[u]=sp;
    p[u]=pos++;
    fp[p[u]]=u;
    if(son[u] == -1) return;
    getpos(son[u],sp);
    for(int i=head[u];i!=-1;i=edge[i].nxt)
    {
        int v=edge[i].to;
        if(v!=son[u]&&v!=fa[u])
            getpos(v,v);
    }
}
struct Node
{
    int l,r;
    int maxx;
    int summ;
}segtree[maxn*3];
void push_up(int i)
{
    segtree[i].summ = segtree[i<<1].summ + segtree[(i<<1)|1].summ;
    segtree[i].maxx = max(segtree[i<<1].maxx,segtree[(i<<1)|1].maxx);
}
void build(int i,int l,int r)
{
    segtree[i].l = l;
    segtree[i].r = r;
    if(l == r)
    {
        segtree[i].summ = segtree[i].maxx = ee[fp[l]];
        return ;
    }
    int mid = (l + r)/2;
    build(i<<1,l,mid);
    build((i<<1)|1,mid+1,r);
    push_up(i);
}

void update(int i,int k,int val)//???????k???val
{
    if(segtree[i].l == k && segtree[i].r == k)
    {
        segtree[i].summ = segtree[i].maxx = val;
        return;
    }
    int mid = (segtree[i].l + segtree[i].r)/2;
    if(k <= mid)update(i<<1,k,val);
    else update((i<<1)|1,k,val);
    push_up(i);
}
int query(int i,int l,int r)
{
    if(segtree[i].l==l&&segtree[i].r==r)
        return segtree[i].maxx;
    int mid=(segtree[i].l+segtree[i].r)/2;
    if(r<=mid)
        return query(i<<1,l,r);
    else if(l>mid)
        return query(i<<1|1,l,r);
    else
        return max(query(i<<1,l,mid),query(i<<1|1,mid+1,r));
}
int query2(int i,int l,int r)
{
    if(segtree[i].l==l&&segtree[i].r==r)
        return segtree[i].summ;
    int mid=(segtree[i].l+segtree[i].r)/2;
    if(r<=mid)
        return query2(i<<1,l,r);
    else if(l>mid)
        return query2((i<<1)|1,l,r);
    else
        return query2(i<<1,l,mid)+query2((i<<1)|1,mid+1,r);
}
int findd(int u,int v)
{
    int f1=top[u],f2=top[v];
    int tmp=-1000000000;
    while(f1!=f2)
    {
        if(deep[f1]<deep[f2])
        {
            swap(f1,f2);
            swap(u,v);
        }
        tmp=max(tmp,query(1,p[f1],p[u]));
        u=fa[f1];
        f1=top[u];
    }
    if(deep[u]>deep[v])
        swap(u,v);
    return max(tmp,query(1,p[u],p[v]));
}
int findd2(int u,int v) //??u->v??????????
{
    int f1 = top[u], f2 = top[v];
    int tmp = 0;
    while(f1 != f2)
    {
        if(deep[f1] < deep[f2])
        {
            swap(f1,f2);
            swap(u,v);
        }
        tmp += query2(1,p[f1],p[u]);
        u = fa[f1];
        f1 = top[u];
    }
    if(deep[u] > deep[v]) swap(u,v);
    return tmp + query2(1,p[u],p[v]);
}
int main ()
{
    int T;
    int n;
    while(~scanf("%d",&n))
    {
        init();
        for(int i=0;i<n-1;i++)
        {
            int xx,yy;
            scanf("%d %d",&xx,&yy);
            add_edge(xx,yy);
            add_edge(yy,xx);
        }
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&ee[i]);
        }
        dfs1(1,0,0);
        getpos(1,1);
        build(1,0,pos-1);
        char op[20];
        int u,v;
        int opp;
        cin>>opp;
        while(opp--)
        {
            scanf("%s",op);
            scanf("%d%d",&u,&v);
            if(op[1]=='M')
                printf("%d\n",findd(u,v));
            else if(op[1]=='S')
                printf("%d\n",findd2(u,v));
            else
                update(1,p[u],v);
        }
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值