139. 回文子串的最大长度（哈希）

如果一个字符串正着读和倒着读是一样的，则称它是回文的。

给定一个长度为N的字符串S，求他的最长回文子串的长度是多少。

输入格式
输入将包含最多30个测试用例，每个测试用例占一行，以最多1000000个小写字符的形式给出。

输入以一个以字符串“END”（不包括引号）开头的行表示输入终止。

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

每个输出占一行。

输入样例：
abcbabcbabcba
abacacbaaaab
END
输出样例：
Case 1: 13
Case 2: 6
难度： 中等
时/空限制： 1s / 64MB
总通过数： 398
总尝试数： 931
来源： 《算法竞赛进阶指南》

题意： 正方向反方向哈希，然后枚举回文中点，二分判断长度。

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

using namespace std;
typedef unsigned long long ll;
const int maxn = 1000000 + 7;
const int MOD = 131;
char s[maxn];
ll f1[maxn],f2[maxn],p[maxn];

ll H1(int i,int j)
{
    return f1[j] - f1[i - 1] * p[j - i + 1];
}

ll H2(int i,int j)
{
    return f2[i] - f2[j + 1] * p[j - i + 1];
}

int main()
{
    p[0] = 1;
    for(int i = 1;i <= maxn;++i)
    {
        p[i] = p[i - 1] * 131;
    }
    int kase = 0;
    while(1)
    {
        scanf("%s",s + 1);
        int n = (int)strlen(s + 1);
        if(s[1] == 'E' && s[2] == 'N' && s[3] == 'D')
        {
            break;
        }
        for(int i = 1;i <= n;i++)f1[i] = f1[i - 1] * 131 + s[i] - 'a' + 1;
        for(int i = n;i >= 1;i--)f2[i] = f2[i + 1] * 131 + s[i] - 'a' + 1;
        int ans = 0;
        for(int pos = 1;pos <= n;pos++)
        {
            int l = 0,r = min(pos - 1,n - pos),mid;
            while(l < r)
            {
                mid = (l + r + 1) / 2;
                if(H1(pos - mid,pos - 1) == H2(pos + 1,pos + mid))
                {
                    l = mid;
                }
                else
                    r = mid - 1;
            }
            ans = max(ans,l * 2 + 1);
            
            l = 0;r = min(pos - 1,n - pos + 1);
            while(l < r)
            {
                mid = (l + r + 1) / 2;
                if(H1(pos - mid,pos - 1) == H2(pos,pos + mid - 1))
                    l = mid;
                else
                    r = mid - 1;
            }
            ans = max(ans,l * 2);
        }
        printf("Case %d: ",++kase);
        cout << ans << endl;
    }
    return 0;
}

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
#include <iostream>
#include <map>
#include <string>

using namespace std;

typedef long long ll;
typedef unsigned long long ull;

const int maxn = 1e6 + 7;
const int mod = 99991;

char s[maxn];
ull f1[maxn],f2[maxn],p[maxn];

int n;

ull H1(int i,int j) {
    return f1[j] - f1[i - 1] * p[j - i + 1];
}

ull H2(int i,int j) {
    return f2[i] - f2[j + 1] * p[j - i + 1];
}

int main() {
    p[0] = 1;
    for(int i = 1;i < maxn;i++) {
        p[i] = p[i - 1] * mod;
    }
    int kase = 0;
    while(~scanf("%s",s + 1)) {
        if(s[1] == 'E') break;
        n = strlen(s + 1);
        memset(f1,0,sizeof(f1));
        memset(f2,0,sizeof(f2));
        for(int i = 1;i <= n;i++) {
            f1[i] = f1[i - 1] * mod + s[i];
        }
        for(int i = n;i >= 1;i--) {
            f2[i] = f2[i + 1] * mod + s[i];
        }
        
        int ans = 1;
        for(int i = 1;i <= n;i++) { //枚举中心点位置
            int l = 0,r = min(i - 1,n - i);
            int res = 0;
            while(l <= r) {
                int mid = (l + r) >> 1;
                if(H1(i - mid,i - 1) == H2(i + 1,i + mid)) {
                    res = mid;
                    l = mid + 1;
                }
                else r = mid - 1;
            }
            ans = max(ans,res * 2 + 1);
            
            l = 0,r = min(i,n - i);
            res = 0;
            while(l <= r) { //右半部分长度
                int mid = (l + r) >> 1;
                if(H1(i - mid + 1,i) == H2(i + 1,i + mid)) {
                    res = mid;
                    l = mid + 1;
                }
                else {
                    r = mid - 1;
                }
            }
            
            ans = max(ans,res * 2);
        }
        
        printf("Case %d: %d\n",++kase,ans);
    }
    return 0;
}