Description
给你一个有n个点的树,每个点的包括一个位运算opt和一个权值x,位运算有&,l,^三种,分别用1,2,3表示。
每次询问包含三个数x,y,z,初始选定一个数v。然后v依次经过从x到y的所有节点,每经过一个点i,v就变成v opti
xi,所以他想问你,最后到y时,希望得到的值尽可能大,求最大值?给定的初始值v必须是在[0,z]之间。每次修
改包含三个数x,y,z,意思是把x点的操作修改为y,数值改为z。
Input
第一行三个数n,m,k。k的意义是每个点上的数,以及询问中的数值z都 <2^k。之后n行
每行两个数x,y表示该点的位运算编号以及数值
之后n - 1行,每行两个数x,y表示x和y之间有边相连
之后m行,每行四个数,Q,x,y,z表示这次操作为Q(1位询问,2为修改),x,y,z意义如题所述
0 <= n , m <= 100000 , k <= 64
Output
对于每个操作1,输出到最后可以造成的最大刺激度v
Sample Input
5 5 3
1 7
2 6
3 7
3 6
3 1
1 2
2 3
3 4
1 5
1 1 4 7
1 1 3 5
2 1 1 3
2 3 3 3
1 1 3 2
Sample Output
7
1
5
Sol
首先这题的核心思想就是把0和全集都过一遍,然后贪心选取,那么这道题就是在树上动态维护这样的操作。
树链剖分和LCT均可。
对于一个节点,我们维护这个点用全集经过以及用0经过的值。
显然可以开64个进行维护,但是这样会TLE。
设v0是0经过的结果,v1是全集经过的结果。
考虑所有位合并到一起进行,如果合并的话,那么公式就是:
(~a.v0&b.v0)|(a.v0&b.v1),(~a.v1&b.v0)|(a.v1&b.v1)
意义就是:要成为0的话,那么对于前面的a,a0中0的部分要经过b的0的部分并结合,已经是1的部分要经过b的1的部分并结合(这样才能保证合法),之后把这两种走法合并起来。
要成为1的话就要用a1去更新,同理。
因为位运算不满足交换律,所以还要维护两种方向的结果。
LCT的话在rev时还要记得交换l->r,r->l的值。
由于本人自带大常数,LCT自带大常数,bzoj还是老爷机,所以在bzojTLE了,luogu可以ac。
splay的时候记得注意合并顺序,坚决不能弄反。
初值就是\(0\quad opt\quad x\)以及\(全集全集\quad opt\quad x\)的结果啦,这样后面就不用考虑opt对合并的影响啦。
Code
#include <bits/stdc++.h>
#define ull unsigned long long
using namespace std;
int n,m,op,x,y,rev[100005],fa[100005],c[100005][2];ull gg=0,z;
struct val{ull v0,v1;val(ull v0=0,ull v1=0):v0(v0),v1(v1){}}lr[100005],rl[100005],v[100005];
bool rt(int x){return c[fa[x]][0]!=x&&c[fa[x]][1]!=x;}
bool R(int x,int y){return c[x][1]==y;}
val merge(val a,val b){val c=val((~a.v0&b.v0)|(a.v0&b.v1),(~a.v1&b.v0)|(a.v1&b.v1));return c;}
void up(int x)
{
lr[x]=rl[x]=v[x];
if(c[x][0]) lr[x]=merge(lr[c[x][0]],lr[x]),rl[x]=merge(rl[x],rl[c[x][0]]);
if(c[x][1]) lr[x]=merge(lr[x],lr[c[x][1]]),rl[x]=merge(rl[c[x][1]],rl[x]);
}
void rever(int x){rev[x]^=1;swap(c[x][0],c[x][1]);swap(lr[x],rl[x]);}
void dn(int x){if(!rev[x]) return;rev[x]=0;rever(c[x][0]);rever(c[x][1]);}
void ud(int x){if(!rt(x)) ud(fa[x]);dn(x);}
void ro(int x)
{
int y=fa[x],z=fa[y],a=R(y,x),b=!a;
if(!rt(y)) c[z][R(z,y)]=x;
fa[x]=z;fa[y]=x;fa[c[x][b]]=y;c[y][a]=c[x][b];c[x][b]=y;up(y);up(x);
}
void spl(int x){ud(x);for(int y=fa[x],z=fa[y];!rt(x);ro(x),y=fa[x],z=fa[y])if(!rt(y))ro(R(z,y)^R(y,x)?y:x);}
void acs(int x){for(int y=0;x;y=x,x=fa[x]) spl(x),c[x][1]=y,up(x);}
void mrt(int x){acs(x);spl(x);rever(x);}
void lnk(int x,int y){mrt(x);fa[x]=y;}
ull query(int x,int y,ull z)
{
mrt(x);acs(y);spl(y);
ull ans=0,a0=lr[y].v0,a1=lr[y].v1,i;int k;
for(k=63,i=(1ull<<k);k>=0;i>>=1,k--){if(a0&i) ans+=i;else if((a1&i)&&z>=i) ans+=i,z-=i;}
return ans;
}
int main()
{
scanf("%d%d%*d",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%d%llu",&op,&z);
if(op==1) v[i]=val(gg,z);if(op==2) v[i]=val(z,~gg);if(op==3) v[i]=val(z,~z);
}
for(int i=1;i<n;i++) scanf("%d%d",&x,&y),lnk(x,y);
for(int i=1;i<=m;i++)
{
scanf("%d%d%d%llu",&op,&x,&y,&z);
if(op==1){printf("%llu\n",query(x,y,z));continue;}else mrt(x);
if(y==1) v[x]=val(gg,z);if(y==2) v[x]=val(z,~gg);if(y==3) v[x]=val(z,~z);
up(x);
}
}