SPOJ 1811 Longest Common Substring 后缀自动机
模板来源:http://www.neroysq.com/?p=76
思路:http://blog.sina.com.cn/s/blog_7812e98601012dfv.html
题意就是求两个字符串的最长公共子串,串长最大250000。
以串A构建一个后缀自动机,用串B来匹配。枚举串B的每一位B[i]即考虑串B中所有以B[i]为结尾的子串,维护的值为以B[i]为末尾能匹配的最大长度tmpL。
假设走到B[i]时已经匹配好的串为str,如果当前节点有B[i]这个儿子,直接向下走,++tmpL。
如果没有,沿着fail指针向前回退,直到找到一个有B[i]儿子的节点。
#include <cstdio> #include <cstring> #include <cstdlib> #include <algorithm> using namespace std; const int MAXN = 250010; const int SigmaSize = 26; struct sanode { sanode *f, *ch[SigmaSize]; int l; }; struct Suffix_Automaton { sanode pool[ ( MAXN << 1 ) + 20 ]; sanode *init; sanode *tail; int tot; void Clear() { tot = 0; init = pool; tail = init; return; } void Insert( int c ) { sanode *q = pool + ( ++tot ), *p = tail; q->l = p->l + 1; for ( ; p && !p->ch[c]; p = p->f ) p->ch[c] = q; tail = q; if ( !p ) q->f = init; else { if ( p->ch[c]->l == p->l + 1 ) q->f = p->ch[c]; else { sanode *r = pool + ( ++tot ), *u = p->ch[c]; *r = *u; r->l = p->l + 1; u->f = q->f = r; for ( ; p && p->ch[c] == u; p = p->f ) p->ch[c] = r; } } } }; char str[MAXN + 20]; Suffix_Automaton SAM; int main() { int len; scanf( "%s", str + 1 ); SAM.Clear(); len = strlen( str + 1 ); for ( int i = 1; i <= len; ++i ) SAM.Insert( str[i] - 'a' ); sanode *p = SAM.init; scanf( "%s", str + 1 ); len = strlen( str + 1 ); int tmpL = 0, ans = 0; for ( int i = 1; i <= len; ++i ) { if ( p->ch[ str[i] - 'a' ] ) //可以向下匹配的时候就继续向下匹配 p = p->ch[ str[i] - 'a' ], ++tmpL; else //如果当前p没有str[i]这个孩子 { while ( p && !p->ch[ str[i] - 'a' ] )
p = p->f; //沿着fail指针向前找,直到找到有str[i]儿子的结点,或者到根节点 if( p ) //如果能找到一个有str[i]儿子的节点 { tmpL = p->l + 1; p = p->ch[ str[i] - 'a' ]; } else //直到回到根也没有找到 { p = SAM.init; tmpL = 0; } } ans = max( ans, tmpL ); } printf( "%d\n", ans ); return 0; }