用后缀自动机求两个长串的最长公共子串,效果拔群。多样例的时候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; } |
原文:http://www.cnblogs.com/chanme/p/3551560.html