lightoj1348(树链剖分&线段树)

常规树链剖分题,把点权哈希到数组中,用线段树维护。

#include<cstdio>
#include<vector>
#include<cstring>
using namespace std;
const int N = 30000+10;
int dep[N],siz[N],fa[N],id[N],son[N],val[N],top[N],num;
vector<int>v[N];
void dfs1(int u, int f, int d) {
    dep[u] = d;
    siz[u] = 1;
    son[u] = 0;
    fa[u] = f;
    for (int i = 0; i < v[u].size(); i++) {
        int ff = v[u][i];
        if (ff == f) continue;
        dfs1(ff, u, d + 1);
        siz[u] += siz[ff];
        if (siz[son[u]] < siz[ff])
            son[u] = ff;
    }
}
void dfs2(int u, int tp) {
    top[u] = tp;
    id[u] = ++num;
    if (son[u]) dfs2(son[u], tp);
    for (int i = 0; i < v[u].size(); i++) {
        int ff = v[u][i];
        if (ff == fa[u] || ff == son[u]) continue;
        dfs2(ff, ff);
    }
}
int tree[N<<2];
void update(int pos,int w,int l,int r,int root)
{
    if(l==r) {tree[root] = w; return ; }
    int mid = (l+r)/2; int tmp = root*2;
    if(pos<=mid) update(pos,w,l,mid,tmp);
    else if(pos>mid) update(pos,w,mid+1,r,tmp+1);
    tree[root] = tree[tmp] + tree[tmp+1];
}
int query(int L,int R,int l,int r,int root)
{
    if(L<=l&&R>=r){return tree[root]; }
    int mid = (l+r)/2; int tmp = root*2;
    int ans = 0;
    if(mid>=L) ans += query(L,R,l,mid,tmp);
    if(mid<R) ans += query(L,R,mid+1,r,tmp+1);
    return ans ;
}
int _find(int u,int v){
    int f1=top[u],f2=top[v],tmp=0;
    while(f1!=f2){
        if(dep[f1]<dep[f2]){
            swap(f1,f2);
            swap(u,v);
        }
        tmp+=query(id[f1],id[u],1,num,1);
        u=fa[f1];f1=top[u];
    }
    if(dep[u]>dep[v]) swap(u,v);
    tmp+=query(id[u],id[v],1,num,1);
    return tmp;
}
int main()
{
    int cases,n,a,b,c,t=1;
    scanf("%d",&cases);
    while(cases--)
    {
        scanf("%d",&n);
        for(int i=0;i<=n;i++) v[i].clear();
        memset(tree,0,sizeof(tree));
        for(int i=1;i<=n;i++) scanf("%d",&val[i]);
        for(int i=1;i<n;i++)
        {
            scanf("%d%d",&a,&b); a++; b++;
            v[a].push_back(b);
            v[b].push_back(a);
        }
        num = 0; dfs1(1,0,1); dfs2(1,1);
        for(int i=1;i<=n;i++) update(id[i],val[i],1,num,1);
        int q ; scanf("%d",&q);
        printf("Case %d:\n",t++);
        while(q--)
        {
            scanf("%d%d%d",&a,&b,&c); b++;
            if(a==1) update(id[b],c,1,num,1);
            else c++,printf("%d\n",_find(b,c));
        }
    }
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值