[BZOJ1264][AHOI2006]基因匹配Match(DP + 树状数组)

Posted 蒟蒻zht的博客

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[BZOJ1264][AHOI2006]基因匹配Match(DP + 树状数组)相关的知识,希望对你有一定的参考价值。

传送门

 

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

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

 

——代码

技术分享
 1 #include <cstdio>
 2 #include <cstring>
 3 #include <iostream>
 4 
 5 using namespace std;
 6 
 7 const int MAXN = 200001;
 8 int n, cnt, ans;
 9 int a[MAXN], c[MAXN];
10 int head[MAXN], to[MAXN << 4], next[MAXN << 4];
11 
12 inline int query(int x)
13 {
14     int ret = 0;
15     for(; x; x -= x & -x) ret = max(ret, c[x]);
16     return ret;
17 }
18 
19 inline void update(int x, int d)
20 {
21     for(; x <= n; x += x & -x) c[x] = max(c[x], d);
22 }
23 
24 inline void add(int x, int y)
25 {
26     to[cnt] = y;
27     next[cnt] = head[x];
28     head[x] = cnt++;
29 }
30 
31 int main()
32 {
33     int i, j, x;
34     scanf("%d", &n);
35     n *= 5;
36     memset(head, -1, sizeof(head));
37     for(i = 1; i <= n; i++) scanf("%d", &a[i]);
38     for(i = 1; i <= n; i++) scanf("%d", &x), add(x, i);
39     for(i = 1; i <= n; i++)
40         for(j = head[a[i]]; j != -1; j = next[j])
41         {
42             x = query(to[j] - 1) + 1;
43             update(to[j], x);
44             ans = max(ans, x);
45         }
46     printf("%d", ans);
47     return 0;
48 }
View Code

 

以上是关于[BZOJ1264][AHOI2006]基因匹配Match(DP + 树状数组)的主要内容,如果未能解决你的问题,请参考以下文章

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

BZOJ 1264 [AHOI2006]基因匹配Match

bzoj1264: [AHOI2006]基因匹配Match

BZOJ1264 [AHOI2006]基因匹配Match 动态规划 树状数组

bzoj 1264: [AHOI2006]基因匹配Match (树状数组优化dp)

[BZOJ1264][AHOI2006]基因匹配Match(DP + 树状数组)