回文子串的最大长度(哈希,二分)

回文子串的最大长度

题目链接
如果一个字符串正着读和倒着读是一样的，则称它是回文的。
给定一个长度为 N 的字符串 S，求他的最长回文子串的长度是多少。

输入格式

输入将包含最多 30 个测试用例，每个测试用例占一行，以最多 1000000 个小写字符的形式给出。
输入以一个以字符串 END 开头的行表示输入终止。

输出格式

对于输入中的每个测试用例，输出测试用例编号和最大回文子串的长度（参考样例格式）。
每个输出占一行。

输入样例：

abcbabcbabcba
abacacbaaaab
END

输出样例：

Case 1: 13
Case 2: 6

算法分析

子串为 字符串连续的一段,判断回文我们采用哈希的方式判断,选取子串的情况我们用,枚举中间点(用马拉车方法将字符串变为奇数个) 之后二分操作的半径,找到最大的满足半径,之后判断马拉车后的字符个数.

代码分析

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
using namespace std;
const int maxn=2000010,P=131;
typedef unsigned long long ull;
ull hl[maxn],hr[maxn],p[maxn];
char str[maxn];
ull get(ull h[],int l,int r )
{
    return h[r]-h[l-1]*p[r-l+1];
}
int main()
{
    int t=1;
    while(cin>>(str+1),strcmp(str+1,"END"))//strcmp相等时返回0
    {
        int len=strlen(str+1);
        for(int i=len*2;i;i-=2)		//马拉车算法处理,将字符都变为奇数情况
        {
            str[i]=str[i/2];
            str[i-1]='z'+1;
        }
        len*=2;
        p[0]=1;
        for(int i=1,j=len;i<=len;i++,j--)
        {
            hl[i]=hl[i-1]*P+str[i]-'a'+1;
            hr[i]=hr[i-1]*P+str[j]-'a'+1; //反求后端前缀和
            p[i]=p[i-1]*P;
        }
        int res=0;
        for(int i=1;i<=len;i++)   //选取中间点
        {
            int l=0,r=min(i-1,len-i);
            while(l<r)      //二分选取满足的子串半径
            {
                int mid=l+r+1>>1;
                if(get(hl,i-mid,i-1)!=get(hr,len-(i+mid)+1,len-(i+1)+1))
                    r=mid-1;
                else
                    l=mid;
            }
            if(str[i+l]<='z')		//两端由字符包围说明字符多一个
                res=max(res,l+1);
            else					//两端由非字母包围说明非字符多一个
                res=max(res,l);
        }
        printf("Case %d: %d\n",t++,res);
    }
}