首页 > 代码库 > 【莫队算法】【权值分块】bzoj2223 [Coci 2009]PATULJCI
【莫队算法】【权值分块】bzoj2223 [Coci 2009]PATULJCI
不带修改主席树裸题<=>莫队+权值分块裸题。
复杂度O(m*sqrt(n))。
P.S.题目描述坑爹,第二个数是权值的范围。
#include<cstdio>#include<algorithm>#include<cmath>using namespace std;#define N 300001#define M 10001int f,c;inline void R(int &x){ c=0;f=1; for(;c<‘0‘||c>‘9‘;c=getchar())if(c==‘-‘)f=-1; for(x=0;c>=‘0‘&&c<=‘9‘;c=getchar())(x*=10)+=(c-‘0‘); x*=f;}void P(int x){ if(x<10)putchar(x+‘0‘); else{P(x/10);putchar(x%10+‘0‘);}}int sum=1,Lim,l[105],n,m,num[10001],r[105],cnts[105],a[N],b[10001],anss[M],num2[N];struct Ask{int l,r,p;}Q[M];bool operator < (const Ask &a,const Ask &b){return num2[a.l]!=num2[b.l] ? num2[a.l]<num2[b.l] : a.r<b.r;}void val_mb(){ int sz=sqrt(Lim); if(!sz) sz=1; for(;sum*sz<Lim;++sum) { l[sum]=r[sum-1]+1; r[sum]=sum*sz; for(int i=l[sum];i<=r[sum];++i) num[i]=sum; } l[sum]=r[sum-1]+1; r[sum]=Lim; for(int i=l[sum];i<=Lim;++i) num[i]=sum;}void mo_mb(){ int tot=1,sz=sqrt(n); if(!sz) sz=1; for(;tot*sz<n;++tot) { int r=tot*sz; for(int i=(tot-1)*sz+1;i<=r;++i) num2[i]=tot; } for(int i=(tot-1)*sz+1;i<=n;++i) num2[i]=tot;}void Insert(const int &x){++b[x]; ++cnts[num[x]];}void Delete(const int &x){--b[x]; --cnts[num[x]];}int Query(const int &L,const int &R){ int goal=(R-L+1>>1); for(int j=1;j<=sum;++j) if(cnts[j]>goal) for(int i=l[j];i<=r[j];++i) if(b[i]>goal) return i; return -1;}int main(){ R(n); R(Lim); for(int i=1;i<=n;++i) R(a[i]); val_mb(); R(m); for(int i=1;i<=m;++i) R(Q[i].l),R(Q[i].r),Q[i].p=i; mo_mb(); sort(Q+1,Q+m+1); for(int i=Q[1].l;i<=Q[1].r;++i) Insert(a[i]); anss[Q[1].p]=Query(Q[1].l,Q[1].r); for(int i=2;i<=m;++i) { if(Q[i].l<Q[i-1].l) for(int j=Q[i-1].l-1;j>=Q[i].l;--j) Insert(a[j]); else for(int j=Q[i-1].l;j<Q[i].l;++j) Delete(a[j]); if(Q[i].r<Q[i-1].r) for(int j=Q[i-1].r;j>Q[i].r;--j) Delete(a[j]); else for(int j=Q[i-1].r+1;j<=Q[i].r;++j) Insert(a[j]); anss[Q[i].p]=Query(Q[i].l,Q[i].r); } for(int i=1;i<=m;++i) if(anss[i]==-1) puts("no"); else {printf("yes "); P(anss[i]); puts("");} return 0;}
【莫队算法】【权值分块】bzoj2223 [Coci 2009]PATULJCI
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。