最长回文子串

问题描述

从标准输入读取一个整数N（N<=30)，代表字符串的个数，接下来是N个字符串（字符串长度<=10^6)

输出每个字符串的最长回文子串长度。 

输入：

3
abababa
aaaabaa
acacdas

输出

7
5
3

解题方法：manacher算法

• 不妨假设给定的字符串为S。
• 统一回文长度为奇数、偶数的情况的预处理。处理过程不赘述（abba）---> （$#a#b#b#a#）
• 数组 P[i] 来记录以字符S[i]为中心的最长回文子串向左/右扩张的长度（包括S[i]，也就是把该回文串“对折”以后的长度），比如S和P的对应关系：（巧妙的是P[i]-1正好是原字符串中回文串的总长度）

  S  #  1  #  2  #  2  #  1  #  2  #  3  #  2  #  1  #
  P  1  2  1  2  5  2  1  4  1  2  1  6  1  2  1  2  1

  计算P[i]，从左向右，这里需要引入以S[i]之前字符S[id]为中心的回文所能到达的最远字符的编号mostFar。 在计算P[i]时，为避免重复计算，分以下情况讨论：
• 1. 当mostFar>i 时（j为i关于id 的镜像位置），若P[j]<mostFar-i，则P[i] = P[j].  否则，P[i]=mostFar-i.(这时，P[i]的值可能大于mostFar-i，因此后面还得继续找)；
• 2. 当mostFar<=i 时，我们就无法从前面的P值推断P[i]，所以P[i]=1，后面再继续查找。
• 计算玩P数组后，max(P[0],...P[len(S')-1])-1即为所求。

代码

#include <iostream>
#include <string>
#include<algorithm>
using namespace std;
// 将字符串处理，比如abcd变成：$#a#b#c#d#. 
// 这样可以将奇数长度的回文和偶数长度回文一起处理。 
string preProcess(string & str) {
<span style="white-space:pre">	</span>string ret="$#";
	for (int i=0; i<str.length(); i++) {
		ret += str[i];
		ret += "#";
	} 
	return ret;
}
int p[2000009]={0};
void getMaxLength(string str) {
	// mostFar记录了i之前回文到达的最右坐标
	// id是与之对应的中心坐标。 p记录的是以i为中心的回文半价，单个字母为1. 
	int id=0, mostFar=0, maxL=0;
	p[0] = 0;
	for (int i=1; i<str.length(); i++) {
		if (mostFar > i) {
			int j=2*id-i; // j and i are symmetric by id. 
			if (p[j] < mostFar-i) // extension needed
				p[i] = p[j];
			else
		        p[i] = mostFar-i; // extension needed 
		} else 
			p[i] = 1;
		// extension
		while ((i+p[i]<str.length()) && (i-p[i]>=0) && str.at(i+p[i]) == str.at(i-p[i]))
			++p[i];
		if (p[i]+i > mostFar) { // update mostFar and id.
			mostFar = (p[i]+i);
			id = i;
		}
	}
	for (int i=0; i<str.length(); i++)
		maxL = max(maxL, p[i]);
	cout<<maxL-1<<endl;
} 
int main() {
	int n=0;
	cin>>n;
	string str="";
	for (int i=0; i<n; i++) {
		cin>>str;
		if (str.length() < 2)
			cout<<str.length()<<endl;
		else {
			str = preProcess(str);
			getMaxLength(str);
		}
	}
	return 0;
}

原文：http://www.felix021.com/blog/read.php?2040