导弹拦截

题目描述

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

输入导弹依次飞来的高度（雷达给出的高度数据是 \le 50000

$\le 50000$ 的正整数），计算这套系统最多能拦截多少导弹，如果要拦截所有导弹最少要配备多少套这种导弹拦截系统。

思路：Dilworth定理：求一个序列里面最少有多少最长不上升序列等于求这个序列里最长上升序列的长度（不会证明）所以这道题变成求最长上升序列和最长不降序列。

树状数组做法（注意在求两个值的时候树状数组里存的是0~a[i] 的最大值 和 a[i]~top的最大值）

#include<bits/stdc++.h>
#define read(x) scanf("%d",&x)
using namespace std;
const int maxn=1e6+6;
typedef long long LL;
int c[maxn];
int f[maxn];
int a[maxn],top;
int lowbit(int x) {
    return x&(-x);
}
int get1(int x) {
    int ret=0;
    while(x>0) {
        ret=max(f[x],ret);
        x-=lowbit(x);
    }
    return ret;
}
void add1(int x,int v) {
    while(x<=top) {
        f[x]=max(f[x],v);
        x+=lowbit(x);
    }
}
int get2(int x) {
    int ret=0;
    while(x<=top) {
        ret=max(f[x],ret);
        x+=lowbit(x);
    }
    return ret;
}
void add2(int x,int v) {
    while(x>0) {
        f[x]=max(f[x],v);
        x-=lowbit(x);
    }
}
int cnt=1;
int main() {
    while(read(a[cnt]) != EOF) {
        top=max(top,a[cnt]);
        c[cnt]=a[cnt];
        cnt++;
    }
    cnt--;
    int ans=0,num=0;
    for(int i=1;i<=cnt;i++) {
        int x=get1(a[i])+1;
        ans=max(ans,x);
        add1(a[i]+1,x);
    }
    memset(f,0,sizeof(f));
    for(int i=1;i<=cnt;i++) {
        int x=get2(a[i])+1;
        num=max(num,x);
        add2(a[i],x);
    }
    printf("%d\n%d\n",num,ans);
    return 0;
}