• 【[TJOI2018]异或】


    写板子了,可持久化(Trie)的板子了

    其实和主席树写法类似,还是存好左右儿子之后存好权值

    之后差分去查询就好了

    这道题第一问我们直接(dfs)序转化成区间

    第二问搞成(x,y,lca(x,y),fa[lca])之后一起差分就好了

    代码

    #include<algorithm>
    #include<iostream>
    #include<cstring>
    #include<cstdio>
    #define maxn 100005
    #define re register
    #define LL long long
    #define max(a,b) ((a)>(b)?(a):(b))
    #define min(a,b) ((a)<(b)?(a):(b))
    inline int read() {
    	char c=getchar();int x=0;while(c<'0'||c>'9') c=getchar();
    	while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();
    	return x;
    }
    struct E{int v,nxt;}e[maxn<<1];
    int n,m,__,num,bit[35],cnt;
    int rt[maxn],ch[maxn*35*2][2],v[maxn*35*2],a[maxn],Rt[maxn];
    int head[maxn],top[maxn],son[maxn],sum[maxn],fa[maxn],to[maxn],_to[maxn],deep[maxn];
    inline void add(int x,int y) {e[++num].v=y;e[num].nxt=head[x];head[x]=num;}
    int change(int pre,int t) {
    	int root=++cnt;
    	int val=bit[t];
    	v[root]=v[pre]+1;
    	if(!t) return root;
    	ch[root][val^1]=ch[pre][val^1];
    	ch[root][val]=change(ch[pre][val],t-1);
    	return root;
    }
    int query(int l,int r,int t) {
    	if(!t) return 0;
    	int val=bit[t];
    	int now=v[ch[r][val^1]]-v[ch[l][val^1]];
    	if(now) return (1<<(t-1))+query(ch[l][val^1],ch[r][val^1],t-1);
    	return query(ch[l][val],ch[r][val],t-1);
    }
    int ask(int a,int b,int c,int d,int t) {
    	if(!t) return 0;
    	int val=bit[t];
    	int now=v[ch[a][val^1]]+v[ch[b][val^1]]-v[ch[c][val^1]]-v[ch[d][val^1]];
    	if(now) return (1<<(t-1))+ask(ch[a][val^1],ch[b][val^1],ch[c][val^1],ch[d][val^1],t-1);
    	return ask(ch[a][val],ch[b][val],ch[c][val],ch[d][val],t-1);
    }
    void dfs1(int x) {
    	int maxx=-1;sum[x]=1;
    	for(re int i=head[x];i;i=e[i].nxt)
    	if(!deep[e[i].v])
    	{
    		deep[e[i].v]=deep[x]+1;
    		int now=0,t=a[e[i].v];
    		memset(bit,0,sizeof(bit));
    		while(t) bit[++now]=(t&1),t>>=1;
    		Rt[e[i].v]=change(Rt[x],32);
    		fa[e[i].v]=x;dfs1(e[i].v);
    		sum[x]+=sum[e[i].v];
    		if(sum[e[i].v]>maxx) maxx=sum[e[i].v],son[x]=e[i].v;
    	}
    }
    void dfs2(int x,int topf) {
    	top[x]=topf;to[++__]=x;_to[x]=__;
    	if(!son[x]) return;
    	dfs2(son[x],topf);
    	for(re int i=head[x];i;i=e[i].nxt) if(!top[e[i].v]) dfs2(e[i].v,e[i].v);
    }
    inline int LCA(int x,int y) {
    	while(top[x]!=top[y])
    	{if(deep[top[x]]<deep[top[y]]) std::swap(x,y); x=fa[top[x]];}
    	if(deep[x]<deep[y]) return x;return y;
    }
    int main()
    {
    	n=read(),m=read();
    	for(re int i=1;i<=n;i++) a[i]=read();
    	int x,y,opt,val;
    	for(re int i=1;i<n;i++) x=read(),y=read(),add(x,y),add(y,x);
    	deep[1]=1,dfs1(1),dfs2(1,1);
    	for(re int i=1;i<=n;i++) 
    	{
    		int now=0,t=a[to[i]];
    		memset(bit,0,sizeof(bit));
    		while(t) bit[++now]=(t&1),t>>=1;
    		rt[i]=change(rt[i-1],32);
    	}
    	while(m--)
    	{
    		opt=read();x=read();
    		if(opt==1) val=read();
    		if(opt==2) y=read(),val=read();
    		int now=0;
    		memset(bit,0,sizeof(bit));
    		while(val) bit[++now]=(val&1),val>>=1;
    		if(opt==1) printf("%d
    ",query(rt[_to[x]-1],rt[_to[x]+sum[x]-1],32));
    		if(opt==2) {int lca=LCA(x,y);printf("%d
    ",ask(Rt[x],Rt[y],Rt[lca],Rt[fa[lca]],32));}
    	}
    	return 0;
    }
    
  • 相关阅读:
    LuoguP1126 机器人搬重物(BFS)
    POJ1950----DFS
    C
    B
    A
    C
    B
    A
    A
    B
  • 原文地址:https://www.cnblogs.com/asuldb/p/10252615.html
Copyright © 2020-2023  润新知