【NOIP1999普及组】导弹拦截题解（DP+单调性优化+二分）

题目：luogu1020.
题目大意：第一问求一个最长不上升子序列，第二问求方案数.
设序列长度为$n$，$1\le n\le 10^5$.

我们现在知道LIS问题有一个$O(n^2)$的算法，但是这个算法在这个问题中并不能拿到满分.

我们现在考虑第一个问题的优化，我们考虑到枚举到点$i$的时候，最长的长度为$f[i]$，然后我们再记录一个值$g[i]$，表示长度为$i$的子序列的结尾的最大值.

显然，$g[i]$具有单调性，也就是说$g[i]$会越来越小，所以直接每次转移的时候二分找到$g[i]$刚好比$a[i]$小或等于的$i+1$.时间复杂度$O(n\log n)$.

根据Dilworth定理，第二问答案就是最长不下降子序列的长度.

代码如下：

#include<bits/stdc++.h>
  using namespace std;
#define Abigail inline void
const int N=100000;
const int INF=1<<31-1;
int n,a[N+9],f[N+9],g[N+9],ans1,ans2;
Abigail into(){
  while (~scanf("%d",&a[++n]));
}
Abigail work(){
  --n;
  int l,r,mid,now;
  f[1]=1;g[1]=a[1];now=1;g[0]=INF;
  for (int i=2;i<=n;i++){
    l=0;r=now;
    while (l+1<r){
      mid=l+r>>1;
      if (g[mid]>=a[i]) l=mid;
      else r=mid;
    }
    if (g[r]>=a[i]) mid=r;
    else mid=l;
    f[i]=mid+1;
    if (now<f[i]) now++,g[f[i]]=a[i];
    else g[f[i]]=max(g[f[i]],a[i]);
    ans1=max(ans1,f[i]);
  }
  f[1]=1;g[1]=a[1];now=1;g[0]=-INF;
  for (int i=2;i<=n;i++){
    l=0;r=now;
    while (l+1<r){
      mid=l+r>>1;
      if (g[mid]<a[i]) l=mid;
      else r=mid;
    }
    if (g[r]<a[i]) mid=r;
    else mid=l;
    f[i]=mid+1;
    if (now<f[i]) now++,g[f[i]]=a[i];
    else g[f[i]]=min(g[f[i]],a[i]);
    ans2=max(ans2,f[i]);
  }
}
Abigail outo(){
  printf("%d\n%d\n",ans1,ans2);
}
int main(){
  into();
  work();
  outo();
  return 0;
}