参考资料:Palindromic Tree——回文树【处理一类回文串问题的强力工具】(请注意,其中似乎有一些错误)
回文自动机似乎和回文树是同一个东西qwq?
回文自动机(PAM)是一种处理回文串的工具。它的每个结点表示一个本质不同的回文串,转移边\(c\)表示在当前字符串的首尾分别加一个字符\(c\)。
回文自动机由两棵树组成,根结点分别称为\(odd\)和\(even\)。\(even\)表示空串,长度为\(0\),长度为偶数的回文串在它的子树上;\(odd\)表示一个“虚拟”的串,长度为\(-1\),长度为奇数的回文串在它的子树上。\(odd\)的直接儿子表示只有一个字母的回文串。沿着转移边\(c\)走一步就在当前串首尾各加上一个字符\(c\)。和AC自动机类似,一个结点的\(fail\)指针指向它的最长回文真后缀(定义\(fail[even]=odd\))。比如wqwqqwqwq的回文自动机长这样(数字表示结点编号,红箭头表示\(fail\)指针):
我画着画着发现这个字符串里回文串比想象的多
和后缀自动机类似,构造回文自动机也采用每次插入一个字符的方法。设原串是\(S\),当前位置是\(pos\),要加入的字符是\(c\),则可能会多一些以字符\(c\)结尾的回文串。而多的这些字符串可以看成是一个回文串\([a,pos-1]\)满足\(S_{a-1}=c\)后面加一个字符\(c\)。于是要找到最长的这样的回文串\([a,pos-1]\),即从\(pos-1\)这个结点开始爬\(fail\)链,直到\(p\)点满足\(S_{pos-len[p]-1}=S_{pos}\)。爬\(fail\)链最终会到长度为\(-1\)的\(even\),由于\(pos-(-1)-1=pos\),所以这个式子最终一定会成立。这个过程即代码中的\(get\_fail\)函数。
设第一个满足如上条件的点是\(p\)。如果\(p\)已经有了\(c\)字符的转移,则直接增加它的\(cnt\)(该字符串出现次数)即可;如果没有,则新建结点\(q\),\(q\)的长度显然是\(p\)的长度加\(2\),\(q\)的\(fail\)是从\(p\)的\(fail\)往上爬,找到第一个在后面加字符\(c\)仍为回文串的地方(方法同上述\(get\_fail\)),把它加字符\(c\)后转移到的点作为\(q\)的\(fail\)。
注意,如果要统计每个回文串的出现次数(即\(cnt\)),建完后要在\(fail\)树上做一遍树上递推(因为插入的时候只在当前点结尾的最长回文串的结点\(cnt\)上加\(1\)。如果一个串出现了,它的最长回文真后缀一定也出现了)。由于回文自动机是两棵树,所以不需要像后缀自动机求\(Right\)集合大小一样拓扑排序,只要按标号从大到小做即可。
题目:洛谷3649
把每个结点的长度\(len\)乘上出现次数\(cnt\)然后加起来就好了。
代码:
#include <cstdio>
#include <algorithm>
#include <cctype>
#include <cstring>
#include <string>
#define _ 0
using namespace std;
namespace zyt
{
const int N = 3e5 + 10;
template<typename T>
inline bool read(T &x)
{
char c;
bool f = false;
x = 0;
do
c = getchar();
while (c != EOF && c != '-' && !isdigit(c));
if (c == EOF)
return false;
if (c == '-')
f = true, c = getchar();
do
c = getchar();
while (isdigit(c));
if (f)
x = -x;
return true;
}
inline bool read(string &s)
{
static char buf[N];
if (scanf("%s", buf) == -1)
return false;
else
{
s = buf;
return true;
}
}
template<typename T>
inline void write(T x)
{
static char buf[20];
char *pos = buf;
if (x < 0)
putchar('-'), x = -x;
do
*pos++ = x % 10 + '0';
while (x /= 10);
while (pos > buf)
putchar(*--pos);
}
const int CH = 26;
typedef long long ll;
string s;
namespace Palindrome_Auto_Machine
{
struct node
{
int len, cnt, fail, s[CH];
}tree[N];
int cnt, last, odd, even, pos;
char s[N];
void init()
{
last = even = 0, odd = 1, cnt = 1, pos = 0;
s[0] = '#';
tree[odd].len = -1, tree[even].len = 0;
tree[odd].fail = tree[even].fail = odd;
}
int get_fail(int p)
{
while (s[pos - tree[p].len - 1] != s[pos])
p = tree[p].fail;
return p;
}
void insert(const char c)
{
s[++pos] = c;
int x = c - 'a', p = get_fail(last);
if (!tree[p].s[x])
{
tree[++cnt].len = tree[p].len + 2;
tree[cnt].fail = tree[get_fail(tree[p].fail)].s[x];
tree[p].s[x] = cnt;
}
last = tree[p].s[x];
++tree[last].cnt;
}
void build(const string &str)
{
for (int i = 0; i < str.size(); i++)
insert(str[i]);
for (int i = cnt; i > 0; i--)
tree[tree[i].fail].cnt += tree[i].cnt;
}
inline ll solve()
{
ll ans = 0;
for (int i = 0; i <= cnt; i++)
ans = max(ans, (ll)tree[i].cnt * tree[i].len);
return ans;
}
}
int work()
{
using Palindrome_Auto_Machine::init;
using Palindrome_Auto_Machine::build;
using Palindrome_Auto_Machine::solve;
read(s);
init();
build(s);
write(solve());
return (0^_^0);
}
}
int main()
{
return zyt::work();
}