【BZOJ 3729】3729: Gty的游戏 (Splay维护dfs序+博弈)

未经博主同意不得转载

 

 

3729: Gty的游戏

Time Limit: 20 Sec  Memory Limit: 128 MB
Submit: 448  Solved: 150

Description

某一天gty在与他的妹子玩游戏。
妹子提出一个游戏,给定一棵有根树,每个节点有一些石子,每次可以将不多于L的石子移动到父节点,询问
将某个节点的子树中的石子移动到这个节点先手是否有必胜策略。
gty很快计算出了策略。
但gty的妹子十分机智,她决定修改某个节点的石子或加入某个新节点。
gty不忍心打击妹子,所以他将这个问题交给了你。
另外由于gty十分绅士,所以他将先手让给了妹子。

Input

第一行两个数字,n和L,n<=5*10^4,L<=10^9
第二行n个数字,表示每个节点初始石子数。
接下来n-1行,每行两个整数u和v,表示有一条从u到v的边。
接下来一行一个数m,表示m组操作。
接下来m行,每行第一个数字表示操作类型
若为1,后跟一个数字v,表示询问在v的子树中做游戏先手是否必胜。
若为2,后跟两个数字x,y表示将节点x的石子数修改为y。
若为3,后跟三个数字u,v,x,表示为u节点添加一个儿子v,初始石子数为x。
在任意时刻,节点数不超过5*10^4。

Output

对于每个询问,若先手必胜,输出"MeiZ",否则输出"GTY"。
另,数据进行了强制在线处理,对于m组操作,除了类型名以外,都需要异或之前回答为"MeiZ"的个数。

Sample Input

2 1000
0 0
1 2
1
1 1

Sample Output

GTY

HINT

Source

 

【分析】

  膜奥爷爷啦~~

  好吧,就是首先,阶梯尼姆,很明显吧。就是把距离根节点为奇数层的异或起来就好了。

  那就是差不多维护子树的异或和,但是树是动态的。【表示动态树我真的很垃圾,splay忘光,LCT不会,前面做的题都是离线的【如果离线就很快就会做了

     所以接下来要学学动态树了。。

  用splay维护dfs序,和之前差不多嘛,一棵树的子树的dfs序记录了st和ed之后,查询区间[st,ed]就行了。

  splay就是能求出键值在某范围里面的东西【具体维护什么,和啊,最大值啊,都是你决定的】,但实际splay树上维护的是splay树子树上的东西。

  这时,你只要把st splay到跟,ed splay到根的右儿子,那么ed的做儿子表示的区间就是[st+1,ed-1]你直接问它的子树就好了。

  实际上,并不需要实际的dfs序,只要你按顺序插入的,那splay树就是有序的,你记录每个点的st和ed,到时找前驱后继就好了。

  可以选择看图:

  

  左图是原树,右图是splay树。红色是值。splay树的中序遍历就是dfs序,splay树上红色的两个维护的是子树内dep为奇和偶的点的异或和。

  有意义的点我打在左端点,右端点的值均为0。

  splay的时候upd一下把左右儿子的值跟自己异或一下就能维护那个异或和了。

 

  1 #include<cstdio>
  2 #include<cstdlib>
  3 #include<cstring>
  4 #include<iostream>
  5 #include<algorithm>
  6 #include<map>
  7 using namespace std;
  8 #define Maxn 400010
  9 
 10 int n,Mod;
 11 
 12 int w[Maxn];
 13 int nl[Maxn],nr[Maxn];
 14 map<int,int> z;
 15 
 16 struct node
 17 {
 18     int x,y,next;
 19 }t[Maxn*2];
 20 int len,first[Maxn];
 21 void ins(int x,int y)
 22 {
 23     t[++len].x=x;t[len].y=y;
 24     t[len].next=first[x];first[x]=len;
 25 }
 26 
 27 struct sp
 28 {
 29     int son[2],fa,d,a1,a2;
 30     bool dep;
 31 }tr[Maxn*2];int tot;
 32 
 33 void upd(int x)
 34 {
 35     int lc=tr[x].son[0],rc=tr[x].son[1];
 36     tr[x].a1=tr[lc].a1^tr[rc].a1;
 37     tr[x].a2=tr[lc].a2^tr[rc].a2;
 38     if(!tr[x].dep) tr[x].a2^=tr[x].d;
 39     else tr[x].a1^=tr[x].d;
 40 }
 41 void rot(int x)
 42 {
 43     int fa=tr[x].fa,yy=tr[fa].fa;
 44     int w=tr[fa].son[0]==x?1:0;
 45     tr[fa].son[1-w]=tr[x].son[w];
 46     if(tr[x].son[w]) tr[tr[x].son[w]].fa=fa;
 47     if(yy)
 48     {
 49         if(tr[yy].son[0]==fa) tr[yy].son[0]=x;
 50         else tr[yy].son[1]=x;
 51     }tr[x].fa=yy;
 52     tr[x].son[w]=fa;tr[fa].fa=x;
 53     upd(fa);upd(x);
 54 }
 55 void splay(int x,int nf)
 56 {
 57     while(tr[x].fa!=nf)
 58     {
 59         int fa=tr[x].fa,yy=tr[fa].fa;
 60         if(yy==nf) rot(x);
 61         else
 62         {
 63             if((tr[yy].son[0]==fa)==(tr[fa].son[0]==x)) {rot(fa);rot(x);}
 64             else {rot(x);rot(x);}
 65         }
 66     }
 67 }
 68 int Lower(int x)
 69 {
 70     splay(x,0);
 71     x=tr[x].son[1];
 72     while(tr[x].son[0]) x=tr[x].son[0];
 73     return x;
 74 }
 75 void insert(int fa,int x)
 76 {
 77     int r=Lower(fa);
 78     splay(fa,0);splay(r,fa);
 79     tr[r].son[0]=x;
 80     tr[x].fa=r;
 81     upd(r);upd(fa);
 82 }
 83 void dfs(int x,int fa)
 84 {
 85     tr[x].dep=!tr[fa].dep;
 86     if(fa)
 87     {
 88         tr[nl[x]].son[1]=nr[x];tr[nr[x]].fa=nl[x];
 89         tr[nl[x]].d=w[x];upd(nr[x]);upd(nl[x]);
 90         insert(fa,x);
 91     }
 92     for(int i=first[x];i;i=t[i].next) if(t[i].y!=fa)
 93     {
 94         int y=t[i].y;
 95         dfs(y,x);
 96     }
 97 }
 98 
 99 int main()
100 {
101     scanf("%d%d",&n,&Mod);Mod++;
102     for(int i=1;i<=n;i++) {scanf("%d",&w[i]);z[i]=i;w[i]%=Mod;}
103     for(int i=1;i<n;i++)
104     {
105         int x,y;
106         scanf("%d%d",&x,&y);
107         ins(x,y);ins(y,x);
108     }
109     for(int i=1;i<=n;i++) nl[i]=i,nr[i]=i+n;tot=n*2;
110     tr[0].dep=0;tr[1].son[1]=n+1;
111     tr[1].d=w[1];tr[n+1].fa=1;upd(n+1);upd(1);
112     dfs(1,0);
113     int q,nw=0;
114     scanf("%d",&q);
115     while(q--)
116     {
117         int opt,ans;
118         int x,y,k;
119         scanf("%d",&opt);
120         if(opt==1)
121         {
122             scanf("%d",&x);
123             x^=nw;
124             x=z[x];
125             splay(nl[x],0);splay(nr[x],nl[x]);
126             if(!tr[x].dep) ans=tr[tr[nr[x]].son[0]].a1;
127             else ans=tr[tr[nr[x]].son[0]].a2;
128             if(ans==0) printf("GTY\n");
129             else {printf("MeiZ\n");nw++;}
130         }
131         else if(opt==2)
132         {
133             scanf("%d%d",&x,&y);
134             x^=nw;y^=nw;
135             y%=Mod;
136             x=z[x];
137             splay(nl[x],0);
138             if(!tr[nl[x]].dep) tr[nl[x]].a2^=y^tr[nl[x]].d;
139             else tr[nl[x]].a1^=y^tr[nl[x]].d;
140             tr[nl[x]].d=y;
141         }
142         else
143         {
144             scanf("%d%d%d",&x,&y,&k);
145             x^=nw;y^=nw;k^=nw;
146             k%=Mod;x=z[x];z[y]=y=++n;
147             nl[y]=++tot;nr[y]=++tot;
148             tr[nl[y]].dep=!tr[nl[x]].dep;
149             tr[nl[y]].d=k;tr[nr[y]].d=0;
150             tr[nr[y]].fa=nl[y];
151             tr[nl[y]].fa=nl[x];upd(nl[y]);upd(nr[y]);
152             insert(nl[x],nl[y]);
153         }
154     }
155     return 0;
156 }
View Code

 

【AC了还是很兴奋的。。。以后要多做点dfs序,splay之类的。。

 

2017-03-30 16:39:42

posted @ 2017-03-30 15:12  konjak魔芋  阅读(266)  评论(0编辑  收藏  举报