[PAT-A 1085]Perfect Sequence

题目大意：
给定一个正整数数列和正整数p，设这个数列中最大值是M，最小值是m，如果M<m*p,则称这个数列是完美数列
给定参数p和一些正整数，从中选择尽可能多的数构成一个完美数列。求选择出的数的个数.
即从N个正整数中选择若干个数，使得选出的这些数中的最大值不超过最小值的p倍，问满足条件的选择方案中，选出的数的最大个数。

思路：
1）将数列从小到大排序。
2）从左至右扫描序列，对其中的每一个数a[i]，在a[i+1]~a[n-1]内二分查找第一个超过a[i]*p的数，位置为j，则j-i记为完美数列的个数，
取所有j-i的最大值即为所求答案

AC代码1：（暴力解，两次循环）

#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn = 100010;
int main() {
	int n, p, a[maxn];
	scanf("%d%d", &n, &p);
	for (int i = 0; i < n; i++) {
		scanf("%d", &a[i]);
	}
	sort(a, a + n);
	int i = 0, j = 0, count = 1;
	while (i < n && j < n) {
		while (j < n && a[j] <= (long long)a[i] * p) {
			count = max(count, j - i + 1);
			j++;
		}
		i++;
	}
	printf("%d\n", count);
	return 0;
}

AC代码2（二分查找）：

#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn = 100010;
int n, p, a[maxn];
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;//说明第一个大于x的数在mid后面
		else right = mid;//说明第一个大于x的数在mid之前(含mid)
	}
	return left;
}
int main() {
	(void)scanf("%d %d", &n, &p);
	for (int i = 0; i < n; i++) {
		(void)scanf("%d", &a[i]);
	}
	sort(a, a + n);
	int ans = 1;
	for (int i = 0; i < n; i++) {
		int j = binarySearch(i, (long long)a[i] * p);
		ans = max(ans, j - i);
	}
	printf("%d", ans);
	return 0;
}

AC代码3(upper_bound()函数)：
（其实代码2就是自己手写upper_bound()函数）

#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn = 100010;
int n, p, a[maxn];
int main() {
	(void)scanf("%d %d", &n, &p);
	for (int i = 0; i < n; i++) {
		(void)scanf("%d", &a[i]);
	}
	sort(a, a + n);
	int ans = 1;
	for (int i = 0; i < n; i++) {
		//upper_bound()返回从begin到end-1位置第一个大于num的数，返回地址，减去数组首地址，得到对应的下标
		//lower_bound()返回从begin到end-1位置第一个大于等于num的数，返回地址，减去数组首地址，得到对应的下标
		int j = upper_bound(a + i + 1, a + n, (long long)a[i] * p) - a;
		ans = max(ans, j - i);
	}
	printf("%d", ans);
	return 0;
}

AC代码4:two pointers

#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn = 100010;
int main() {
	int n, p, a[maxn];
	(void)scanf("%d %d", &n, &p);
	for (int i = 0; i < n; i++) {
		(void)scanf("%d", &a[i]);
	}
	sort(a, a + n);
	int i = 0, j = 0, count = 1;
	while (i < n && j < n) {
		while (j < n && a[j] <= (long long)a[i] * p) {
			count = max(count, j - i + 1);
			j++;
		}
		i++;
	}
	printf("%d\n", count);
	return 0;
}