首页 > 其他 > 详细

SPOJ LCS 后缀自动机

时间:2014-02-17 04:22:33      阅读:523      评论:0      收藏:0      [点我收藏+]

用后缀自动机求两个长串的最长公共子串,效果拔群。多样例的时候memset要去掉。

解题思路就是跟CLJ的一模一样啦。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
#pragma warning(disable:4996)
#include<cstring>
#include<string>
#include<iostream>
#include<cmath>
#include<vector>
#include<algorithm>
#define maxn 250050
using namespace std;
 
struct State{
    State *suf, *go[26];
    int val;
    State() :suf(0), val(0){
        memset(go, 0, sizeof(go));
    }
}*root, *last;
 
State statePool[maxn * 2], *cur;
 
void init()
{
    cur = statePool;
    root = last = cur++;
}
 
void extend(int w)
{
    State *p = last, *np = cur++;
    np->val = p->val + 1;
 
    while (p&&!p->go[w]) p->go[w] = np, p = p->suf;
    if (!p) np->suf = root;
    else{
        State *q = p->go[w];
        if (p->val + 1 == q->val){
            np->suf = q;
        }
        else{
            State *nq = cur++;
            memcpy(nq->go, q->go, sizeof q->go);
            nq->val = p->val + 1;
            nq->suf = q->suf;
            q->suf = nq;
            np->suf = nq;
            while (p&&p->go[w] == q){
                p->go[w] = nq, p = p->suf;
            }
        }
    }
    last = np;
}
 
char stra[maxn], strb[maxn];
 
int main()
{
    while (~scanf("%s%s", stra, strb))
    {
        init();
        //memset(statePool, 0, sizeof(statePool));
        int lena = strlen(stra);
        for (int i = 0; i < lena; i++){
            extend(stra[i] - ‘a‘);
        }
        int ans = 0;
        int lenb = strlen(strb);
        State *p = root;
        int len = 0;
        for (int i = 0; i < lenb; i++){
            if (p->go[strb[i] - ‘a‘]){
                len++; ans = max(ans, len);
                p = p->go[strb[i] - ‘a‘];
            }
            else{
                while (p&&!p->go[strb[i] - ‘a‘]){
                    p = p->suf;
                }
                if (!p) {
                    p = root;
                    len = 0;
                }
                else{
                    len = p->val + 1;
                    ans = max(len, ans);
                    p = p->go[strb[i] - ‘a‘];
                }
            }
        }
        printf("%d\n", ans);
    }
    return 0;
}

SPOJ LCS 后缀自动机

原文:http://www.cnblogs.com/chanme/p/3551560.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!