BZOJ3489: A simple rmq problem
阿新 • • 發佈:2018-08-04
dig con tor 更新 for 數據結構 屬於 pac acc
隊友丟過來的毒瘤數據結構 ...KDtree模板題?告辭 不會KDtree ....那我只能考慮別的做法了 很顯然 每個a[i] 都有上一個出現的位置last和下一個出現的位置next 對於查詢的區間而言 當且僅當 [l,r]屬於(last,next)的子集時並且了l<=i<=r這點的值才會產生貢獻 因此 我們考慮主席樹維護以pos為位置的歷史版本 對於主席樹中的[next,n]區間更新對應的動態開點的線段樹中插入[i,a[i]] 查詢就是普通的樹套樹的查詢即可 空間(nlog^2n) 時間(nlog^2n)
做法:主席樹+動態開點線段樹
/************************************************************** Problem: 3489 User: c20161007 Language: C++ Result: Accepted Time:23948 ms Memory:476340 kb ****************************************************************/ #include <bits/stdc++.h> const int MAXN=1e5+10; #define ll long long using namespace std; ll read(){ ll x=0,f=1;char ch=getchar(); while(!isdigit(ch)){if(ch==‘-‘)f=-1;ch=getchar();} while(isdigit(ch))x=x*10+ch-‘0‘,ch=getchar(); return f*x; } int vis[MAXN]; typedef struct node{ int l,r,vul,id; friend bool operator<(node aa,node bb){ return aa.l<bb.l; } }node; node que[MAXN];int rt[MAXN]; typedef struct Node{ int l,r,id; }Node; Node d[32*MAXN];int a[MAXN]; int cnt,cnt1,cnt2,n,m; typedef struct nOde{ int l,r,max; }nOde; nOde D[371*MAXN]; void update1(int &x,int y,int l,int r,int t,int vul){ x=++cnt2;D[x]=D[y];D[x].max=max(vul,D[x].max); //cout<<l<<" "<<r<<" "<<t<<" "<<vul<<endl; if(l==r)return ; int mid=(l+r)>>1; if(t<=mid)update1(D[x].l,D[y].l,l,mid,t,vul); else update1(D[x].r,D[y].r,mid+1,r,t,vul); } int t1; void update(int &x,int y,int l,int r,int t,int vul){ x=++cnt;d[x]=d[y];update1(d[x].id,d[y].id,1,n,t1,vul); // cout<<l<<":::"<<r<<" "<<d[x].id<<endl; if(l==r)return ; int mid=(l+r)>>1; if(t<=mid)update(d[x].l,d[y].l,l,mid,t,vul); else update(d[x].r,d[y].r,mid+1,r,t,vul); } int ans=0; int Ql,Qr; void querty1(int x,int l,int r,int ql,int qr){ if(ql<=l&&r<=qr){ans=max(ans,D[x].max);return ;} int mid=(l+r)>>1; if(ql<=mid)querty1(D[x].l,l,mid,ql,qr); if(qr>mid)querty1(D[x].r,mid+1,r,ql,qr); } void querty(int x,int l,int r,int ql,int qr){ // cout<<l<<" "<<r<<" "<<d[x].id<<endl; if(ql<=l&&r<=qr){querty1(d[x].id,1,n,Ql,Qr);return ;} int mid=(l+r)>>1; if(ql<=mid)querty(d[x].l,l,mid,ql,qr); if(qr>mid)querty(d[x].r,mid+1,r,ql,qr); } int main(){ n=read();m=read(); for(int i=1;i<=n;i++)que[i].vul=read(),a[i]=que[i].vul,que[i].id=i; for(int i=1;i<=n;i++)que[i].l=vis[a[i]]+1,vis[a[i]]=i; for(int i=1;i<=n;i++)vis[i]=n+1; for(int i=n;i>=1;i--)que[i].r=vis[a[i]]-1,vis[a[i]]=i; sort(que+1,que+n+1); // for(int i=1;i<=n;i++)cout<<que[i].l<<" "<<que[i].r<<" "<<que[i].vul<<endl; int L=1; for(int i=1;i<=n;i++){ // cout<<que[i].l<<" "<<L<<endl; while(L<que[i].l)rt[L]=rt[L-1],L++; t1=que[i].id; if(que[i].l!=que[i-1].l)update(rt[que[i].l],rt[que[i].l-1],1,n,que[i].r,que[i].vul),L++; else update(rt[que[i].l],rt[que[i].l],1,n,que[i].r,que[i].vul); // cout<<que[i].l<<" "<<rt[que[i].l]<<" "<<L<<endl; } for(;L<=n;L++)rt[L]=rt[L-1]; // for(int i=1;i<=n;i++)cout<<rt[i]<<" "; // cout<<endl; //for(int i=1;i<=n;i++)vis[i]=0; //for(int i=1;i<=n;i++)L[i]=R[vis[a[i]]]+1,vis[a[i]]=i; int l,r,lx,rx,res=0; for(int i=1;i<=m;i++){ l=read();r=read();lx=l;rx=r; //cout<<l<<" "<<r<<endl; lx=min((l+res)%n+1,(r+res)%n+1); rx=max((l+res)%n+1,(r+res)%n+1); // cout<<lx<<":::"<<rx<<" "<<rt[lx]<<endl; Ql=lx;Qr=rx; ans=0;querty(rt[lx],1,n,rx,n);res=ans; //ans=0;querty(1,1,n,r,n);res printf("%d\n",res); } return 0; }
BZOJ3489: A simple rmq problem