洛谷 P4592 [TJOI2018]异或 解题报告

P4592 [TJOI2018]异或

题目描述

现在有一颗以\(1\)为根节点的由\(n\)个节点组成的树，树上每个节点上都有一个权值\(v_i\)。现在有\(Q\)次操作，操作如下：

• 1 x y：查询节点\(x\)的子树中与\(y\)异或结果的最大值

• 2 x y：查询路径\(x\)到\(y\)上点与\(z\)异或结果最大值

输入输出格式

输入格式：

第一行是两个数字\(n,Q\);

第二行是\(n\)个数字用空格隔开，第\(i\)个数字\(v_i\)表示点\(i\)上的权值

接下来\(n-1\)行，每行两个数，\(x,y\)，表示节点\(x\)与\(y\)之间有边

接下来\(Q\)行，每一行为一个查询，格式如上所述.

输出格式：

对于每一个查询,输出一行，表示满足条件的最大值。

说明

对于\(10\%\)的数据，有\(1<n,Q\leq100\)

对于\(20\%\)的数据，有\(1<n,Q\leq1000\)

对于\(40\%\)的数据，有\(1<n,Q\leq10000\)

对于\(100\%\)的数据，有\(1<n,Q\leq100000\)

对于\(100\%\)的数据，有查询\(1\)中的\(y\leq2^{30}\),查询\(2\)中的\(z\leq2^{30}\)。

---

区间异或最大值可以用可持久化字典树实现。

这个题建对DFS序建一颗，对自根向下的链建，分别处理两种询问就可以了。

---

Code:

#include <cstdio>
#include <cctype>
const int N=1e5+10;
int read()
{
    int x=0;char c=getchar();
    while(!isdigit(c)) c=getchar();
    while(isdigit(c)) {x=x*10+c-'0';c=getchar();}
    return x;
}
int max(int x,int y){return x>y?x:y;}
#define ls ch[now][0]
#define rs ch[now][1]
#define rep(i,a,b) for(int i=a;i<=b;i++)
int head[N],to[N<<1],Next[N<<1],cnt;
void add(int u,int v)
{
    to[++cnt]=v,Next[cnt]=head[u],head[u]=cnt;
}
int poi[N],n,Q;
namespace work1//子树
{
    int ch[N*32][2],mx[N*32],tot,dfn[N],low[N],ha[N],dfsclock,root[N];
    void dfs(int now,int fa)
    {
        dfn[now]=++dfsclock;
        ha[dfsclock]=now;
        for(int i=head[now];i;i=Next[i])
            if(to[i]!=fa)
                dfs(to[i],now);
        low[now]=dfsclock;
    }
    void Insert(int las,int &now,int dep,int id)
    {
        if(!now) now=++tot;
        if(dep<0){mx[now]=id;return;}
        int bit=poi[ha[id]]>>dep&1;
        Insert(ch[las][bit],ch[now][bit],dep-1,id);
        ch[now][bit^1]=ch[las][bit^1];
        mx[now]=max(mx[ls],mx[rs]);
    }
    void init()
    {
        dfs(1,0);
        rep(i,1,n) Insert(root[i-1],root[i],30,i);
    }
    int query(int now,int les,int dep,int x)
    {
        if(dep<0) return x^poi[ha[mx[now]]];
        int bit=x>>dep&1;
        if(mx[ch[now][bit^1]]>=les) return query(ch[now][bit^1],les,dep-1,x);
        return query(ch[now][bit],les,dep-1,x);
    }
    void work()
    {
        int x=read(),y=read();
        printf("%d\n",query(root[low[x]],dfn[x],30,y));
    }
}
namespace work2
{
    int f[N][20],ch[N*32][2],mx[N*32],root[N],dep[N],en[N*32],tot,tmp;
    void swap(int &x,int &y){tmp=x,x=y,y=tmp;}
    void Insert(int las,int &now,int de,int id)
    {
        if(!now) now=++tot;
        if(de<0){mx[now]=dep[id],en[now]=id;return;}
        int bit=poi[id]>>de&1;
        Insert(ch[las][bit],ch[now][bit],de-1,id);
        ch[now][bit^1]=ch[las][bit^1];
        mx[now]=max(mx[ls],mx[rs]);
    }
    void dfs(int now)
    {
        for(int i=1;f[now][i-1];i++) f[now][i]=f[f[now][i-1]][i-1];
        for(int i=head[now];i;i=Next[i])
        {
            int v=to[i];
            if(v==f[now][0]) continue;
            dep[v]=dep[now]+1;
            f[v][0]=now;
            Insert(root[now],root[v],30,v);
            dfs(v);
        }
    }
    void init()
    {
        dep[1]=1;
        Insert(0,root[1],30,1);
        dfs(1);
    }
    int LCA(int x,int y)
    {
        if(dep[x]<dep[y]) swap(x,y);
        for(int i=18;~i;i--)
            if(dep[f[x][i]]>=dep[y])
                x=f[x][i];
        if(x==y) return x;
        for(int i=18;~i;i--)
            if(f[x][i]!=f[y][i])
                x=f[x][i],y=f[y][i];
        return f[x][0];
    }
    int query(int now,int les,int de,int x)
    {
        if(de<0) return poi[en[now]]^x;
        int bit=x>>de&1;
        if(mx[ch[now][bit^1]]>=les) return query(ch[now][bit^1],les,de-1,x);
        return query(ch[now][bit],les,de-1,x);
    }
    void work()
    {
        int x=read(),y=read(),z=read(),lca=LCA(x,y);
        printf("%d\n",max(query(root[x],dep[lca],30,z),query(root[y],dep[lca],30,z)));
    }
}
int main()
{
    n=read(),Q=read();
    rep(i,1,n) scanf("%d",poi+i);
    for(int u,v,i=1;i<n;i++)
        u=read(),v=read(),add(u,v),add(v,u);
    work1::init();
    work2::init();
    rep(i,1,Q)
    {
        if(read()==1) work1::work();
        else work2::work();
    }
    return 0;
}

---

2018.11.3