PAT乙级 1030 完美数列 (25分)

1030 完美数列 (25分)

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

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

输入格式：

输入第一行给出两个正整数 N 和 p，其中 N（≤10​⁵​​ ）是输入的正整数的个数，p（≤10^​9​​ ）是给定的参数。第二行给出 N 个正整数，每个数不超过 10^​9​​ 。

输出格式：

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

输入样例：

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

输出样例：

8

思路：

首先一个坑点是，这个题目涉及两个10⁹级别大小的整型相乘，所以用int的话可能会有溢出的风险。
可不可以写成 M/p<=m 呢？不可以，因为 9/8<=1 结果是True(因为整型数字相除，余数被省略了)


在《算法笔记》（胡凡，曾磊著）中提到过一种“算法”，叫做 two pointers 算法。书上说 two pointers 更像是一种编程技巧，而不像算法，所以很少有教材会写到。two pointers 相关知识，可以百度搜索一下。
对本题而言，具体做法如下：


1.对所有数据从小到大排序。
2.设置p1 p2两个变量，这两个变量都是从0开始，并且假定p2 大于p1。就像一个人追着一个人赶路一样。
  2.1 如果 nums[p2] 小于 nums[p1]*p，这时候p2就说：“p1你不要再动了，你已经够大了，我要往前走一步！”，然后p2加一。
  2.2 如果 nums[p2] 大于 nums[p1]*p，这时候p2就说：“p1你太小了，我站着不动，你往前走一步吧！”，然后p1加一。
2.3更新一下最大值
3.输出最大值


可以看出在这个循环中，不是p1增加了1就是p2增加了1，而这两个变量都不能超过n，所以这个算法的时间复杂度在最坏情况下是2n.


关于下面的统计数据，出现了多个峰值，是非常有趣的，其实主要分布在两个区间内，一个是30毫秒左右的，一个是75毫秒左右的，75ms主要是因为cin cout相对比较慢一些，左右的小峰值和判定系统的状态有关，比如我同一段代码多次提交有时候能20毫秒有时候30毫秒。我浏览了一下其他同学写的博客，还有一种方法是双层for循环，效率也很高而且也好理解。O(n²)的方法肯定是要超时的。

统计数据：

题解 方法一：

#include <cstdio>
#include <algorithm>
using namespace std;

long long n,p,nums[110000],maxx=-1,p1=0,p2=0;

int main(){
    scanf("%lld %lld",&n,&p);
    for(int i=0;i<n;scanf("%lld",&nums[i++]));
    sort(nums,nums+n);
    while(p1<n && p2<n){
        if(nums[p2]<=nums[p1]*p)
            maxx = maxx>(p2-p1+1)?maxx:p2-p1+1; //更新最大值
        nums[p2]<=nums[p1]*p ? p2++ : p1++;
    }
    printf("%lld",maxx);
    return 0;
}

题解 方法二：

#include <cstdio>
#include <algorithm>
using namespace std;

long long n,p,nums[110000],maxx=-1,p1=0,p2=0;

int main(){
    scanf("%lld %lld",&n,&p);
    for(int i=0;i<n;scanf("%lld",&nums[i++]));
    sort(nums,nums+n);
    while(p2<n && p1<n){
        maxx = maxx>(p2-p1+1)?maxx:p2-p1+1;
        nums[p2]<=nums[p1]*p ? p2++ : 0;
        nums[p2]<=nums[p1]*p ? 0 : p1++;
    }
    printf("%lld",maxx);
    return 0;
}

题解 方法三

双for循环结构。
-_-||| 我怎么觉得这个方法更好一些呢。。。

#include <cstdio>
#include <algorithm>
using namespace std;

long long n,nums[110000],p,res;

int main(){
    scanf("%lld %lld",&n,&p);
    for(int i=0;i<n;scanf("%lld",&nums[i++]));
    sort(nums,nums+n);
    for(int i=0;i<n;++i)
        for(int j=i+res;nums[j] <= nums[i]*p && j<n;++j)
            res = j-i+1;
    printf("%lld",res);
    return 0;
}

题解 方法四：

二分搜索，从时间复杂度上看，肯定是不如前面几种方法，但是也是一种思路，而且速度还可以。

#include <cstdio>
#include <algorithm>
using namespace std;

long long n,M,m,p,a[100005],ma=1,ta;

int main(){
    scanf("%lld %lld",&n,&p);
    for(int i=0;i<n;scanf("%lld",&a[i++]));
    sort(a,a+n);
    for(int i=0;i+ma<n;++i){
        ta = upper_bound(a+i,a+n,a[i]*p)-a-i; //下限是a+i，算是一个优化
        ma = max(ma,ta);                      //用upper_bound而不是lower_bound
    }
    printf("%lld",ma);
    return 0;
}