BZOJ4134 ljw和lzr的hack比赛

今年清华集训D1T1-_-

首先我们可以想到n^2的算法，我们对每个点记录一个sg[x]，代表以x为根的子树的游戏的sg函数值，然后我们再对每个点记录rem[x]，代表删掉从x到根的路径后剩余的游戏的sg值的异或和。

容易发现在我们计算x的sg值的时候，x的rem是他所有儿子的sg值的异或和，而对于x的每一个儿子y，y的子树里的所有rem值都会异或rem[x]^sg[y]

然后我们求第一个不出现在子树内的rem集合里的非负整数即是x的sg值（不考虑已经被hack的点）

考虑优化这个算法，我们可以对每个点建立一个权值线段树，维护子树内所有的rem值是否出现过，然后在线段树上爬即可查出这个点的sg值

然后对于一个点，我们如果他没有预先被hack，那么我们把他的rem插入这个点的线段树中，然后把这棵线段树和他所有儿子的线段树进行线段树合并即可

由于我们要有把某个子树里的所有值都异或一个数的操作，那么我们可以把权值线段树的区间设为0~2的幂-1，然后在线段树上打异或标记，如果对应的位是1的话就交换左右儿子即可

最后再dfs一遍树处理每个点的rem判断是否为0即可输出答案

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<map>
#include<set>
#include<bitset>
#include<queue>
#include<stack>
using namespace std;
#define MAXN 200010
#define MAXM 4000010
#define INF 1000000000
#define MOD 1000000007
#define eps 1e-8
#define ll long long
struct vec{
	int to;
	int fro;
};
vec mp[MAXN*2];
int tai[MAXN],cnt;
int n,m;
int sg[MAXN],rem[MAXN];
int N,D;
int son[MAXM][2];
int siz[MAXM],rev[MAXM];
int rt[MAXN];
int tot;
bool fucked[MAXN],okay[MAXN];
inline void be(int x,int y){
	mp[++cnt].to=y;
	mp[cnt].fro=tai[x];
	tai[x]=cnt;
}
inline void bde(int x,int y){
	be(x,y);
	be(y,x);
}
inline void torev(int x,int y,int d){
	rev[x]^=y;
	if(rev[x]&(1<<d)){
		swap(son[x][0],son[x][1]);
	}
}
inline void ud(int x){
	siz[x]=siz[son[x][0]]+siz[son[x][1]];
}
inline void pd(int x,int d){
	if(rev[x]){
		torev(son[x][0],rev[x],d-1);
		torev(son[x][1],rev[x],d-1);
		rev[x]=0;
	}
}
void ins(int &x,int y,int z,int p,int d){
	if(!x){
		x=++tot;
		son[x][0]=son[x][1]=siz[x]=rev[x]=0;
	}
	if(y==z){
		siz[x]=1;
		return ;
	}
	pd(x,d);
	int mid=y+z>>1;
	if(p<=mid){
		ins(son[x][0],y,mid,p,d-1);
	}else{
		ins(son[x][1],mid+1,z,p,d-1);
	}
}
int merge(int x,int y,int l,int r,int d){
	if(!x||!y){
		return x+y;
	}
	if(l==r){
		siz[x]|=siz[y];
		return x;
	}
	int mid=l+r>>1;
	pd(x,d);
	pd(y,d);
	son[x][0]=merge(son[x][0],son[y][0],l,mid,d-1);
	son[x][1]=merge(son[x][1],son[y][1],mid+1,r,d-1);
	ud(x);
	return x;
}
int ask(int x,int y,int z,int d){
	if(y==z){
		return y;
	}
	pd(x,d);
	int mid=y+z>>1;
	if(siz[son[x][0]]==mid-y+1){
		return ask(son[x][1],mid+1,z,d-1);
	}else{
		return ask(son[x][0],y,mid,d-1);
	}
}
void dfs(int x,int f){
	int i,y;
	rem[x]=0;
	for(i=tai[x];i;i=mp[i].fro){
		y=mp[i].to;
		if(y!=f){
			dfs(y,x);
			rem[x]^=sg[y];
		}
	}
	if(!fucked[x]){
		ins(rt[x],0,N,rem[x],D);
	}
	for(i=tai[x];i;i=mp[i].fro){
		y=mp[i].to;
		if(y!=f){
			torev(rt[y],rem[x]^sg[y],D);
			rt[x]=merge(rt[x],rt[y],0,N,D);
		}
	}
	sg[x]=ask(rt[x],0,N,D);
}
void find(int x,int f){
	int i,y;
	if(!fucked[x]&&!rem[x]){
		okay[x]=1;
	}
	for(i=tai[x];i;i=mp[i].fro){
		y=mp[i].to;
		if(y!=f){
			rem[y]^=rem[x]^sg[y];
			find(y,x);
		}
	}
}
int main(){
	int i,x,y;
	scanf("%d",&n);
	for(i=1;i<=n;i++){
		scanf("%d",&fucked[i]);
	}
	for(i=1;i<n;i++){
		scanf("%d%d",&x,&y);
		bde(x,y);
	}
	for(N=1,D=0;N<=n;N<<=1,D++){
		
	}
	N--;
	D--;
	dfs(1,0);
	if(!sg[1]){
		printf("-1\n");
		return 0;
	}
	find(1,0);
	
	for(i=1;i<=n;i++){
		if(okay[i]){
			printf("%d\n",i);
		}
	}
	return 0;
}

/*

*/