3730: 震波 Time Limit:? 15 Sec??Memory Limit:? 256 MBSubmit:? 1202??Solved:? 288Description 在一片土地上有N個城市，通過N-1條無向邊互相連接，形成一棵樹的結構，相鄰兩個城市的距離為1，其中第i個城市的價值為value[i]。
Input 第一行包含兩個正整數N和M。
Output 包含若干行，對于每個詢問輸出一行一個正整數表示答案。
Sample Input 8 1 Sample Output 11100101 HINT
1<=N,M<=100000
我們要詢問和一個點距離<=k的所有點權和 并且支持修改
又要求強制在線
怎么搞呢？
動態點分治
搞出點分樹 利用點分樹log層的優勢來搞事
考慮我們要維護什么
對于每個分治根，與其子樹中的距離就是深度
所以考慮對于每個分治根，在子樹中，利用其深度進行信息的維護
怎么搞呢~開線段樹
每個分治根都開 以深度為下標 權值插入就行了
還有一個東西要注意下
每次由分治根向上爬 去更新答案
但是會出現算重的情況啊
所以每個點開兩個線段樹，一個表示以它為分治根，其子樹的信息
一個表示它的分治父親中，它的子樹中的信息
搞的時候剪掉就行了
這里有圖和樣例解釋
#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>
using namespace std;typedef long long ll;inline int read()
{int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}return x*f;
}
void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}const int N=500100,inf=0X3f3f3f3f;int last[N],ecnt;
struct EDGE{int to,nt;}e[N<<1];
inline void add(int u,int v)
{e[++ecnt]=(EDGE){v,last[u]};last[u]=ecnt;}int root,n,sum,ans,f[N],size[N];
int V[N];//ancestor
int anc[N][20],dep[N];void dfs(int u)
{for(int i=1;(1<<i)<=dep[u];++i)anc[u][i]=anc[anc[u][i-1]][i-1];for(int i=last[u];i;i=e[i].nt)if(e[i].to^anc[u][0]){dep[e[i].to]=dep[u]+1;anc[e[i].to][0]=u;dfs(e[i].to);}
}inline int getlca(int u,int v)
{if(dep[u]<dep[v])swap(u,v);int len=dep[u]-dep[v];for(int i=0;(1<<i)<=len;++i)if(len&(1<<i))u=anc[u][i];if(u==v)return u;for(int i=17;~i;i--)if(anc[u][i]^anc[v][i])u=anc[u][i],v=anc[v][i];return anc[u][0];
}inline int getdis(int u,int v)
{int lca=getlca(u,v);return dep[u]+dep[v]-(dep[lca]<<1);}bool vis[N];void getroot(int u,int fa)
{size[u]=1;f[u]=0;for(int i=last[u];i;i=e[i].nt)if(e[i].to^fa&&!vis[e[i].to]){getroot(e[i].to,u);size[u]+=size[e[i].to];f[u]=max(f[u],size[e[i].to]);}f[u]=max(f[u],sum-size[u]);if(f[u]<f[root])root=u;
}int rt[N][2];struct seg_tree{int ls,rs,sum;}tr[N<<4];
int sz=0;void insert(int &k,int l,int r,int x,int val)
{if(!k)k=++sz;tr[k].sum+=val;if(l==r)return ;int mid=(l+r)>>1;x<=mid?insert(tr[k].ls,l,mid,x,val):insert(tr[k].rs,mid+1,r,x,val);
}int query_seg(int k,int l,int r,int x,int y)
{if(!k)return 0;if(l>=x&&r<=y)return tr[k].sum;int mid=(l+r)>>1;if(y<=mid)return query_seg(tr[k].ls,l,mid,x,y);if(x>mid)return query_seg(tr[k].rs,mid+1,r,x,y);return query_seg(tr[k].ls,l,mid,x,y)+query_seg(tr[k].rs,mid+1,r,x,y);
}// source destination
void seek(int u,int ffa,int d,int dest,int p)
{insert(rt[dest][p],0,n,d,V[u]);for(int i=last[u];i;i=e[i].nt)if(ffa^e[i].to&&!vis[e[i].to])seek(e[i].to,u,d+1,dest,p);
}int fa[N];void init(int u,int ffa)
{vis[u]=1;fa[u]=ffa;seek(u,0,0,u,0);for(int i=last[u];i;i=e[i].nt)if(!vis[e[i].to]){sum=size[e[i].to];root=0;getroot(e[i].to,u);seek(e[i].to,u,1,root,1);init(root,u);}
}void modify(int u,int val)
{val=val-query_seg(rt[u][0],0,n,0,0);insert(rt[u][0],0,n,0,val);for(int dis,v=u;fa[v];v=fa[v]){dis=getdis(fa[v],u);insert(rt[v][1],0,n,dis,val);insert(rt[fa[v]][0],0,n,dis,val);}
}void query(int u,int k)
{int res=query_seg(rt[u][0],0,n,0,k);for(int dis,v=u;fa[v];v=fa[v]){dis=getdis(fa[v],u);res-=query_seg(rt[v][1],0,n,0,k-dis);res+=query_seg(rt[fa[v]][0],0,n,0,k-dis);}print(res);ans=res;puts("");
}int main()
{n=sum=read();int Q=read();register int i,opt,u,v;for(i=1;i<=n;++i)V[i]=read();for(i=1;i<n;++i){u=read();v=read();add(u,v);add(v,u);}dfs(1);f[0]=inf;getroot(1,0);init(root,0);while(Q--){opt=read();u=read()^ans;v=read()^ans;opt?modify(u,v):query(u,v);}return 0;
}
/*
8 1
1 10 100 1000 10000 100000 1000000 10000000
1 2
1 3
2 4
2 5
3 6
3 7
3 8
0 3 111100101
*/
總結
以上是生活随笔 為你收集整理的BZOJ 3730: 震波 动态树分治 线段树 lca 的全部內容，希望文章能夠幫你解決所遇到的問題。
如果覺得生活随笔 網站內容還不錯，歡迎將生活随笔 推薦給好友。
