【bzoj2724】【蒲公英】【分块】
来源:互联网 发布:网络推广新渠道 编辑:程序博客网 时间:2024/05/17 03:23
Description
Input
修正一下
l = (l_0 + x - 1) mod n + 1, r = (r_0 + x - 1) mod n + 1
Output
Sample Input
6 3
1 2 3 2 1 2
1 5
3 6
1 5
1 2 3 2 1 2
1 5
3 6
1 5
Sample Output
1
2
1
2
1
HINT
修正下:
n <= 40000, m <= 50000
题解:
区间众数.
对序列分块,t[i][j]表示第i块到第j块的众数.
把序列按权值和位置排序之后就可以二分查找一个数在一段区间内的出现次数.
所以块之外的部分暴力二分查找即可.
代码:
#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#include<cmath> #define N 40010#define M 210using namespace std;int a[N],bl[N],n,m,bk,x,y,cnt,ans,h[N],l[N],r[N],sum[N],t[M][M];struct block{int l,r;}st[M];struct xu{int x,p;}q[N];void getbk(){ bk=sqrt(n); if (n%bk) cnt=n/bk+1;else cnt=n/bk; for (int i=1;i<=n;i++) bl[i]=(i-1)/bk+1; for (int i=1;i<=cnt;i++) st[i].l=(i-1)*bk+1,st[i].r=i*bk; }int find(int x){ int l=1,r=n; while(l<=r){ int mid=(l+r)>>1; if (x<h[mid]) r=mid-1; else if (x>h[mid]) l=mid+1; else return mid; }}bool cmp(xu a,xu b){ if (a.x==b.x) return a.p<b.p; else return a.x<b.x;}void pre(){ for (int i=1;i<=n;i++) q[i].x=a[i],q[i].p=i; sort(q+1,q+n+1,cmp); for (int i=1;i<=n;i++){ if (!l[q[i].x]) l[q[i].x]=i; r[q[i].x]=i; } int mx(0),v; for (int i=1;i<=cnt;i++){ memset(sum,0,sizeof(sum));mx=0; for (int j=st[i].l;j<=n;j++){ sum[a[j]]++; if (sum[a[j]]>mx||(sum[a[j]]==mx)&&a[j]<v){ mx=sum[a[j]];v=a[j]; } t[i][bl[j]]=v; } }}int up(int k,int x,int y){ int l=x,r=y,ans; while(l<=r){ int mid=(l+r)>>1; if (q[mid].p>k) r=mid-1; else {ans=mid;l=mid+1;} } return ans;}int down(int k,int x,int y){ int l=x,r=y,ans; while (l<=r){ int mid=(l+r)>>1; if (q[mid].p<k) l=mid+1; else {ans=mid;r=mid-1;} } return ans;}int query(int x,int y){ int mx(0),v; if (bl[x]==bl[y]){ for (int i=x;i<=y;i++){ int temp=up(y,l[a[i]],r[a[i]])-down(x,l[a[i]],r[a[i]])+1; if (temp>mx||(temp==mx)&&(a[i]<v)){ v=a[i];mx=temp; } } } else{ if (bl[x]+1<bl[y]){ v=t[bl[x]+1][bl[y]-1]; mx=up(y,l[v],r[v])-down(x,l[v],r[v])+1; } for (int i=x;i<=st[bl[x]].r;i++){ int temp=up(y,l[a[i]],r[a[i]])-down(x,l[a[i]],r[a[i]])+1; if (temp>mx||(temp==mx)&&(a[i]<v)){ v=a[i];mx=temp; } } for (int i=st[bl[y]].l;i<=y;i++){ int temp=up(y,l[a[i]],r[a[i]])-down(x,l[a[i]],r[a[i]])+1; if (temp>mx||(temp==mx)&&(a[i]<v)){ v=a[i];mx=temp; } } } return v;}int main(){ scanf("%d%d",&n,&m); for (int i=1;i<=n;i++) scanf("%d",&a[i]); getbk(); for (int i=1;i<=n;i++) h[i]=a[i]; sort(h+1,h+n+1); for (int i=1;i<=n;i++) a[i]=find(a[i]); pre(); for (int i=1;i<=m;i++){ scanf("%d%d",&x,&y); x=(x+ans-1)%n+1;y=(y+ans-1)%n+1; if (x>y) swap(x,y); ans=h[query(x,y)]; printf("%d\n",ans); }}
0 0
- 【bzoj2724】【蒲公英】【分块】
- [分块] [BZOJ2724] [Violet 6] 蒲公英
- bzoj2724 [Violet 6]蒲公英 分块
- [BZOJ2724][Violet 6]蒲公英(分块)
- BZOJ2724: [Violet 6]蒲公英(分块)
- [bzoj2724]蒲公英
- BZOJ2724 [Violet 6]蒲公英 解题报告【数据结构】【分块】
- 【Violet 6】【BZOJ2724】蒲公英
- 【bzoj2724】[Violet 6]蒲公英
- bzoj2724: [Violet 6]蒲公英
- BZOJ2724: [Violet 6]蒲公英
- 【bzoj2724】[Violet 6]蒲公英
- bzoj2724 蒲公英【解法一】
- bzoj2724 蒲公英【解法二】
- BZOJ2724: [Violet 6]蒲公英
- bzoj2724: [Violet 6]蒲公英
- bzoj2724[Violet 6][蒲公英]
- 【bzoj2724】[Violet 6]蒲公英
- Java 高效读取大数据文件—最优方法
- 使用 Linux 文件恢复工具
- 教你一步一步实现图标无缝变形切换
- [RaspberryPi]树莓派初体验 - 系统安装与SSH连接
- C#,VB.Net,VB6调用Codesoft打印标签
- 【bzoj2724】【蒲公英】【分块】
- 注释方法解析
- Java多线程线程安全实现方式
- 目前系统系统用到的优化手段
- 【android】:android知道包名打开另一个应用
- Cloudera的CDH和Apache的Hadoop的区别
- Android App 版本更新实现
- tikz-定义坐标点
- 框架 day56 BOS项目练习(activiti,历史数据操作,流程变量,任务组,监听器,网关,spring整合)