查询区间内等于x的数的个数(分块)

问题:有一个有n个数的数组,q次查询,每次查询l, r ,x,查询区间[l,r]内等于x的数的个数

思路:分块。


就把这题当成是分块的入门题来讲解一下分块。分块其实就是一种比较优美的暴力(我觉得),一般的分块都是把长度为n的数组分成每一块为sqrt(n)个数的多个块。然后对于区间的操作就可以不再是一个一个数进行处理,而是可以一个块一个块的处理,每次查询最多会涉及到sqrt(n)个块,这样复杂度就降了下来。

(先进行几点说明,本文中的数组下标从1开始,块的编号也是从1开始,每个块的大小k = sqrt(n))

分块过后的对于区间[l,r]的操作就分为两种情况:

设cl为区间左端点l所在的块号,cr为区间右端点r所在的块号。那么怎么求一个下标所属的块号呢,很简单:(x-1)/k + 1。

第一种情况是cl == cr,也就是左右端点在同一个块内,这样直接从l循环到r处理操作就可以了,复杂度很显然最多是sqrt(n)

第二种情况我们可以统一的假定为 cl != cr,这样我们可以把区间分为三部分,一是l到cl*k这部分可以直接循环操作,后面就是第cl+1块到cr-1块这部分,最后就是(cr-1)*k+1到r这部分。


这一题要求的是一个区间内等于x的有多少个数,我们可以对每一个块进行排序,然后就可以用二分求出等于x的个数。


#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>

using namespace std;

const int maxn = 300010;

int a[maxn];
int b[maxn];

int n,k;
int li(int x) //求x的块号
{
    return (x-1)/k + 1;
}

int query(int l,int r,int x)
{
    int i;
    int cl = li(l);
    int cr = li(r);
    int cnt = 0;
    if(cl == cr)
    {
        for(i=l;i<=r;i++)
            if(a[i] == x)
                cnt++;
        return cnt;
    }
    for(i=l;i<=cl*k;i++)
    {
        if(a[i] == x)
            cnt++;
    }
    for(i=cl+1;i<cr;i++)
    {
        cnt += upper_bound(b+(i-1)*k+1,b+i*k+1,x) - lower_bound(b+(i-1)*k+1,b+i*k+1,x);
    }
    for(i=(cr-1)*k+1;i<=r;i++)
    {
        if(a[i] == x)
            cnt++;
    }

    return cnt;
}
int main(void)
{
    int i,j;
    scanf("%d",&n);
    for(i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        b[i] = a[i];
    }
    k = sqrt(n);
    for(i=1;i<=k+1;i++)
    {
        if((i-1)*k+1 > n)
            break;
        sort(b+(i-1)*k+1,b+min(i*k,n)+1);
    }
    int q;
    scanf("%d",&q);
    while(q--)
    {
        int l,r,x;
        scanf("%d%d%d",&l,&r,&x);
        int ans = query(l,r,x);
        printf("%d\n",ans);
    }

    return 0;
}

  • 3
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
powell方法是计算最值的有效方法,一般情况下是无约束的,目前Powell已发展了一种称为BOBYQA的带区间约束算法。 压缩文件里面optimization.h是无约束powell寻优的一个类,以前在baidu上找到的。powell.h和powell.cpp是带区间约束的powell算法BOBYQA的C++包装,由于BOBYQA算法只能计算二维以上,一维算法是俺自己写的。调用时,一维算法采用ExecuteBrent,二维以上采用BoundedPowell。 BOBYQA目录下是powell方法祖师爷M. J. D. Powell写的FORTRAN代码,powell.lib是俺把Fortran代码编译成的静态库。由于这个lib是Fortran Power Station 4.0做的,msfrt40.dll就成了一个甩不掉的累赘。 如果使用Fortran编程的话,就不需要俺的包装了,直接用Fortran代码吧。 使用方法请参考Fortran代码中main.f,有问题请联系[email protected]。 另外链接时好像还需要fps4.0中的msfrt.lib,原来没有加进去,现在加不进去了。需要的话请邮件联系。 鉴于很多同志不清楚用法,下面略加说明。 BoundedPowell函的参说明。 第一个参,待优化函的指针,这个函必须定义成_stdcall类型,其本身带有三个参,第一个参是优化变量个数,第二个参是优化变量组,第三个参就是函的值,由于使用Fortran的关系,这几个参都必须传递地址。 举例说,假定待优化函为f(x)=x^2,那么只有一个参,这个函应该定义为 void _stdcall objfun(int *n, double *para, double *f) { *f=para[0]*para[0]; } 第二个参int n,待优化计算变量的个数 第三个参double *x, 待优化变量,这是一个数组,长度为n 第四个参double *xlb,变量的下界,组 第五个参double *xub,变量的上界,组 第五个参double rhobeg,第六个参double rhoend,这两个参是Powell同志定义的两个半径,具体我也说不清楚,你想搞清楚的话可以搜Powerll同志的文献。一般你把rhobeg设置为1,rhoend设置为计算的精度,比如你希望精度为万分之一,rhoend就为1e-4。 第七个参int maxfun,最大迭代次 第八个参,BOBYQA计算时的返回代码,具体如下 // //BOBYQA iflag返回值的含义: //IFLAG=1, Return from BOBYQA because NPT is not in the required interval //IFLAG=2, Return from BOBYQA because one of the differences XU(I)-XL(I)/6X is less than 2*RHOBEG. //IFLAG=3, Return from BOBYQA because FCN has been called MAXFUN times. //IFLAG=4, Return from BOBYQA because of much cancellation in a denominator. //IFLAG=5, Return from BOBYQA because a trust region step has failed to reduce Q. //

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值