博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
bzoj 3196 Tyvj 1730 二逼平衡树
阅读量:5871 次
发布时间:2019-06-19

本文共 3067 字,大约阅读时间需要 10 分钟。

题目:

第一道的树套树。

线段树套平衡树。找了一个当模板。感觉还好啦。但跑得很慢。卡时过了。

1.好神奇呀,竟然都要加-INF和INF。别忘了求siz时-1;

2.去右儿子找的时候别忘了减去左儿子和自己的siz;

3.凡是有splay的地方就别忘了 & ;

4.想要 < 的就把它放在return y?y:cr里,想要 <= 的也就把它放在那个地方,大约就是这样就行;

5.线段树的 tt 别和Splay的 tot 混用,不然线段树节点的数组就得……

6.2e6是怎么算的呀?1e6为什么不行?

#include
#include
#include
using namespace std;const int N=5e4+5,M=2e6+5,INF=1e9;int n,m,t[N],tot,tt,fa[M],c[M][2],val[M],siz[M];//struct Node{ int l,r,ls,rs,rt;}a[N<<1];void updt(int cr){ siz[cr]=siz[c[cr][0]]+siz[c[cr][1]]+1;}void rotate(int x,int &rt){ int y=fa[x],z=fa[y],d=(x==c[y][1]); if(y==rt)rt=x; else c[z][y==c[z][1]]=x; fa[y]=x;fa[x]=z;fa[c[x][!d]]=y; c[y][d]=c[x][!d];c[x][!d]=y; updt(y);updt(x);}void splay(int x,int &rt){ while(x!=rt) { int y=fa[x],z=fa[y]; if(y!=rt) { if((x==c[y][0])^(y==c[z][0]))rotate(x,rt); else rotate(y,rt); } rotate(x,rt); }}int find_k(int cr,int x){ if(siz[c[cr][0]]==x-1)return cr; if(siz[c[cr][0]]>=x)return find_k(c[cr][0],x); return find_k(c[cr][1],x-c[cr][0]-1);}int find_pre(int cr,int x){ if(!cr)return 0; if(val[cr]>=x)return find_pre(c[cr][0],x); int y=find_pre(c[cr][1],x); return y?y:cr;}void insert(int &rt,int x){ int y=find_pre(rt,x);splay(y,rt); y=find_k(c[rt][1],1);splay(y,c[rt][1]); fa[++tot]=y;c[y][0]=tot;siz[tot]=1;val[tot]=x; updt(y);updt(rt);}void del(int &rt,int x){ int y=find_pre(rt,x);splay(y,rt); y=find_k(c[rt][1],2);splay(y,c[rt][1]); fa[c[y][0]]=0;c[y][0]=0; updt(y);updt(rt);}void build(int l,int r,int cr){ a[cr].l=l;a[cr].r=r;a[cr].rt=++tot; fa[tot]=0;c[tot][1]=tot+1;c[tot][0]=0;val[tot]=-INF;siz[tot]=2; fa[++tot]=tot-1;c[tot][1]=c[tot][0]=0;val[tot]=INF;siz[tot]=1; for(int i=l;i<=r;i++)insert(a[cr].rt,t[i]); if(l==r)return; int mid=((l+r)>>1); a[cr].ls=++tt;build(l,mid,tt); a[cr].rs=++tt;build(mid+1,r,tt);}int find_Suc(int rt,int x) //Suc表示找的不是正常后继(>=){ if(!rt)return 0; if(val[rt]
=L&&a[cr].r<=R)return query_rnk(a[cr].rt,x); int mid=((a[cr].l+a[cr].r)>>1),ret=0; if(mid>=L)ret+=query_rnk(L,R,a[cr].ls,x); if(mid
=L&&a[cr].r<=R)return query_Rnk(a[cr].rt,x); int mid=((a[cr].l+a[cr].r)>>1),ret=0; if(mid>=L)ret+=query_Rnk(L,R,a[cr].ls,x); if(mid
>1); if(P<=mid)mdfy(P,a[cr].ls,x); else mdfy(P,a[cr].rs,x);}int query_pre(int L,int R,int cr,int x){ if(a[cr].l>=L&&a[cr].r<=R)return val[find_pre(a[cr].rt,x)]; int ret=-INF,mid=((a[cr].l+a[cr].r)>>1); if(mid>=L)ret=max(ret,query_pre(L,R,a[cr].ls,x)); if(mid
=L&&a[cr].r<=R)return val[find_suc(a[cr].rt,x)]; int ret=INF,mid=((a[cr].l+a[cr].r)>>1); if(mid>=L)ret=min(ret,query_suc(L,R,a[cr].ls,x)); if(mid
>1); if(query_Rnk(x,y,1,mid)>=k)ans=mid,r=mid-1;//不+1:包含了全部==自己的 else l=mid+1; } printf("%d\n",ans); } if(op==3)mdfy(x,1,y),t[x]=y; if(op==4)printf("%d\n",query_pre(x,y,1,k)); if(op==5)printf("%d\n",query_suc(x,y,1,k)); } return 0;}

 

转载于:https://www.cnblogs.com/Narh/p/9221951.html

你可能感兴趣的文章
继续上一篇剩下的例子
查看>>
检查脚本
查看>>
为什么PHP能够受到大家追捧和喜爱,又为什么饱受嘲讽?
查看>>
网线(水晶头)制作流程 请附件下载
查看>>
基于APPIUM的移动自动化测试
查看>>
使用PowerShel导入和导出Hyper-v虚拟机
查看>>
NodeJs+Qunit的使用方式
查看>>
求两个数的最大公约数
查看>>
for循环语句的用法
查看>>
ERROR 1820 (HY000): You must reset your password
查看>>
解决QT5 VS2010调试时不能显示字符串的内容
查看>>
快速构建Windows 8风格应用34-构建Toast通知
查看>>
Ubuntu下部署zabbix 开源监控系统
查看>>
M283-bsp包问题
查看>>
Android的手机震动
查看>>
mysql 数据库的维护,优化
查看>>
pssh安装使用
查看>>
Erlang资料
查看>>
Linux查询命令man手册各章节解释
查看>>
linux下面实现执行rm命令,显示do not use rm command
查看>>