有的算法只会模板题真不行,有的算法没有理解深刻照样不会做题,今天又看了一遍回文树,加深了一下理解,顺便整理一下模板
超好理解
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;
const int MAXN = 100050, SIZE = 26;
/*1.len[i]表示编号为i的节点表示的回文串的长度(一个节点表示一个回文串)
2.next[i][c]表示编号为i的节点表示的回文串在两边添加字符c以后变成的回文串的编号(和字典树类似)。
3.fail[i]表示节点i失配以后跳转不等于自身的节点i表示的回文串的最长后缀回文串(和AC自动机类似)。
4.cnt[i]表示节点i表示的本质不同的串的个数(建树时求出的不是完全的,最后count()函数跑一遍以后才是正确的)
5.num[i]表示以节点i表示的最长回文串的最右端点为回文串结尾的回文串个数。
6.last指向新添加一个字母后所形成的最长回文串表示的节点。
7.S[i]表示第i次添加的字符(一开始设S[0] = -1(可以是任意一个在串S中不会出现的字符))。
8.p表示添加的节点个数。
9.n表示添加的字符个数。
*/
struct Palindromic_Tree
{
int next[MAXN][SIZE];//next指针
int fail[MAXN];//fail指针
int cnt[MAXN];//表示节点i表示的回文串的个数(建树时求出的不是完全的,最后Count()函数跑一遍以后才是正确的)
int num[MAXN];//表示以节点i表示的最长回文串的最右端点为回文串结尾的回文串个数
int len[MAXN];//len[i]表示节点i表示的回文串的长度(一个节点表示一个回文串)
int S[MAXN];//存放添加的字符
int last;//指向上一个字符所在的节点,方便下一次add(指向新添加一个字母后所形成的最长回文串表示的节点)
int n;//字符数组指针(表示添加的字符个数。)
int p;//节点指针(表示添加的节点个数)
int NewNode(int L) {//新建节点
for(int i = 0; i < SIZE; ++i)
next[p][i] = 0;
cnt[p] = num[p] = 0;
len[p] = L;
return p++;
}
void Init() {//初始化
p = n = 0;
NewNode(0); NewNode(-1);
last = 0;
S[n] = -1;//开头放一个字符集中没有的字符,减少特判
fail[0] = 1;
}
int GetFail(int x) {//和KMP一样,失配后找一个尽量最长的
while(S[n - len[x] - 1] != S[n]) x = fail[x];
return x;
}
void Add(int c) {
S[++n] = c;
int cur = GetFail(last);//通过上一个回文串找这个回文串的匹配位置
if(!next[cur][c]) {//如果这个回文串没有出现过,说明出现了一个新的本质不同的回文串
int now = NewNode(len[cur] + 2);//新建节点
fail[now] = next[GetFail(fail[cur])][c];//和AC自动机一样建立fail指针,以便失配后跳转
next[cur][c] = now;
num[now] = num[fail[now]] + 1;
}
last = next[cur][c];
cnt[last]++;
}
void Count() {
for(int i = p - 1; i >= 0; --i)
cnt[fail[i]] += cnt[i];
//父亲累加儿子的cnt,因为如果fail[v]=u,则u一定是v的子回文串!
}
}tree;
char str[MAXN];
int main(){
while(~scanf("%s",str))
{
tree.Init();
int len=strlen(str);
for(int i=0;i<len;i++)
{
tree.Add(str[i]-'a');
printf("%d",tree.p-2);
if(i==len-1) printf("\n");
else printf(" ");
}
}
return 0;
}