loj #2013. 「SCOI2016」幸运数字

A 国共有

在经过每一座城市时，游览者就会有机会与这座城市的幸运数字拍照，从而将这份幸运保存到自己身上。然而，幸运是不能简单叠加的，这一点游览者也十分清楚。他们迷信着幸运数字是以异或的方式保留在自己身上的。例如，游览者拍了

有些聪明的游览者发现，只要选择性地进行拍照，便能获得更大的幸运值。例如在上述三个幸运值中，只选择

第一行包含两个正整数 $G_iGi 表示 i ii 号城市的幸运值。随后 n−1 n - 1n−1 行，每行包含两个正整数 x xx、y yy，表示 x xx 号城市和 y yy 号城市之间有一条道路相连。随后 q qq 行，每行包含两个正整数 x xx、y yy，表示这名旅行者的旅行计划是从 x xx 号城市到 y yy 号城市。$

输出需要包含

样例输入

样例输出

14
11

$G_i leq 2 ^ {60}N≤20000,Q≤200000,Gi≤260$

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#define maxn 20010
using namespace std;
int n,q,num,head[maxn],fa[maxn][15],dep[maxn];
long long st[maxn][15][62],t[62],lb[62],ans;
struct node{int to,pre;}e[maxn*2];
long long qread(){
    long long i=0,j=1;
    char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')j=-1;ch=getchar();}
    while(ch<='9'&&ch>='0'){i=i*10+ch-'0';ch=getchar();}
    return i*j;
}
void Insert(int from,int to){
    e[++num].to=to;
    e[num].pre=head[from];
    head[from]=num;
}
void dfs(int x){
    for(int i=1;dep[x]-(1<<i)>0;i++)
        fa[x][i]=fa[fa[x][i-1]][i-1];
    for(int i=head[x];i;i=e[i].pre){
        int to=e[i].to;
        if(to==fa[x][0])continue;
        dep[to]=dep[x]+1;fa[to][0]=x;
        dfs(to);
    }
}
void merge(long long a[],long long c[]){
    for(int i=0;i<=60;i++)t[i]=a[i];
    for(int i=0;i<=60;i++)
        for(int j=60;j>=0;j--){
            if(!t[i])break;
            if(t[i]>>j&1){
                if(!c[j]){c[j]=t[i];break;}
                else t[i]^=c[j];
            }
        }
}
int lca(int x,int y){
    if(dep[x]<dep[y])swap(x,y);
    for(int i=14;i>=0;i--)
        if(fa[x][i]&&dep[fa[x][i]]>dep[y])
        x=fa[x][i];
    if(dep[x]>dep[y])x=fa[x][0];
    if(x==y)return x;
    for(int i=14;i>=0;i--)
        if(fa[x][i]&&fa[y][i]&&fa[x][i]!=fa[y][i])
        x=fa[x][i],y=fa[y][i];
    return fa[x][0];
}
int cal(int x,int goal){
    if(dep[x]==goal)return x;
    for(int i=14;i>=0;i--)
        if(fa[x][i]&&dep[fa[x][i]]>goal)
        x=fa[x][i];
    return fa[x][0];
}
int main(){
    scanf("%d%d",&n,&q);
    for(int i=1;i<=n;i++){
        long long x=qread();
        for(int j=60;j>=0;j--)
            if(x>>j&1){
                st[i][0][j]=x;
                break;
            }
    }
    for(int i=1;i<n;i++){
        int x=qread(),y=qread();
        Insert(x,y);Insert(y,x);
    }
    dep[1]=1;
    dfs(1);
    for(int i=1;1<<i<=n;i++)
        for(int j=1;j<=n;j++)
        if(dep[j]-(1<<i)>=0)
        merge(st[j][i-1],st[j][i]),
        merge(st[cal(j,dep[j]-(1<<i-1))][i-1],st[j][i]);
    int x,y,z,p1,p2;
    while(q--){
        memset(lb,0,sizeof(lb));
        scanf("%d%d",&x,&y);
        ans=0;z=lca(x,y);
        p1=log2(dep[x]-dep[z]+1);
        merge(st[x][p1],lb);
        merge(st[cal(x,dep[z]+(1<<p1)-1)][p1],lb);
        p2=log2(dep[y]-dep[z]+1);
        merge(st[y][p2],lb);
        merge(st[cal(y,dep[z]+(1<<p2)-1)][p2],lb);
        for(int i=60;i>=0;i--)
            if(!(ans>>i&1))ans^=lb[i];
        printf("%lld
",ans);
    }
    return 0;
}

loj #2013. 「SCOI2016」幸运数字

#2013. 「SCOI2016」幸运数字

题目描述

输入格式

输出格式

样例

样例输入

样例输出

数据范围与提示