牛客竞赛&&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 ...
随机推荐
- sql server 中DateName()函数及DatePart()函数
Datepart():返回代表指定日期的指定日期部分的整数 语法:Datepart(datepart,date) 返回类型:int DateName():返回代表指定日期的指定日期部分的字符串 语法 ...
- word论文文献引用上标括号
参考 http://jingyan.baidu.com/article/c45ad29c310734051753e20d.html 在插入参考文献引用的尾注时,默认为上标数据且没有中括号.现在要统一加 ...
- field_automation源码分析
field_automation主要通过uvm_field_utils_begin和uvm_field_utils_end来加到自己写的class中. uvm_field_utils_begin,在o ...
- html5-css列表和表格
td{ /*width: 150px; height: 60px;*/ padding: 10px; text-align: center;} table{ width ...
- Another Meaning (KMP + DP)
先用KMP重叠匹配求出各个匹配成功的尾串位置.然后利用DP去求,那转移方程应该是等于 上一个状态 (无法匹配新尾巴) 上一个状态 + 以本次匹配起点为结尾的状态(就是说有了新的位置) + 1 (单单一 ...
- Swift之关键字使用(I)
static和class的使用 static 使用 在非class的类型(包括enum和struct)中,一般使用static来描述类型作用域.在这个类型中,我们可以在类型范围中声明并使用存储属性,计 ...
- HCatalog 学习之路
最近在使用sqoop把数据从hive数仓导出到mysql数据库中接触到了hcatalog,所以特意学习了解一下相关知识,据悉hcatalog还是apache顶级项目. 学习参考: HCatalog 介 ...
- android使用ARouter跳转activity(阿里巴巴开源的)
android使用ARouter跳转activity(阿里巴巴开源的) 使用ARouter方式,点击按钮跳转到其他activitypublic void buyOrSell(String str){ ...
- 2018年Android面试题含答案--适合中高级
1.java中==和equals和hashCode的区别 基本数据类型的==比较的值相等. 类的==比较的内存的地址,即是否是同一个对象,在不覆盖equals的情况下,同比较内存地址,原实现也为 = ...
- vue 去掉路由中的#
在router.js中修改, const router = new VueRouter({ mode: 'history', routes: [...] })