GDUT - 专题学习2 A - 导弹

A - 导弹

题目

某国为了防御敌国的导弹袭击，发展出一种导弹拦截系统。但是这种导弹拦截系统有一个缺陷：虽然它的第一发炮弹能够到达任意的高度，但是以后每一发炮弹都不能高于前一发的高度。某天，雷达捕捉到敌国的导弹来袭。由于该系统还在试用阶段，所以只有一套系统，因此有可能不能拦截所有的导弹。

输入导弹依次飞来的高度（雷达给出的高度数据是不大于30000的正整数），计算这套系统最多能拦截多少导弹，如果要拦截所有导弹最少要配备多少套这种导弹拦截系统。

输入格式

1行，若干个整数（个数≤100000）

输出格式

2行，每行一个整数，第一个数字表示这套系统最多能拦截多少导弹，第二个数字表示如果要拦截所有导弹最少要配备多少套这种导弹拦截系统。

输入输出样例

输入

389 207 155 300 299 170 158 65

输出

6
2

思路

首先要解决的是输入问题，这里的输入并没有明确终止输入的条件，所以直接用

while (cin >> a[++n])

就可以了

其实题目要求的最多拦截数和最少拦截系统意思就是要求序列的最长下降子序列的长度和最长上升序列的长度，简单易懂。

就以样例数据为例

389 207 155 300 299 170 158 65

其中最长下降子序列应该为

389 300 299 170 158 65

共6个元素

最长上升子序列其一应该为

155 300

共2个元素

那么如何求得此结果呢？

首先需要一个用来储存所有导弹高度的数组a，然后还需要一个储存下降子序列的数组b，且b的长度为len

如果a[i] <= b[len]，就说明a[i]元素是可以接在b后面的，且接完后b仍然有序。

b[++len] = a[i]

如果a[i] > b[len]，就说明a[i]接不上，这时候我们就硬接，如何硬接？用二分法找到数组b中第一个小于a[i]的元素，然后用a[i]代替之，因为len已经储存了上一个下降子序列的长度，所以就相当于储存了从a[1]到a[i]的最长下降子序列的长度了，这段话请细品。

此处因为是严格大于的，所以用upper_bound而不是lower_bound，但是这里有个非常重要的点，就是二分查找函数默认只能查找有序序列，而这里的有序指的是升序，所以我们这个时候就要加一个比较函数greater <int> ()，否则就会出现奇怪的错误。

b[upper_bound(b + 1, b + len + 1, a[i], greater<int>()) - b] = a[i];

最后只需输出长度len即可。

反之亦然，只需将比较符号反过来，再将upper_bound改成lower_bound即可。

代码

#include<iostream>
#include<algorithm>
using namespace std;
int a[100005], ans1[100005], t1 = 1, ans2[100005], t2 = 1, n = 0;//t1是降序子序列长度,t2反之
int main()
{
	while (cin >> a[++n]);//输入
	ans1[1] = ans2[1] = a[1];//先输入一个方便比较
	for (int i = 2; i <= n - 1; ++i)
	{
        //求最大下降序列长度
		if (ans1[t1] >= a[i])
		{
			ans1[++t1] = a[i];
		}
		else
		{
			ans1[upper_bound(ans1 + 1, ans1 + t1 + 1, a[i], greater<int>()) - ans1] = a[i];
		}
        //求最大上升序列长度
		if (ans2[t2] < a[i])
		{
			ans2[++t2] = a[i];
		}
		else
		{
			ans2[lower_bound(ans2 + 1, ans2 + t2 + 1, a[i]) - ans2] = a[i];
		}
	}
	cout << t1 << endl << t2;
	return 0;
}