Palindromes and Super Abilities URAL - 1960 (回文树板子)

Palindromes and Super Abilities URAL - 1960 (回文树板子)

After solving seven problems on Timus Online Judge with a word “palindrome” in the problem name, Misha has got an unusual ability. Now, when he reads a word, he can mentally count the number of unique nonempty substrings of this word that are palindromes.
Dima wants to test Misha’s new ability. He adds letters s 1, …, s n to a word, letter by letter, and after every letter asks Misha, how many different nonempty palindromes current word contains as substrings. Which n numbers will Misha say, if he will never be wrong?
Input
The only line of input contains the string s 1… s n, where s i are small English letters (1 ≤ n ≤ 10 5).
Output
Output n numbers separated by whitespaces, i-th of these numbers must be the number of different nonempty substrings of prefix s 1… s i that are palindromes.
Example
input output
aba
1 2 3

题目大意:

给定一个字符串
输出前i个字符串的不同本质的回文串的个数

解题思路：

回文树板子
回文树介绍
首先我们定义一些变量。

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表示添加的字符个数。

AC代码:

#include <bits/stdc++.h>
using namespace std; 
const int maxn= 1e5+10 ;
const int N = 26 ; 
struct Palindromic_Tree {
	int Next[maxn][N] ;//Next指针，Next指针和字典树类似，指向的串为当前串两端加上同一个字符构成
	int fail[maxn] ;//fail指针，失配后跳转到fail指针指向的节点
	int cnt[maxn] ;//cnt[i]表示节点i表示的本质不同的串的个数（建树时求出的不是完全的，最后count()函数跑一遍以后才是正确的）
	int num[maxn] ;//
	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 < N ; ++ i ) Next[p][i] = 0 ;
		cnt[p] = 0 ;
		num[p] = 0 ;
		len[p] = l ;
		return p ++ ;
	}
 
	void init () {//初始化
		p = 0 ;
		newnode (  0 ) ;
		newnode ( -1 ) ;
		last = 0 ;
		n = 0 ;
		S[n] = -1 ;//开头放一个字符集中没有的字符，减少特判
		fail[0] = 1 ;
	}
 
	int get_fail ( int x ) {//和KMP一样，失配后找一个尽量最长的
		while ( S[n - len[x] - 1] != S[n] ) x = fail[x] ;
		return x ;
	}
 
	void add ( int c ) {
		c -= 'a' ;
		S[++ n] = c ;
		int cur = get_fail ( last ) ;//通过上一个回文串找这个回文串的匹配位置
		if ( !Next[cur][c] ) {//如果这个回文串没有出现过，说明出现了一个新的本质不同的回文串
			int now = newnode ( len[cur] + 2 ) ;//新建节点
			fail[now] = Next[get_fail ( 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;
string s;
int ans[maxn];
int main(){
	cin>>s;
	tree.init();
	int len = s.length();
	for(int i=0;i<len;i++){
		tree.add(s[i]);
		ans[i] = tree.p-2;
	}
	for(int i=0;i<len;i++){
		cout<<ans[i]<<" ";
	}
	cout<<endl;
	return 0;
}