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);
}
}

【bzoj4811】[Ynoi2017]由乃的OJ 树链剖分/LCT+贪心的更多相关文章

  1. Luogu3613 睡觉困难综合征/BZOJ4811 Ynoi2017 由乃的OJ 树链剖分、贪心

    传送门 题意:给出一个$N$个点的树,树上每个点有一个位运算符号和一个数值.需要支持以下操作:修改一个点的位运算符号和数值,或者给出两个点$x,y$并给出一个上界$a$,可以选取一个$[0,a]$内的 ...

  2. BZOJ4811 [Ynoi2017]由乃的OJ 树链剖分

    原文链接http://www.cnblogs.com/zhouzhendong/p/8085286.html 题目传送门 - BZOJ4811 题意概括 是BZOJ3668长在树上并加上修改和区间询问 ...

  3. [BZOJ4811][YNOI2017]由乃的OJ(树链剖分+线段树)

    起床困难综合症那题,只要从高往低贪心,每次暴力跑一边看这一位输入0和1分别得到什么结果即可. 放到序列上且带修改,只要对每位维护一个线段树,每个节点分别记录0和1从左往右和从右往左走完这段区间后变成的 ...

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

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4811 因为位运算的结果有可合并性,所以可以树链剖分,线段树维护: 细节很多,特别要注意从左往 ...

  5. bzoj4811 [Ynoi2017]由乃的OJ 树链剖分+贪心+二进制

    题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=4811 题解 我现在为什么都写一题,调一天啊,马上真的退役不花一分钱了. 考虑这道题的弱化版 N ...

  6. 【bzoj4811】[Ynoi2017]由乃的OJ 树链剖分+线段树区间合并

    题解: 好像和noi那题并没有什么区别 只是加上了修改和变成树上 比较显然我们可以用树链剖分来维护

  7. 【BZOJ4811】[Ynoi2017]由乃的OJ 树链剖分+线段树

    [BZOJ4811][Ynoi2017]由乃的OJ Description 由乃正在做她的OJ.现在她在处理OJ上的用户排名问题.OJ上注册了n个用户,编号为1-",一开始他们按照编号排名. ...

  8. Bzoj 2243: [SDOI2011]染色 树链剖分,LCT,动态树

    2243: [SDOI2011]染色 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 5020  Solved: 1872[Submit][Status ...

  9. Bzoj 1036: [ZJOI2008]树的统计Count 树链剖分,LCT

    1036: [ZJOI2008]树的统计Count Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 11102  Solved: 4490[Submit ...

随机推荐

  1. maven 安装 过程

    maven 安装 过程 1 下载: apache-maven-3.0.3-bin.zip 压缩包 2 将压缩包解压到需要安装的目录文件中. 比如解压到: Z:\zr_anzhungwenjian\ap ...

  2. Aggregate可以做字符串累加、数值累加、最大最小值

    Aggregate("", (m, n) => m + n + ".").TrimEnd('.'); .Aggregate(0, (m, n) => ...

  3. logistic回归和线性回归

    1.输出: 线性回归输出是连续的.具体的值(如具体房价123万元) 回归 逻辑回归的输出是0~1之间的概率,但可以把它理解成回答“是”或者“否”(即离散的二分类)的问题 分类 2.假设函数 线性回归: ...

  4. spring mvc default-servlet mvc:resources mvc:default-servlet-handler区别

    mvc:default-servlet-handler其实就是default-servlet 交由web容器自己处理 mvc:resources spring来处理 没有被映射的url web容器来处 ...

  5. Eclipse中建立Maven项目后,Java Resources资源文件下没有src/main/java文件夹

    当建立好一个Maven项目后,在Java Resources资源文件夹下没有看到src/main/java文件夹,然后手动去创建Source Folder时,提示该文件已存在,如图: 有一个解决办法: ...

  6. MySQL5.7 在CentOS 下的安装

    尝试了在版本的CentOS6.6 和CentOS7.2 下安装,在6.6下比较复杂些.特地做下记录 在CentOS7.2 下安装,需要在官网下载 mysql-5.7.16-1.el7.x86_64.r ...

  7. xml和configparser模块

    一.xml模块 xml是实现不同语言或程序之间进行数据交换的协议,跟json差不多,但json使用起来更简单, 但至今很多传统公司如金融行业的很多系统的接口还主要是xml. xml的格式如下,就是通过 ...

  8. Android适配器Adapter的学习

    Android中有很多的适配器,首先看看这些适配器的继承结构 这些适配器中,BaseAdapter用的最多,也用的最熟,先放过他,从ArrayAdapter开始 一个listAdapter用来管理一个 ...

  9. Tornado抽象方法抽象类

    #!/usr/bin/env python #抽象方法抽象类 import abc class Foo(metaclass=abc.ABCMeta): def f1(self): raise Exce ...

  10. 关于Android使用Instrumentation做功能测试的时候遇到的一个问题

    最近在看测试方面的东西,看到官网上有一个使用Instrumentation做功能测试的例子,自己敲了敲,但是在自己的手机上就是测不过. 经过调试,我发现是我手机上的输入法把输入事件拦截了,需要多输入一 ...