bzoj 4811 由乃的OJ
bzoj 4811 由乃的OJ
- 考虑树链剖分.
- 树剖后用一颗线段树维护一段连续区间,类似于一个函数,各位上进入 \(0/1\) ,输出的数字分别是什么.注意到最多只有 \(64\) 位,可以用一个 \(unsigned\ long\ long\) 的大数状压表示,合并两段区间时推导一下可以做到 \(O(1)\) .
- 注意 \(3 种\)位运算混在一起,满足交换律,却不满足结合律,所以从区间左/右侧进入同一个数,最后得到的数不同.需要维护两个方向的函数.
- 修改时在线段树上单点修改就好,查询用树剖的基本操作,先找出 \(x\to y\) 这一段的函数,再贪心构造解.
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define mp make_pair
#define pii pair<int,int>
typedef unsigned long long ull;
ull read(){
ull nm=0,fh=1;char cw=getchar();
for(;!isdigit(cw);cw=getchar()) ;
for(;isdigit(cw);cw=getchar()) nm=nm*10+(cw-'0');
return nm*fh;
}
const int MAXN=2e5+10;
ull cnt=0,head[MAXN],to[MAXN<<1],nx[MAXN<<1];
inline void addedge(ull u,ull v)
{
++cnt;
to[cnt]=v;
nx[cnt]=head[u];
head[u]=cnt;
}
inline void ins(ull u,ull v)
{
addedge(u,v);
addedge(v,u);
}
const unsigned long long U1=1;
ull fa[MAXN],siz[MAXN],mxson[MAXN],top[MAXN],dfn[MAXN],rnk[MAXN],dep[MAXN],idx=0;
void dfs1(ull u)
{
siz[u]=U1;
dep[u]=dep[fa[u]]+U1;
for(ull i=head[u];i>0;i=nx[i])
{
ull v=to[i];
if(v==fa[u])
continue;
fa[v]=u;
dfs1(v);
siz[u]+=siz[v];
if(siz[v]>siz[mxson[u]])
mxson[u]=v;
}
}
void dfs2(ull u,ull Tp)
{
dfn[u]=++idx;
rnk[idx]=u;
top[u]=Tp;
if(mxson[u]>0)
dfs2(mxson[u],Tp);
for(ull i=head[u];i>0;i=nx[i])
{
ull v=to[i];
if(v==fa[u] || v==mxson[u])
continue;
dfs2(v,v);
}
}
inline ull calc(ull x,ull bas,ull op)
{
if(op==1)
return x&bas;
if(op==2)
return x|bas;
return x^bas;
}
ull val[MAXN],opt[MAXN];
unsigned long long maxv=0;
struct node{
ull l,r;
ull d[2][2];//d[左入/右入][全0/全1]
void init(ull v,ull op)
{
d[0][0]=d[1][0]=calc(0,v,op);
d[0][1]=d[1][1]=calc(maxv,v,op);
}
void merge(node L,node R)
{
d[0][0]=(L.d[0][0]&R.d[0][1])|((~L.d[0][0])&R.d[0][0]);
d[1][0]=(R.d[1][0]&L.d[1][1])|((~R.d[1][0])&L.d[1][0]);
d[0][1]=(L.d[0][1]&R.d[0][1])|((~L.d[0][1])&R.d[0][0]);
d[1][1]=(R.d[1][1]&L.d[1][1])|((~R.d[1][1])&L.d[1][0]);
}
void inverse()
{
swap(d[0][0],d[1][0]);
swap(d[0][1],d[1][1]);
}
};
struct SegmentTree{
#define root Tree[o]
#define lson Tree[o<<1]
#define rson Tree[o<<1|1]
node Tree[MAXN<<2];
void BuildTree(ull o,ull l,ull r)
{
root.l=l,root.r=r;
if(l==r)
{
root.init(val[rnk[l]],opt[rnk[l]]);
return;
}
ull mid=(l+r)>>1;
BuildTree(o<<1,l,mid);
BuildTree(o<<1|1,mid+1,r);
root.merge(lson,rson);
}
void update(ull o,ull pos,ull newv,ull newopt)
{
ull l=root.l,r=root.r;
if(l==r)
{
root.init(newv,newopt);
return;
}
ull mid=(l+r)>>1;
if(pos<=mid)
update(o<<1,pos,newv,newopt);
else
update(o<<1|1,pos,newv,newopt);
root.merge(lson,rson);
}
node query(ull o,ull L,ull R)
{
ull l=root.l,r=root.r;
if(L<=l && r<=R)
return root;
node res;
res.init(0,3);
if(l>R || r<L)
return res;
res.merge(query(o<<1,L,R),query(o<<1|1,L,R));
return res;
}
}T;
ull n,m,k;
void init()
{
ull rt=(n+1)>>1;
fa[rt]=0;
for(ull i=0;i<k;i++)
maxv+=(U1<<i);
dfs1(rt);
dfs2(rt,rt);
T.BuildTree(1,1,n);
}
ull solve(ull x,ull y,ull lim)
{
ull tot=0;
node res1,res2,cur;
res1.init(0,3);
res2.init(0,3);
while(top[x]!=top[y])
{
if(dep[top[x]]>dep[top[y]])
{
res1.merge(T.query(1,dfn[top[x]],dfn[x]),res1);
x=fa[top[x]];
}
else
{
res2.merge(T.query(1,dfn[top[y]],dfn[y]),res2);
y=fa[top[y]];
}
}
if(dep[x]>dep[y])
{
node c=T.query(1,dfn[y],dfn[x]);
res1.merge(c,res1);
}
else
{
node c=T.query(1,dfn[x],dfn[y]);
res2.merge(c,res2);
}
res1.inverse();
cur.merge(res1,res2);
for(ull i=k;i-->0;)
{
if(cur.d[0][0] & (U1<<i))
tot+=(U1<<i);
else if( ( (cur.d[0][1]) & (U1<<i) ) >0 && lim>=(U1<<i) )
lim-=(U1<<i),tot+=(U1<<i);
}
return tot;
}
int main()
{
n=read(),m=read(),k=read();
for(ull i=1;i<=n;++i)
opt[i]=read(),val[i]=read();
for(ull i=1;i<n;++i)
{
ull u=read(),v=read();
ins(u,v);
}
init();
while(m--)
{
ull op=read();
ull x=read(),y=read(),z=read();
if(op==1)
cout<<solve(x,y,z)<<endl;
else
T.update(1,dfn[x],z,y);
}
return 0;
}
bzoj 4811 由乃的OJ的更多相关文章
- bzoj 4811: [Ynoi2017]由乃的OJ
树链剖分,用zkw线段树维护每条链两个方向上对每一位的变换情况,由于位数较少,可以用两个unsigned long long表示 #include<cstdio> typedef unsi ...
- BZOJ 4811 [Ynoi2017]由乃的OJ ——Link-Cut Tree
直接维护按照顺序经过每一段,初始的1可以变成什么,初始为0可以变成什么. 然后答案就可以和起床困难综合征一样贪心处理了. 写起来并不好写. 发现交换左右子树之后答案会改变,GG 调了一天,最后还是T掉 ...
- BZOJ 4811 树链剖分+线段树
思路: 感觉这题也可神了.. (还是我太弱) 首先发现每一位不会互相影响,可以把每一位分开考虑,然后用树链剖分或者LCT维护这个树 修改直接修改,询问的时候算出来每一位填0,1经过这条链的变换之后得到 ...
- bzoj 1565 [NOI2009]植物大战僵尸 解题报告
1565: [NOI2009]植物大战僵尸 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 2161 Solved: 1000[Submit][Stat ...
- 博主自传——蒟蒻的OI之路
博主来自河北石家庄市第二中学,现在读高二,主攻信息学竞赛(其实并没有学习其他学科竞赛). NOIP中人品大爆发,使劲挤进河北省一等奖队伍,侥幸留在竞赛团队中(差点就淘汰出局啦). 关于我的ID,YOU ...
- 【BZOJ】3224: Tyvj 1728 普通平衡树(某不科学的oj)
http://www.lydsy.com/JudgeOnline/problem.php?id=3224 无力吐槽,无力吐槽,无力吐槽....... bzoj竟然不能用time(0)我竟然不造!!re ...
- BZOJ 3595: [Scoi2014]方伯伯的Oj SBT+可持久化Treap
3595: [Scoi2014]方伯伯的Oj Time Limit: 6 Sec Memory Limit: 256 MBSubmit: 102 Solved: 54[Submit][Status ...
- bzoj 3678 wangxz与OJ
3678: wangxz与OJ Time Limit: 10 Sec Memory Limit: 128 MBhttp://www.lydsy.com/JudgeOnline/problem.php ...
- BZOJ 3595: [Scoi2014]方伯伯的Oj Splay + 动态裂点 + 卡常
Description 方伯伯正在做他的Oj.现在他在处理Oj上的用户排名问题. Oj上注册了n个用户,编号为1-”,一开始他们按照编号排名.方伯伯会按照心情对这些用户做以下四种操作,修改用户的排名和 ...
随机推荐
- [小问题笔记(七)] JS和 jQuery常用语句笔记(隐藏/显示/禁用标签 日期操作 跳转等)
隐藏/显示标签 $("#div1").css("display", "none");$("#div2").css(&qu ...
- WinCE数据通讯之Web Service篇
准备写个WinCE平台与数据库服务器数据通讯交互方面的专题文章,今天先整理个Web Service通讯方式. 公司目前的硬件产品平台是WinCE5.0,数据通讯是连接服务器与终端的桥梁,关系着终端的数 ...
- 修饰器Decorator
类的修饰 许多面向对象的语言都有修饰器(Decorator)函数,用来修改类的行为.目前,有一个提案将这项功能,引入了 ECMAScript. @testable class MyTestableCl ...
- vue削笔机
// 1.vue 是异步更新dom // 2.vue修改多次状态,vue 只会渲染一次 // 3.vue变化侦查机制,每次状态的变化都会发出一个渲染信号. // 检查队列中是否存在,不存在将渲染操作添 ...
- CSP(Content Security Policy) 入门教程
参考: http://www.ruanyifeng.com/blog/2016/09/csp.html https://developer.mozilla.org/en-US/docs/Web/HTT ...
- VMWare虚拟机网络配置
Bridged(桥接模式) 桥接模式相当于虚拟机和主机在同一个真实网段,VMWare充当一个集线器功能(一根网线连到主机相连的路由器上),所以如果电脑换了内网,静态分配的ip要更改.图如下: NAT( ...
- Codeforces Round #349 (Div. 2)
第一题直接算就行了为了追求手速忘了输出yes导致wa了一发... 第二题技巧题,直接sort,然后把最大的和其他的相减就是构成一条直线,为了满足条件就+1 #include<map> #i ...
- 第八天 1-7 实战:创建一个root无法删除的文件
实战:创建一个root无法删除的文件 简介:Linux文件的最底层(内核级别)属性的查看与修改 命令:lsattr.chattr Linux文件除了具有基本权限rwx,及特殊权限(SUID.SGID. ...
- [嵌入式培训 笔记]-----Vim编辑器使用简介
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 第一讲小结 1. 光标在屏幕文本中的移动既 ...
- 河南省多校联盟二-C
1281: 邪能炸弹 时间限制: 1 秒 内存限制: 128 MB提交: 222 解决: 80 题目描述 正在入侵艾泽拉斯的古尔丹偶然间得到了一颗邪能炸弹,经过研究,他发现这是一颗威力极其巨大且难 ...