P1020 [NOIP1999 普及组] 导弹拦截

题目描述

题目描述
 某国为了防御敌国的导弹袭击，发展出一种导弹拦截系统。但是这种导弹拦截系统有一个缺陷：虽然它的第一发炮弹能够到达任意的高度，但是以后每一发炮弹都不能高于前一发的高度。某天，雷达捕捉到敌国的导弹来袭。由于该系统还在试用阶段，所以只有一套系统，因此有可能不能拦截所有的导弹。
 输入导弹依次飞来的高度（雷达给出的高度数据是\le 50000≤50000的正整数），计算这套系统最多能拦截多少导弹，如果要拦截所有导弹最少要配备多少套这种导弹拦截系统。
 题目来源

解题思路

 狄尔沃斯定理(Dilworth’s theorem)亦称偏序集分解定理，是关于偏序集的极大极小的定理，该定理断言：对于任意有限偏序集，其最大反链中元素的数目必等于最小链划分中链的数目。此定理的对偶形式亦真，它断言：对于任意有限偏序集，其最长链中元素的数目必等于其最小反链划分中反链的数目，由偏序集P按如下方式产生的图G称为偏序集的可比图：G的节点集由P的元素组成，而e为G中的边，仅当e的两端点在P中是可比较的，有限全序集的可比图为完全图 。
 其实就是：最少的下降序列个数就等于整个序列最长上升子序列的长度。所以求最少系统数目就是求最少下降序列个数，就是求最长上升子序列的长度。而对于一个系统的最多拦截导弹数，就是求最长非上升子序列（两个问题无关）。

求最长非递增序列

定义数组d1,数组长度len1，求数组a的最长非递增序列(数组都从1开始)
初始状态： d1[1]=a[1]，len1=1。
第i步： 如果a[i]<=d1[len1],则d1[++len1]=a[i].
否则a[i]>=d1[len1],则找到在一个下降序列里寻找一个小于x的数，替换掉它.
结果： 最后i=n(数组a的长度)结束，最长非递增序列的长度就为len1

解释：a[i]<=d1[len1]，可以直接加入序列则直接加入。而对于a[i]>=d1[len1]，直接替换掉第一个小于a[i]的数，无疑换为一个更大的数更有利于递增序列的发展（变长），那么对以后a[i+1]难道不会有影响，两者不会矛盾吗？不会，因为每一次比较都是比较d1数组的最后一个元素，所以a[i]替换为更大的数没有影响，我们要求数组长度，那么这个地方是谁有那么重要吗？反正下一个已经确定，比较也是和下一个甚至之后的比。那解释完不会有负面影响，讲讲为什么要这样做？因为这样会导致这个数变大，以后如果又找到一个a[x]大于d1[len]就又要比较，这个更新为更大的a[i]会影响a[x]的插入，才能调整正确的递增序列。

AC代码

#include<iostream>
#include<algorithm>
using namespace std;
int a[100000];//记录依次飞来的导弹
int d1[100000];//记录最长上升子序列
int d2[100000];//记录最长非上升子序列
int main()
{
	int k = 0;
	while (cin >> a[++k]); k--;//输入导弹高度，存储于数组0到k
	d1[1] = a[1];//初始都不用比较直接纳入
	d2[1] = a[1];
	int len1 = 1, len2 = 1;//分别记录两序列长度
	for (int i = 2; i <= k; i++)//从a[2]开始求最长上升子序列的长度和最长非上升子序列长度
	{
		if (a[i] <= d2[len2])d2[++len2] = a[i];//如果小于最长非上升子序列最小（后）的元素则直接插入
		else//否则更新
		{
			int p2 = upper_bound(d2 + 1, d2 + len2, a[i], greater<int>())-d2;//找到非递增数组d2中第一个（最大的）小于a[i]的下标
			d2[p2] = a[i];//更新大于a[i]的得不偿失，更新等于a[i]的没有意义，所以只能更新第一个小于a[i]的（保证非递增非上升且利于最长）
		}
		if (a[i] > d1[len1])d1[++len1] = a[i];//同理但要上升，所以等于就不能直接插入
		else
		{
			int p1 = lower_bound(d1 + 1, d1 + len1, a[i]) - d1;//上升子序列，更新第一个大于等于（因为可能有等于的情况）a[i]的
			d1[p1] = a[i];//更新小于a[i]的不能保证上升，所以更新第一个大于a[i]的使其更小更有利于长度增长
		}
	}
	cout << len2 << endl << len1;//输出最长非上升序列长度（单系统最多拦截导弹数），输出最长上升子序列长度（最少降序子序列数，系统数目）
	return 0;
}