Description

首先知道A nand B=not(A and B) (运算操作限制了数位位数为K)比如2 nand 3,K=3,则2 nand 3=not (2 and 3)=not 2=5。
给出一棵树,树上每个点都有点权,定义树上从a到b的费用为0与路径上的点的权值顺次nand的结果,例如:从2号点到5号点顺次经过2->3->5,权值分别为5、7、2,K=3,那么最终结果为0 nand 5 nand 7 nand 2=7 nand 7 nand 2=0 nand 2=7,现在这棵树需要支持以下操作。
①    Replace a b:将点a(1≤a≤N)的权值改为b。
②    Query a b:输出点a到点b的费用。
请众神给出一个程序支持这些操作。

Input

第一行N,M,K,树的节点数量、总操作个数和运算位数。
    接下来一行N个数字,依次表示节点i的权值。
接下来N-1行,每行两个数字a,b(1≤a,b≤N)表示a,b间有一条树边。
接下来M行,每行一个操作,为以上2类操作之一。
 

Output

对于操作②每个输出一行,如题目所述。

Sample Input

3 3 3
2 7 3
1 2
2 3
Query 2 3
Replace 1 3
Query 1 1

Sample Output

4
7

HINT

100%的数据N、M≤100000,K≤32

我们发现nand这个东西啊,excited。。。

既没有结合律又没有交换律。
考虑一般算法,按位拆分,分别维护0依次nand的结果,1依次nand的结构,然后就可以合并了。
因为要考虑方向所以线段树从前到后的nand和从后往前的nand都需要维护。
然后讨论讨论套上树链剖分就行了。
时间复杂度每次O(log^3N)。
#include<cstdio>
#include<cctype>
#include<queue>
#include<cstring>
#include<algorithm>
#define rep(i,s,t) for(int i=s;i<=t;i++)
#define dwn(i,s,t) for(int i=s;i>=t;i--)
#define ren for(int i=first[x];i;i=next[i])
using namespace std;
const int BufferSize=1<<16;
char buffer[BufferSize],*head,*tail;
inline char Getchar() {
if(head==tail) {
int l=fread(buffer,1,BufferSize,stdin);
tail=(head=buffer)+l;
}
return *head++;
}
typedef unsigned int uint;
inline uint read() {
uint x=0,f=1;char c=Getchar();
for(;!isdigit(c);c=Getchar()) if(c=='-') f=-1;
for(;isdigit(c);c=Getchar()) x=x*10+c-'0';
return x*f;
}
const int maxn=100010;
int n,m,k,first[maxn],next[maxn<<1],to[maxn<<1],e;
void AddEdge(int u,int v) {
to[++e]=v;next[e]=first[u];first[u]=e;
to[++e]=u;next[e]=first[v];first[v]=e;
}
int fa[maxn],son[maxn],siz[maxn],dep[maxn];
void dfs(int x) {
dep[x]=dep[fa[x]]+1;siz[x]=1;
ren if(to[i]!=fa[x]) {
fa[to[i]]=x;dfs(to[i]);
siz[x]+=siz[to[i]];if(siz[to[i]]>siz[son[x]]) son[x]=to[i];
}
}
int top[maxn],pos[maxn],cnt;
void build(int x,int tp) {
top[x]=tp;pos[x]=++cnt;
if(son[x]) build(son[x],tp);
ren if(to[i]!=fa[x]&&to[i]!=son[x]) build(to[i],to[i]);
}
int lca(int x,int y) {
int f1=top[x],f2=top[y];
while(f1!=f2) {
if(dep[f1]<dep[f2]) swap(f1,f2),swap(x,y);
x=fa[f1];f1=top[x];
}
return dep[x]<dep[y]?x:y;
}
struct Node {
uint sum0,sum1;
}T[maxn<<2],T2[maxn<<2];
uint all,val[maxn];
uint nand(uint x,uint y) {return (x&y)^all;}
Node operator + (Node A,Node B) {
Node C;C.sum0=C.sum1=0;
rep(i,0,k-1) {
if(A.sum0>>i&1) {
if(B.sum1>>i&1) C.sum0|=1u<<i;
}
else {
if(B.sum0>>i&1) C.sum0|=1u<<i;
}
if(A.sum1>>i&1) {
if(B.sum1>>i&1) C.sum1|=1u<<i;
}
else {
if(B.sum0>>i&1) C.sum1|=1u<<i;
}
}
return C;
}
void maintain(int o) {
int lc=o<<1,rc=lc|1;
T[o]=T[lc]+T[rc];
T2[o]=T2[rc]+T2[lc];
}
void update(int o,int l,int r,int p,uint v) {
if(l==r) T[o]=T2[o]=(Node){nand(0,v),nand(all,v)};
else {
int mid=l+r>>1,lc=o<<1,rc=lc|1;
if(p<=mid) update(lc,l,mid,p,v);
else update(rc,mid+1,r,p,v);
maintain(o);
}
}
Node ans;
int flag;
void query(int o,int l,int r,int ql,int qr,int tp) {
if(ql<=l&&r<=qr) {
if(!flag) ans=(!tp?T[o]:T2[o]),flag=1;
else ans=ans+(!tp?T[o]:T2[o]);
}
else {
int mid=l+r>>1,lc=o<<1,rc=lc|1;
if(!tp) {
if(ql<=mid) query(lc,l,mid,ql,qr,tp);
if(qr>mid) query(rc,mid+1,r,ql,qr,tp);
}
else {
if(qr>mid) query(rc,mid+1,r,ql,qr,tp);
if(ql<=mid) query(lc,l,mid,ql,qr,tp);
}
}
}
int Ql[maxn],Qr[maxn],Top;
void query(int x,int y) {
int z=lca(x,y),f,ql,qr;
f=top[x];flag=0;
while(f!=top[z]) {
ql=pos[f];qr=pos[x];
query(1,1,n,ql,qr,1);
x=fa[f];f=top[x];
}
ql=pos[z];qr=pos[x];
query(1,1,n,ql,qr,1);
f=top[y];
while(f!=top[z]) {
ql=pos[f];qr=pos[y];
Ql[++Top]=ql;Qr[Top]=qr;
y=fa[f];f=top[y];
}
ql=pos[z]+1;qr=pos[y];
if(ql<=qr) Ql[++Top]=ql,Qr[Top]=qr;
while(Top) {
query(1,1,n,Ql[Top],Qr[Top],0);
Top--;
}
printf("%u\n",ans.sum0);
}
int main() {
n=read();m=read();k=read();
rep(i,0,k-1) all|=1ll<<i;
rep(i,1,n) val[i]=read();
rep(i,2,n) AddEdge(read(),read());
dfs(1);build(1,1);
rep(i,1,n) update(1,1,n,pos[i],val[i]);
rep(i,1,m) {
char c=Getchar();while(!isalpha(c)) c=Getchar();
int x=read(),y=read();
if(c=='Q') query(x,y);
else update(1,1,n,pos[x],y);
}
return 0;
}

  

BZOJ2908: 又是nand的更多相关文章

  1. 【BZOJ2908】又是nand 树链剖分+线段树

    [BZOJ2908]又是nand escription 首先知道A nand B=not(A and B) (运算操作限制了数位位数为K)比如2 nand 3,K=3,则2 nand 3=not (2 ...

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

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

  3. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

  4. NAND Flash的基本操作——读、写、擦除

    基本操作 这里将会简要介绍一下NAND Flash的基本操作在NAND Flash内部是如何进行的,基本操作包括:读.写和擦除.   读:     当我们读取一个存储单元中的数据时(如图2.4),是使 ...

  5. 怎么看时序图--nand flash的读操作详解(转载)

    出处:http://blog.chinaunix.net/uid-28852942-id-3992727.html这篇文章不是介绍 nand flash的物理结构和关于nand flash的一些基本知 ...

  6. MTD下的Nand驱动

    目录 MTD下的Nand驱动 引入 平台设备资源文件 关键数据结构 平台框架 s3c24xx_nand_probe nand_scan s3c2410_nand_add_partition add_m ...

  7. 怎么看时序图--nand flash的读操作详解

    出处:http://blog.chinaunix.net/uid-28852942-id-3992727.html 这篇文章不是介绍 nand flash的物理结构和关于nand flash的一些基本 ...

  8. 转载:百为STM32开发板教程之十二——NAND FLASH

    http://bbs.21ic.com/icview-586200-1-1.html 百为STM32开发板教程之十二——NAND FLASH 参考资料:百为stm32开发板光盘V3\百为stm32开发 ...

  9. 3D NAND闪存是个啥?让国内如此疯狂

    Repost: https://news.mydrivers.com/1/477/477251.htm 上个月底武汉新芯科技主导的国家级存储器产业基地正式动工,在大基金的支持下该项目将投资240亿美元 ...

随机推荐

  1. Python 网络爬虫(新闻采集脚本)

    =====================爬虫原理===================== 通过Python访问新闻首页,获取首页所有新闻链接,并存放至URL集合中. 逐一取出集合中的URL,并访问 ...

  2. ASP.NET MVC 5 使用autofac实现DI

    使用Nuget添加Autofac.MVC的引用 启动项设置 注册Controller 注册ModelBinder 注册相关的web abstraction 为View层启用属性注入 为Action F ...

  3. dreamweaver cs6 mac破解版

    http://www.sdifenzhou.com/dreamweaver-cs6-mac.html

  4. SQL中循环和条件语句

    .if语句使用示例: declare @a int begin print @a end else begin print 'no' end .while语句使用示例: declare @i int ...

  5. SQL表值函数(上月添加1-28)

    ALTER function [dbo].[fn_getdate3] ( ) ) RETURNS @Table_Var TABLE ( LastTime datetime ) as begin Dec ...

  6. CentOS 安装Paramiko模块

    转自:http://www.cnblogs.com/hyli/p/3910585.html 1.下载安装包: https://pypi.python.org/packages/source/p/par ...

  7. c# 保存数据到txt (追加)

    StringBuilder sb = new StringBuilder(); sb.AppendLine(DateTime.Now.ToString("yyyy-MM-dd hh:mm:s ...

  8. CISCO VLAN ACL

    对于cisco VLAN ACL 首先得定义 standard ACL或 extented ACL用于抓取流量 注意这里的抓取流量不是最终的对流量的操作,而是决定什么样的流量用VLAN ACL 来处理 ...

  9. vue.js 请求数据

    VUE.JS var vm = new Vue({ el:"#list", data:{ gridData: "", }, mounted: function( ...

  10. js 函数

    函数:封装了某一块功能 四要素: 1.返回类型 2.函数名 3.参数列表4.函数体强类型语言 返回类型 函数名 首字母大写 参数列表string(字符串) Show (int a){ 函数体 }弱类型 ...