牛客竞赛&&mjt的毒瘤赛
题目链接
https://ac.nowcoder.com/acm/contest/368/F
思路
询问可以离线。
然后每个节点上建32个权值线段树(权值不大,其实只要20颗)
记录每一位权值为x(如果是根节点的话)的01和
然后从根节点向上合并。
访问到需要访问的就查询。
大体这样,不过细节和处理极其恶心(毒瘤mjt或者我码力太差)
合并的时候。
因为两个数字相同要,合并到底的时候要这样.
if(!siz[a]&&!siz[b]) return a;
if(siz[a]==siz[b]) {
for(int i=0;i<=20;++i) xo[a][i]=0;
siz[a]=0;siz[b]=0;
} else {
for(int i=0;i<=20;++i) xo[a][i]=(bool)((1<<i)&l);
siz[a]=1,siz[b]=0;
}
return a;
因为要异或上x,所以就要讨论子树中0的个数,因为出现两个的树的删去了,还要维护siz
空间\时间复杂度是\(nlog^2\)
代码
/*
询问可以离线。
然后每个节点上建32个权值线段树(权值不大,其实只要20颗)
记录每一位权值为x(如果是根节点的话)的01和
然后从根节点向上合并。
访问到需要访问的就查询。
大体这样,不过细节和处理极其恶心(毒瘤mjt或者我码力太差)
合并的时候。
因为两个数字相同要,合并到底的时候要这样
if(!siz[a]&&!siz[b]) return a;
if(siz[a]==siz[b]) {
for(int i=0;i<=20;++i) xo[a][i]=0;
siz[a]=0;siz[b]=0;
} else {
for(int i=0;i<=20;++i) xo[a][i]=(bool)((1<<i)&l);
siz[a]=1,siz[b]=0;
}
return a;
因为要异或上x,所以就要讨论子树中0的个数,因为出现两个的树的删去了,还要维护siz
空间\时间复杂度是$nlog^2$
*/
#include<bits/stdc++.h>
using namespace std;
const int N=2e5+7;
vector<int> G[N];
int n,m,w[N],cnt;
long long ans[N];
vector<pair<int,int> > Q[N];
int ls[N*20],rs[N*20],xo[N*20][21],rt[N*20],siz[N*20];
void pushup(int rt) {
for(int i=0;i<=20;++i)
xo[rt][i]=xo[ls[rt]][i]+xo[rs[rt]][i];
siz[rt]=siz[ls[rt]]+siz[rs[rt]];
}
int zibi;
int merge(int l,int r,int a,int b) {
if(!a||!b) return a+b;
if(l==r) {
if(!siz[a]&&!siz[b]) return a;
if(siz[a]==siz[b]) {
for(int i=0;i<=20;++i) xo[a][i]=0;
siz[a]=0;siz[b]=0;
} else {
for(int i=0;i<=20;++i) xo[a][i]=(bool)((1<<i)&l);
siz[a]=1,siz[b]=0;
}
return a;
}
int mid=(l+r)>>1;
ls[a]=merge(l,mid,ls[a],ls[b]);
rs[a]=merge(mid+1,r,rs[a],rs[b]);
pushup(a);
return a;
}
void insert(int &now,int l,int r,int k) {
now=++cnt;
if(l==r) {
for(int i=0;i<=20;++i) xo[now][i]=(bool)((1<<i)&k);
siz[now]=1;
return;
}
int mid=(l+r)>>1;
if(k<=mid) insert(ls[now],l,mid,k);
else insert(rs[now],mid+1,r,k);
pushup(now);
}
int query(int l,int r,int k,int dep,int rt) {
if(!rt||r<=k) return 0;
if(l>k) return xo[rt][dep];
int mid=(l+r)>>1;
// int ans=0;
if(k<=mid) return query(mid+1,r,k,dep,rs[rt])+query(l,mid,k,dep,ls[rt]);
// if(k>mid)
return query(mid+1,r,k,dep,rs[rt]);
}
int query_siz(int l,int r,int k,int rt) {
if(!rt||r<=k) return 0;
if(l>k) return siz[rt];
int mid=(l+r)>>1;
// int ans=0;
if(k<=mid) return query_siz(mid+1,r,k,rs[rt])+query_siz(l,mid,k,ls[rt]);
// if(k>mid)
return query_siz(mid+1,r,k,rs[rt]);
}
void dfs(int u,int f) {
insert(rt[u],1,n,w[u]);
for(vector<int>::iterator it=G[u].begin();it!=G[u].end();++it) {
if(*it==f) continue;
dfs(*it,u);
rt[u]=merge(1,n,rt[u],rt[*it]);
}
for(vector<pair<int,int> >::iterator it=Q[u].begin();it!=Q[u].end();++it) {
for(int i=0;i<=20;++i) {
int js=query_siz(1,n,it->first,rt[u]);
if(it->first&(1<<i)) {
ans[it->second]+=(1LL<<i)*(long long)(js-query(1,n,it->first,i,rt[u]));
} else {
ans[it->second]+=(1LL<<i)*(long long)query(1,n,it->first,i,rt[u]);
}
}
}
}
int main(){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;++i) scanf("%d",&w[i]);
for(int i=2;i<=n;++i) {
int a,b;
scanf("%d%d",&a,&b);
G[a].push_back(b);
G[b].push_back(a);
}
for(int i=1;i<=m;++i) {
int u,x;
scanf("%d%d",&u,&x);
Q[u].push_back(make_pair(x,i));
}
dfs(1,0);
for(int i=1;i<=m;++i) printf("%lld\n",ans[i]);
return 0;
}
牛客竞赛&&mjt的毒瘤赛的更多相关文章
- 2020牛客竞赛 DP F 碎碎念
作者:儒生雄才1链接:https://ac.nowcoder.com/discuss/366644来源:牛客网 题目连接:https://ac.nowcoder.com/acm/contest/300 ...
- nowcoder(牛客网)提高组模拟赛第一场 解题报告
T1 中位数(二分) 这个题是一个二分(听说是上周atcoder beginner contest的D题???) 我们可以开一个数组b存a,sort然后二分b进行check(从后往前直接遍历check ...
- 题解——牛客网OI赛制测试赛2
T1 规律题 考虑先全部选中再去重即可 #include <cstdio> #include <algorithm> #include <cstring> #inc ...
- nowcoder(牛客网)提高组模拟赛第四场 解题报告
T1 动态点分治 就是模拟..... 但是没有过!! 看了题解之后发现.... 坑点:有可能 \(x<=r\),但是
- 牛客网提高组模拟赛第七场 T3 洞穴(附bitset介绍)
就是DP. 我们可以很简单的想到要枚举中间点,进行边数的转移. 但是因为边长数据范围很大,所以我们考虑log的倍增. 状态设计为\(dp[i][j][k]\),为从节点\(i\)走\(2^k\)步能否 ...
- 牛客网提高组模拟赛第七场 T2 随机生成树
其实看懂题就很水啦qwq,就是求\(1-N\)的约数啦. 暴力求的话时间复杂度是\(O(NlogN)\)的,其实正解是枚举每个数的倍数......这样的时间复杂度是\(\frac{N}{1}+\fra ...
- 牛客网提高组模拟赛第五场 T1同余方程(异或)(位运算)
区间不好做,但是我们可以转化成前缀来做.转化为前缀之后之后就是二维前缀和. 但是我还是不怎么会做.所以只能去看吉老师的题解 (确定写的那么简单真的是题解???). 我们要求模一个数余0,就等于找它的倍 ...
- nowcoder(牛客网)普及组模拟赛第一场 解题报告
蒟蒻我可能考了一场假试 T1 绩点 这题没什么好说的,应该是只要会语言的就会做. T2 巨大的棋盘 一个模拟题吧qwq,但是要注意取模的时候先加上n或者m再取模,要不然会错的. #include< ...
- 牛客竞赛-Who killed Cock Robin
Who killed Cock Robin? I, said the Sparrow, With my bow and arrow,I killed Cock Robin. Who saw him d ...
随机推荐
- CSS常用样式属性
1.CSS字体和文本相关属性 属性 font-family 规定文本的字体系列,比如:“serif” ''sans-serif" font-size 规定文本的字体尺寸 font-style ...
- Lua 可变参数之arg与select
function TestFunc(...) local arg = { ... } --Lua .2以后不再支持默认arg参数,{}与...之间要有空格 print("输入的参数个数:&q ...
- CentOS下nginx+php的配置及nginx开机启动配置
关闭防火墙 (不然外链接是访问不了 apache) service iptables stop 关闭安全系统 SELinux( 不然报403 访问页面错误 ) 1.Nginx安装主要在于配置文件的修改 ...
- AmiGO2:在线浏览和查询GO信息的利器
GO数据库的信息是非常庞大的,为了有效的检索和浏览GO数据库的信息,官方提供了AmiGO, 可以方便的浏览,查询和下载对应信息,官网如下 http://amigo.geneontology.org/a ...
- Vector集合——单列集合的“祖宗”类
是实现可增长的对象数组:所以底层也是数组: 与collection集合不同的是,vector是同步的,意味着是单线程的,意味着效率低,速度慢, 所以在jdk1.2版本之后被ArrayList集合所取代 ...
- Linux基础命令---修改用户密码
passwd 更改用户密码,超级用户可以修改所有用户密码,普通用户只能修改自己的密码.这个任务是通过调用LinuxPAM和LibuserAPI来完成的.本质上,它使用LinuxPAM将自己初始化为一个 ...
- websocket 原理
自己写一个websocket import socket, hashlib, base64 sock = socket.socket() sock.bind(('127.0.0.1', 9000)) ...
- Linux文件管理和编辑常用命令
Linux文件管理和编辑常用命令 mkdir 命令 功能说明 mkdir 命令用于创建一个目录,mkdir是make directory的缩写 格式: mkdir [选项] 目录名 mkdir 命令的 ...
- [转载]FlipClock.js时钟,计数,3D翻转插件
1.FlipClock.js能够自动定义计数,时钟的翻牌效果,调用简单,下面简单记录下用法 2.官网地址:http://www.flipclockjs.com/ 3.调用2个文件 <link h ...
- n个整数中1出现的次数
整数中1出现的次数(从1到n整数中1出现的次数) (两种方法:1.规律.2暴力求解) 题目描述 求出1 ~ 13的整数中1出现的次数,并算出100 ~ 1300的整数中1出现的次数?为此他特别数了一下 ...