B1030 完美数列 (25分)/A1085 Perfect Sequence (25分)

给定一个正整数数列,和正整数 p,设这个数列中的最大值是 M,最小值是 m,如果 M≤mp,则称这个数列是完美数列。

现在给定参数 p 和一些正整数,请你从中选择尽可能多的数构成一个完美数列。

输入格式:
输入第一行给出两个正整数 N 和 p,其中 N(≤10^​5)是输入的正整数的个数,p(≤10​e9)是给定的参数。第二行给出 N 个正整数,每个数不超过 10^9
​​ 。

输出格式:
在一行中输出最多可以选择多少个数可以用它们组成一个完美数列。

输入样例:

10 8
2 3 20 4 5 1 6 7 8 9

输出样例:

8

思路:
开一个数组a[],从小到大排序,确定一个左端点a[i]和右端点a[j],使得a[j]<=a[i]*p成立,且j-i最大。
从左至右扫描,对于每一个a[i],在a[i]~a[n-1]内二分查找第一个超过a[i]*p的数的位置j,这样 j-i 就是对a[i]来说满足a[i]*p的最大长度,然后所有 j-i 的最大值

注意:
题目所给数据规模来看,O(n2)肯定会超时(这也是使用二分查找的原因)。同时a[i]*p超出了int的数据范围,所以使用long long数据类型。

样例解释:
10 8
2 3 20 4 5 1 6 7 8 9

在这里插入图片描述
排完序后,初始时刻,i==0,对于a[i] (也就是1),在a[1]~a[9]中寻找第一个大于8的数(也就是9)的位置 j,此时j == 8,j - i == 8

#include <cstdio>
#include <algorithm>
using namespace std;
const int maxn = 100010;
int n, p, a[maxn];

//从a[i+1]~a[n-1]中寻找第一个大于x的数的位置
int binarySearch(int i , long long x)
{
    if(a[n - 1] <= x) return n; //所有数都不大于x,返回n
    int left = i + 1, right = n - 1, mid;
    while(left < right)
    {
        mid = (left + right) / 2;
        if(a[mid] <= x)     
            left = mid + 1;
        else
            right = mid - 1;
    }
    return left;    //由于结束时left==right,因此返回left和right都可以
}

int main()
{
    scanf("%d%d", &n, &p);
    for(int i = 0; i < n; i++)
    {
        scanf("%d", &a[i]);
    }
    sort(a, a + n);
    int ans = 1;    //最大长度,初值为1(表示至少有一个数)
    for(int i = 0; i < n; i++)
    {
        int j = binarySearch(i, (long long)a[i] * p);
        ans = max(ans, j - i);  //更新最大长度
    }
    printf("%d\n",ans);
    return 0;
}

使用upper_bound函数,会更加简洁。

lower_bound( )和upper_bound( )都是利用二分查找的方法在一个排好序的数组中进行查找的。

在从小到大的排序数组中,
lower_bound( begin,end,num):从数组的begin位置到end-1位置二分查找第一个大于或等于num的数字,找到返回该数字的地址,不存在则返回end。通过返回的地址减去起始地址begin,得到找到数字在数组中的下标。

upper_bound( begin,end,num):从数组的begin位置到end-1位置二分查找第一个大于num的数字,找到返回该数字的地址,不存在则返回end。通过返回的地址减去起始地址begin,得到找到数字在数组中的下标。

在这里插入图片描述

#include <cstdio>
#include <algorithm>
using namespace std;
const int maxn = 100010;
int n, p, a[maxn];

int main()
{
    scanf("%d%d", &n, &p);
    for(int i = 0; i < n; i++)
        scanf("%d", &a[i]);
    sort(a, a + n);
    int ans = 1;
    for(int i = 0; i < n; i++)
    {
        
        int j = upper_bound(a + i + 1, a + n, (long long)a[i] * p) - a; //注意upper_bound函数是左闭右开,返回的是地址,两个地址直接相减获取的是两个地址之间的距离
        ans = max(ans, j - i);
    }
    printf("%d", ans);
    return 0;
}

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值