标签:
区间求众数,强制在线。
分块处理,求出每一块区间的众数。对于一次询问的区间,其中众数要么是中间大块的众数,要么是两边小块中的数。
扫描所求区间中的边缘小块,求其中每个数在区间内的出现次数,找出次数最多的就是区间众数了。
求数的出现次数,可以预先用vector或者普通数组存这个数在序列中的每一个出现位置,然后二分查找即可。
不管什么算法,逢二分必挂的被动技能又强行发动了,调二分调了半小时……
用时比别人多了10000+ms,也是迷。
1 /*by SilverN*/ 2 #include<iostream> 3 #include<algorithm> 4 #include<cstring> 5 #include<cstdio> 6 #include<cmath> 7 #include<vector> 8 #include<map> 9 using namespace std; 10 const int mxn=50010; 11 const int block=180; 12 int read(){ 13 int x=0,f=1;char ch=getchar(); 14 while(ch<‘0‘ || ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();} 15 while(ch>=‘0‘ && ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();} 16 return x*f; 17 } 18 int n,m; 19 // 20 vector<int>pos[mxn];//数字出现位置 21 int v[mxn];//离散化标号对应的数字 22 map<int,int>mp;int mct=0;//数字对应的离散化标号 23 // 24 int f[510][510]; 25 int b[mxn];//分块 26 // 27 int a[mxn]; 28 int cnt[mxn]; 29 void init(int x){ 30 memset(cnt,0,sizeof cnt); 31 int i,j; 32 int mx=0,num=0; 33 for(i=(x-1)*block+1;i<=n;i++){ 34 cnt[v[i]]++; 35 if(cnt[v[i]]>mx || (cnt[v[i]]==mx && a[i]<num)){ 36 mx=cnt[v[i]]; 37 num=a[i]; 38 } 39 f[x][b[i]]=num; 40 } 41 return; 42 } 43 int find(int x,int L,int R){ 44 /* x=mp[x]; 45 int t=upper_bound(pos[x].begin(),pos[x].end(),R)-lower_bound(pos[x].begin(),pos[x].end(),L); 46 return t;*/ 47 x=mp[x]; 48 int l=0,r=pos[x].size()-1; 49 int ansl,ansr; 50 while(l<=r){ 51 int mid=(l+r)>>1; 52 if(pos[x][mid]<=R)l=mid+1; 53 else r=mid-1; 54 } 55 ansr=l; 56 l=0,r=pos[x].size()-1; 57 while(l<=r){ 58 int mid=(l+r)>>1; 59 if(pos[x][mid]<L)l=mid+1; 60 else r=mid-1; 61 } 62 ansl=l; 63 // printf("ask:%d res:%d %d\n",x,ansl,ansr); 64 return ansr-ansl; 65 } 66 int query(int l,int r){ 67 int mxnum=f[b[l]+1][b[r]-1]; 68 int mx=find(mxnum,l,r); 69 int i,j; 70 int ed=min(b[l]*block,r); 71 for(i=l;i<=ed;i++){ 72 int tmp=find(a[i],l,r); 73 if(tmp>mx || (tmp==mx && a[i]<mxnum)){ 74 mx=tmp; mxnum=a[i]; 75 } 76 } 77 if(b[l]!=b[r]) 78 for(i=(b[r]-1)*block+1;i<=r;i++){ 79 int tmp=find(a[i],l,r); 80 if(tmp>mx || (tmp==mx && a[i]<mxnum)){ 81 mx=tmp; mxnum=a[i]; 82 } 83 } 84 return mxnum; 85 } 86 int main(){ 87 n=read();m=read(); 88 int i,j; 89 for(i=1;i<=n;i++){ 90 a[i]=read(); 91 if(!mp[a[i]])mp[a[i]]=++mct; 92 v[i]=mp[a[i]]; 93 pos[v[i]].push_back(i); 94 } 95 for(i=1;i<=n;++i){b[i]=(i-1)/block+1;} 96 for(i=1;i<=b[n];i++){init(i);} 97 int x=0; 98 int ql,qr; 99 while(m--){ 100 ql=read();qr=read(); 101 ql=(ql+x-1)%n+1; 102 qr=(qr+x-1)%n+1; 103 if(ql>qr)swap(ql,qr); 104 x=query(ql,qr); 105 printf("%d\n",x); 106 } 107 return 0; 108 }
标签:
原文地址:http://www.cnblogs.com/SilverNebula/p/5937144.html