最长回文子串【Manacher】

Manacher

很巧妙地回文串匹配方法。

这个算法就是通过利用现在已经知道的回文串匹配情况进行之后的回文串判断

巧妙跳过了很多曾经匹配过的子串

看了两篇博文，帮助很大~~~

https://www.felix021.com/blog/read.php?2040

http://blog.csdn.net/ggggiqnypgjg/article/details/6645824

本方法中的一些定义：

s：输入的原始字符串数组

str：更新s之后的数组

p[i]，以i为中心的最长回文子串.

而且呢，神奇的是：原字符串中的最长回文子串长度就是p[i]-1

（因为p[i]时是加上了#的半径，而原字符串里面有一半是#嘛，再把算了两边的重点去掉，所以是p[i]-1）

记录最长回文子串：

id：目前匹配的最长回文串的中心位置

mx：目前匹配的最长回文串的扩展到右侧的边界坐标，也就是中心位置+半径。

因为不知道回文串是奇数还是偶数，所以需要实现处理一下，在每个字母中间加上#，并在开头，即下标为0的位置加上$

这样也就保证了字符串在数组中从下标1开始，并且是偶数个

a b c b c==> $ # a # b # c # b # c #

最重要的是这句：

p[i]=min{p[2*id-i],  mx-i };

id-(i-id)是i关于id的对称点。

mx-i就是i根据已有的匹配结果所能匹配的长度。也就是p[id]+id-i

附上前面提到的博文中的图片

hduoj3068

最长回文

Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 14875    Accepted Submission(s): 5446

Problem Description

给出一个只由小写英文字符a,b,c...y,z组成的字符串S,求S中最长回文串的长度.
回文就是正反读都是一样的字符串,如aba, abba等

 

Input

输入有多组case,不超过120组,每组输入为一行小写英文字符a,b,c...y,z组成的字符串S
两组case之间由空行隔开(该空行不用处理)
字符串长度len <= 110000

 

Output

每一行一个整数x,对应一组case,表示该组case的字符串中所包含的最长回文长度.

 

Sample Input

  
  

   
   aaaa

abab
  
  

 

Sample Output

  
  

   
   4
3
  
  

 

代码：

#include <cstdio>
#include <cstring> 
#define maxn 220010
#define min(x,y) x<y?x:y
using namespace std;
int n,p[maxn];
char s[maxn],str[maxn];

void init(){
	str[0]='$';
	str[1]='#'; 
	for(int i=0;i<n;i++){		//必须从0开始遍历 
		str[2*i+2]=s[i];
		str[2*i+3]='#';
	} 
	n=n*2+2;
	str[n]='\0';
	return; 
}


void kp(){
	
	int mx=0,id=0; 
	//循环
	//1.根据现有情况更新p[i]
	//2.根据后面字符串情况更新p[i] 
	//3.照最长 
	for(int i=1;i<n;i++){
		//if(p[i]<mx-i)
		if(mx>i)
			p[i]=min(p[2*id-i],mx-i);
		else
			p[i]=1;
		
		while(str[i+p[i]]==str[i-p[i]])	p[i]++; 
		
		//更新最长 
		if(p[i]+i>mx){
			mx=p[i]+i;
			id=i;
		}
	}	
}


int main(){
	while(scanf("%s",s)!=EOF){
		n=strlen(s);
		init();
		kp();
		int ans=0;		
		for(int i=0;i<n;i++){
			if(ans<p[i]) ans=p[i];
		} 
		printf("%d\n",ans-1);
	}
	return 0;
} 

UVA10817