树状数组求最长上升子序列

先说最简单的作法：

一种是最多见的dp方法，令f[i]表示以A[i]元素结尾的LIS长度，那么,F[i]=max{F[j]+1) 其中1<=j<i，A[j]<A[i]，边界是初始化F[i]=1，复杂度O(n^2)。

#include <cstdio>
#include <algorithm>
using namespace std;
#define FOR(i,a,b) for(int i=(a);i<=(b);++i)

int n, a[1005], f[1005];
int main () {
  scanf("%d", &n);
  FOR(i, 1, n) scanf("%d", &a[i]), f[i] = 1;
  FOR(i, 1, n)
    FOR(j, 1, i - 1)
      if (a[j] < a[i])
        f[i] = max(f[i], f[j] + 1);
  int ans = 0;
  FOR(i, 1, n)
    ans = max(ans, f[i]);
  printf("%d\n", ans);
  return 0;
}

另外一种是有些贪心思想，利用二分：

若是子序列长度相同，那么最末尾的元素越小，就越有优点，因而对于长度相同的子序列，咱们老是用更小的来替换。

#include <cstdio>
#include <algorithm>
using namespace std;
const int inf=0x3f3f3f3f;
int a[1005], n;

int main() {
    for (int i = 0; i < 1005; ++i)
        a[i] = inf;
    scanf("%d", &n);
    for (int i = 0, t; i < 1005; ++i) {
        scanf("%d", &t);
        *lower_bound(a, a + n, t) = t;
    }
    printf("%d\n", lower_bound(a, a + n, inf) - a);
    return 0;
}

 

另外一种，用树状数组，这种方法和第一种方法思路是一致的，只不过在寻找比本身小的元素中，最长的那个LIS的过程，利用了树状数组来加速。以前的思想是：当前元素的LIS是位置在本身以前的，且LIS是最大的那个值+1。先不论树状数组是怎么工做的，对于原始数组v，复制一个数组a来排序，以后遍历数组v。这里因为是按顺序遍历v的，因此对于v[i]，先操做的元素确定位置上在v[i]以前。

而后，对于v[i]，获得其在a中的位置p，想办法获得在p以前的LIS长度是多少(就是操做query(p-1))，假设是q，那么ans=max(ans, q+1)。这里，因为对于v[i]实际操做的是p，p是v[i]值大小的相对位置，因而又知足：在位置p前面的那些元素，必然大小是<v[i]的，综上，的确能够赞成：这和以前LIS解法一的思想是一致的了。

如今的要讲树状数组怎么工做，即如何完成上面的”想办法“，无非就是要完成查询位置p以前的最大值（区间最大值），以及对于当前元素v[i]和其在a中位置的p，获得LIS后更新位置p以前的LIS。以往BIT老是用来求区间和，其实这只是它的一个应用，说到底它是一个数据结构，方法就是将sum和add的加法操做变成最值操做。  理解到这儿：其实，这彻底能够用线段数来实现，只不过这里要处理的区间刚好是[1, p]，即刚好必定是从1开始的，用树状数据能简单解决，就没必要用线段数了。

#include <cstdio>
#include <algorithm>
using namespace std;

#define FOR(i,a,b) for(int i=(a);i<=(b);++i)
#define maxN 10005
int n, L, v[maxN], a[maxN], T[maxN];

void update(int i, int x) { 
  for (; i <= L; i += i & -i)
    T[i] = max(T[i], x);
}
int query(int i) {
  int ans = 0;
  for (; i; i -= i & -i)
    ans = max(ans, T[i]);
  return ans;
}
int main () {
  scanf("%d", &n);
  FOR(i, 0, n - 1) scanf("%d", &v[i]), a[i] = v[i];
  sort(a, a + n);
  L = unique(a, a + n) - a;
  int ans = 1, t;
  FOR(i, 0, n - 1) {
    int p = lower_bound(a, a + L, v[i]) - a + 1;
    t = query(p - 1) + 1;
    ans = max(ans, t);
    update(p, t);
  }
  printf("%d\n", ans);
  return 0;
}