「Luogu P1210」回文检测 解题报告

题面

这是一道诡异的黄题

居然让你求一串吧啦吧啦的东西中 字母（大小写）最长的回文串的长度，还要输出完整的串 吐血

思路：

保持淡定，我们啥都不会，就会Manacher，那就用Manacher大法！

1、把字符串处理成只有字母（去皮）

2、把字符串中插入#号（掺假）

3、Manacher跑一跑最长回文串长度（市场检查）

4、求出最长长度后再重新代入原来的字符串（打回原形！妖孽！）

然后，就完了~

就是去皮的时候，顺便记录一下字符的位置，这样重新代入时就比较方便

并且由于字符串有很多行，读入时就要用getchar()，直到EOF为止

似乎Noip2018普及T1有人因为这个WA？

Code：

#include<bits/stdc++.h>
#define M 20010
using namespace std;
struct node{
    char c;
    int id;
}a[M<<1];
char st[M],s[M<<1];
int l1,l2,len,ans,res;
int p[M<<1];
bool check(char c)
{
    return (c>='a'&&c<='z')||(c>='A'&&c<='Z');
}
void init()
{
    int i;
    for(i=0;i<l1;i++)
        if(check(st[i]))
        {
            a[l2].c=(st[i]>='a'&&st[i]<='z')? st[i]:st[i]-'A'+'a';
            a[l2++].id=i;//记录编号（位置）
        }
    s[0]='!';
    s[1]='#';
    for(i=0;i<l2;i++)
    {
        s[i*2+2]=a[i].c;
        s[i*2+3]='#';
    }
    len=l2*2+2;
    s[len]='?';
    return;
}
void manacher()//大法不解释
{
    int id=0,mx=0;
    for(int i=0;i<len;i++)
    {
        if(i<mx)
            p[i]=min(p[id*2-i],mx-i);
        else
            p[i]=1;
        while(s[i-p[i]]==s[i+p[i]])
                p[i]++;
        if(i+p[i]>mx)
        {
            id=i;
            mx=i+p[i];
        }
    }
    return;
}
void Print(int l,int r)
{
    if(l&1)//由扩展后的s数组变回a数组
        l=(l-3)/2;
    else
        l=(l-2)/2;
    if(r&1)
        r=(r-3)/2;
    else
        r=(r-2)/2;
    for(int k=a[l].id;k<=a[r].id;k++)
        printf("%c",st[k]);
    return;
}
int main()
{
    int i;
    char c=getchar();//读入
    while(c!=EOF)
    {
        st[l1++]=c;
        c=getchar();
    }
    init();//初始化
    manacher();//大法
    for(i=0;i<len;i++)//求最长
        ans=max(ans,p[i]);
    for(i=0;i<len;i++)//求第一次出现的位置
        if(p[i]==ans)
            break;
    ans--;
    printf("%d\n",ans);
    Print(i-ans+1,i+ans-1);//重新代入
    return 0;
}