BZOJ2453: 维护队列
n<=10000个数m<=10000次操作:1、L~R中有多少种数字,2、单点修改,修改次数<=1000。
一种通法是记pre[i]表示i前最近的和a[i]相等的数字,然后问题变成:[L,R]中有多少个pre[i]<L。
方法一:pre数组的动态维护可以用set。对每个数开一个set存所有出现位置,然后每次改变一个数的时候就是set的单点查询、单点删除、单点插入、查前驱后继,以此变成若干个pre的单点修改。然后树状数组套主席树即可。未写。
方法二:分块,每个块里的pre排个序,完整块查询时直接在排序好的数组里二分即可。
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<stdio.h> 2 #include<string.h> 3 #include<stdlib.h> 4 #include<algorithm> 5 #include<math.h> 6 //#include<iostream> 7 using namespace std; 8 9 int n,m,q; 10 #define maxn 10011 11 #define maxm 111 12 int a[maxn],b[maxn],v[maxn],last[maxn*maxm],bel[maxn],tot; 13 void erase() {for (int i=1;i<=n;i++) last[v[i]]=0;} 14 void reset() 15 { 16 for (int i=1;i<=n;i++) a[i]=last[v[i]],b[i]=a[i],last[v[i]]=i; 17 for (int i=1;i<=n;i+=m) sort(b+i,b+i+min(m,n-i+1)); 18 } 19 int query(int x,int y) 20 { 21 int ans=0; 22 if (bel[x]==bel[y]) for (int i=x;i<=y;i++) ans+=(a[i]<x); 23 else 24 { 25 for (int i=bel[x]+1;i<bel[y];i++) 26 { 27 int L=(i-1)*m+1,R=i*m+1; 28 while (L<R) 29 { 30 const int mid=(L+R)>>1; 31 if (b[mid]>=x) R=mid; 32 else L=mid+1; 33 } 34 ans+=L-(i-1)*m-1; 35 } 36 for (int i=x;bel[i]==bel[x];i++) ans+=(a[i]<x); 37 for (int i=y;bel[i]==bel[y];i--) ans+=(a[i]<x); 38 } 39 return ans; 40 } 41 int main() 42 { 43 scanf("%d%d",&n,&q); 44 m=(int)sqrt(n); 45 for (int i=1;i<=n;i++) bel[i]=(i-1)/m+1; 46 tot=bel[n]; 47 48 for (int i=1;i<=n;i++) scanf("%d",&v[i]); 49 reset(); 50 51 char c;int x,y; 52 while (q--) 53 { 54 while ((c=getchar())!='R' && c!='Q'); 55 scanf("%d%d",&x,&y); 56 if (c=='R') erase(),v[x]=y,reset(); 57 else printf("%d\n",query(x,y)); 58 } 59 return 0; 60 }