L2-008. 最长对称子串----暴力+dp+Manacher

L2-008. 最长对称子串

时间限制

100 ms

内存限制

65536 kB

代码长度限制

8000 B

判题程序

Standard

作者

陈越

对给定的字符串，本题要求你输出最长对称子串的长度。例如，给定"Is PAT&TAP symmetric?"，最长对称子串为"s PAT&TAP s"，于是你应该输出11。

输入格式：

输入在一行中给出长度不超过1000的非空字符串。

输出格式：

在一行中输出最长对称子串的长度。

输入样例：

Is PAT&TAP symmetric?

输出样例：

11

题目链接：https://www.patest.cn/contests/gplt/L2-008

去年做过这个题，明明记得我是反转后用最长公共子串做的，结果今年就不对了，加测试数据了？

可以暴力，以每个i为对称轴，然后暴力。

暴力代码：

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
int main(){
    string s;
    getline(cin,s);
    int len=s.length();
    int mix=0;
    for(int i=0;i<len;i++){
        int x=1;
        int r=i+1;
        int l=i-1;
        while(r<len&&l>=0){//回文串长度为奇数
            if(s[r]==s[l]){
                r++;
                l--;
                x+=2;
            }
            else{
                break;
            }
        }
        mix=max(mix,x);
        r=i+1;
        l=i;
        x=0;
        while(r<len&&l>=0){//回文串长度为偶数
            if(s[r]==s[l]){
                x+=2;
                r++;
                l--;
            }
            else{
                break;
            }
        }
        mix=max(mix,x);
    }
    cout<<mix<<endl;
    return 0;
}

上面就是个水代码，祥巨告诉我可以dp，我就强行dp了一波，发现还是暴力，而且，怎么越看和最长公共子串就越像呢？

dp代码，未验证正确性：

for(int i=0;i<n;i++){
    int j=i;
    while(j>=0){
        if(s[j]==s[i]&&(i-j<2||dp[j+1][i-1])){
            dp[j][i]=true;
            Max=max(Max,i-j+1);
        }
        j--;
    }
}

接下来就是Manacher算法了，今晚刚学的，祥巨告诉我的，然而真的很好使，但是他的时间复杂度我是真的没看懂。

先看

http://www.cnblogs.com/biyeymyhjob/archive/2012/10/04/2711527.html
再看进阶：

https://segmentfault.com/a/1190000003914228

代码：

#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
using namespace std;
int p[10000];
void Manacher(string s){
    memset(p,0,sizeof(p));
    int mx=0,id=0;
    for(int i=1;i<=(int)s.size();i++){
        if(mx>i){
            p[i]=(p[2*id-i])<(mx-i)?p[2*id-i]:(mx-i);
        }
        else{
            p[i]=1;
        }
        while(s[i-p[i]]==s[i+p[i]]){
            p[i]++;
        }
        if(i+p[i]>mx){
            mx=i+p[i];
            id=i;
        }
    }
    int Max=0,ii;
    for(int i=1;i<(int)s.size();i++){
        if(p[i]>Max){
            ii=i;
            Max=p[i];
        }
    }
    Max--;
    int ans=0;
    int start=ii-Max;
    int End=ii+Max;
    for(int i=start;i<=End;i++){
        if(s[i]!='#'){
            ans++;
        }
    }
    cout<<ans<<endl;
}
int main(){
    string s;
    getline(cin,s);
    int len=s.length();
    string str;
    str+="$#";
    for(int i=0;i<len;i++){
        str+=s[i];
        str+="#";
    }
    Manacher(str);
    return 0;
}