震波

在一片土地上有 $N$ 个城市，通过 $N-1$ 条无向边互相连接，形成一棵树的结构，相邻两个城市的距离为 $1$，其中第 $i$ 个城市的价值为 $value[i]$。
不幸的是，这片土地常常发生地震，并且随着时代的发展，城市的价值也往往会发生变动。

接下来你需要在线处理 $M$ 次操作：
- $0~x~k$ 表示发生了一次地震，震中城市为 $x$ ，影响范围为 $k$ ，所有与 $x$ 距离不超过 $k$ 的城市都将受到影响，该次地震造成的经济损失为所有受影响城市的价值和。
- $1~x~y$ 表示第$x$个城市的价值变成了 $y$ 。

为了体现程序的在线性，操作中的 $x$ 、$y$ 、$k$ 都需要异或你程序上一次的输出来解密，如果之前没有输出，则默认上一次的输出为 $0$ 。

 

 

---

 

 

Solution
考虑建出点分树，每一个点开一棵以长度为下标的线段树。
查询时可以一步步往上跳。假设他到父亲的权值为t,每此查询与x距离为k k+t ...的点。
可是这样会算重。
那么可以在每个点再开一棵线段树，下标为点到他点分树上父亲的权值。
每次减去第二棵线段树中k+t的值。也就是自己与父亲相同的贡献

卡常毒题
orz邵神树状数组常树优秀

 

#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define maxn 200006
using namespace std;
int n,m,head[maxn],a[maxn],val[maxn],siz,tot,deep[maxn];
int f[maxn],vis[maxn],par[maxn],sz[maxn],rt;
int cnt,root[maxn],st[maxn][20],sc,dfn[maxn],l2[maxn];
struct node{
    int v,nex;
}e[maxn];
struct no{
    int ls,rs,v;
}tr[maxn*140];
inline void lca(int k,int fa){
    st[dfn[k]=++sc][0]=deep[k]=deep[fa]+1;
    for(int i=head[k];i;i=e[i].nex)if(e[i].v!=fa)lca(e[i].v,k),st[++sc][0]=deep[k];
}
inline void add(int &k,int l,int r,int pl,int v){
    if(!k)k=++cnt;
    if(l==r){tr[k].v+=v;return;}
    int mid=l+r>>1;
    if(pl<=mid)add(tr[k].ls,l,mid,pl,v);
    else add(tr[k].rs,mid+1,r,pl,v);
    tr[k].v=tr[tr[k].ls].v+tr[tr[k].rs].v;
}
inline int ask(int k,int l,int r,int pl){
    if(!k||pl<0)return 0;
    if(l==r)return tr[k].v;
    int mid=l+r>>1;
    if(pl<=mid)return ask(tr[k].ls,l,mid,pl);
    else return tr[tr[k].ls].v+ask(tr[k].rs,mid+1,r,pl);
    
}
inline int dis(int x,int y){
    if(dfn[y]<dfn[x])swap(x,y);
    int i=l2[dfn[y]-dfn[x]+1];
    int dl=min(st[dfn[x]][i],st[dfn[y]-(1<<i)+1][i]);
    return deep[x]+deep[y]-(dl<<1);
}
inline void dfs(int o,int la,int k,int fa){
    add(root[o],0,n,dis(o,k),val[k]);
    if(la)add(root[o+n],0,n,dis(k,la),val[k]);
    sz[k]=1;
    for(int i=head[k];i;i=e[i].nex){
        if(e[i].v==fa||vis[e[i].v])continue;
        dfs(o,la,e[i].v,k);
        sz[k]+=sz[e[i].v];
    }
}
inline void findr(int k,int fa){
    sz[k]=1;f[k]=0;
    for(int i=head[k];i;i=e[i].nex){
        if(e[i].v==fa||vis[e[i].v])continue;
        findr(e[i].v,k);
        sz[k]+=sz[e[i].v];
        f[k]=max(f[k],sz[e[i].v]);
    }
    f[k]=max(f[k],siz-sz[k]);
    if(f[k]<f[rt])rt=k;
}
inline void work(int k,int la){
    vis[k]=1;
    dfs(k,la,k,0);
    for(int i=head[k];i;i=e[i].nex){
        if(!vis[e[i].v]){
            siz=sz[e[i].v];rt=0;findr(e[i].v,k);
            par[rt]=k;
            work(rt,k);
        }
    }
}
inline int R(){
    int v=0;char ch;
    while(!isdigit(ch=getchar()));v=ch-48;
    while(isdigit(ch=getchar()))v=(v<<1)+(v<<3)+ch-48;
    return v;
}
int main()
{
    n=R();m=R();
    for(int i=1;i<=n;i++)val[i]=R();
    for(int i=1,t1,t2;i<n;i++){
        t1=R();t2=R();
        e[++tot].v=t2;e[tot].nex=head[t1];head[t1]=tot;
        e[++tot].v=t1;e[tot].nex=head[t2];head[t2]=tot;
    }
    lca(1,0);
    for(int i=2;i<=sc;i++)l2[i]=l2[i/2]+1;
    for(int j=1;j<=l2[sc];j++)
        for(int i=1;i+(1<<j)<=sc+1;i++)st[i][j]=min(st[i][j-1],st[i+(1<<j-1)][j-1]);
    rt=0;f[0]=1e9;siz=n;findr(1,0);
    work(rt,0);
    int ans=0;
    for(int i=1,op,y,x,k;i<=m;i++){
        op=R();
        if(op==0){
            x=R();k=R();
            x^=ans,k^=ans;
            ans=0;int fi=x;
            for(int d=0;x;){
                ans+=ask(root[x],0,n,k-d);d=dis(par[x],fi);
                if(par[x])ans-=ask(root[n+x],0,n,k-d);
                x=par[x];
            }
            printf("%d\n",ans);
        }
        if(op==1){
            x=R();y=R();
            x^=ans;y^=ans;
            int fi=x;
            for(int d=0;x;){
                add(root[x],0,n,d,y-val[fi]);d=dis(fi,par[x]);
                if(par[x])add(root[n+x],0,n,d,y-val[fi]);
                x=par[x];
            }
            val[fi]=y;
        }
    }
    return 0;
}

 

 

---