题目链接:POJ3974【Palindrome】
• 题目描述:
• Andy是一个聪明夫人的计算机专业的学生,他正在上算法课,教授问学生一个简单的问题:“你们能不能提出一个有效的算法,来找出一个字符串中最长回文的长度?”
• 如果一个字符串向前读和向后读是相同的,则称其为回文,例如"madam"是回文,而"acm"则不是回文。
• 同学们认识到这是一个经典的问题,但是他们无法找到一个比遍历所有子串并检查它们是否是回文更好的解决方案,但显然这样的算法根本没有效率。过了一会儿,Andy举手说:“OK!我有一个更好的算法。”在Andy开始解释他的想法之前,停了一会儿,然后说:“好吧,我有一个更好的算法!”
• 如果您认为您知道Andy的最终解决方案,就证明它。给定一个最多1000000个字符的字符串,请您查找并输出该字符串中最长回文的长度。
• 输入
• 您的程序将对最多达30个的测试用例进行测试,每个测试用例在一行中以最多1000000个小写字符的字符串形式给出。输入以字符串"END"开头的一行结束(为了清楚起见,用引号)。
• 输出
• 对于输入中的每个测试用例,输出测试用例编号和最长回文的长度。
输入样例:
abcbabcbabcba
abacacbaaaab
END
输出样例:
Case 1: 13
Case 2: 6
题意:
输入多行字符串,输入END代表结束
要求输出每行字符串内最大回文的长度。
数据范围:
最多30行,每行最多1000000个字符
思路:
马拉车模板题
找到p数组求大值
最大值减一就是答案
马拉车具体讲解可以参考这篇博客Manacher(马拉车)算法总结
代码:
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
using namespace std;
const int maxn = 1000000+5;
char s[maxn],s_new[maxn*2];
int p[maxn*2],ans,len;
int init()
{
len = strlen(s);
s_new[0] = '$';
s_new[1] = '#';
int m = 2;
for(int i = 0;i < len;i++)
{
s_new[m++] = s[i];
s_new[m++] = '#';
}
s_new[m++] = '\0';
return m;
}
void Manacher()
{
int l = init();
memset(p,0,sizeof(p));
int mx = 0,di,ans = 0;
for(int i = 0;i < l;i++)
{
p[i] = mx > i ? min(p[2 * di - i], mx - i) : 1;
while(s_new[i-p[i]] == s_new[i+p[i]]) p[i]++;
if(i+p[i] > mx)
{
mx = i+p[i];
di = i;
ans = max(ans,p[i]);
}
}
printf("%d\n",ans-1);
}
int main()
{
int k = 1;
while(~scanf("%s",s))
{
if(!strcmp(s,"END")) break;
printf("Case %d: ",k++);
Manacher();
}
return 0;
}
注意:
新的s_new和p数组的长度要是s的两倍
新字符串的头尾要记得加上特殊字符串