CF1055F Tree and XOR
就是选择两个数找第k大对儿
第k大?二分+trie上验证
O(nlognlogn)
直接按位贪心
维护可能的决策点(a,b)表示可能答案的对儿在a和b的子树中选择
所有可能决策点都贡献这一位是0,看是否<=k
然后更新出下一层的决策点
但是空间太小,
所以要滚动
我的方法:
维护trie节点和控制区间,维护每个区间的元素,维护决策点
注意,(a,b)(b,a)算两个。考虑a!=b时候贡献*=2
#include<bits/stdc++.h>
#define reg register int
#define il inline
#define fi first
#define se second
#define mk(a,b) make_pair(a,b)
#define numb (ch^'0')
#define pb push_back
#define solid const auto &
#define enter cout<<endl
#define pii pair<int,int>
// #define int long long
using namespace std;
typedef long long ll;
template<class T>il void rd(T &x){
char ch;x=;bool fl=false;
while(!isdigit(ch=getchar()))(ch=='-')&&(fl=true);
for(x=numb;isdigit(ch=getchar());x=x*+numb);
(fl==true)&&(x=-x);
}
template<class T>il void output(T x){if(x/)output(x/);putchar(x%+'');}
template<class T>il void ot(T x){if(x<) putchar('-'),x=-x;output(x);putchar(' ');}
template<class T>il void prt(T a[],int st,int nd){for(reg i=st;i<=nd;++i) ot(a[i]);putchar('\n');} namespace Miracle{
const int N=1e6+;
ll k,a[N];
int b[][N],num[];
ll ans;
struct tr{
int ls,rs;
int l,r;
void init(){
ls=rs=l=r=;
}
void up(int ll,int rr){
l=ll;r=rr;
}
int sz(){
if(l&&r) return r-l+;
return ;
}
}t[][N];
int tot[];
int n;
pii s[][N];
int q[];
struct node{
int nxt,to;
ll val;
}e[N];
int hd[N],cnt;
void add(int x,int y,ll z){
e[++cnt].nxt=hd[x];
e[cnt].to=y;e[cnt].val=z;
hd[x]=cnt;
}
void dfs(int x,ll dis){
b[][++num[]]=x;a[x]=dis;
for(reg i=hd[x];i;i=e[i].nxt){
int y=e[i].to;
dfs(y,dis^e[i].val);
}
}
int tmp;
void build(int x,int d){
int las=tmp,now=tmp^;
int bc=t[las][x].r;
for(reg i=t[las][x].l;i<=t[las][x].r;++i){
if(!(a[b[las][i]]&(1LL<<d))){
b[now][++num[now]]=b[las][i];
}else{
b[now][bc--]=b[las][i];
}
}
if(num[now]>=t[las][x].l){
++tot[now];
t[las][x].ls=tot[now];
t[now][tot[now]].init();
t[now][tot[now]].up(t[las][x].l,num[now]);
}
if(num[now]<t[las][x].r){
++tot[now];
t[las][x].rs=tot[now];
t[now][tot[now]].init();
t[now][tot[now]].up(num[now]+,t[las][x].r);
}
num[now]=t[las][x].r;
}
int main(){
rd(n);rd(k);
int y;ll w;
for(reg i=;i<=n;++i){
rd(y);rd(w);add(y,i,w);
}
dfs(,);
t[][++tot[]].up(,num[]);
num[]=;tot[]=;
build(tot[],);
s[][++q[]]=mk(,);
for(reg d=;d>=;--d){
ll con=;
int las=tmp;
int now=tmp^;
for(reg i=;i<=q[las];++i){
if(s[las][i].fi==s[las][i].se){
con+=(ll)t[now][t[las][s[las][i].fi].ls].sz()*t[now][t[las][s[las][i].se].ls].sz()+
(ll)t[now][t[las][s[las][i].fi].rs].sz()*t[now][t[las][s[las][i].se].rs].sz();
}else{
con+=(ll)t[now][t[las][s[las][i].fi].ls].sz()*t[now][t[las][s[las][i].se].ls].sz()*+
(ll)t[now][t[las][s[las][i].fi].rs].sz()*t[now][t[las][s[las][i].se].rs].sz()*;
}
}
// cout<<" con "<<con<<endl;
if(k>con){//
k-=con;
ans+=(1LL<<d);
q[now]=;
for(reg i=;i<=q[las];++i){
if(t[las][s[las][i].fi].ls&&t[las][s[las][i].se].rs){
s[now][++q[now]]=mk(t[las][s[las][i].fi].ls,t[las][s[las][i].se].rs);
}
if(s[las][i].fi!=s[las][i].se){
if(t[las][s[las][i].fi].rs&&t[las][s[las][i].se].ls){
s[now][++q[now]]=mk(t[las][s[las][i].fi].rs,t[las][s[las][i].se].ls);
}
}
}
}else{
q[now]=;
for(reg i=;i<=q[las];++i){
if(t[las][s[las][i].fi].ls&&t[las][s[las][i].se].ls){
s[now][++q[now]]=mk(t[las][s[las][i].fi].ls,t[las][s[las][i].se].ls);
}
if(t[las][s[las][i].fi].rs&&t[las][s[las][i].se].rs){
s[now][++q[now]]=mk(t[las][s[las][i].fi].rs,t[las][s[las][i].se].rs);
}
}
}
tmp^=;
if(d){
swap(now,las);
tot[now]=;
num[now]=;
for(reg i=;i<=tot[las];++i){
build(i,d-);
}
}
}
ot(ans);
return ;
} }
signed main(){
Miracle::main();
return ;
} /*
Author: *Miracle*
*/
然后这个代码又臭又长
一个很nb的写法:
1.首先不用建树,p<=i,直接v[i]=v[p]^w,一行搞定
2.不用维护决策两个点,只用维护每个元素,可能匹配的子树节点位置!以及自己的权值属于的位置
3.每次更新节点的size和a,s和k比较大小
4.更新b,
也根本不用0/1滚动,b和a数组已经区分了位置。
#include<bits/stdc++.h>
#define reg register int
#define il inline
#define fi first
#define se second
#define mk(a,b) make_pair(a,b)
#define numb (ch^'0')
#define pb push_back
#define solid const auto &
#define enter cout<<endl
#define pii pair<int,int>
using namespace std;
typedef long long ll;
template<class T>il void rd(T &x){
char ch;x=;bool fl=false;
while(!isdigit(ch=getchar()))(ch=='-')&&(fl=true);
for(x=numb;isdigit(ch=getchar());x=x*+numb);
(fl==true)&&(x=-x);
}
template<class T>il void output(T x){if(x/)output(x/);putchar(x%+'');}
template<class T>il void ot(T x){if(x<) putchar('-'),x=-x;output(x);putchar(' ');}
template<class T>il void prt(T a[],int st,int nd){for(reg i=st;i<=nd;++i) ot(a[i]);putchar('\n');} namespace Miracle{
const int N=1e6+;
ll k,s,t,v[N];
int n;
int a[N],b[N],ch[N][],sz[N];
int tot;
ll ans;
int get(int x,int c){
return ch[x][c]?ch[x][c]:ch[x][c]=++tot;
}
int main(){
rd(n);rd(k);int p;ll w;
for(reg i=;i<=n;++i) rd(p),rd(w),v[i]=v[p]^w;
for(reg i=;i<=n;++i) a[i]=b[i]=;
for(reg d=;d>=;--d){
for(reg i=;i<=tot;++i) ch[i][]=ch[i][]=sz[i]=;
tot=s=t=;
for(reg i=;i<=n;++i) sz[a[i]=get(a[i],v[i]>>d&)]++;
for(reg i=;i<=n;++i) s+=sz[ch[b[i]][v[i]>>d&]];
if(s<k) k-=s,t=,ans|=(1LL<<d);
for(reg i=;i<=n;++i) b[i]=ch[b[i]][(v[i]>>d&)^t];
}
ot(ans);
return ;
} }
signed main(){
Miracle::main();
return ;
} /*
Author: *Miracle*
*/
tql!
CF1055F Tree and XOR的更多相关文章
- 解题:CF1055F Tree and XOR
题面 树上路径是可以通过到根的路径和LCA差出来的,所以建立一棵Trie树按位贪心即可......吗? 发现空间并不够,需要我们每层现建,要记录每个数和它异或答案之后在这一层插进去的编号 #inclu ...
- [atAGC052F]Tree Vertices XOR
结论 注意到如果$x$周围有偶数个1,对$x$操作显然不会改变$a_{x}$,因此不妨强制操作的点周围要有奇数个1,不难发现此时恰好会改变该点,即令$a_{x}=a_{x}\oplus 1$ 称$\{ ...
- [atAGC052B]Tree Edges XOR
定义两点的距离$d(x,y)$为$x$到$y$路径上边权异或和,则两棵树相同当且仅当$\forall 1\le i\le n$,$d(1,i)$相同 新建一个节点0,连边$(0,1)$,初始权值为0, ...
- bzoj2006 [NOI2010]超级钢琴 (及其拓展)
bzoj2006 [NOI2010]超级钢琴 给定一个序列,求长度在 \([L,\ R]\) 之间的区间和的前 \(k\) 大之和 \(n\leq5\times10^5,\ k\leq2\times1 ...
- CF241B Friends
CF241B Friends 和Tree and Xor思路一样CF1055F Tree and XOR 直接找到第k大val,可以直接建出trie,然后按位贪心 考虑比val大的数的和 还是用b[i ...
- BZOJ3282: Tree
传送门 又是权限题= =,过了NOIp我就要去当一只权限狗! LCT裸题,get到了两个小姿势. 1.LCA操作应该在access中随时updata 2.Link操作可以更简单 void Link(i ...
- luogu P2574 XOR的艺术 (线段树)
luogu P2574 XOR的艺术 (线段树) 算是比较简单的线段树. 当区间修改时.\(1 xor 1 = 0,0 xor 1 = 1\)所以就是区间元素个数减去以前的\(1\)的个数就是现在\( ...
- 1014: [JSOI2008]火星人prefix - BZOJ
Description 火星人最近研究了一种操作:求一个字串两个后缀的公共前缀.比方说,有这样一个字符串:madamimadam,我们将这个字符串的各个字符予以标号:序号: 1 2 3 4 5 6 7 ...
- Mail.Ru Cup 2018 Round 2 Solution
A. Metro Solved. 题意: 有两条铁轨,都是单向的,一条是从左往右,一条是从右往左,Bob要从第一条轨道的第一个位置出发,Alice的位置处于第s个位置,有火车会行驶在铁轨上,一共有n个 ...
随机推荐
- Sublime Text3安装教程,配置教程,常用插件安装等方法
前言: sublimeText3的特点: 1.Sublime Text 是一款跨平台代码编辑器,在Linux.OS X和Windows下均可使用. 2.Sublime Text 是可扩展的,并包含大量 ...
- PHP学习(运算符)
PHP运算符一般分为算术运算符.赋值运算符.比较运算符.三元运算符.逻辑运算符.字符串连接运算符.错误控制运算符. 算术运算符 主要是用于进行算术运算的,例如:加法运算.减法运算.乘法运算.除法运算 ...
- Java版阿里云通信短信发送API接口实例(新)
阿里云通信(原名阿里大于)的短信服务(Short Message Service)是阿里云为用户提供的一种通信服务的能力,支持快速发送短信验证码.短信通知等. 完美支撑双11期间2亿用户,发送6亿短信 ...
- 使用multiprocessing的问题总结
Python2中的multiprocessing模块,规避了GIL(Global Interpreter Lock)带来的多线程无法实际并发的缺点,设计了几乎与threading模块一模一样的API, ...
- Hdu 4810
2014-05-02 15:53:50 题目连接 2013年南京现场赛的题目,现场的时候,排在我们前面的队伍基本都过了这题,我们后面的队伍也有不少过了这题,唯独我们没有.. 后来是Qingyu Sha ...
- 判断字符s是否为正整数和正小数
下面一段代码是判断是否为正整数和正小数的:
- Win10系统使用Docker安装oracle并通过Navicat for oracle进行登录
一.安装Docker Linux系统可以直接采用命令进行Docker安装: Win7系统安装Dokcer实际通过Boot2Docker在Windows下安装一个VirtualBox来实现: Boot2 ...
- 2018-7-5-dotnet-设计规范-·-抽象定义
title author date CreateTime categories dotnet 设计规范 · 抽象定义 lindexi 2018-07-05 15:48:20 +0800 2018-2- ...
- 洛谷P3366 【模板】最小生成树(kuskal)
#include<bits/stdc++.h> using namespace std; ; ; struct node{ int cnt,fa; }f[maxn]; inline voi ...
- part10.2-字符设备驱动模型