【算法】【C++】初学算法，在此记录

0.预留位置

1.使用Manacher算法求最长回文子串

给出一个字符串，要求计算出这一字符串的最长回文子串的长度。
如果遍历每一个字符，并以该字符为中心向两边查找，则其时间
复杂度为O(n2)。
Manacher算法，又被戏称为“马拉车”算法，可以在时间复杂度
为O(n)的情况下求解一个字符串的最长回文子串的长度。

1.

回文串分为奇回文（abcba）和偶回文（baccba），为了避免区分奇偶的麻烦，把字符串变一下形式，在字符串之间加入特殊符号，并且首位再加不同的特殊符号以防越界，于是字符串变成奇回文（上述偶回文得到的新字符串：@#a#b#c#c#b#a!）

2.

给得到的新字符串str_new定义一个辅助数组int p[]，p[i]表示在str_new中以第i个位置为中心得到的最长回文串的半径

3.

可想而知，若新字符串str_new的最大半径为p[i]，则原字符串的最长回文子串长度为p[i]-1

4.

对于p[i]，如果i<mx，设j是i关于id对称点，如图所示，则基于以下三
种情况，可以求出p[i]的值：
（1）以j为中心的回文串有一部分在以id为中心的回文串之外。因为mx
是以id为中心的最长回文的右边界，所以以i为中心的回文串不可能会
有字符在以id为中心的回文串之外；否则mx就不是以id为中心的最长回
文的右边界。所以，在这种情况下，p[i]=mx–i。
（2）以j为中心的回文串全部在以id为中心的回文串的内部，则
p[i]=p[j]，而且p[i]不可能再增加。
（3）以j为中心的回文串的左端正好与以id为中心的回文串的左端重合。
则p[i]=p[j]或p[i]=mx–i，并且p[i]还有可能会继续增加，即while
(str_new[i-p[i]]==str_new[i+p[i]]) p[i]++;

5.

从str[1]开始遍历到末尾
先通过

if(i<mx)    p[i]=std::min(p[2*id-i],mx-i);
else        p[i]=1;

确定好id-mx范围内的以i为中心的最长子串半径，然后再通过

while(str[i+p[i]]==str[i-p[i]])     p[i]+=1;

确定最长半径

若得到以i为中心的子串最右位置超过mx，则更新范围

if(p[i]+i>mx)           //超出最右边就更新范围
        {
   
            id=i;
            mx=p[i]+i;
            mxsub=std::max(mxsub,p[i]);     //更新最长子串长度
        }

代码：输入一个字符串，用Manacher算法求其中最长回文子串的长度

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
const int max_len = 1e6+3;      //输入字符串的最大长度
const int inPos=max_len+2;      //字符串输入位置，从str数组中间输入字符串
int p[max_len*2];               //p[i]:以i为中心的回文串的最大长度
int len;
char str[max_len*2];            //用于储存输入的字符串

//把字符串各字符之间插入字符，首位再各插入不同字符防止越界
void change_str()
{
   
    int index=0;
    len=strlen(str);            //获得输入字符串的长度
    str[index++]='@';
    for(int i=inPos; i<len; i++)
    {
   
        str[index++]='#';
        str[index++]=str[i];
    }
    str[index++]='#';
    str[index]='!';
    len=index;                  //获得改造后字符串的长度
}

int manacher()
{
   
    int mx=0,id=1;              //mx表示最右端位置 id表示中心点位置
    int mxsub=0;                //最长子串长度
    for(int i=1; i<len; i++)