[BZOJ1264][AHOI2006]基因匹配Match（DP + 树状数组）

传送门

有点类似LCS，可以把 a[i] 在 b 串中的位置用一个链式前向星串起来，由于链式前向星是从后往前遍历，所以可以直接搞。

状态转移方程 f[i] = max(f[j]) + 1 ( 1 <= j  < i && a[i] == b[j] )

——代码

 #include <cstdio>
 #include <cstring>
 #include <iostream>

 using namespace std;

 const int MAXN = ;
 int n, cnt, ans;
 int a[MAXN], c[MAXN];
 int head[MAXN], to[MAXN << ], next[MAXN << ];

 inline int query(int x)
 {
     int ret = ;
     for(; x; x -= x & -x) ret = max(ret, c[x]);
     return ret;
 }

 inline void update(int x, int d)
 {
     for(; x <= n; x += x & -x) c[x] = max(c[x], d);
 }

 inline void add(int x, int y)
 {
     to[cnt] = y;
     next[cnt] = head[x];
     head[x] = cnt++;
 }

 int main()
 {
     int i, j, x;
     scanf("%d", &n);
     n *= ;
     memset(head, -, sizeof(head));
     for(i = ; i <= n; i++) scanf("%d", &a[i]);
     for(i = ; i <= n; i++) scanf("%d", &x), add(x, i);
     for(i = ; i <= n; i++)
         for(j = head[a[i]]; j != -; j = next[j])
         {
             x = query(to[j] - ) + ;
             update(to[j], x);
             ans = max(ans, x);
         }
     printf("%d", ans);
     return ;
 }