区间第k小

题目大意

给你一个序列元素均在[0,n)内,并给定常数w。
每次在线询问区间第k小,要求忽略区间出现次数>w的数。

根号算法

我们不妨考虑如何离线。
可以使用莫队+线段树做到 n根号n log n,非常菜鸡。
可以发现我们只有n次询问却有n根号n次插入删除,如果我们使用线段树是很亏的。
考虑经典套路即平衡复杂度:
kth是有办法O(1)-O(√n)的,具体做法是维护值域分块,以及一个桶。
值域分块就是把值域做分块,维护每一块有多少个数。
那么可以用√n的时间确定答案在哪一块,然后暴力扫这个块,通过桶定位答案。
现在它要求我们在线,我们思考我们为何要使用莫队,因为我们想要得到一段区间的值域分块和桶。
如何在线得到一段区间的值域分块?
我们可以把序列分块,并预处理第i块到第j块的值域分块,这个空间和时间都是n根号的。
得到一段区间的值域分块,先利用预处理的,再将剩余部分暴力加入,这个非常可行。
麻烦的是桶不能这样做,但是我们并不需要知道一个区间的桶,具体的我们只需要知道一个数在一个区间的出现次数。
这个可以预处理每个数在前i块中出现了多少次,然后只需要对剩余部分暴力维护出桶。
于是这题就做完了。

#pragma GCC optimize(2)
#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define min(a,b) (a<b?a:b)
#define max(a,b) (a>b?a:b)
using namespace std;
const int maxn=100000+10,B=320;
int belong[maxn],ans[maxn/B+10][maxn/B+10][maxn/B+10],sum[maxn/B+10][maxn],a[maxn],cnt[maxn],b[maxn/B+10],sta[80];
int i,j,k,l,r,t,n,m,x,y,tot,top,w,now,wdc;
bool czy;
int read(){
    int x=0,f=1;
    char ch=getchar();
    while (ch<'0'||ch>'9'){
        if (ch=='-') f=-1;
        ch=getchar();
    }
    while (ch>='0'&&ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}
inline int query(int l,int r,int k,int f){
    int i,j,t;
    int ll=belong[l],rr=belong[r];
    fo(i,1,belong[n])
        if (k>b[i]) k-=b[i];
        else break;
    j=i;
    fo(i,(j-1)*B+1,min(j*B,n)+1){
        if (i>min(j*B,n)) break;
        if (f) t=sum[rr-1][i]-sum[ll][i]+cnt[i];
        else t=cnt[i];
        if (t>w) t=0;
        if (k>t) k-=t;else break;
    }
    if (i>min(j*B,n)) return n+1;
    return i;
}
void write(int x){
    if (!x){
        putchar('0');
        putchar('\n');
        return;
    }
    top=0;
    while (x){
        sta[++top]=x%10;
        x/=10;
    }
    while (top) putchar('0'+sta[top--]);
    putchar('\n');
}
int main(){
    freopen("kth.in","r",stdin);freopen("kth.out","w",stdout);
    n=read();w=read();m=read();czy=read();
    fo(i,1,n) a[i]=read(),a[i]++;
    fo(i,1,n) belong[i]=(i-1)/B+1;
    fo(i,1,belong[n]){
        fo(j,1,n) cnt[j]=0;
        fo(j,1,belong[n]) b[j]=0;
        fo(j,(i-1)*B+1,n){
            x=a[j];
            if (cnt[x]==w) b[belong[x]]-=w;
            else if (cnt[x]<w) b[belong[x]]++;
            cnt[x]++;
            if (j==n||belong[j]!=belong[j+1])
                fo(k,1,belong[n]) ans[i][belong[j]][k]=b[k];
        }
    }
    fo(i,1,n) sum[belong[i]][a[i]]++;
    fo(x,1,n)
        fo(i,1,belong[n]) sum[i][x]+=sum[i-1][x];
    fo(i,1,n) cnt[i]=0;
    while (m--){
        wdc++;
        fo(i,1,belong[n]) b[i]=0;
        j=read();k=read();y=read();
        if (czy){
            j^=now;
            k^=now;
            y^=now;
        }
        l=belong[j];r=belong[k];
        if (r-l<=0){
            fo(i,j,k){
                x=a[i];
                if (cnt[x]==w) b[belong[x]]-=w;
                else if (cnt[x]<w) b[belong[x]]++;
                cnt[x]++;
            }
            now=query(j,k,y,0)-1;
            write(now);
            fo(i,j,k) if (cnt[a[i]]) cnt[a[i]]=0;
            continue;
        }
        fo(i,1,belong[n]) b[i]=ans[l+1][r-1][i];
        fo(i,j,l*B){
            x=a[i];
            t=cnt[x]+sum[r-1][x]-sum[l][x];
            if (t==w) b[belong[x]]-=w;
            else if (t<w) b[belong[x]]++;
            cnt[x]++;
        }
        fo(i,(r-1)*B+1,k){
            x=a[i];
            t=cnt[x]+sum[r-1][x]-sum[l][x];
            if (t==w) b[belong[x]]-=w;
            else if (t<w) b[belong[x]]++;
            cnt[x]++;
        }
        now=query(j,k,y,1)-1;
        write(now);
        fo(i,j,l*B) if (cnt[a[i]]) cnt[a[i]]=0;
        fo(i,(r-1)*B+1,k) if (cnt[a[i]]) cnt[a[i]]=0;
    }
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值