「POJ-3974」Palindrome-manacher

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
const int MAXN = 2010000;
int p[MAXN];
char s1[MAXN], s2[MAXN];
inline int manacher(const int len, const char *s2, int *p) {
    register int ans = 0, max = -1, id = -1;
    for (register int i = 1; i <= len; i++) {
        if (id + max >= i) p[i] = std::min(p[id * 2 - i], id + max - i);
        while (i - p[i] - 1 >= 0 && i + p[i] + 1 <= len && s2[i - p[i] - 1] == s2[i + p[i] + 1]) p[i]++;
        if (i + p[i] > id + max) max = p[i], id = i;
        ans = std::max(ans, p[i]);
    }
    return ans;
}
inline int solve(const int len, const char *s1, char *s2) {
    for (register int i = 0; i < len; i++) s2[i << 1 | 1] = s1[i], s2[(i << 1) + 2] = '#';
    return s2[0] = '#', s2[len << 1 | 1] = '\0', memset(p, 0, sizeof(p)), manacher(len << 1 | 1, s2, p);
}
int main() {
    register int cnt = 0;
    while (scanf("%s", s1)) {
        if (strcmp(s1, "END") == 0)break;
        printf("Case %d: %d\n", ++cnt, solve(strlen(s1), s1, s2));
    }
}