51Nod 1088 最长回文子串 【Manacher算法】带模板详解

51Nod 1088 最长回文子串

输入一个字符串Str，输出Str里最长回文子串的长度。

回文串：指aba、abba、cccbccc、aaaa这种左右对称的字符串。

串的子串：一个串的子串指此（字符）串中连续的一部分字符构成的子（字符）串
例如 abc 这个串的子串：空串、a、b、c、ab、ac、bc、abc

 收起

输入

输入Str（Str的长度 <= 1000)

输出

输出最长回文子串的长度L。

输入样例

daabaac

输出样例

5

提到最长回文串就肯定是Manacher算法没问题了  这个算法刚开始我也有点迷 后来才领会了其中真意 咱们一步一步来看

首先  是该算法对字符串的处理   在原字符串的每两个字符串中加入分隔符 这样做有什么用呢  就是 不必在考虑回文串到底是奇数串还是偶数串   奇数串有回文中心   例如 aba 中 b就是回文中心  可以以b为中心向两边扩散  而偶数串 例如 abba 就没有回文中心了 或者说是以 bb 两个字母为回文中心  这样奇数串和偶数串还得分别处理  无论是时间还是代码量都会增长  Manacher算法第一步就是插入分隔符    aba 成为 $#a#b#a# 而 abba $#a#b#b#a# 这样做无论是  奇数串还是偶数串都变成了奇数串 

	str+='$';
	str+='#';
	for(int i = 0; i < s.length();i++){
		str+=s[i];
		str+='#';
	}

接着 就是匹配了   这一点有点复杂 我们就先用暴力算法 然后慢慢改进成Manacher算法 匹配很简单 就是先找到回文中心然后向两边扩散

for(int i = 0; i < s.length;i++){
    int len = 1;
    while(s[i]-len == s[i]+len){
        len++;
    }
}

这样我们就能得到以每个字符为中心的回文串长度 但是时间复杂度是O(n^3) 就很复杂了 那么 怎样才能降低时间复杂度呢？Manacher算法的精髓就在这里 只有没匹配过的字符才进行匹配 匹配过的就不再匹配  理论上 时间复杂度就为O(n) 那么 Manacher算法是怎样做到的呢 

第一 把已经匹配过的字符对应的回文子串长度记录下来  这样再用到的时候无需再匹配 直接查表

第二 如果新串的回文中心(i)是在旧串(p)的范围之内的话 当匹配以i为中心 时 存在 以j为中心的子串（i+j = 2p） 关于 p 对称 已经 匹配好 所以len[j] 就可以在len[i]的基础上扩充  如果不在的话 只能老老实实匹配了  完了我们还要把新匹配的串更新为子串p方便后面的回文中心的匹配  由于我们插入了分隔符 所以记录子串的半径就是原来子串的长度+1 毕竟我们匹配的时候会多混进来一个#

这样我们就完成字符串的匹配了  当然每次得到一个回文串的长度的时候都和ans取最大值更新ans  因为我们的得到

这样Manacher算法就算完成了  

ps: 第一次写的向外扩张用的while循环  后来做 1089 TLE 了 改成for以后 就AC了  个人是这样理解的  while 循环做++运算 每次都要进行操作 如果用for 运算属于循环一部分 不知道对不对

void Manacher(string s){
	//插入分隔符
	string str;
	str+='$';
	str+='#';
	for(int i = 0; i < s.length();i++){
		str+=s[i];
		str+='#';
	}
	
	int MaxR = 0; // 最大回文半径
	int po = 0;   //回文串的中间字符
	int ans = 0;
	//Manacher匹配 
	for(int i = 1; i <= str.length();i++){
		
		//新的中心在 之前的最大半径之内 匹配对应段 
		if(MaxR>i){
			Len[i] = min(MaxR-i,Len[2*po-i]);	
		}else{
			Len[i] = 1;
		}
		//以po为中心开始向两边匹配 
        for(; str[i+Len[i]] == str[i-Len[i]]; Len[i]++);
		
		//若新串长度大于旧串长度 更新最大长度 更新回文串中心 
		if(Len[i]+i > MaxR){
			MaxR = Len[i]+i;
			po = i; 
		}
		ans=max(ans,Len[i]);
	} 
	
	//输出最长字符串长度
	printf("%d\n",ans-1);
}

然后我们的题目就直接套模板就可以了

#include<iostream>
#include<cmath>
#include<string> 

using namespace std;

const int N = 1000+5;

int Len[N];

void Manacher(string s){
	//插入分隔符
	string str;
	str+='$';
	str+='#';
	for(int i = 0; i < s.length();i++){
		str+=s[i];
		str+='#';
	}
	
	int MaxR = 0; // 最大回文半径
	int po = 0;   //回文串的中间字符
	int ans = 0;
	//Manacher匹配 
	for(int i = 1; i <= str.length();i++){
		
		//新的中心在 之前的最大半径之内 匹配对应段 
		if(MaxR>i){
			Len[i] = min(MaxR-i,Len[2*po-i]);	
		}else{
			Len[i] = 1;
		}
		//以po为中心开始向两边匹配 
        for(; str[i+Len[i]] == str[i-Len[i]]; Len[i]++);
		
		//若新串长度大于旧串长度 更新最大长度 更新回文串中心 
		if(Len[i]+i > MaxR){
			MaxR = Len[i]+i;
			po = i; 
		}
		ans=max(ans,Len[i]);
	} 
	
	//输出最长字符串长度
	printf("%d\n",ans-1);
}
int main()
{
	string s;
	cin >> s;
	Manacher(s);
}