manacher

马拉车

今天做到一道字符串找最大回文子串的题，蠢蠢的我直接暴力，两种暴力都做了，第一个O(n^3)，从左到右枚举，很显然会tle，然后进行了一下优化，从中间开始进行枚举，本来以为可以了，
结果O(n^2)，还是TLE，由于第一次接触到这样的算法题，直接面向csdn搜了一些题解，都用的一个方法，叫做manacher，也就是马拉车算法，为了做这个题我也是花了四五个小时的时间来看这个算法，大概有了一些思路后，看着其他大佬的代码自己也是写了出来
先讲一下我自己对manacher的理解吧
这个算法第一个让人惊艳的地方，就是通过“插空”操作，把奇数和偶数长度的字符串均改为长度为奇数的字符串，比如字符串abcde，操作后为
!#a#b#c#d#e#@其中操作后的字符串第一个和最后一个元素为“哨兵”，就不需要考虑边界问题，然后插空操作的字符就要选择和输入字符串不一样的字符，尽量选择不常用的吧
然后这个算法的需要的一些数据
p[],p数组是来判断某一个位置最大回文串的半径比如说i的位置最大的回文串长度是5那么半径就是3，即p[i]=3（以某一个位置的字符为中心的回文串长度一定是奇数，可以自己写几个试一下；
还有坐标参数，mid和rt，mid的含义是现在已知的最长的回文串的最中间位置的下表，rt则是已知的最长的回文串最远能到达的地方，可以rt与mid是绑定的，二者一定是同时发生变化的。
还需要一个ans来记录每一次更新mid与rt求得的回文串的长度，取最大
接下来的一步也是重点（难，）怎样通过1到i-1的p[]来计算p[i],首先分两种情况，第一个是i要小于rt，也就是说i仍在之前所求的最长回文串之间，这个时候就要充分利用回文串对称的性质，一定有一个j和i是关于mid对称的，然而p[j]是我们已经求出来的，这时候就要判断一下，如果i加上p[j]是不是会大于rt，此时的p[i]应该是从p[j]与rt-i中选取最小值。第二种情况就是i大于rt，这个时候无法通过前边所求来求p[i],只能枚举。如果说此时的i+p[i]大于rt，那么就要更新rt与mid，然后进行后边的回文串的长度的求解
代码如下：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string.h>
using namespace std;
const int N=2000000;
char s[N],str[N];
int p[N];
int n;
int m;
void manacher()
{
    int mid=0;
    int rt=0;
    int res=0;
    for(int i=1;i<=m;i++)
    {
       p[i]=i<rt?min(p[2*mid-i],rt-i):1;
       while(str[i+p[i]]==str[i-p[i]])p[i]++;
       if(i+p[i]>rt)
       {
           rt=i+p[i];
           mid=i;
       }
       res=max(res,p[i]-1);
    }
     cout<<res<<endl;   
}
int main()
{
   str[0]='!';
   str[1]='#';
   scanf("%s",s);
   n=strlen(s);
   for(int i=0;i<n;i++)
   {
       str[i*2+2]=s[i];
       str[i*2+3]='#';
   }
   m=2*n+1;
   str[m+1]='#';
   manacher();
   return 0;

}

欢迎指正