P1020 [NOIP1999 提高组] 导弹拦截

题目：

题目描述

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

输入导弹依次飞来的高度，计算这套系统最多能拦截多少导弹，如果要拦截所有导弹最少要配备多少套这种导弹拦截系统。

输入格式

一行，若干个整数，中间由空格隔开。

输出格式

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

输入输出样例

输入 #1复制

389 207 155 300 299 170 158 65

输出 #1复制

6
2

说明/提示

对于前 50% 数据（NOIP 原题数据），满足导弹的个数不超过 104 个。该部分数据总分共 100 分。可使用O(n^2) 做法通过。
对于后 50% 的数据，满足导弹的个数不超过 105 个。该部分数据总分也为 100分。请使用 O(nlogn) 做法通过。

对于全部数据，满足导弹的高度为正整数，且不超过 5×10^4。

此外本题开启 spj，每点两问，按问给分。

第一问：

利用dp方法解决，状态转移方程为：

        

由于需要进行双层循环，时间复杂度为O(n^2)

第二问：用贪心的思想，从左到右变量导弹，维护每个系统当前拦截的最低炮弹数组，当前炮弹的高度由不低于该高度的最低系统进行更新。

O(n^2)代码：

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10;
int a[N];
int dp[N];
int ans=1;
int f[N];
int main(){
	int n;
	int cnt=0;
	int ma=0;
	while(~scanf("%d",&n))
	{
		a[cnt++]=n;
	}
	for(int j=0;j<cnt;j++)
		{for(int k=0;k<j;k++)
		if(a[j]<=a[k])dp[j]=max(dp[j],dp[k]+1);
		ma=max(ma,dp[j]);
		}
	f[ans]=a[0];
	for(int i=1;i<cnt;i++)
	{
		int k1=0;
		int k2=-1;
		for(int j=1;j<=ans;j++)
		if(f[j]>=a[i]){
			if(k1==0){
				k1=f[j];
				k2=j;
			}else if(f[j]<k1){
				k1=f[j];
				k2=j;
			}
		}
		if(k2==-1){
			ans++;
			f[ans]=a[i];
		}else {
			f[k2]=a[i];
		}
	}
	cout<<ma+1<<endl;
	cout<<ans;
	return 0;
}

对于第一问的优化：

dp数组的含义进行更改：dp[i]=j；表示长度为i的上升子序列当前状态下以j结尾；

那么每次更新当前a[j]时，利用二分法找出dp[1]到dp[len]中大于a[j]的最小值，因为dp是降序序列，所以用upper_bound方法，将找出的下标k的dp[k+1]=a[j];

对于第二问的优化：

利用第一问的思路，会发现f数组是一个非递减序列，每次更新都不会改变这种性质。所以每次跟新不再遍历，用二分法查找。

总结：对于优化成O（n*logn）的方法都是使用贪心法（或者有什么其他的方法，我想不到/(ㄒoㄒ)/~~），将查找转换成二分查找。

AC代码如下：

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10;
int a[N];
int dp[N];
int ans=1;
int f[N];
int find_2(int x,int l)
{
	int r=1;
	while(r<=l){
		int mid=(r+l)>>1;
		if(f[mid]>x)l=mid-1;
		else if(f[mid]<x)r=mid+1;
		else return mid;
	}
	return r;
}
int main(){
	int n;
	int cnt=0;
	int ma=0;
	while(~scanf("%d",&n))
	{
		a[cnt++]=n;
	}
	dp[1]=a[0];
	int len=1;
	for(int j=1;j<cnt;j++)
		{
			if(a[j]<=dp[len]){
				len++;
				dp[len]=a[j];
				continue;
			}
			int k=upper_bound(dp+1,dp+len+1,a[j],greater<int>())-dp;
			dp[k]=a[j];
		}
	f[ans]=a[0];
	for(int i=1;i<cnt;i++)
	{
		int k=find_2(a[i],ans);
		if(k<=ans)f[k]=a[i];else f[++ans]=a[i];
	}
	cout<<len<<endl;
	cout<<ans;
	return 0;
}

（我想将lower_bound代替二分查找但是失败了，就直接改成写过的二分查找/(ㄒoㄒ)/~~）