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的更多相关文章

  1. 解题:CF1055F Tree and XOR

    题面 树上路径是可以通过到根的路径和LCA差出来的,所以建立一棵Trie树按位贪心即可......吗? 发现空间并不够,需要我们每层现建,要记录每个数和它异或答案之后在这一层插进去的编号 #inclu ...

  2. [atAGC052F]Tree Vertices XOR

    结论 注意到如果$x$周围有偶数个1,对$x$操作显然不会改变$a_{x}$,因此不妨强制操作的点周围要有奇数个1,不难发现此时恰好会改变该点,即令$a_{x}=a_{x}\oplus 1$ 称$\{ ...

  3. [atAGC052B]Tree Edges XOR

    定义两点的距离$d(x,y)$为$x$到$y$路径上边权异或和,则两棵树相同当且仅当$\forall 1\le i\le n$,$d(1,i)$相同 新建一个节点0,连边$(0,1)$,初始权值为0, ...

  4. bzoj2006 [NOI2010]超级钢琴 (及其拓展)

    bzoj2006 [NOI2010]超级钢琴 给定一个序列,求长度在 \([L,\ R]\) 之间的区间和的前 \(k\) 大之和 \(n\leq5\times10^5,\ k\leq2\times1 ...

  5. CF241B Friends

    CF241B Friends 和Tree and Xor思路一样CF1055F Tree and XOR 直接找到第k大val,可以直接建出trie,然后按位贪心 考虑比val大的数的和 还是用b[i ...

  6. BZOJ3282: Tree

    传送门 又是权限题= =,过了NOIp我就要去当一只权限狗! LCT裸题,get到了两个小姿势. 1.LCA操作应该在access中随时updata 2.Link操作可以更简单 void Link(i ...

  7. luogu P2574 XOR的艺术 (线段树)

    luogu P2574 XOR的艺术 (线段树) 算是比较简单的线段树. 当区间修改时.\(1 xor 1 = 0,0 xor 1 = 1\)所以就是区间元素个数减去以前的\(1\)的个数就是现在\( ...

  8. 1014: [JSOI2008]火星人prefix - BZOJ

    Description 火星人最近研究了一种操作:求一个字串两个后缀的公共前缀.比方说,有这样一个字符串:madamimadam,我们将这个字符串的各个字符予以标号:序号: 1 2 3 4 5 6 7 ...

  9. Mail.Ru Cup 2018 Round 2 Solution

    A. Metro Solved. 题意: 有两条铁轨,都是单向的,一条是从左往右,一条是从右往左,Bob要从第一条轨道的第一个位置出发,Alice的位置处于第s个位置,有火车会行驶在铁轨上,一共有n个 ...

随机推荐

  1. phonegap支付宝2.0移动快捷支付插件IOS版

    坑爹的支付宝,一两年都没有更新sdk了,这两天突然更新sdk,而且更新的变化特别大,所以只能对之前的支付宝快捷支付插件重新写了一遍. 这样既顺应了支付宝的更新,同时也支持了ios8. 废话少说,集成过 ...

  2. 洛谷 P1027 Car的旅行路线 最短路+Dijkstra算法

    目录 题面 题目链接 题目描述 输入输出格式 输入格式 输出格式 输入输出样例 输入样例 输出样例 说明 思路 AC代码 总结 题面 题目链接 P1027 Car的旅行路线 题目描述 又到暑假了,住在 ...

  3. Libevent:4event loop

    一:运行loop       一旦一些events在event_base注册之后(下一节会讨论如何创建和注册events),就可以使Libevent等待events,并且在events准备好时能够通知 ...

  4. React 从零搭建项目 使用 create-react-app脚手架

    一.安装 npm install -g create-react-app 版本校验:create-react-app --version 二.创建项目 create-react-app指令默认调用np ...

  5. 数组的查找,删除 Day07

    package com.sxt.arraytest2; /* * 形参列表的作用:1.接受方法调用处传来的实参 * 2.规定了实参传入数据的类型 */ import java.util.Arrays; ...

  6. @codeforces - 1214H@ Tiles Placement

    目录 @description@ @solution@ @part - 1@ @part - 2@ @accepted code@ @details@ @description@ 给定一个 n 点的树 ...

  7. 06多次查询某区间内topk问题

            题目描述:给定一个数组,需要多次查找不同区间内的,第k大或者第k小的元素.         考虑题目是多次查找,如果采用只对查询区间内的元素进行排序的思路,然后输出第k大的数的策略,那 ...

  8. vue 后期追回的属性不更新视图问题

    this.$set(this.problemList[index], 'sub', [])   因为原始数组是有set,get而追加的没有,所以需要用这种方式   // 添加 this.$set(th ...

  9. 解决ViewState过于庞大的问题

    这里是我将ViewState持久化保持在服务器端的代码,这样ViewState不占用网络带宽,因此其存取只是服务器的磁盘读取时间.并且它很 小,可以说是磁盘随便转一圈就能同时读取好多ViewState ...

  10. EC Round 33 F. Subtree Minimum Query 主席树/线段树合并

    这题非常好!!! 主席树版本 很简单的题目,给一个按照指定节点的树,树上有点权,你需要回答给定节点的子树中,和其距离不超过k的节点中,权值最小的. 肯定首先一想,按照dfs序列建树,然后按照深度为下标 ...