https://vjudge.net/problem/URAL-1297
题意:
给出一个字符串求最长回文子串
代码:
//论文题,把字符串反过来复制一遍到后边,中间用一个没出现的字符隔开,然后就是枚举当以i位置为中间位时的最长回文串
//是多大,就是求这中间的heigh数组的最小值。用rmq预处理。枚举位置时分回文串长度是奇数还是偶数。
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int MAXN=;
int sa[MAXN+],ra[MAXN+],he[MAXN+],buc[MAXN+],xx[MAXN+],yy[MAXN+],f[MAXN+][];
int len,m;
char s[MAXN+];
void get_suf()
{
int *x=xx,*y=yy;
for(int i=;i<m;i++) buc[i]=;
for(int i=;i<len;i++) buc[x[i]=s[i]]++;
for(int i=;i<m;i++) buc[i]+=buc[i-];
for(int i=len-;i>=;i--) sa[--buc[x[i]]]=i;
for(int k=;k<=len;k<<=){
int p=;
for(int i=len-;i>=len-k;i--) y[p++]=i;
for(int i=;i<len;i++) if(sa[i]>=k) y[p++]=sa[i]-k;
for(int i=;i<m;i++) buc[i]=;
for(int i=;i<len;i++) buc[x[y[i]]]++;
for(int i=;i<m;i++) buc[i]+=buc[i-];
for(int i=len-;i>=;i--) sa[--buc[x[y[i]]]]=y[i];
swap(x,y);
p=;x[sa[]]=;
for(int i=;i<len;i++){
if(y[sa[i-]]==y[sa[i]]&&y[sa[i-]+k]==y[sa[i]+k])
x[sa[i]]=p-;
else x[sa[i]]=p++;
}
if(p>=len) break;
m=p;
}
for(int i=;i<len;i++) ra[sa[i]]=i;
int k=;
for(int i=;i<len;i++){
if(ra[i]==) { he[]=; continue; }
if(k) k--;
int j=sa[ra[i]-];
while(s[i+k]==s[j+k]&&i+k<len&&j+k<len) k++;
he[ra[i]]=k;
}
}
void rmq1(int n)
{
for(int i=;i<n;i++)
f[i][]=he[i];
for(int j=;(<<j)<n;j++){
for(int i=;i+(<<j)-<n;i++)
f[i][j]=min(f[i][j-],f[i+(<<(j-))][j-]);
}
}
int rmq2(int l,int r)
{
if(r<l) swap(r,l);
l++;
int k=;
while(<<(k+)<=r-l+) k++;
return min(f[l][k],f[r-(<<k)+][k]);
}
int main()
{
while(scanf("%s",s)==){
int x=strlen(s);
s[x]='';
len=x+;
for(int i=x-;i>=;i--) s[len++]=s[i];
m=;
get_suf();
rmq1(len);
int id=,max_len=;
for(int i=;i<x;i++){
int tmp;
if(i!=){
tmp=rmq2(ra[i],ra[*x-i+]);
if(tmp*>max_len){
max_len=tmp*;
id=i-tmp;
}
}
if(i!=x&&i!=){
tmp=rmq2(ra[i],ra[*x-i]);
if(tmp*->max_len){
max_len=tmp*-;
id=i-tmp+;
}
}
}
for(int i=id;i<id+max_len;i++) printf("%c",s[i]);
printf("\n");
}
return ;
}