BZOJ 3956: Count 主席树 可持久化线段树 单调栈

https://www.lydsy.com/JudgeOnline/problem.php?id=3956

从描述可以得到性质: 每个好点对 ( 除了差值为1的好点对 ) 中间的数 ( i , j ) 一定有一个最大值, 视这个值控制这个好点对, 那么每个值最多只控制一个好点对.

然后我们就可以通过单调栈 ( 求出lp[i]左边第一个大于i的数的坐标, rp[i]右边第一个大于i的数的坐标, lf[i]左边第一个大于等于i的数的坐标 (为避免重复只统计有重复区间的值里的最左边一个，即lp=lf则不统计) )  找到每个点控制的好点对 ( 也就是找到所有好点对 ) .

然后用主席树维护询问(对每个l建r作为查找值的树).

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 using namespace std;
 const int maxn=;
 int n,m,T;
 int a[maxn]={},sta[maxn]={},tail=;
 int lp[maxn]={},rp[maxn]={},lf[maxn]={},rt[maxn]={};
 int lc[maxn*]={},rc[maxn*]={},siz[maxn*]={},tot=;
 struct nod{
     int x,y;
 }e[maxn];int cnt=;
 bool mcmp(nod aa,nod bb){ return aa.x<bb.x; }
 int read(){
     int w=,f=;char ch=getchar();
     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
     while(ch>=''&&ch<=''){w=w*+ch-'';ch=getchar();}
     return w*f;
 }
 void build(int &x,int y,int l,int r,int z){
     x=++tot;siz[x]=siz[y]+;lc[x]=lc[y];rc[x]=rc[y];
     if(l==r)return;
     int mid=(l+r)/;
     if(z<=mid)build(lc[x],lc[y],l,mid,z);
     else build(rc[x],rc[y],mid+,r,z);
 }
 int getsum(int x,int y,int l,int r,int z,int t){
     if(z<=l&&r<=t)return siz[y]-siz[x];
     int mid=(l+r)/,ans=;
     if(z<=mid)ans=getsum(lc[x],lc[y],l,mid,z,t);
     if(t>mid)ans+=getsum(rc[x],rc[y],mid+,r,z,t);
     return ans;
 }
 inline void fir(){
     int i;
     a[]=a[n+]=(<<);
     for(i=;i<=n;++i){
         while(a[sta[tail]]<=a[i])tail--;
         lp[i]=sta[tail];sta[++tail]=i;
     }
     sta[]=n+;tail=;
     for(i=n;i>;--i){
         while(a[sta[tail]]<=a[i])tail--;
         rp[i]=sta[tail];sta[++tail]=i;
     }
     tail=;sta[]=;
     for(i=;i<=n;++i){
         while(a[sta[tail]]<a[i])tail--;
         lf[i]=sta[tail];sta[++tail]=i;
     }
 }
 inline void init(){
     int i,j;
     for(i=;i<=n;++i){
         if(lp[i]&&rp[i]<=n&&lp[i]==lf[i]){//最近的比a[i]大的数与a[i]间没有与a[i]相等的数，避免重复记录
             e[++cnt].x=lp[i];e[cnt].y=rp[i];//lp[i]和rp[i]必须合法
         }
     }sort(e+,e+cnt+,mcmp);
     for(i=j=;i<=n;++i){
         rt[i]=rt[i-];
         while(j<=cnt&&e[j].x==i){build(rt[i],rt[i],,n,e[j].y);j++;}
     }
 }
 int main(){
     int i,x,y,ans=;
     n=read();m=read();T=read();
     for(i=;i<=n;++i)a[i]=read();
     fir(); init();
     for(i=;i<=m;++i){
         x=read();y=read();
         if(T){x=(x+ans-)%n+;y=(y+ans-)%n+;}
         if(x>y)swap(x,y);
         ans=getsum(rt[x-],rt[y],,n,x,y)+y-x;
         printf("%d\n",ans);
     }
     return ;
 }