【题解】Casting Spells LA 4975 UVa 1470 双倍回文 SDOI 2011 BZOJ 2342 Manacher

首先要吐槽LRJ，书上给的算法标签是“有难度，需要结合其他数据结构”，学完Manacher才发现几乎一裸题

题目的意思是问原串中有多少个ww^Rww^R这样的子串，其中w^R表示w的反串

比较容易看出来，ww^Rww^R本身是一个回文串，ww^R也是一个回文串

最裸的暴力是，我们枚举每一个回文串，然后判断这个回文串的左半边是不是也是个回文串

然后我们考虑用Manacher

我们考虑Manacher的工作原理，是在充分利用原先的信息的前提下，不重复，不遗漏的枚举每个回文串

也就是说，在Manacher的运算过程当中，每个回文串我们都会考虑到

也就是说，我们可以在Manacher的运作过程当中，顺便完成答案的计算，具体操作见下

我们可以画出来这样一张图

# a # b # b # a # a # b # b # a #

0 1 2 3 4 5 6 7 8

我们要更新答案的时候，处在最中间的8号位置，不失一般性，我们处在 i 号位置，要求这个位置是字符#，这个要求的原因很显然

当 i 号位置的回文半径到达4的倍数的时候，说明我们左半边的串的长度是偶数，设当前的回文半径为 r

这时候左半边的串的中心位置就是 i - r/2 ，可以自己手算一下

如果 i - r/2 处的回文半径大于等于 i - r/2 的话，那么显然左半边是回文串，我们找到了一个满足要求子串，就可以ans = max( ans, r )

代码如下：

 #include <cstring>
 #include <cstdio>
 #include <algorithm>

 using namespace std;
 const int MAXN = ;

 int n;
 char str[MAXN], s[MAXN<<];

 void input() {
     scanf( "%s", str ), n = strlen(str);
     int p = ;
     for( int i = ; i < n; ++i )
         s[p++] = '#', s[p++] = str[i];
     s[p++] = '#';
 }

 int rd[MAXN<<]; // 回文半径
 void manacher() {
     int mx = , p = , len = *n+, ans = ;
     for( int i = ; i < len; ++i ) {
         if( i < mx ) rd[i] = min( rd[*p-i], mx-i );
         else rd[i] = ;
         while( i+rd[i] < len && i-rd[i] >=  && s[i+rd[i]] == s[i-rd[i]] ) {
             if( s[i] == '#' && rd[i] %  ==  && rd[i-rd[i]/] >= rd[i]/ )
                 ans = max( ans, rd[i] );
             ++rd[i];
         }
         if( i+rd[i] > mx ) mx = i+rd[i], p = i;
     }
     printf( "%d\n", ans );
 }

 int main() {
     int T; scanf( "%d", &T );
     while( T-- ) input(), manacher();
     return ;
 }