bzoj4811 [Ynoi2017]由乃的OJ 树链剖分+位运算

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=4811

因为位运算的结果有可合并性，所以可以树链剖分，线段树维护；

细节很多，特别要注意从左往右运算和从右往左计算是不同的，在不同条件下一定要区分！！！

这篇博客写得很好(我就是模仿它写的)：https://blog.csdn.net/a1799342217/article/details/78818480

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef unsigned long long ull;
int const maxn=;
int n,m,p,in[maxn],id[maxn],to[maxn],tp[maxn],head[maxn],ct,dep[maxn],fa[maxn],siz[maxn],tim;
ull op[maxn][];
struct N{
    int to,next;
    N(int t=,int n=):to(t),next(n) {}
}edge[maxn<<];
struct data{ull p0,p1;}tr[maxn<<],tl[maxn<<];
struct T{int l,r;}t[maxn<<];
void add(int x,int y){edge[++ct]=N(y,head[x]);head[x]=ct;}
void dfs1(int x,int f)
{
    fa[x]=f;dep[x]=dep[f]+;siz[x]=;
    for(int i=head[x];i;i=edge[i].next)
    {
        int u=edge[i].to;
        if(u==f)continue;
        dfs1(u,x);siz[x]+=siz[u];
        if(siz[u]>siz[to[x]])to[x]=u;
    }
}
void dfs2(int x)
{
    in[id[x]=++tim]=x;
    if(to[x])tp[to[x]]=tp[x],dfs2(to[x]);
    for(int i=head[x];i;i=edge[i].next)
    {
        int u=edge[i].to;
        if(u!=fa[x]&&u!=to[x])tp[u]=u,dfs2(u);
    }
}
data update(ull op,ull w)
{
    data ret;
    if(op==)ret.p0=(&w),ret.p1=((~)&w);
    if(op==)ret.p0=(|w),ret.p1=((~)|w);
    if(op==)ret.p0=(^w),ret.p1=((~)^w);
    return ret;
}
data pushup(data x,data y)
{
    data ret;
    ret.p0=(x.p0&y.p1) | ((~x.p0)&y.p0);
    ret.p1=(x.p1&y.p1) | ((~x.p1)&y.p0);
    return ret;
}
void build(int x,int l,int r)
{
    t[x].l=l;t[x].r=r;
    if(l==r)
    {
        tl[x]=tr[x]=update(op[in[l]][],op[in[l]][]);//in[l]而非x!
        return;
    }
    int mid=((l+r)>>);
    build(x<<,l,mid); build(x<<|,mid+,r);
    tl[x]=pushup(tl[x<<],tl[x<<|]);
//    tr[x]=pushup(tr[x<<1],tr[x<<1|1]);
    tr[x]=pushup(tr[x<<|],tr[x<<]);//!!!!!!!!!!

}
data sch(int x,int l,int r,int fl)
{
    if(t[x].l>=l&&t[x].r<=r)
        return fl?tr[x]:tl[x];
//      int mid=((l+r)>>1);
    int mid=((t[x].l+t[x].r)>>);
    if(r<=mid)return sch(x<<,l,r,fl);
    else if(l>mid)return sch(x<<|,l,r,fl);
    else return pushup(sch((x<<)+fl,l,r,fl),sch((x<<|)-fl,l,r,fl));//从左往右或从右往左pushup有不同
}
data find(int x,int y)
{
    data ans1=update(,),ans2=update(,);//ans1为x到lca，ans2为y到lca
    while(tp[x]!=tp[y])
    {
//        if(dep[x]>dep[y])
        if(dep[tp[x]]>dep[tp[y]])
        {
            ans1=pushup(ans1,sch(,id[tp[x]],id[x],));
            x=fa[tp[x]];
        }
        else
        {
            ans2=pushup(sch(,id[tp[y]],id[y],),ans2);
            y=fa[tp[y]];
        }
    }
    if(dep[x]<dep[y])return pushup(pushup(ans1,sch(,id[x],id[y],)),ans2);
    else return pushup(pushup(ans1,sch(,id[y],id[x],)),ans2);
}
void query(int x,int y,ull z)//ull
{
    data ans=find(x,y);ull s=,ret=;
    for(int i=p-;i>=;i--)
    {
        if((1ull<<i)&ans.p0)ret+=(1ull<<i);
        else if(((1ull<<i)&ans.p1)&&s+(1ull<<i)<=z)
            ret+=(1ull<<i),s+=(1ull<<i);
    }
    printf("%llu\n",ret);
}
void modify(int x,int p,ull op,ull w)//ull
{
    if(t[x].l==t[x].r)
    {
        tl[x]=tr[x]=update(op,w);
        return;//
    }
    int mid=((t[x].l+t[x].r)>>);
    if(p<=mid)modify(x<<,p,op,w);
    else modify(x<<|,p,op,w);
    tl[x]=pushup(tl[x<<],tl[x<<|]);
//  tr[x]=pushup(tr[x<<1],tr[x<<1|1]);
    tr[x]=pushup(tr[x<<|],tr[x<<]);
}
int main()
{
    scanf("%d%d%d",&n,&m,&p);
    for(int i=;i<=n;i++)scanf("%llu%llu",&op[i][],&op[i][]);
    for(int i=,x,y;i<n;i++)scanf("%d%d",&x,&y),add(x,y),add(y,x);
    dfs1(,);tp[]=;dfs2();build(,,n);
    for(int i=,q,x,y;i<=m;i++)
    {
        ull z;//
        scanf("%d%d%d%llu",&q,&x,&y,&z);
        if(q==)query(x,y,z);
        else modify(,id[x],y,z);//
    }
    return ;
}