manacher算法第一次使用,看了半天才勉强理解算法的原理,真心觉得发明此算法的是大牛,佩服佩服。
找了一个链接,对于manacher算法有一个很好的介绍:
http://www.cnblogs.com/xiaoningmeng/p/5861154.html
看完算法就看看题目吧,这是HDOJ3294题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3294
比较裸的manacher算法应用,最后再做好字符串的转换即可。
代码如下:
#include <stdio.h>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int maxn = 2e6+10;
char str[maxn];//初始字符串
char s[maxn<<1];//构造的新字符串
int p[maxn<<1]; //以s[i]为中心的最长回文子串的最右端到s[i]的位置的长度
int pos;//记录回文字符串中心位置
int mannacher(char *ss){
int len=strlen(ss);
int t=0;
s[t++]='@';//防止数组越界
//构造新字符串
s[t++]='#';
for (int i=0; i<len; i++){
s[t++]=ss[i];
s[t++]='#';
}
s[t]=0;
int mx=0,id=0,ans=0;
for (int i=0;i<t;i++){
if (mx>i) p[i]=min(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;
}
if(p[i]-1>ans){
ans=p[i]-1;
pos=i;
}
}
return ans; //最长回文串长度
}
int main(){
char c,a[2];
while(~scanf("%s %s",a,str)){
c=a[0];
int len=mannacher(str);//字符串s的最长回文子串长度
if (len<2) printf("No solution!\n");
else {
int start=(pos-len+1)/2-1; //最大回文子串开始的位置
int e=start+len-1; //最大回文子串结束的位置
printf("%d %d\n",start,e);
int diff=c-'a';
for (int i=start;i<=e;i++)
(str[i]-'a'>=diff)?printf("%c",str[i]-diff):printf("%c",str[i]+26-diff);
printf("\n");
}
}
return 0;
}