bzoj1878: [SDOI2009]HH的项链

Description
HH有一串由各种漂亮的贝壳组成的项链。HH相信不同的贝壳会带来好运,所以每次散步 完后,他都会随意取出一段贝壳,思考它们所表达的含义。HH不断地收集新的贝壳,因此他的项链变得越来越长。有一天,他突然提出了一个问题:某一段贝壳中,包含了多少种不同的贝壳?这个问题很难回答。。。因为项链实在是太长了。于是,他只好求助睿智的你,来解决这个问题。

Input
第一行:一个整数N,表示项链的长度。
第二行:N个整数,表示依次表示项链中贝壳的编号(编号为0到1000000之间的整数)。
第三行:一个整数M,表示HH询问的个数。
接下来M行:每行两个整数,L和R(1 ≤ L ≤ R ≤ N),表示询问的区间。
N ≤ 50000,M ≤ 200000

Output
M行,每行一个整数,依次表示询问对应的答案。

Sample Input
6
1 2 3 4 3 5
3
1 2
3 5
2 6

Sample Output
2
2
4

学习莫队算法之后手写的第一道题,,感觉莫队其实不是很难,,,思路非常清晰暴力,只是时间复杂度计算的很玄学。
这题大致题意就是在给你M个询问,让你计算L~R之间有多少个不相同的数,
我记得某位大佬说过,只要是可以离线的区间查询问题,用莫队算法都会产生奇效,
那么这道题,有多组区间查询,还不能用线段树ST表什么乱七八糟的东西,,那我们就用莫队分块来搞一搞
大体思路就是对查询进行分块,以左端点所在的块为关键字排序,左端点在同一块上的查询就按照右端点进行排序,
那么我们固定一个当前区间,枚举排好序的查询中每组查询的左右端点,在移动的过程中便更新完了答案
看起来非常暴力是不是,,,莫队算法最核心的部分就是最开始对查询进行的排序,
我们把数组分为k块;根据我们的排序可知,我们的左端点最多移动(n*q)/k,右端点最多移动 n*k次,
所以说当我们把块的大小设为sqrt(n)时能保证时间复杂度稳定为n*sqrt(n),
其他就没什么好说的了,,看代码吧

#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
struct jgt{int l,r,k,ans;}qst[1000005];
int len,n,q,num[1000005],cot[1000005],bel[1000005];
int ll,rr,ansl;
bool cmp1(jgt a,jgt b){return bel[a.l]<bel[b.l]||(bel[a.l]==bel[b.l]&&a.r<b.r);}//以左端点所在的块为第一关键字,右端点的位置为第二关键字排序, 
bool cmp2(jgt a,jgt b){return a.k<b.k;}
void del(int a)//从当前区间删除元素 
{
    cot[a]--;
    if(cot[a]==0)ansl--;
}
void add(int a)//向当前区间添加元素 
{
    if(cot[a]==0)ansl++;
    cot[a]++;
}
int main()
{
    scanf("%d",&n);
    len=ceil(sqrt(1.0*n));//确定区间大小

    for(int i=1;i<=n;i++)
    {
        scanf("%d",&num[i]);
        bel[i]=(i-1)/len+1;//确定所属的块 
    }
    scanf("%d",&q);
    for(int i=1;i<=q;i++)//输入查询 
    {
        qst[i].k=i;
        scanf("%d%d",&qst[i].l,&qst[i].r);
    }
    sort(qst+1,qst+q+1,cmp1);//排序

    ll=qst[1].l;rr=ll;add(num[ll]);//确定当前区间 

    for(int i=1;i<=q;i++)
    {
        while(ll>qst[i].l) {ll--;add(num[ll]);}//移动左端点 
        while(ll<qst[i].l) {del(num[ll]);ll++;}

        while(rr>qst[i].r) {del(num[rr]);rr--;}//移动右端点 
        while(rr<qst[i].r) {rr++;add(num[rr]);}

        qst[i].ans=ansl;//储存答案 
    }

    sort(qst+1,qst+q+1,cmp2);

    for(int i=1;i<=q;i++)   printf("%d\n",qst[i].ans);

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值