震波

Time Limit: 15 Sec  Memory Limit: 256 MB
Submit: 1573  Solved: 358
[Submit][Status][Discuss]

Description

在一片土地上有N个城市,通过N-1条无向边互相连接,形成一棵树的结构,相邻两个城市的距离为1,其中第i个城市的价值为value[i]。
不幸的是,这片土地常常发生地震,并且随着时代的发展,城市的价值也往往会发生变动。
接下来你需要在线处理M次操作:
0 x k 表示发生了一次地震,震中城市为x,影响范围为k,所有与x距离不超过k的城市都将受到影响,该次地震造成的经济损失为所有受影响城市的价值和。
1 x y 表示第x个城市的价值变成了y。
为了体现程序的在线性,操作中的x、y、k都需要异或你程序上一次的输出来解密,如果之前没有输出,则默认上一次的输出为0。

Input

第一行包含两个正整数N和M。
第二行包含N个正整数,第i个数表示value[i]。
接下来N-1行,每行包含两个正整数u、v,表示u和v之间有一条无向边。
接下来M行,每行包含三个数,表示M次操作。

Output

包含若干行,对于每个询问输出一行一个正整数表示答案。

Sample Input

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 1

Sample Output

11100101

HINT

1<=N,M<=100000

1<=u,v,x<=N

1<=value[i],y<=10000

0<=k<=N-1

Source

题解:先建立点分树

ra :表示该节点点分树子树内到其距离为某值的所有节点的权值和

rb:示该节点点分树子树内到其父亲节点距离为某值的所有节点的权值和

另一个问题:如何快速求出两点间距离?使用RMQLCA,预处理后 O(1) 查询LCA深度即可。

然后就是类似是容斥吧,维护即可。

修改的话就是计算delta,来解决。

每次不管怎么样,log层,每层的线段树上解决,

n log n是极限,所以可以认为是 n log^2n的时间复杂度。

然后就可以解决了,线段树的基本操作吧。

 #include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<cstdio> #define N 100007
using namespace std;
inline int read()
{
int x=,f=;char ch=getchar();
while(ch<''||ch>''){if (ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<=''){x=(x<<)+(x<<)+ch-'';ch=getchar();}
return x*f;
} int n,m,lastans,S,tot,xj;
int cnt,hed[N],rea[N*],nxt[N*];
int pos[N],mn[][N*],dep[N];
int value[N],vis[N],sz[N],fa[N],ra[N],rb[N],rt,f[N],Log[N*];
int ls[N*],rs[N*],sum[N*],tp; void add(int u,int v)
{
nxt[++cnt]=hed[u];
hed[u]=cnt;
rea[cnt]=v;
}
void get_root(int u,int fa)
{
sz[u]=,f[u]=;
for (int i=hed[u];i!=-;i=nxt[i])
{
int v=rea[i];
if (vis[v]||v==fa)continue;
get_root(v,u),sz[u]+=sz[v],f[u]=max(f[u],sz[v]);
}
f[u]=max(f[u],S-sz[u]);
if (f[u]<f[rt])rt=u;
}
inline int dis(int x,int y)
{
int t=dep[x]+dep[y],k;
x=pos[x],y=pos[y];
if (x>y)swap(x,y);
k=Log[y-x+];
return t-(min(mn[k][x],mn[k][y-(<<k)+])<<);
}
void update(int p,int a,int l,int r,int &x)
{
if (!x)x=++xj;
sum[x]+=a;
if (l==r)return;
int mid=(l+r)>>;
if (p<=mid)update(p,a,l,mid,ls[x]);
else update(p,a,mid+,r,rs[x]);
}
int query(int p,int l,int r,int x)
{
if (!x) return ;
if (l==r) return sum[x];
int mid=(l+r)>>;
if (p<=mid)return query(p,l,mid,ls[x]);
else return query(p,mid+,r,rs[x])+sum[ls[x]];
}
void calc(int u,int par,int p)
{
update(dis(u,p),value[u],,n-,ra[p]);
if (fa[p])update(dis(u,fa[p]),value[u],,n-,rb[p]);
for (int i=hed[u];i!=-;i=nxt[i])
{
int v=rea[i];
if (vis[v]||v==par) continue;
calc(v,u,p);
}
}
void divide(int u)
{
calc(u,,u),vis[u]=;int yS=S;
for (int i=hed[u];i!=-;i=nxt[i])
{
int v=rea[i];
if (vis[v])continue;
rt=;
if (sz[v]>sz[u])S=yS-sz[u];
else S=sz[v];
get_root(v,),fa[rt]=u;
divide(rt);
}
}
inline void modify(int x,int y)
{
for (int i=x;i!=;i=fa[i])
{
update(dis(x,i),y-value[x],,n-,ra[i]);
if (fa[i])update(dis(x,fa[i]),y-value[x],,n-,rb[i]);
}
value[x]=y;
}
inline int solve(int x,int p)
{
int res=;
for (int i=x;i!=;i=fa[i])
{
if (dis(x,i)<=p)res+=query(p-dis(x,i),,n-,ra[i]);
if (fa[i]&&dis(x,fa[i])<=p)res-=query(p-dis(x,fa[i]),,n-,rb[i]);
}
return res;
}
void dfs(int u,int par)
{
pos[u]=++tot,mn[][tot]=dep[u];
for (int i=hed[u];i!=-;i=nxt[i])
{
int v=rea[i];
if (v==par) continue;
dep[v]=dep[u]+,dfs(v,u),mn[][++tot]=dep[u];
}
}
int main()
{
memset(hed,-,sizeof(hed)); n=read(),m=read();
for (int i=;i<=n;i++)
value[i]=read();
for (int i=;i<n;i++)
{
int x=read(),y=read();
add(x,y),add(y,x);
}
dfs(,);
for (int i=;i<=tot;i++)Log[i]=Log[i>>]+;
for (int i=;(<<i)<=tot;i++)
for (int j=;j<=tot-(<<i)+;j++)
mn[i][j]=min(mn[i-][j],mn[i-][j+(<<(i-))]);
f[]=<<,S=n;
get_root(,),divide(rt);
while(m--)
{
int flag=read(),x=read()^lastans,y=read()^lastans;
if (flag)modify(x,y);
else printf("%d\n",lastans=solve(x,y));
}
}

最新文章

  1. track message forwards, avoiding request loops, and identifying the protocol capabilities of all senders along the request/response chain
  2. struts2 配置拦截器
  3. SSAS Cube 维度成员关系Rigid 和 Flexible
  4. hdu 4026 2011上海赛区网络赛F TSP ****
  5. 点击文字label同时选中checkbox radio
  6. Excel01-不同的单元格输入同一数据
  7. CentOS6.X 系统安装图文教程,超详细
  8. java学习之Java中JDK,JRE和JVM之间的关系(转载)
  9. Aprior算法
  10. 第一节: dingo/API 最新版 V2.0 之安装讲解(连载)
  11. 【学习笔记】JS知识点整理
  12. 关于SASS
  13. java第六次课后作业
  14. SharePoint 2013 新特性 (三) 破改式 &mdash;&mdash; 设计管理器的使用 [2.HTML变身模板页]
  15. Chrome 开发者工具
  16. Generative Adversarial Nets[Vanilla]
  17. Ubuntu16.04下使用rdesktop命令远程连接windows机器
  18. JavaScript RegExp Object 正则表达式入门
  19. dubbo2.5.3升级到dobbo2.8.4(dubbox) jar
  20. Java网络编程一:基础知识详解

热门文章

  1. DaemonSet 典型应用场景【转】
  2. Maven搭建Struts2+Spring3+Hibernate4框架
  3. Bootstrap历练实例:导航内的下拉菜单
  4. mysql crash cource 书中实例
  5. tensorflow目标检测API之训练自己的数据集
  6. 【思维题 线段树】cf446C. DZY Loves Fibonacci Numbers
  7. 分割catalina.out 每天生成一个文件
  8. paper:synthesizable finit state machine design techniques using the new systemverilog 3.0 enhancements之output encoded style with registered outputs(Good style)
  9. phpExcel使用方法二
  10. Golang map并发 读写锁