SPOJ LCS2
Posted mjtcn
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了SPOJ LCS2相关的知识,希望对你有一定的参考价值。
LCS2 - Longest Common Substring II
题意:
求N(N<=10)个串的最长公共子串。
分析:
poj2774上那道题,对一个串建立后缀自动机,另一个在上面匹配。
这道题是对多个串求。那么同样,让每个串在后缀自动机上匹配,然后记录在后缀自动机的每个节点上记录,当前串在这个位置的最大匹配数,h数组。
然后mn数组,每次对于这所有的节点的h取小,为从第2个串到现在所有的串,都能在这个节点上匹配的长度。
因为一旦某个节点匹配上了,那么它的父节点(parent树)的父节点都会匹配上(因为父节点是当前点的后缀),所以按拓扑倒序,更新父节点的h,为父节点的len,(即最大长度)。
代码:
1 #include<bits/stdc++.h> 2 using namespace std; 3 typedef long long LL; 4 5 inline int read() { 6 int x=0,f=1;char ch=getchar();for(;!isdigit(ch);ch=getchar())if(ch==‘-‘)f=-1; 7 for (;isdigit(ch);ch=getchar())x=x*10+ch-‘0‘;return x*f; 8 } 9 10 const int N = 100010; 11 12 struct Suffix_Automaton{ 13 int fa[N<<1], trans[N<<1][26], len[N<<1], Last, Index; 14 int v[N], sa[N<<1], mn[N<<1], h[N<<1]; 15 char s[N]; 16 17 void extend(int c) { 18 int P = Last,NP = ++Index; 19 len[NP] = len[P] + 1; 20 for (; P&&!trans[P][c]; P=fa[P]) trans[P][c] = NP; 21 if (!P) fa[NP] = 1; //- 22 else { 23 int Q = trans[P][c]; 24 if (len[P] + 1 == len[Q]) fa[NP] = Q; 25 else { 26 int NQ = ++Index; 27 fa[NQ] = fa[Q]; 28 len[NQ] = len[P] + 1; 29 memcpy(trans[NQ], trans[Q], sizeof trans[Q]); 30 fa[Q] = NQ; 31 fa[NP] = NQ; 32 for (; P&&trans[P][c]==Q; P=fa[P]) trans[P][c] = NQ; 33 } 34 } 35 Last = NP; 36 } 37 void build() { 38 Last = Index = 1; 39 scanf("%s",s+1); 40 int n = strlen(s+1); 41 for (int i=1; i<=n; ++i) extend(s[i] - ‘a‘); 42 for (int i=1; i<=Index; ++i) v[len[i]] ++; 43 for (int i=1; i<=n; ++i) v[i] += v[i-1]; 44 for (int i=1; i<=Index; ++i) sa[ v[len[i]]-- ] = i; // sa[i] 排名为i的节点。按深度排名(拓扑用) 45 } 46 void calcc() { 47 int n = strlen(s+1), now = 0, p = 1; 48 memset(h, 0, sizeof(h)); 49 for (int i=1; i<=n; ++i) { 50 int c = s[i] - ‘a‘; 51 if (trans[p][c]) p = trans[p][c], now ++; 52 else { 53 for (; p&&!trans[p][c]; p=fa[p]); 54 if (!p) now = 0, p = 1; 55 else now = len[p] + 1, p = trans[p][c]; 56 } 57 h[p] = max(h[p], now); 58 } 59 for (int i=Index; i>=1; --i) { // 拓扑倒序,parent树中从深度深的到浅的 60 int t = sa[i]; 61 mn[t] = min(mn[t], h[t]); 62 if (h[t] && fa[t]) h[fa[t]] = len[fa[t]]; 63 } 64 } 65 void solve() { 66 build(); 67 memset(mn, 0x3f, sizeof(mn)); 68 while (scanf("%s",s+1) != EOF) calcc(); 69 int ans = 0; 70 for (int i=1; i<=Index; ++i) ans = max(ans, mn[i]); 71 printf("%d",ans); 72 } 73 }sam; 74 75 int main() { 76 sam.solve(); 77 return 0; 78 }
以上是关于SPOJ LCS2的主要内容,如果未能解决你的问题,请参考以下文章