构建：像线性的莫队那样，依旧是按sqrt(n)为一块分块。

 int dfs(int x){

     int size=;

     dfn[x]=++ind;

     for (int i=;i<=;i++)

      if (bin[i]<=deep[x])

      fa[x][i]=fa[fa[x][i-]][i-];

      else break;

     for (int i=first[x];i;i=next[i]){

         int pur=go[i];

         if (pur!=fa[x][]){

             fa[pur][]=x;

             deep[pur]=deep[x]+;

             size+=dfs(pur);

             if (size>=blo){

                 blonum++;

                 for (int k=;k<=size;k++)

                  belong[st[top--]]=blonum;

                 size=;

             }

         }

     }

     st[++top]=x;

     return size+;

 }

然后呢，我们可以发现一些树上莫队的性质：

用S(v, u)代表 v到u的路径上的结点的集合。
用root来代表根结点，用lca(v, u)来代表v、u的最近公共祖先。
那么
S(v, u) = S(root, v) xor S(root, u) xor lca(v, u)
其中xor是集合的对称差。
简单来说就是节点出现两次消掉。
lca很讨厌，于是再定义
T(v, u) = S(root, v) xor S(root, u)
观察将curV移动到targetV前后T(curV, curU)变化：
T(curV, curU) = S(root, curV) xor S(root, curU)
T(targetV, curU) = S(root, targetV) xor S(root, curU)
取对称差：
T(curV, curU) xor T(targetV, curU)= (S(root, curV) xor S(root, curU)) xor (S(root, targetV) xor S(root, curU))
由于对称差的交换律、结合律：
T(curV, curU) xor T(targetV, curU)= S(root, curV) xorS(root, targetV)
两边同时xor T(curV, curU)：
T(targetV, curU)= T(curV, curU) xor S(root, curV) xor S(root, targetV)
即
T(targetV, curU)= T(curV, curU) xor T(curV, targetV)
也就是说，更新的时候，xor T(curV, targetV)就行了。
即，对curV到targetV路径(除开lca(curV, targetV))上的结点，将它们的存在性取反即可。
 
按照以上方式，我们就可以每次从u1走到u2,v1走到v2，然后要求query的时候对他们的公共祖先的存在性取反，然后求完答案后再取反回去。
 
bzoj3757 苹果树
 #include <cstdio>

 #include <iostream>

 #include <cstring>

 #include <algorithm>

 #include <string>

 #include <queue>

 #include <stack>

 #include <cmath>

 #include <map>

 #include <vector>

 #include <functional>

 #include <ctime>

 #include <cstdlib>

 #include <sstream>

 #include <set>

 #include <deque>

 #define Rep(i, l, r) for (int i = l; i <= r; ++i)

 #define Req(i, l, r) for (int i = l; i >= r; --i)

 #define N 100005

 int tot,go[N],next[N],first[N],dfn[N],belong[N];

 int st[N],res[N],n,m,x,y,fa[N][],bin[],col[N],ans,ind,deep[N],blo,blonum;

 int top,pd[N],root;

 struct op{

     int u,v,id,a,b;

 }q[N];

 int p[N];

 void insert(int x,int y){

     tot++;

     go[tot]=y;

     next[tot]=first[x];

     first[x]=tot;

 }

 void add(int x,int y){

     insert(x,y);

     insert(y,x);

 }

 bool cmp(op q,op w){

     if (belong[q.u]==belong[w.u]) return dfn[q.v]<dfn[w.v];

     else return belong[q.u]<belong[w.u];

 }

 int dfs(int x){

     int size=;

     dfn[x]=++ind;

     for (int i=;i<=;i++)

      if (bin[i]<=deep[x])

      fa[x][i]=fa[fa[x][i-]][i-];

      else break;

     for (int i=first[x];i;i=next[i]){

         int pur=go[i];

         if (pur!=fa[x][]){

             fa[pur][]=x;

             deep[pur]=deep[x]+;

             size+=dfs(pur);

             if (size>=blo){

                 blonum++;

                 for (int k=;k<=size;k++)

                  belong[st[top--]]=blonum;

                 size=;

             }

         }

     }

     st[++top]=x;

     return size+;

 }

 void reverse(int x){

     if (pd[x]) {

         p[col[x]]--;

         pd[x]=;

         if (p[col[x]]==) ans--;

     }

     else{

         p[col[x]]++;

         pd[x]=;

         if (p[col[x]]==) ans++;

     }

 }

 void solve(int u,int v){

     while (u!=v){

         if (deep[u]>deep[v]) reverse(u),u=fa[u][];

         else reverse(v),v=fa[v][];

     }

 }

 int lca(int x,int y){

     if (deep[x]<deep[y]) std::swap(x,y);

     int t=deep[x]-deep[y];

     for (int i=;i>=;i--)

      if (t&bin[i]){

         x=fa[x][i];

      }

     for (int i=;i>=;i--)

      if (fa[x][i]!=fa[y][i]){

         x=fa[x][i];

         y=fa[y][i];

      }

     if (x==y) return x;

     else return fa[x][];

 }

 int main(){

     bin[]=;

     for (int i=;i<;i++) bin[i]=bin[i-]*;

     scanf("%d%d",&n,&m);

     blo=sqrt(n);

     for (int i=;i<=n;i++)

      scanf("%d",&col[i]);

     for (int i=;i<=n;i++){

         scanf("%d%d",&x,&y);

         if (x==)

             root=y;

         else

         if (y==)

             root=x;

         else

         add(x,y);

     }

     dfs(root);

     blonum++;

     while (top) belong[st[top--]]=blonum;

     for (int i=;i<=m;i++){

         scanf("%d%d",&q[i].u,&q[i].v);

         if (dfn[q[i].u]>dfn[q[i].v]) std::swap(q[i].v,q[i].u);

         scanf("%d%d",&q[i].a,&q[i].b);

         q[i].id=i;

     }

     std::sort(q+,q++m,cmp);

     int t=lca(q[].u,q[].v);

     solve(q[].u,q[].v);

     reverse(t);

     res[q[].id]=ans;

     if (p[q[].a]&&p[q[].b]&&q[].a!=q[].b) res[q[].id]--;

     reverse(t);

     for (int i=;i<=m;i++){

         solve(q[i-].u,q[i].u);

         solve(q[i-].v,q[i].v);

         t=lca(q[i].u,q[i].v);

         reverse(t);

         res[q[i].id]=ans;

         if (p[q[i].a]&&p[q[i].b]&&q[i].a!=q[i].b) res[q[i].id]--;

         reverse(t);

     }

     for (int i=;i<=m;i++)

      printf("%d\n",res[i]);

 }
bzoj3052 WC糖果公园
 #include <cstdio>

 #include <iostream>

 #include <cstring>

 #include <algorithm>

 #include <string>

 #include <queue>

 #include <stack>

 #include <cmath>

 #include <map>

 #include <vector>

 #include <functional>

 #include <ctime>

 #include <cstdlib>

 #include <sstream>

 #include <set>

 #include <deque>

 #define Rep(i, l, r) for (int i = l; i <= r; ++i)

 #define Req(i, l, r) for (int i = l; i >= r; --i)

 #define N 200005

 long long res[N];

 long long pre[N],col[N],v[N],w[N],ans;

 int belong[N],go[N],tot,first[N],next[N];

 int deep[N],bin[],fa[N][],st[N],top,x,y,n,m,c1,c2,dfn[N],ind,blo,blonum;

 int pd[N],p[N],qe,ty;

 struct op{

     int x,y,id,t;

     long long pre;

 }c[N],b[N];

 bool cmp(op a,op b){

     if (belong[a.x]==belong[b.x]&&belong[a.y]==belong[b.y]) return a.t<b.t;

     else

     if (belong[a.x]==belong[b.x]) return belong[a.y]<belong[b.y];

     return belong[a.x]<belong[b.x];

 }

 void insert(int x,int y){

     tot++;

     go[tot]=y;

     next[tot]=first[x];

     first[x]=tot;

 }

 void add(int x,int y){

     insert(x,y);

     insert(y,x);

 }

 int dfs(int x){

     int size=;

     dfn[x]=++ind;

     for (int i=;i<=;i++)

      if (deep[x]>=bin[i])

       fa[x][i]=fa[fa[x][i-]][i-];

      else break;

     for (int i=first[x];i;i=next[i]){

         int pur=go[i];

         if (pur!=fa[x][]){

             fa[pur][]=x;

             deep[pur]=deep[x]+;

             size+=dfs(pur);

             if (size>=blo){

                 blonum++;

                 for (int j=;j<=size;j++)

                  belong[st[top--]]=blonum;

                 size=;

             }

         }

     }

     st[++top]=x;

     return size+;

 }

 void reverse(int x){

     if (pd[x]){

         ans-=w[p[col[x]]]*v[col[x]];

         p[col[x]]--;

         pd[x]=;

     }

     else{

         pd[x]=;

         p[col[x]]++;

         ans+=w[p[col[x]]]*v[col[x]];

     }

 }

 void change(int x,int y){

     if (pd[x]){

         reverse(x);

         col[x]=y;

         reverse(x);

     }

     else{

         col[x]=y;

     }

 }

 void solve(int x,int y){

     while (x!=y){

         if (deep[x]>deep[y]) reverse(x),x=fa[x][];

         else reverse(y),y=fa[y][];

     }

 }

 int lca(int x,int y){

     if (deep[x]<deep[y]) std::swap(x,y);

     int t=deep[x]-deep[y];

     for (int i=;i>=;i--)

      if (t&bin[i])

       x=fa[x][i];

     for (int i=;i>=;i--)

      if (fa[x][i]!=fa[y][i])

       x=fa[x][i],y=fa[y][i];

     if (x==y) return x;

     else return fa[x][];

 }

 int main(){

     //freopen("tx.txt","r",stdin);

     bin[]=;

     for (int i=;i<=;i++) bin[i]=bin[i-]*;

     scanf("%d%d%d",&n,&m,&qe);

     blo=pow(n,2.0/)*0.5;

     for (int i=;i<=m;i++)

      scanf("%lld",&v[i]);

     for (int i=;i<=n;i++)

      scanf("%lld",&w[i]);

     for (int i=;i<n;i++){

         scanf("%d%d",&x,&y);

         add(x,y);

     }

     for (int i=;i<=n;i++){

      scanf("%lld",&col[i]);

      pre[i]=col[i];

      }

     dfs();

     //blonum++;

     while (top) belong[st[top--]]=blonum;

     for (int i=;i<=qe;i++){

         scanf("%d",&ty);

         if (!ty){

             c1++;

             scanf("%d%d",&c[c1].x,&c[c1].y);

             c[c1].pre=pre[c[c1].x];

             pre[c[c1].x]=c[c1].y;

         }

         else{

             c2++;

             scanf("%d%d",&b[c2].x,&b[c2].y);

             b[c2].t=c1;

             if (dfn[b[c2].x]>dfn[b[c2].y]) std::swap(b[c2].x,b[c2].y);

             b[c2].id=c2;

         }

     }

     std::sort(b+,b++c2,cmp);

     for (int i=;i<=b[].t;i++)

      change(c[i].x,c[i].y);

     solve(b[].x,b[].y);

     int t=lca(b[].x,b[].y);

     reverse(t);

     res[b[].id]=ans;

     reverse(t);

     for (int i=;i<=c2;i++){

         for (int j=b[i-].t+;j<=b[i].t;j++)

          change(c[j].x,c[j].y);

         for (int j=b[i-].t;j>b[i].t;j--)

          change(c[j].x,c[j].pre);

         solve(b[i-].x,b[i].x);

         solve(b[i-].y,b[i].y);

         int t=lca(b[i].x,b[i].y);

         reverse(t);

         res[b[i].id]=ans;

         reverse(t);

     }

     for (int i=;i<=c2;i++)

      printf("%lld\n",res[i]);

 }

树上莫队 wowow的更多相关文章

  1. 【BZOJ 3735】苹果树 树上莫队(树分块+离线莫队+鬼畜的压行)

    2016-05-09 UPD:学习了新的DFS序列分块,然后发现这个东西是战术核导弹?反正比下面的树分块不知道要快到哪里去了 #include<cmath> #include<cst ...

  2. 【BZOJ-3757】苹果树 块状树 + 树上莫队

    3757: 苹果树 Time Limit: 20 Sec  Memory Limit: 256 MBSubmit: 1305  Solved: 503[Submit][Status][Discuss] ...

  3. [BZOJ 3052] [wc2013] 糖果公园 【树上莫队】

    题目链接:BZOJ - 3052 题目分析 这道题就是非常经典的树上莫队了,并且是带修改的莫队. 带修改的莫队:将询问按照 左端点所在的块编号为第一关键字,右端点所在的块为第二关键字,位于第几次修改之 ...

  4. spoj COT2 - Count on a tree II 树上莫队

    题目链接 http://codeforces.com/blog/entry/43230树上莫队从这里学的,  受益匪浅.. #include <iostream> #include < ...

  5. BZOJ 4129: Haruna’s Breakfast [树上莫队 分块]

    传送门 题意: 单点修改,求一条链的mex 分块维护权值,$O(1)$修改$O(S)$求mex...... 带修改树上莫队 #include <iostream> #include < ...

  6. 【WC2013】糖果公园 [树上莫队]

    题意: 一棵树,修改一个点的颜色,询问两点路径上每种颜色的权值$val[c]$*出现次数的权值$cou[w[c]]$的和 sro VFK 树上莫队 按照王室联邦的方法分块,块的大小直径个数有保证,并不 ...

  7. Codeforces 852I Dating 树上莫队

    Dating 随便树上莫队搞一搞就好啦. #include<bits/stdc++.h> #define LL long long #define LD long double #defi ...

  8. 【XSY1642】Another Boring Problem 树上莫队

    题目大意 给你一棵\(n\)个点的树,每个点有一个颜色\(c_i\),每次给你\(x,y,k\),求从\(x\)到\(y\)的路径上出现次数第\(k\)多的颜色的出现次数 \(n,q\leq 1000 ...

  9. SP10707 COT2 - Count on a tree II (树上莫队)

    大概学了下树上莫队, 其实就是在欧拉序上跑莫队, 特判lca即可. #include <iostream> #include <algorithm> #include < ...

随机推荐

  1. 微信、QQ和手机号之间不得不说的故事！

    发文字,发图片,发心情,视频聊天,查看附近的人,微信能干的事情QQ都可以,那么它们有什么区别,我QQ用得好好的为什么要我联系人都导到微信去?我们很早就有了QQ,但是在QQ时代,我们虽然用QQ发消息聊天 ...

  2. 探讨c#中的unchecked是什么意思,起什么作用?

    Checked与Unchecked     对于因为整数类型参与算术操作和类型转换时产生的“溢出异常”——System.OverflowException,在某些算法来讲不算真正的“异常”,相反这种溢 ...

  3. HTML5迷你游戏作验证码

    验证码最常见的是各种变形的字符,因为识别程序进化得越来越聪明,验证码也变得越来越难以识别,给用户造成了很多的麻烦和反感. 已经有很多人尝试过各种改进的验证码,比如动画的验证码,做题目的验证码,要回答问 ...

  4. C++基础学习笔记----第十三课（操作符重载-下）

    本节主要讲使用成员函数重载操作符,包括[],=,(),->四种操作符的重载以及&&和||的问题. 类的成员函数进行操作符重载 基本概念 类的成员函数也可以进行操作符的重载.类的普 ...

  5. FilterDispatcher 的作用(struts2.1.3以前，新版本改了)

    org.apache.struts2.dispatcher.FilterDispatcher是Struts2的主要的Filter,负责四个方面的功能: (1)执行Actions (2)清除Action ...

  6. 企业版IDP的申请及“In House”发布

    企业版IDP,即iOS Development Enterprise Program.注意是$299/Year那种,并不是$99/Year的那种. 这种方式的IDP其最大的好处在于:可以发布“In H ...

  7. CentOS6.3(64位)下安装Oracle11gR2(64)服务器

    安装环境 Linux服务器:Centos6.3 64位 Oracle服务器:Oracle11gR2 64位 系统要求 1.Linux安装Oracle系统要求 系统要求 说明 内存 必须高于1G的物理内 ...

  8. gcc -L -l的使用

    -l参数就是用来指定程序要链接的库,-l参数紧接着就是库名,那么库名跟真正的库文件名有什么关系呢?就拿数学库来说,他的库名是m,他的库文件名是libm.so,很容易看出,把库文件名的头lib和尾.so ...

  9. (转)Facade模式

    Facade模式要求一个子系统的外部与其内部的通信必须通过一个统一的Facade对象进行.Facade模式提供一个高层次的接口,使得子系统更易于使用. 就如同医院的接待员一样,Facade模式的Fac ...

  10. （转）VS无法启动调试：“生成下面的模块时,启用了优化或没有调试信息“

    中调试项目遇到错误提示,Visual Studio 2010(或VS2008或VS2005)启动调试的时候,弹出提示信息: 生成下面的模块时,启用了优化或没有调试信息: C:\WINDOWS\Micr ...