【DP】拦截导弹（C++）

题目描述

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

输入格式

一行，空格隔开的若干个正整数。

输出格式

两行，每行一个整数。 第一个整数表示最多能拦截的导弹数; 第二个整数表示要拦截所有导弹最少要配备的系统数。

数据范围与提示

2<=n<=500

输入样例

389 207 155 300 299 170 158 65

输出样例

6
2

这是一道经典的动态规划题，分成两个部分，一个部分是求一套最多能拦截的导弹数，另一个部分是求系统数，自然，我们分成两部分来求。

第一部分：当$a_j>=a_i$时，状态转移方程为$f_i=max(f_i,f_j+1)$
第二部分：死循环暴力枚举，从第一个开始扫，扫到这一套系统无法拦截为止，接着再增加一套，如果全部没了，则代表扫完了，那么当前的系统数就是最优的。


参考代码_{（本代码在NFLSOJ提交通过）}：

#include <bits/stdc++.h>
using namespace std;
int a[10001],n=1,f[10001];
bool ok[10001];
int main() {
	while(scanf("%d",&a[n])!=EOF)
		n++;
	n--;
	for(int i=1; i<=n; i++) {
		for(int j=i; j>=1; j--) {
			if(a[i]<=a[j]) {
				f[i]=max(f[i],f[j]+1);
			}
		}
	}
	int mx=0;
	for(int i=1; i<=n; i++) {
		//cout<<"### "<<f[i]<<"\n";
		mx=max(mx,f[i]);
	}
	printf("%d\n",mx);
	int mn=INT_MAX,ans=0,sum=0;
	while(1) {
		sum++;
		mn=INT_MAX;
		for(int i=1; i<=n; i++) {
			if(!ok[i]&&a[i]<=mn) {
				mn=min(a[i],mn);
				ok[i]=1;
				ans++;
				//cout<<"### "<<a[i]<<"   ans="<<ans<<"   sum="<<sum<<endl;
			}
		}
		if(ans==n) {
			printf("%d\n",sum);
			return 0;
		}
	}

	return 0;
}