CodeChef - PALPROB Contest Page

题目：一个长度为len的回文串的回文值定义1+这个字符串的前len/2个字符组成的字符串的回文值，如果一个字符串不是回文串，那么它的回文值为0。

给你一个字符串，求它的回文子串的回文值之和

思路：用half来保存每个字符串的前len/2个字符组成的回文串的节点

代码：

#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<iostream>
#include<algorithm>
#include<ctime>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<string>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<stack>
#include<list>
#include<numeric>
using namespace std;
#define LL long long
#define ULL unsigned long long
#define INF 0x3f3f3f3f3f3f3f3f
#define mm(a,b) memset(a,b,sizeof(a))
#define PP puts("*********************");
template<class T> T f_abs(T a){ return a > 0 ? a : -a; }
template<class T> T gcd(T a, T b){ return b ? gcd(b, a%b) : a; }
template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}
// 0x3f3f3f3f3f3f3f3f

const int MAXN = 100005, SIZE = 26;
struct Palindromic_Tree {
	int next[MAXN][SIZE];//next指针
	int fail[MAXN];//fail指针
	LL cnt[MAXN];//表示节点i表示的本质不同的串的个数(建树时求出的不是完全的,最后Count()函数跑一遍以后才是正确的)
	int sc[MAXN];
	int half[MAXN];
	int len[MAXN];//len[i]表示节点i表示的回文串的长度
	int S[MAXN];//存放添加的字符 S[MAXN << 1]
	int last;//指向上一个字符所在的节点,方便下一次add
	int n;//字符数组指针 //lbd, rbd
	int p;//节点指针
	/*int lbd,rbd;
	int last[2];//0表示在字符串的首部加上一个字符，1表示在尾部加上一个字符*/
	int NewNode(int l) {//新建节点
		for(int i = 0; i < SIZE; ++i)
			next[p][i] = 0;
		cnt[p] = sc[p] = 0;
		len[p] = l;
		return p++;
	}
	void Init(/*int x*/) {//初始化
		p = n = 0;
		NewNode(0); NewNode(-1);
		last = 0;
		S[n] = -1;
		fail[0] = 1;
		half[0]=half[1]=1;
	}
	int GetFail(int x/*, int d*/) {//和KMP一样,失配后找一个尽量最长的
		while(S[n - len[x] - 1] != S[n]) x = fail[x];
		return x;
	}
	void Add(int c/*, int d*/) {//d=0表示在首部添加字符,d=1表示在尾部添加字符
		S[++n] = c;
		int cur = GetFail(last);//通过上一个回文串找这个回文串的匹配位置
		if(!next[cur][c]) {//如果这个回文串没有出现过,说明出现了一个新的本质不同的回文串
			int now = NewNode(len[cur] + 2);//新建节点
			fail[now] = next[GetFail(fail[cur]/*, d*/)][c];//和AC自动机一样建立fail指针,以便失配后跳转
			next[cur][c] = now;
			if(len[now]==1) half[now]=0;
			else{
                 int pos=half[cur];
                 while(S[n-len[pos]-1]!=S[n]||2*len[pos]+4>len[now])
                    pos=fail[pos];
                 half[now]=next[pos][c];
			}
			sc[now]=1+(len[now]/2==len[half[now]]?sc[half[now]]:0);
		}
		last = next[cur][c];
		cnt[last]++;
	}
	LL Count() {
	    LL sum=0;
		for(int i = p - 1; i >= 0; --i){
			cnt[fail[i]] += cnt[i];
			sum+=cnt[i]*sc[i];
		}
		return sum;
	}
}tree;
char str[MAXN];
int main(){

    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%s",str);
        tree.Init();
        for(int i=0;str[i]!='\0';i++){
            tree.Add(str[i]-'a');
        }
        printf("%lld\n",tree.Count());
    }
    return 0;
}