[pta随笔]patB1030 完美数列——双指针算法

原题：

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

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

输入格式：

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

输出格式：

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

输入样例：

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

输出样例：

8

另外：本题卡时间，200ms。 

关于测试点：

1、2两个主测试点应该和sample大差不差

3测试点是单一数据，不做特判会输出0

4测试点是海量数据，如果是普通的双循环会炸

5测试点是大数据，不用long long会炸

 

 这题目卡了好久，主要卡了大家都在卡的测试点4，然而最后去csdn看了大佬的代码想来想去才意识到原来是我一开始的算法就有漏洞，不过没想到的是居然除了带海量数据的测试点4过不去以外其他的都能秒，还不需要开特判：

//测试点4不能过，算法本身有漏洞，仅供学习参考
#include<bits/stdc++.h>
using namespace std;
int main()
{
	deque<long long>li;
	long long save=0,x;
	long long n, p;
	cin >> n >> p;
	for (long long i = 0; i < n; i++)
	{
		cin >> x;
		li.emplace_front(x);
	}
	sort(li.begin(), li.end());
	x = li.size();
	for (long long i = 0; i < x; i++)//头处理
	{
		if (li.back() <= li[i] * p) 
		{
			if (li.size() - i > save) save = li.size() - i;
			break;
		}
	}
	for (long long i = x - 1; i >= 0; i--)//尾处理
	{
		if (li[i] <= li.front() * p)
		{
			if (i+1 > save) save = i+1;
			break;
		}
	}
	cout << save << endl;
	return 0;
}

这个算法本身就是排序完毕之后，从头砍到不能砍，再从尾砍到不能砍，时间复杂度很低，乍一看很好测，然而：

input:
5 8
0 1 2 16 10000
result:
1
true answer:
2

 也就是说，如果最佳数据在中间，而不在两边，就无法得到正确答案。

那么这么方便的混分算法不能用的话，我考虑的就是两种算法：

一种，排序完开双重循环，i在头，j在尾，双方各自遍历一遍，如果是最大完美数列，就保存下来。时间复杂度是O(n^2)，如果要削减就要剪枝，但是这玩意剪枝挺麻烦的，即使ij相遇终止依然会超时。(当然可能是因为我没有在多花点时间处理而是直接选了第二种吧)

另一种，双指针，ij放在开头，j先走，直到不再是完美数列后让i一步一步走，如此往复直到ij都走到末尾，感觉怪怪的，但是时间复杂度是O(n)-O(n^2)，而且几乎不会到n^2那么大。本题的测试点4就只要42ms。

其实看题解的时候就发现那些算dalao的代码其实本质就是双指针，不过却没明说这一点...双指针算法确实是一个很经典很有效的好东西，不过我也就是个在acwing看了点基础课的萌新罢了...

代码如下：

#include<bits/stdc++.h>
using namespace std;
int main()
{
	deque<long long>li;
	long long save = 0, x;
	long long n, p;
	cin >> n >> p;
	for (long long i = 0; i < n; i++)
	{
		cin >> x;
		li.emplace_front(x);
	}
	sort(li.begin(), li.end());
	x = li.size();
	for (long long i = 0, j = 0; i < x; i++)
	{
		while (j < x&&li[j] <= li[i] * p ) //记得注意越界问题
        {
			j++;
		}
		if (j - i > save) save = j - i;
	}
	if (x == 1) save = 1;
	cout << save << endl;
	return 0;
}