Chino的数列
题解：
一道练代码能力的题目。。
首先很显然他是一道平衡树裸题
第5个操作是势能分析维护最大值最小值就可以了
另外设置虚点和noip2017队列那题一样（不过我只写过线段树）
具体细节：
1.内存池，要直接判断(!x) 因为可能进去就是0
2.输出的时候有重复的要都输出
3.search的时候要down
4.为了简单我写的down是更新自己的，而这个标记会对updata产生影响，所以updata的时候要down两个儿子
5.splay虚点变变三个点要判断点数，分三种情况 1个点，2个点，>2个点
这个洛谷评测是什么鬼。。数据测出来是对的显示wa
另外本题也可以用分块做。。以后写下试试
代码：
#include <bits/stdc++.h> using namespace std; #define rint register int #define IL inline #define rep(i,h,t) for (rint i=h;i<=t;i++) #define dep(i,t,h) for (rint i=t;i>=h;i--) #define me(x) memset(x,0,sizeof(x)) #define ll long long const int INF=1e9; const int N=4e6; queue<int> q; int s[N],count2[N],count3[N],data[N],maxa[N],mina[N],rs[N],ls[N],lazy[N],fa[N],root; #define max(a,b) (a)>(b)?(a):(b) #define min(a,b) (a)<(b)?(a):(b) int cnt,minxx=INF; struct phs{ IL void down(rint x) { if (lazy[x]==INF) return; maxa[x]=mina[x]=data[x]=lazy[x]; if (ls[x]) lazy[ls[x]]=lazy[x]; if (rs[x]) lazy[rs[x]]=lazy[x]; lazy[x]=INF; } IL void updata(rint x) { rint l=ls[x],r=rs[x],tmp=data[x]; down(l); down(r); count2[x]=count2[l]+count2[r]+count3[x]; maxa[x]=max(tmp,max(maxa[l],maxa[r])); mina[x]=min(tmp,min(mina[l],mina[r])); } IL void rotate(int x,int y) { int fa1=fa[x]; if (y==1) { if (ls[x]) fa[ls[x]]=fa1; rs[fa1]=ls[x]; } else { if (rs[x]) fa[rs[x]]=fa1; ls[fa1]=rs[x]; } fa[x]=fa[fa1]; if(fa[fa1]) { if (ls[fa[fa1]]==fa1) ls[fa[fa1]]=x; else rs[fa[fa1]]=x; } fa[fa1]=x; if (y==1) ls[x]=fa1; else rs[x]=fa1; updata(fa1); updata(x); } void dfs(int x) { if (fa[x]) dfs(fa[x]); down(x); } IL void splay(int x,int goal) { dfs(x); int fa1=fa[x]; while (fa1!=goal) { if (fa[fa1]==goal) { if (x==ls[fa1]) rotate(x,2); else rotate(x,1); } else if (fa1==ls[fa[fa1]]) if (x==ls[fa1]) rotate(fa1,2),rotate(x,2); else rotate(x,1),rotate(x,2); else if (x==rs[fa1]) rotate(fa1,1),rotate(x,1); else rotate(x,2),rotate(x,1); fa1=fa[x]; } if (goal==0) root=x; } IL int search(rint k) { rint x=root; while (x) { down(x); if (s[x]) { rint l=ls[x],r=rs[x]; rint tmp=count3[x]; count3[x]=1; rint now=ls[x]=q.front(); q.pop(); fa[now]=x; if ((count3[ls[x]]=tmp/2)!=1) s[now]=1; ls[now]=l; fa[l]=now; data[now]=data[x]; updata(now); if (tmp>2) { now=rs[x]=q.front(); q.pop(); data[rs[x]]=data[x]; fa[now]=x; if ((count3[now]=tmp-(tmp/2)-1)!=1) s[now]=1; rs[now]=r; fa[r]=now; updata(now); } down(x); s[x]=0; } if (count2[ls[x]]+1==k) return(x); if (count2[ls[x]]+1<k) k-=count2[ls[x]]+1,x=rs[x]; else x=ls[x]; } } void dfs2(int x) { if (!x) return; cnt++; q.push(x); dfs2(ls[x]); dfs2(rs[x]); ls[x]=rs[x]=count2[x]=s[x]=data[x]=fa[x]=maxa[x]=mina[x]=0; count3[x]=1; lazy[x]=INF; } IL void split(int x,int y) { splay(x,0); splay(y,x); } IL void cl(int x) { down(x); int k1=sqrt(maxa[x]); int k2=sqrt(mina[x]); if (k1==k2) { lazy[x]=k1; down(x); return; } data[x]=sqrt(data[x]); if (ls[x]) cl(ls[x]); if (rs[x]) cl(rs[x]); updata(x); } IL void ycl() { root=1; rs[1]=2; fa[2]=1; count2[1]=2; count2[2]=1; } IL void ou(int x) { down(x); if (ls[x]) ou(ls[x]); if (x!=1&&x!=2) rep(i,1,count3[x]) printf("%d ",data[x]); if (rs[x]) ou(rs[x]); } }S; int main() { freopen("1.in","r",stdin); freopen("1.out","w",stdout); ios::sync_with_stdio(false); int m; cin>>m; maxa[0]=-INF; mina[0]=INF; int len=0; rep(i,3,N-1) q.push(i); rep(i,0,N-1) lazy[i]=INF,count3[i]=1; S.ycl(); rep(i,1,m) { minxx=min(minxx,q.size()); int x,y,z,k; cin>>k; if (k==1) { cin>>x>>y; rint now=S.search(len+1); S.splay(2,0); S.splay(now,2); rint now2=rs[now]=q.front(); q.pop(); count3[now2]=x; data[now2]=y; fa[now2]=now; if (x!=1) s[now2]=1; S.splay(now2,0); len+=x; } if (k==2) { int x; cin>>x; len-=x; x=S.search(x+2); S.split(1,x); // S.dfs2(ls[x]); ls[x]=0; S.splay(x,0); } if (k==3) { int kk; cin>>x>>y>>z>>kk; rint t=y-x+1; rint ans=1; while(z) { if (z&1) ans=(1ll*ans*t)%kk; t=(1ll*t*t)%kk; z>>=1; } int l=S.search(x); int r=S.search(y+2); S.split(l,r); lazy[ls[r]]=ans; S.splay(ls[r],0); } if (k==4) { int x,y; cin>>x>>y; //cout<<x<<" "<<y<<endl; len-=y-x; x=S.search(x),y=S.search(y+2); S.split(x,y); S.down(ls[y]); int now=ls[y]; int mm=maxa[now]; // S.dfs2(ls[now]); S.dfs2(rs[now]); ls[now]=rs[now]=0; lazy[now]=INF; count3[now]=1; s[now]=0; data[now]=mm; S.splay(now,0); printf("%d\n",mm); } if (k==5) { cin>>x>>y; x=S.search(x),y=S.search(y+2); S.split(x,y); int now=ls[y]; S.cl(now); } } S.ou(root); // printf("\n%d\n%d",cnt,minxx); return 0; }
