NOIP1999 提高组] 导弹拦截 题解

题目描述

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

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

输入格式

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

输出格式

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

输入输出样例

输入 #1

389 207 155 300 299 170 158 65

输出 #1

6
2

题目大意

求一个数列的最长不上升子序列的长度和最优经过几次最长不上升子序列的覆盖可以覆盖完整个数列。

思路分析

O(n2)O(n2) 做法这里不提供。仅介绍 O(nlog⁡n)O(nlogn) 做法。

整体思想：贪心+二分

最多能拦截多少很简单。我们用数组维护一个栈。这个栈具有单调性。

然后将原数组从头到尾扫一遍。每次扫到一个数，我们对其进行以下操作：

• 如果这个数 ≤≤ 栈顶：

  将这个数直接入栈。

• 否则：

  找到栈里最小的大于该数的元素，将这个元素替换成这个数。因为栈的单调性，使用二分查找即可。

这样，一共扫 nn 次，每次的最坏时间复杂度 O(log⁡n)O(logn)，总时间复杂度 O(nlog⁡n)O(nlogn)。

那么，如何求要覆盖多少次呢？

先来看一个定理：

Dilworth定理

对于一个偏序集，最少链划分等于最长反链长度。是不是很懵？在这道题中，我可以用人话复述一遍：

最长上升子序列的长度就是能构成的不上升序列的个数。

这样，我们经过简易的推理，就可以得出：

要使用导弹的次数就是最长上升子序列的长度。

如何求最长上升子序列呢？我们把求最长不上升子序列的代码稍微改动一下就好了。

这样，我们就可以在 O(nlog⁡n)O(nlogn) 的时间复杂度内通过本题。

完整代码

#include <bits/stdc++.h>

using namespace std;

const int maxn = 1e6 + 1;

int n;
int a[maxn];
int cnt;
int s[maxn];

void input() {
    do {
        n++;
    } while (cin >> a[n]);
    n--;
}

void flush() {
    memset(s, 0, sizeof(s));
    cnt = 1;
}

int lg_uni_seq() {
    flush();
    s[1] = a[1];
    for (int i = 2; i <= n; i++) {
        if (a[i] <= s[cnt]) {
            s[++cnt] = a[i];
        } else {
            int l = 1, r = cnt;
            while (l < r) {
                int mid = (l + r) >> 1;
                if (a[i] > s[mid]) {
                    r = mid;
                } else {
                    l = mid + 1;
                }
            }
            s[l] = a[i];
        }
    }
    return cnt;
}

int lg_inc_seq() {
    flush();
    s[1] = a[1];
    for (int i = 2; i <= n; i++) {
        if (a[i] > s[cnt]) {
            s[++cnt] = a[i];
        } else {
            int l = 1, r = cnt;
            while (l < r) {
                int mid = (l + r) >> 1;
                if (a[i] <= s[mid]) {
                    r = mid;
                } else {
                    l = mid + 1;
                }
            }
            s[l] = a[i];
        }
    }
    return cnt;
}

int main() {
    input();
    cout << lg_uni_seq() << endl << lg_inc_seq() << endl;
    return 0;
}

 在百忙之中写一篇题解也比较辛苦，别忘了点个赞！