HDU 5919 分块做法

题意和链接前面的blog有,就不再提。。

思路:本来实在是莽不过去。。后来想到其实复杂度均摊下来每块取比sqrt(n)大一些更好。。具体大多少。。嗯。。大概x^2/logx=N  。?我也不会解。。。

所以n取了个2*sqrt(n)结果卡过去了。。。重要的还是思想。。分块后每块排序,和3333套路差不多,时间复杂度O(n*sqrt(n)*log(n))复杂度一看就不小。。不过编码很简单暴力。。毕竟数据结构实在写不了也可以莽一下嘛。。。

 

代码:

 

#include<bits/stdc++.h>

using namespace std;
const int MAXN = 200010;
int n,q,block,num;
int a[MAXN];
int belong[MAXN],L[MAXN],R[MAXN],c[MAXN],d[MAXN],mp[MAXN],tmp[MAXN];

int build(int n){
    block=sqrt(n)*2;
    num=n/block;
    if(n%block!=0) num++;
    for(int i=0;i<num;i++){
        L[i]=i*block+1;
        R[i]=(i+1)*block;
    }
    R[num-1]=n;
    int cnt=0;
    for(int i=1;i<=n;i++){
        if(i<=R[cnt]) belong[i]=cnt;
        else belong[i]=++cnt;
    }
    memset(mp,0,sizeof(mp));
    for(int i = 1; i<= n; i++){
        c[i]=mp[a[i]];
        mp[a[i]] = i;
    }
    for(int i=1;i<=n;i++) d[i]=c[i];
    for(int i=0;i<num;i++){
        sort(c+L[i],c+R[i]+1);
    }
}

int query(int l,int r){
    int ret=0;
    if(belong[l]==belong[r]){
        for(int i=l;i<=r;i++){
            if(d[i]<l) ++ret;
        }
        int x=(ret+1)/2;
        ret=0;
        for(int i=l;i<=r;i++){
            if(d[i]<l) ++ret;
            if(ret==x) return i;
        }
    }
    for(int i=l;i<=R[belong[l]];i++){
        if(d[i]<l) ++ret;
    }
    for(int i=belong[l]+1;i<belong[r];i++){
        tmp[i]=(int)(lower_bound(c+L[i],c+R[i]+1,l)-c)-L[i];
        ret+=tmp[i];
    }
    for(int i=L[belong[r]];i<=r;i++){
        if(d[i]<l) ++ret;
    }
    int x=(ret+1)/2;
    ret=0;
    for(int i=l;i<=R[belong[l]];i++){
        if(d[i]<l) ++ret;
        if(ret==x) return i;
    }
    for(int i=belong[l]+1;i<belong[r];i++){
        int tm=ret;
        ret+=tmp[i];
        if(ret>=x){
            for(int j=L[i];j<=R[i];j++){
                if(d[j]<l) ++tm;
                if(tm==x) return j;
            }
        }
    }
    for(int i=L[belong[r]];i<=r;i++){
        if(d[i]<l) ++ret;
        if(ret==x) return i;
    }
}

int Q[MAXN];

int main(){
    Q[0]=0;
    int t,cas=1;
    scanf("%d",&t);
    while(t--){
        scanf("%d%d",&n,&q);
        for(int i = 1; i <= n; i++) scanf("%d",&a[i]);
        build(n);
        int l,r;
        for(int i=1;i<=q;i++){
            scanf("%d%d",&l,&r);
            l=(l+Q[i-1])%n+1;
            r=(r+Q[i-1])%n+1;
            if(r<l) swap(l,r);
            Q[i]=query(l,r);
        }
        printf("Case #%d:",cas++);
        for(int i=1;i<=q;i++) printf(" %d",Q[i]);
        printf("\n");
    }
    return 0;
}



转载于:https://www.cnblogs.com/zhangxianlong/p/10672532.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值