A1486. 树（王康宁）

题目：http://www.tsinsen.com/A1486

题解：

其实看到和路径有关的就应该想到点分治。

我们找出重心之后遍历每一棵子树得到它的 { x=经过特殊点的个数，y=到rt的异或和}

然后我们按x排序，维护两个头尾指针不断把满足条件的加入trie，然后把左边的放进trie里查询。

但是还有一个问题，所取的两个点不能位于同一棵子树！！！

我yy了一个做法。我们在用三元组来记录{ x=经过特殊点的个数，y=到rt的异或和，ch=所属子树}

然后往trie里插的时候，每条边保留两个ch表示有哪个子树的点从trie往下经过了这里。必须保证这两个ch不同。

然后查询的时候就判断就行了。注意任何时刻往下走的时候都要判断可行性，否则直接返回-1.

代码：

 #include<cstdio>
 #include<cstdlib>
 #include<cmath>
 #include<cstring>
 #include<algorithm>
 #include<iostream>
 #include<vector>
 #include<map>
 #include<set>
 #include<queue>
 #include<string>
 #define inf 1000000000
 #define maxn 250000+5
 #define maxm 8000000+5
 #define eps 1e-10
 #define ll long long
 #define ull unsigned long long
 #define pa pair<int,int>
 #define for0(i,n) for(int i=0;i<=(n);i++)
 #define for1(i,n) for(int i=1;i<=(n);i++)
 #define for2(i,x,y) for(int i=(x);i<=(y);i++)
 #define for3(i,x,y) for(int i=(x);i>=(y);i--)
 #define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go)
 #define for5(n,m) for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)
 #define mod 1000000007
 #define lch k<<1,l,mid
 #define rch k<<1|1,mid+1,r
 #define sqr(x) (x)*(x)
 using namespace std;
 inline int read()
 {
     int x=,f=;char ch=getchar();
     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
     while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}
     return x*f;
 }
 int n,m,k,cnt,tot,head[maxn],v[maxn],w[maxn],t[maxm][][],s[maxn],f[maxn],sum,rt,ans=-;
 bool del[maxn];
 struct edge{int go,next;}e[*maxn];
 struct rec{int x,y,ch;}a[maxn];
 inline bool cmp(rec a,rec b){return a.x<b.x;}
 inline void add(int x,int y)
 {
     e[++tot]=(edge){y,head[x]};head[x]=tot;
     e[++tot]=(edge){x,head[y]};head[y]=tot;
 }
 inline void insert(int y,int ch)
 {
     int x=;
     for3(i,,)
     {
         int j=y>>i&;
         if(!t[x][j][])t[x][j][]=++tot,t[x][j][]=ch;
         else if(t[x][j][]!=ch)t[x][j][]=ch;
         x=t[x][j][];
     }
 }
 inline int query(int y,int ch)
 {
     int x=,ret=;
     for3(i,,)
     {
         int j=y>>i&^;
         if(t[x][j][]&&((t[x][j][]&&t[x][j][]!=ch)||(t[x][j][]&&t[x][j][]!=ch)))ret^=<<i,x=t[x][j][];
         else
         {
             j^=;
             if(t[x][j][]&&((t[x][j][]&&t[x][j][]!=ch)||(t[x][j][]&&t[x][j][]!=ch)))x=t[x][j][];
             else return -;
         }
     }
     return ret;
 }
 inline void getdep(int x,int fa,int w1,int w2,int w3)
 {
     a[++cnt]=(rec){w1,w2,w3};
     for4(i,x)if(!del[y]&&y!=fa)getdep(y,x,w1+v[y],w2^w[y],w3);
 }
 inline void getrt(int x,int fa)
 {
     s[x]=;f[x]=;
     for4(i,x)if(!del[y]&&y!=fa)
     {
         getrt(y,x);
         s[x]+=s[y];
         f[x]=max(f[x],s[y]);
     }
     f[x]=max(f[x],sum-s[x]);
     if(f[x]<f[rt])rt=x;
 }
 inline void work(int x)
 {
     del[x]=;cnt=;
     for4(i,x)if(!del[y])getdep(y,x,v[x]+v[y],w[x]^w[y],y);
     sort(a+,a+cnt+,cmp);
     int tmp=k+v[x],j=cnt;
     for1(i,cnt)
     {
         while(j>i&&a[i].x+a[j].x>=tmp)insert(a[j].y,a[j].ch),j--;
         ans=max(ans,query(a[i].y^w[x],a[i].ch));
     }
     for1(i,cnt)if(a[i].x>=k)ans=max(ans,a[i].y);
     for0(i,tot)t[i][][]=t[i][][]=t[i][][]=t[i][][]=t[i][][]=t[i][][]=;
     tot=;
     for4(i,x)if(!del[y])
     {
         sum=s[y];rt=;
         getrt(y,x);
         work(rt);
     }
 }
 int main()
 {
    freopen("input.txt","r",stdin);
    freopen("output.txt","w",stdout);
    n=read();k=read();
    for1(i,n)v[i]=read();
    for1(i,n)w[i]=read();
    for1(i,n-)add(read(),read());
    tot=;
    sum=n;
    f[rt=]=inf;
    getrt(,);
    work(rt);
    for1(i,n)if(v[i]>=k)ans=max(ans,w[i]);
    cout<<ans<<endl;
    return ;
 }