正题

题目链接:http://www.ybtoj.com.cn/problem/652

题目大意

定义一个元素为一个有序集合包含两个元素C={A,B}C=\{A,B\}C={A,B}
集合C={A,B}C=\{A,B\}C={A,B}的大小以AAA为第一关键字，BBB为第二关键字比较大小。
开始有两个元素S={S,S},T={T,T}S=\{S,S\},T=\{T,T\}S={S,S},T={T,T}且S<TS<TS<T。

然后nnn次加入一个新的由两个之前的元素依次组成的新元素，求出现过的元素小于等于它的有多少个。

1≤n≤5×1041\leq n\leq 5\times 10^41≤n≤5×104

解题思路

如果递归比较是O(n)O(n)O(n)的显然不行，但是我们比较新的元素和旧的元素大小的时候如果我们可以知道以前元素的大小关系就可以快速比较（因为新的元素由旧的元素组成）

所以相当于我们要动态维护大小关系，因为要插入好像只能用平衡树。

然后要在树上查询两个点的大小关系，因为你要在平衡树上边移动边查询，所以不能用查询的时候结构会改变的平衡树（SplaySplaySplay之类的，虽然原本SplaySplaySplay是最方便的），正解是替罪羊的，反正这里用了TreapTreapTreap。

比较大小的时候直接在TreapTreapTreap上暴力跳找LCALCALCA就好了，深度是logloglog级别的，然后不能动态维护深度所以有点麻烦

时间复杂度O(nlog⁡2n)O(n\log^2 n)O(nlog2n)

code

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=51000;
int n,qfn,tot,a[N],b[N],p[N],v[N];
int t[N][2],siz[N],cnt[N],dat[N],fa[N];
int cmp(int x,int y){++qfn;if(x==y)return 2;while(x){if(fa[x]==y)return t[y][0]==x;v[x]=qfn;x=fa[x];}while(y){if(v[fa[y]]==qfn)return t[fa[y]][1]==y;y=fa[y];}return 0;
}
int cap(int x,int y){int tmp=cmp(a[x],a[y]);if(tmp!=2)return tmp;return cmp(b[x],b[y]);
}
void PushUp(int x){siz[x]=siz[t[x][0]]+siz[t[x][1]]+cnt[x];return;
}
void zig(int &x){int y=t[x][0];t[x][0]=t[y][1];fa[t[y][1]]=x;t[y][1]=x;fa[y]=fa[x];fa[x]=y;x=y;PushUp(t[x][1]);PushUp(x);return;
}
void zag(int &x){int y=t[x][1];t[x][1]=t[y][0];fa[t[y][0]]=x;t[y][0]=x;fa[y]=fa[x];fa[x]=y;x=y;PushUp(t[x][0]);PushUp(x);return;
}
void Insert(int &x,int pos){if(!x){x=++tot;p[pos]=x;a[tot]=a[pos];b[tot]=b[pos];cnt[x]=siz[x]=1;dat[x]=rand();return;}int tmp=cap(pos,x),sum=0;if(tmp==2){p[pos]=x;cnt[x]++;PushUp(x);return;}else if(tmp){Insert(t[x][0],pos);fa[t[x][0]]=x;if(dat[t[x][0]]>dat[x])zig(x);}else{Insert(t[x][1],pos);fa[t[x][1]]=x;if(dat[t[x][1]]>dat[x])zag(x);}PushUp(x);return;
}
int Query(int x){int ans=siz[t[x][0]]+cnt[x];while(x){if(t[fa[x]][0]!=x)ans+=siz[t[fa[x]][0]]+cnt[fa[x]];x=fa[x];}return ans;
}
int main()
{//  freopen("comparison.in","r",stdin);
//  freopen("comparison.out","w",stdout);srand(19260817);scanf("%d",&n);n++;t[1][1]=n+1;p[1]=1;a[1]=b[1]=1;p[n+1]=n+1;a[n+1]=b[n+1]=n+1;cnt[1]=cnt[n+1]=1;fa[n+1]=tot=1;PushUp(n+1);PushUp(1);int rt=1;for(int i=2;i<=n;i++){scanf("%d%d",&a[i],&b[i]);a[i]++;b[i]++;a[i]=p[a[i]];b[i]=p[b[i]];Insert(rt,i);printf("%d\n",Query(p[i]));}return 0;
}

YbtOJ#652-集合比较【Treap】相关推荐

【bzoj2555】Substring【后缀平衡树入门】
传送门 (bzoj上不去我也很无奈啊) 题意:维护一个字符串,支持后面加字符串,给定串询问出现次数.强制在线. 数据范围:暴力跑不过前置知识:重量平衡树众所周知,平衡树都用了一些策略保证平衡. 平 ...
数据结构之Treap
1. 概述同splay tree一样,treap也是一个平衡二叉树,不过Treap会记录一个额外的数据,即优先级.Treap在以关键码构成二叉搜索树的同时,还按优先级来满足堆的性质.因而,Treap ...
[POI2008] Poc (原名 Trians) Treap+Hash
这个题和千山鸟飞绝体现出了一种用平衡树解决动态集合问题,主要套路就是蜜汁标记. 这个题我一开始用替罪羊树搞了一下对了28个点,后来我换成了Treap一搞对了14个点,再后来发现被卡了Hash我竟然在自 ...
Java集合—TreeMap底层原理
原文出自:http://cmsblogs.com/?p=1013.尊重作者的成果,转载请注明出处! 个人站点:http://cmsblogs.com TreeMap的实现是红黑树算法的实现,所以要了解 ...
Java集合--TreeMap
转载请注明出处:http://www.cnblogs.com/skywang12345/admin/EditPosts.aspx?postid=3310928 第1部分 TreeMap介绍 TreeM ...
对 Redis 中的有序集合SortedSet的理解
本篇说一下Redis中的有序集合类型,曾几何时,我们想把所有数据存到内存中的数据结构中,但为了多机器共享内存,不得不将这块内存包装成wcf单独部署,同时还要考虑怎么序列化,烦心事太多太多...后 ...
[2021-09-02 contest]CF1251C，可达性统计(bitset优化dp)，Boomerang Tournament(状压dp)，小蓝的好友(mrx)(treap平衡树)
文章目录 CF1251C Minimize The Integer acwing164:可达性统计 Facebook Hacker Cup 2016 Round 1 Boomerang Tournam ...
YbtOJ#791-子集最值【三维偏序】
正题题目链接:http://www.ybtoj.com.cn/contest/123/problem/1 题目大意给出333个长度为nnn的排列A,B,CA,B,CA,B,C.然后一个下标集合SS ...
YbtOJ#883-最大的割【带修线性基】
正题题目链接:http://www.ybtoj.com.cn/contest/118/problem/3 解题思路给出nnn个点,mmm次动态插入一条无向边询问:割掉一些边使得图中至少两点不连通, ...

YbtOJ#652-集合比较【Treap】

正题

题目大意

解题思路

code

YbtOJ#652-集合比较【Treap】相关推荐

最新文章

热门文章