传送门

总算会树形莫队了……

上次听说树形莫队是给树分块,实在看不懂。然后用括号序列的方法做总算能弄明白了

先说一下什么是括号序列,就是在$dfs$的时候,进入的时候记录一下,出去的时候也记录一下

拿样例为例,它的括号序列就是$12443321$

那么我们扩展区间对答案的贡献是可以$O(1)$计算的

假设扩展出的点的颜色是$c$,那么变化量为$val_c*worth_{cnt_c+1}$

因为括号序列它在扩展的时候会把子树里的扫两遍,那么就可以把子树中的答案去掉了

怎么去掉呢?记录一个$vis$然后每次扫到的时候异或一下,看看是否被扫过就行了

然而有几个问题

第一,$LCA$不是路径端点的话不会被记入答案

考虑上面的括号序列,如果要从$4$到$3$那么$2$是不会被计入答案的

第二,起点不是$LCA$的话不会被记入答案

考虑上面,$u$被算了两次,刚好把自己给减掉了

所以上面的两种情况要特判

然后因为是带修莫队,所以再加上时间这一维

还有这题细节挺多的……我因为跳时间轴的时候一个地方写错WA了好久……

 //minamoto
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define ll long long
using namespace std;
#define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
char buf[<<],*p1=buf,*p2=buf;
inline int read(){
#define num ch-'0'
char ch;bool flag=;int res;
while(!isdigit(ch=getc()))
(ch=='-')&&(flag=true);
for(res=num;isdigit(ch=getc());res=res*+num);
(flag)&&(res=-res);
#undef num
return res;
}
char sr[<<],z[];int C=-,Z;
inline void Ot(){fwrite(sr,,C+,stdout),C=-;}
inline void print(ll x){
if(C><<)Ot();if(x<)sr[++C]=,x=-x;
while(z[++Z]=x%+,x/=);
while(sr[++C]=z[Z],--Z);sr[++C]='\n';
}
const int N=;
struct query{
int x,y,l,r,id,t;
}q[N];int q1;
struct change{int x,c;}t[N];int q2;
int n,m,L,R,T,s,Q;ll now,ans[N];
int ver[N<<],edge[N<<],head[N],Next[N<<],top[N],dfn[N],dep[N],sz[N],son[N],lis[N<<],fa[N],tot,num;
int vis[N],col[N],wor[N],cnt[N],val[N];
inline bool cmp(const query a,const query b){
if(a.l^b.l) return a.l<b.l;
if(a.r^b.r) return a.l&?a.r<b.r:a.r>b.r;
return (a.l^a.r)&?a.t<b.t:a.t>b.t;
}
inline void add(int u,int v){
ver[++tot]=v,Next[tot]=head[u],head[u]=tot;
ver[++tot]=u,Next[tot]=head[v],head[v]=tot;
}
void dfs(int u){
dep[u]=dep[fa[u]]+,sz[u]=;
for(int i=head[u];i;i=Next[i]){
int v=ver[i];
if(v!=fa[u]){
fa[v]=u,dfs(v),sz[u]+=sz[v];
if(sz[v]>sz[son[u]]) son[u]=v;
}
}
}
void dfs(int u,int t){
top[u]=t,lis[dfn[u]=++num]=u;
if(son[u]) dfs(son[u],t);
for(int i=head[u];i;i=Next[i]){
if(ver[i]!=fa[u]&&ver[i]!=son[u]) dfs(ver[i],ver[i]);
}
lis[++num]=u;
}
inline int LCA(int u,int v){
while(top[u]!=top[v]){
if(dep[top[u]]<dep[top[v]]) swap(u,v);
u=fa[top[u]];
}
return dep[u]<dep[v]?u:v;
}
inline void sol(int x){
int c=col[x];
(vis[x]^=)?now+=(ll)wor[++cnt[c]]*val[c]:now-=(ll)wor[cnt[c]--]*val[c];
}
inline void mdy(int i){
int u=t[i].x,x=t[i].c,y=col[u];
if(vis[u]) now+=(ll)wor[++cnt[x]]*val[x]-(ll)wor[cnt[y]--]*val[y];
t[i].c=y,col[u]=x;
}
int main(){
n=read(),m=read(),Q=read();
for(int i=;i<=m;++i) val[i]=read();
for(int i=;i<=n;++i) wor[i]=read();
for(int i=;i<n;++i){
int u=read(),v=read();add(u,v);
}
for(int i=;i<=n;++i) col[i]=read();
dfs(),dfs(,);
while(Q--){
int opt=read(),x=read(),y=read();
if(opt){if(dfn[x]>dfn[y]) swap(x,y);
q[++q1]=(query){x,y,dfn[x],dfn[y],q1,q2};
}
else t[++q2]=(change){x,y};
}
s=pow(n,q2?2.0/:1.0/);
for(int i=;i<=n;++i) q[i].l/=s,q[i].r/=s;
sort(q+,q++q1,cmp),L=dfn[q[].l],R=L-,T=;
for(int i=;i<=q1;++i){
int u=q[i].x,v=q[i].y;
int l=dfn[u],r=dfn[v],g=q[i].t;
while(L>l) sol(lis[--L]);
while(R<r) sol(lis[++R]);
while(L<l) sol(lis[L++]);
while(R>r) sol(lis[R--]);
while(T<g) mdy(++T);
while(T>g) mdy(T--);
//LCA不是路径端点的话是不会计算的
//出发点u如果不是LCA也是不会计算的
//这两个都要特判
int p=LCA(u,v);
if(u!=p){sol(u);if(v!=p) sol(p);}
ans[q[i].id]=now;
if(u!=p){sol(u);if(v!=p) sol(p);}
}
for(int i=;i<=q1;++i) print(ans[i]);
Ot();
return ;
}

最新文章

  1. 基于nutz框架理解Ioc容器
  2. SLP alpha 阶段总结
  3. 解读(function($){...})(jQuery)
  4. 使用FFmpeg解码H264-2016.01.14
  5. 饿了么 ---Java面试
  6. Storm集群中执行的各种组件及其并行
  7. 建立Ftp站点
  8. hdu_5213_Lucky(莫队算法+容斥定理)
  9. 查看centos系统位数和强制关闭yum
  10. Unity 案例
  11. pandas处理丢失数据-【老鱼学pandas】
  12. python 爬虫爬取内容时, \xa0 、 \u3000 的含义
  13. 关于php,python,javascript文件或者模块导入引入的区别和联系
  14. IDEA之HttpServletRequest之报错解决方案
  15. 设计模式(二)策略模式(Strategy)
  16. 20170906xlVBA_GetEMailFromDocument
  17. 检查mysql是否运行
  18. python中的&amp;&amp;及||
  19. UsernameToken 【转】
  20. iOS_文件上传进度条的实现思路-AFNettworking

热门文章

  1. Mysql 5.6 MHA (gtid) on Kylin
  2. springmvc 在页面跳转之后 引入文件的路径前面加上了 controller 的映射名
  3. sqlserver 使用维护计划备份
  4. C语言学习笔记--条件编译
  5. DAY11-MYSQL视图、触发器、事务、存储过程、函数
  6. @Value在Controller中取值
  7. springmvc配置式开发下的视图解析器
  8. PLSQL语法
  9. 使用HttpClient进行Get通信
  10. 算法Sedgewick第四版-第1章基础-019一Scanner的用法