USACO 6.2.1 Calf Flac 马拉车算法（Manacher's Algorithm）

http://train.usaco.org/usacoprob2?a=Hfv8bxVT1Ac&S=calfflac

题目大意：一串字符串（不超过200000个字符），求最长回文子串（不考虑非字母字符），并输出这一段（要包含非字母字符，包括换行符）

裸的最大回文子串题，直接用马拉车算法可解：

※马拉车算法讲解可参考：https://www.cnblogs.com/grandyang/p/4475985.html

    但是这篇博客里有些小问题：

        ①变量id指的应该是最右成匹边界对应的字符串中心，而不是最长回文子串的中心

        ②匹配半径P[]时而包含中心（初始为1），时而不包含中心（mx-i），自己写的时候需要统一一下

    

至于非字母字符的问题，最简单粗暴的方法就是直接先把原始字符串过一遍记录下除字符后位置与原始位置的对应关系，顺便把大写字母转变成对应的小写字母

/*
ID: frontie1
TASK: calfflac
LANG: C++
*/
#include <iostream>
#include <cstdio>
#include <string>

using namespace std;

string origin;
char tem[82];
int pos[200000];
int e = 0;
char arr[400010];
int P[400010];

bool ischaracter(char chr)
{
    if((chr >= 'a' && chr <= 'z') || (chr >= 'A' && chr <= 'Z')) return true;
    return false;
}

int Manacher()
{
    int length = 2*e+1;
    int id = 0, mx = 0;
    int outlen = 0, output = 0;
    for(int i = 1; i <= length; ++i){
        P[i] = (mx > i) ? min(mx-i, P[2*id-i]) : 0;
        while(arr[i+P[i]+1] == arr[i-P[i]-1]) ++P[i];
        if(i+P[i] > mx){
            mx = i+P[i];
            id = i;
        }
        if(P[i] > outlen){
            output = i;
            outlen = P[i];
        }
    }
    return output;
}

int main()
{
    freopen("calfflac.in", "r", stdin);
    freopen("calfflac.out", "w", stdout);

    while(cin.getline(tem, 82, '\n') && !cin.eof()){
        origin += tem;
        origin += "\n";
    }

    arr[0] = '$'; arr[1] = '#';

    for(int i = 0; i < origin.length(); ++i){
        if(ischaracter(origin[i])){
            pos[e] = i;
            arr[2*e+2] = (origin[i] <= 'Z' ? 'a'+origin[i]-'A' : origin[i]);
            arr[2*e+3] = '#';
            ++e;
        }
    }
//
//    for(int i = 0; i < e; ++i){
//        cout << pos[i] << ' ';
//        if(i % 80 == 0) cout << endl;
//    }

    arr[2*e+2] = '$';

    int centre = Manacher();

    int st = (centre-P[centre]-1)/2;
    int ed = (centre+P[centre]-3)/2;

    //cout << origin << endl;

    cout << P[centre] << endl;
    for(int i = pos[st]; i <= pos[ed]; ++i){
        cout << origin[i];
    }
    cout << endl;

    return 0;
}

附录：这题在USACO上的测试样例也是很走心hhh，实在不行的话可以先过前两个，然后用第三个的代码试试看（笑）