bzoj2819 Nim

http://www.elijahqi.win/archives/2970
Description
著名游戏设计师vfleaking，最近迷上了Nim。普通的Nim游戏为：两个人进行游戏，N堆石子，每回合可以取其中某一堆的任意多个，可以取完，但不可以不取。谁不能取谁输。这个游戏是有必胜策略的。于是vfleaking决定写一个玩Nim游戏的平台来坑玩家。
为了设计漂亮一点的初始局面，vfleaking用以下方式来找灵感：拿出很多石子，把它们聚成一堆一堆的，对每一堆编号1,2,3,4,…n,在堆与堆间连边，没有自环与重边，从任意堆到任意堆都只有唯一一条路径可到达。然后他不停地进行如下操作：
1.随机选两个堆v,u，询问若在v到u间的路径上的石子堆中玩Nim游戏，是否有必胜策略，如果有，vfleaking将会考虑将这些石子堆作为初始局面之一，用来坑玩家。
2.把堆v中的石子数变为k。
由于vfleaking太懒了，他懒得自己动手了。请写个程序帮帮他吧。
Input
第一行一个数n，表示有多少堆石子。
接下来的一行，第i个数表示第i堆里有多少石子。
接下来n-1行，每行两个数v,u，代表v,u间有一条边直接相连。
接下来一个数q，代表操作的个数。
接下来q行，每行开始有一个字符：
如果是Q，那么后面有两个数v,u，询问若在v到u间的路径上的石子堆中玩Nim游戏，是否有必胜策略。
如果是C，那么后面有两个数v,k，代表把堆v中的石子数变为k。
对于100%的数据：
1≤N≤500000, 1≤Q≤500000, 0≤任何时候每堆石子的个数≤32767
其中有30%的数据：
石子堆组成了一条链，这3个点会导致你DFS时爆栈（也许你不用DFS？）。其它的数据DFS目测不会爆。
注意：石子数的范围是0到INT_MAX
Output
对于每个Q，输出一行Yes或No，代表对询问的回答。
Sample Input
【样例输入】
5
1 3 5 2 5
1 5
3 5
2 5
1 4
6
Q 1 2
Q 3 5
C 3 7
Q 1 2
Q 2 4
Q 5 3
Sample Output
Yes
No
Yes
Yes
Yes
HINT
Source
湖北省队互测
dfs序+树状数组可以解决好多问题
同理 nim游戏的解决套路 然后考虑一个点修改只会影响到他的子树区间即可解决 那么直接树状数组 维护每个点他到根的信息 即可 到时候和lca异或一下就消去了

#include<cstdio>
#include<algorithm>
#define N 500050
#define rg register
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=gc();}
    while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=gc();
    return x*f;
}
struct node{
    int y,next;
}data[N<<1];
int fa[N],mn[N<<1][22],num,tot,n,Log[N<<1],h[N],in[N],pos[N<<1],out[N],dep[N];
int a[N],s[N];
inline void dfs(int x){
    in[x]=++num;mn[++tot][0]=x;pos[x]=tot;
    for (rg int i=h[x];i;i=data[i].next){
        int y=data[i].y;if (y==fa[x]) continue;
        dep[y]=dep[x]+1;fa[y]=x;dfs(y);mn[++tot][0]=x;
    }out[x]=num;
}
inline void init(){
    for (rg int j=1;j<=Log[tot];++j)
        for (rg int i=1;i-1<=tot-(1<<j);++i)
            mn[i][j]=dep[mn[i][j-1]]<dep[mn[i+(1<<j-1)][j-1]]?mn[i][j-1]:mn[i+(1<<j-1)][j-1];
}
inline int lca(int x,int y){
    x=pos[x],y=pos[y];if(x>y) swap(x,y);int t=Log[y-x+1];
    return dep[mn[x][t]]<dep[mn[y-(1<<t)+1][t]]?mn[x][t]:mn[y-(1<<t)+1][t];
}
inline void add(int x,int v){
    for (rg int i=x;i<=n;i+=i&(-i)) s[i]^=v;
}
inline int query(int x){
    int tmp=0;
    for (rg int i=x;i;i-=i&(-i)) tmp^=s[i];return tmp;
}
int main(){
//  freopen("bzoj2819.in","r",stdin);
    //freopen("bzoj2819.out","w",stdout);
    n=read();Log[0]=-1;for (rg int i=1;i<=n<<1;++i) Log[i]=Log[i>>1]+1;
    for (rg int i=1;i<=n;++i) a[i]=read();
    for (rg int i=1;i<n;++i){
        int x=read(),y=read();
        data[++num].y=y;data[num].next=h[x];h[x]=num;
        data[++num].y=x;data[num].next=h[y];h[y]=num;
    }num=0;dfs(1);int Q=read();init();
    for (rg int i=1;i<=n;++i) add(in[i],a[i]),add(out[i]+1,a[i]);
    while(Q--){
        char ch=gc();while(ch!='Q'&&ch!='C') ch=gc();
        if (ch=='Q'){
            int x=read(),y=read();int t=lca(x,y);
            int ans=query(in[x])^query(in[y])^a[t];
            if (ans) puts("Yes");else puts("No");continue;
        }
        int x=read(),k=read();add(in[x],a[x]);add(out[x]+1,a[x]);
        a[x]=k;add(in[x],k);add(out[x]+1,k);
    }
    return 0;
}