构建：像线性的莫队那样，依旧是按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. http调试工具,linux调试工具

    charles Linux下Web性能压力测试工具http_load linux 下的socket 调试工具 netcat Linux下四款Web服务器压力测试工具(http_load.webbenc ...

  2. HDOJ(HDU) 1555 How many days?(水题)

    Problem Description 8600的手机每天消费1元,每消费K元就可以获赠1元,一开始8600有M元,问最多可以用多少天? Input 输入包括多个测试实例.每个测试实例包括2个整数M, ...

  3. cf500A New Year Transportation

    A. New Year Transportation time limit per test 2 seconds memory limit per test 256 megabytes input s ...

  4. ACM2114_S[I](1^3+2^3+3^3）

    #include<iostream> using namespace std; int main() { __int64 n,m,i,j,sum; while(cin>>n) ...

  5. 广州Uber优步司机奖励政策（12月28日到1月3日）

    滴快车单单2.5倍,注册地址:http://www.udache.com/ 如何注册Uber司机(全国版最新最详细注册流程)/月入2万/不用抢单:http://www.cnblogs.com/mfry ...

  6. redis报错

    网站登录异常,redis数据不能写!解决方法汇总! redis---flushdb  ###提示如下错误    ###flushall              清空说有数据,所有库 (error) ...

  7. 读取xml到DataSet中去

    XML如下: <?xml version="1.0" encoding="utf-8" ?> <Config> <System&g ...

  8. EF并发性能文章

    http://www.cnblogs.com/farb/p/ConcurrencyAndTransctionManagement.html

  9. (转)How to renew your Apple Push Notification Push SSL Certificate

    转自:https://blog.serverdensity.com/how-to-renew-your-apple-push-notification-push-ssl-certificate/ It ...

  10. 解析Xcode把应用程序打包成ipa---解决打包完新版本itunes提示不是有效应用程序的问题

    Xcode把应用程序打包成ipa是本文要介绍的内容,不多说,先俩看内容.注意:本方法需要先制作假凭证编译于项目中,否则产生的ipa还是无法于iPhone中运行. 制作方法请参考: http://blo ...