[ural1297]最长回文子串

我也没有题目啦
就是说给一个字符串求出字串中式回文的最长的字串。。好乱啦。

给定一个字符串，求最长回文子串。
算法分析：
穷举每一位，然后计算以这个字符为中心的最长回文子串。注意这里要分两种情况，一是回文子串的长度为奇数，二是长度为偶数。两种情况都可以转化为求一个后缀和一个反过来写的后缀的最长公共前缀。具体的做法是：将整个字符串反过来写在原字符串后面，中间用一个特殊的字符隔开。这样就把问题变为了求这个新的字符串的某两个后缀的最长公共前缀。
这个做法的时间复杂度为O(nlogn)。如果RMQ问题用时间为O(n)的方

#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
char a[210000];
int wr[210000],mc[210000],JS[210000],sa[210000],y[210000],height[210000];
bool cmp(int k1,int k2,int ln){
    return wr[k1]==wr[k2]&&wr[k1+ln]==wr[k2+ln];
}
void get_sa(int n,int m){//构建SA后缀数组
    int i,k,p,ln;
    //memcpy(mc,a,sizeof(a)); 
    for(i=1;i<=n;i++)mc[i]=a[i];
    //a数组:原字符串，mc名次数组
    for(i=0;i<=m;i++)JS[i]=0; 
    for(i=1;i<=n;i++)JS[mc[i]]++;
    for(i=1;i<=m;i++)JS[i]+=JS[i-1];
    for(i=n;i>=1;i--)sa[JS[mc[i]]--]=i;  
    //以上四句为基数排序，不懂的看flash
    ln=1;p=0; 
    //ln为当前子串的长度,p表示有多少不相同的子串
    while(p<n){
        for(k=0,i=n-ln+1;i<=n;i++)y[++k]=i;
        for(i=1;i<=n;i++)if(sa[i]-ln>0)y[++k]=sa[i]-ln;
        for(i=1;i<=n;i++)wr[i]=mc[y[i]];
        //数组y保存的是对第二关键字排序的结果 。
        //数组wr保存的是对第二关键字排序后的mc值
        //以下为对第一关键字排序
        for(i=0;i<=m;i++)JS[i]=0;
        for(i=1;i<=n;i++)JS[wr[i]]++;
        for(i=1;i<=m;i++)JS[i]+=JS[i-1];
        for(i=n;i>=1;i--)sa[JS[wr[i]]--]=y[i];
        memcpy(wr,mc,sizeof(wr));   
        p=1;mc[sa[1]]=1;
        for(i=2;i<=n;i++){
            if(!cmp(sa[i],sa[i-1],ln))p++;
            mc[sa[i]]=p;
        }
        //得到新的mc数组
        m=p;ln*=2;
    }
    a[0]=0;sa[0]=0;
}
void get_he(int n){
    int i,j,k=0;
    for(i=1;i<=n;i++){
        j=sa[mc[i]-1];
        if(k)k--;
        while(a[j+k]==a[i+k]) k++;
        height[mc[i]]=k;
    }
}
int main(){
    scanf("%s",a+1);
    int n=strlen(a+1);
    a[n+1]='^';
    for(int i=n+2;i<=n+n+1;i++){
        a[i]=a[n-(i-n-2)];
    }
    get_sa(n+n+1,300);
    get_he(n+n+1);
    int maxx=0;
    for(int i=1;i<=n+n;i++){
        maxx=max(maxx,height[i]);
    }
    for(int i=1;i<=n-maxx+1;i++){
        bool bk=true;
        for(int j=1;j<=maxx;j++){
            if(a[i+j-1]!=a[i+maxx-j]){
                bk=false;break; 
            }
        }
        if(bk==true){
            for(int j=1;j<=maxx;j++){
                printf("%c",a[i+j-1]);
            }   
        }
    }
    return 0;
}