luogu P1020 导弹拦截

题目描述

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

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

输入格式

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

输出格式

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

输入输出样例

输入
389 207 155 300 299 170 158 65
输出
6
2

这题有两种算法，一种DP算法，时间复杂度O(n^2)，另一种贪心算法，时间复杂度O(nlogn)

思路

第一问就是求最大不上升子序列
第二问 Dilworth定理：偏序集的最少反链划分数等于最长链的长度
我也不太懂，反正就是求最大上升子序列

O(n^2)算法

解释之前写过，这里直接贴代码了
详解点击查看 DP最大上升子序列

#include<iostream>
#include<algorithm>
#include<cstdio>
using namespace std;

int w[100001];
int f[100001];
int d[100001];

int main()
{
	int n=0,ans1=0,ans2=0;
	while(cin>>w[++n]);
	n--;
	for(int i=1;i<=n;i++)
	{
		f[i]=1;
		d[i]=1;
		for(int j=1;j<i;j++)
		{
			if(w[i]<=w[j])
				f[i]=max(f[j]+1,f[i]);
			if(w[i]>w[j])
				d[i]=max(d[j]+1,d[i]);
		}
		ans1=max(ans1,f[i]);
		ans2=max(ans2,d[i]);
	}
	cout<<ans1<<endl<<ans2;
	return 0;
}

O(nlogn)算法

详解点击查看 最长上升子序列优化

#include<iostream>
#include<algorithm>
using namespace std;

int w[100001];
int d[100001];
int u[100001];

int main()
{
	int n=0;
	while(cin>>w[++n]);
	n--;
	int lend=1,lenu=1;
	d[1]=w[1];
	u[1]=w[1];
	for(int i=2;i<=n;i++)
	{
		if(w[i]<=d[lend])
			d[++lend]=w[i];
		else
		{
			int p=upper_bound(d+1,d+lend+1,w[i],greater<int>())-d;
			d[p]=w[i];
		}	
		if(w[i]>u[lenu])
			u[++lenu]=w[i];
		else
		{
			int p=lower_bound(u+1,u+lenu+1,w[i])-u;
			u[p]=w[i];
		}
	}
	cout<<lend<<endl<<lenu;
	return 0;	
}