最長迴文子串——Manacher算法

#include <iostream>
#include <stdio.h>
#include <string>
#include <vector>
using namespace std;
char c[2000010];
void Insert(string s){
    c[0]='@';
    for(int j=0; j<s.length(); j++){
        c[2*(j+1)]= '@';
        c[2*j+1]=s[j];
    }
}

int min(int a, int b){
    return a>b?b:a;
}

int main()
{
    int n;
    int mx;
    int p;
    int max;
    string s;
    vector<int> v;
    cin>>n;
    for(int j=0; j<n; j++){
        cin>>s;
        Insert(s);
        max=0;
        mx=0;
        p=0;
        v.clear();
        int l=2*s.length()+1;
        for(int i=0; i<l; i++){
            //cout<<max<<endl;
            if(mx>i){
                v.push_back(min(mx-i, v[p*2-i]));
            }else{
                v.push_back(1);
            }
            while(c[i+v[i]]==c[i-v[i]] && (i+v[i])<l && (i-v[i])>=0) v[i]++;
            if((v[i]+i)>mx) {
                mx=v[i]+i;
                p=i;
            }
            if(v[i]>max) max=v[i];
        }
        cout<<max-1<<endl;
    }
    return 0;
}

hihoCoder的1032道題的解題代碼

附上最易理解的

算法講解講解

http://m.blog.csdn.net/article/details?id=42061017

以及我理解的Manacher算法:

該算法的思想是從中心向兩邊探索迴文,並標記之前探索的迴文右邊最遠的字符的下標

1)如果是這種情況:

0  1  2  3  4  5  6  7  8  9  10

c  b  c  d   c  b  c   c  e  c    a

比如現在我要探索的是b(下標5),之前探索的以d(下標3)文中心的最長迴文是c  b  c  d  c  b  c,右邊最遠字符是c(下標6),以d爲中心最有兩邊是對稱的,b(下標5)的對稱點是b(下標1),右邊最遠字符c(下標6)關於d的對稱點是下標0的字符c,則b(下標1)到c(下標0)這部分的對應字符是b(下標5)到c(下標0),這部分字符前面已經探索過了,可以跳過匹配過的,從c(下標7)開始探索,因爲這部分是未探索過的。

2)如果是這種情況:

0  1  2  3  4  5  6  7  8  9  10

a  b  c  d  c  b  a   c   e   c    a

比如現在我要探索的是c(下標7),之前探索的以d(下標3)文中心的最長迴文是c  b  c  d  c  b  c,右邊最遠字符是c(下標6),c(下標7)是一個新的字符,自己本身都沒有匹配過,因此只能老老實實的一個一個的兩邊匹配了。

以上是我的理解。

点赞