最长回文子串（后缀数组）

这个也可以利用manacher算法来做

奇数长度的回文串是 从i开始的后缀 与 反转后字符串中从原来的i位子开始的后缀 的公共前缀 
偶数长度类似 
所以和最长公共子串类似的思路, 将原串S与反转后的S中间加入一个不会出现的字符拼接在一起, 对得到的字符串S’求高度数组 
以原来S[i]为中心的奇数长度的回文串长度 = 后缀S’[i…]和后缀 S’[2n-i…]的公共前缀 
以原来S[i-1]和S[i]为中心的回文串长度 = 后缀S’[i…]和后缀S’[2n-i+1…]的公共前缀

由于我这个程序高度数组可能某一个位置上可能为0，也就是排名最后的后缀，但是这样数据是aabaa会报错，进行更改后，现在如果数据是abcd会报成0。

#include<bits/stdc++.h>

using namespace std;

const int maxn=2e5+100;

char s[maxn];
int rk[maxn],tmp[maxn],sa[maxn],n,k,lcp[maxn];

bool cmp_sa(int i,int j)
{
    if(rk[i]!=rk[j]){
        return rk[i]<rk[j];
    }else{
        int ri=i+k<=n?rk[i+k]:-1;
        int rj=j+k<=n?rk[j+k]:-1;
        return ri<rj;
    }
}

void construct_sa()
{
    for(int i=0;i<=n;i++){
        sa[i]=i;
        rk[i]=i<n?s[i]:-1;
    }
    for(k=1;k<=n;k<<=1){
        sort(sa,sa+n+1,cmp_sa);
        tmp[sa[0]]=0;
        for(int i=1;i<=n;i++){
            tmp[sa[i]]=tmp[sa[i-1]]+(cmp_sa(sa[i-1],sa[i])?1:0);
        }
        copy(tmp,tmp+n+1,rk);
    }
}

void construct_lcp()
{
    for(int i=0;i<=n;i++){
        rk[sa[i]]=i;
    }
    int h=0;
    lcp[0]=0;
    for(int i=0;i<n;i++){
        int j=sa[rk[i]-1];
        if(h>0){
            h--;
        }
        for(;j+h<n&&i+h<n;h++){
            if(s[j+h]!=s[i+h]){
                break;
            }
        }
        lcp[rk[i]-1]=h;
    }
    lcp[n]=lcp[n-1];
}
//ST表进行最值查询
int mi[maxn][30];
void init_rmq(int a[],int n)
{
    for(int i=0;i<n;i++){
        mi[i][0]=a[i];
    }
    for(int j=1;(1<<j)<=n;j++){
        for(int i=0;i+(1<<j)<=n;i++){
            mi[i][j]=min(mi[i][j-1],mi[i+(1<<(j-1))][j-1]);
        }
    }
}

int rmq(int l,int r)
{
    int k=0;
    while((1<<(k+1))<=r-l+1){
        k++;
    }
    return min(mi[l][k],mi[r-(1<<k)+1][k]);
}

int main()
{
    scanf("%s",s);
    int len=strlen(s);
    s[len]='$';
    reverse_copy(s,s+len,s+len+1);
    n=len*2+1;
    s[n]=0;
    for(int i=0;i<n;i++){
        cout<<s[i];
    }
    cout<<endl;
    construct_sa();
    construct_lcp();
    for(int i=0;i<=n;i++){
        cout<<sa[i]<<" ";
    }
    cout<<endl;
    for(int i=0;i<=n;i++){
        cout<<rk[i]<<" ";
    }
    cout<<endl;
    init_rmq(lcp,n+1);
    for(int i=0;i<=n;i++){
        cout<<i<<" ";
    }
    cout<<endl;
    for(int i=0;i<=n;i++){
        cout<<lcp[i]<<" ";
    }
    cout<<endl;
    int ans=0;
    if(len%2){
       for(int i=0;i<len;i++){
            int j=len*2-i;
            int l=rmq(min(rk[i],rk[j]),max(rk[i],rk[j]));
            cout<<i<<" "<<rk[i]<<" "<<j<<" "<<rk[j]<<" "<<l<<endl;
            ans=max(ans,2*l-1);
        }
    }else{
        for(int i=0;i<len;i++){
            int j=len*2-i+1;
            int l=rmq(min(rk[i],rk[j]),max(rk[i],rk[j]));
            cout<<i<<" "<<rk[i]<<" "<<j<<" "<<rk[j]<<" "<<l<<endl;
            ans=max(ans,2*l);
        }
    }

    cout<<ans<<endl;
    return 0;
}