【DP】导弹拦截

导弹拦截

导弹拦截

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

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

输入格式
1行，若干个整数（个数≤100000）

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

输入输出样例
输入 #1

389 207 155 300 299 170 158 65

输出 #1

6
2

思路

• 第一问是求最长不上升子序列长度
• 第二问是求最长上升子序列长度（这用到了Dilworth定理）
• 这个定理真的太神奇了，但它具体是什么我不再记录，只是记录一下结论，网上很多讲这个的。
• 这道题算是经典中的经典了吧，因为很多地方都用到了这道题作为例题，我也不讲思路，只是记录一下题解。

代码如下

#include <algorithm>
#include <climits>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <map>
#include <queue>
#include <set>
#include <stack>
#include <string>
#include <unordered_map>
#include <vector>
using namespace std;

vector<int> a, b;//a数组存高度，b数组作为辅助数组
//等等，好像没dp数组，这是因为这篇题解中的dp是隐蔽的
int n;

void init() {
    //初始化数据
    int t = 0;
    while (~scanf("%d", &t)) {
        a.push_back(t);
    }
    n = a.size();
    b.resize(n);
}

void solve() {
    //解决问题
    //求最大拦截数量，即最长不上升子序列长度
    fill(b.begin(), b.end(), INT_MIN);
    for (int i = 0; i < n; i++) {
        *upper_bound(b.begin(), b.end(), a[i], greater<int>()) = a[i];
    }
    int ans = upper_bound(b.begin(), b.end(), INT_MIN + 1, greater<int>()) - b.begin();
    printf("%d\n", ans);
    //最少系统个数，即最长上升子序列长度
    fill(b.begin(), b.end(), INT_MAX);
    for (int i = 0; i < n; i++) {
        *lower_bound(b.begin(), b.end(), a[i]) = a[i];
    }
    ans = lower_bound(b.begin(), b.end(), INT_MAX) - b.begin();
    printf("%d\n", ans);
}

int main(void) {
    init();
    solve();
    return 0;
}

分析一下，这样求最长（上升/不上升/下降/不下降）子序列的时间复杂度是$O(nlogn)$的，应该是最优的算法了吧。
关于求最长（上升/不上升/下降/不下降）子序列的方法，我也有总结，见这篇博文。