【最长回文子串】Manache算法,O(N)时间复杂度

问题描述:


找一个字符串里的最长回文子串



暴力法:定中心,从0长度向两端扩展的方法O(n^2), n >= 10^5还是超时,故只能《O(n^2)

Manacher's 算法:定中心,从p[r],(已能确定以该点为中心两端是回文的长度开始), 向两端扩展,时间复杂度, O(n)



算法如下:


示意图:

【最长回文子串】Manache算法,O(N)时间复杂度_第1张图片


######################

r是当前中心

cen的上一个回文串的中心

######################


定义:

cen:当 前匹配的上一个回文串的中心点
p[i]: 以i为中心两端满足回文的长度
r:从0到len-1遍历
l:r关于cen的对称点
绿色:表示已经确定满足回文的部分



目标: 求P数组


步骤:
0. r 遍历 判断是否小于cen+p[cen],是则下一步,否则p[r]=0, (在不在上一个回文串可达的内部,利用的就是已经匹配的回文串左右相等)
1. 确定p[r] 已经能确定的大小, (在上一个回文串可达的内部,利用的就是已经匹配的回文串左右相等)
    1> 蓝色情况:p[r] = cen+p[cen]-r, 左边的回文长度在边界外,=  r 到边界的距离
    2> 红色情况:p[r] = p[l],在当前cen+-p[cen]两端边界内
2. 在p[r]的基础上,while(str[r-p[r]]==str[r+p[r]])p[r]++,
    This is the key!!!!!!!p[r]不是都从0开始的
    如果是普通的方法,

    以r为中心,向两端扩展则时间复杂度为O(n^2) 


原理实现:

#include
#include
#include
using namespace std;

void solve(string s) {
	int l, r;
	int len = s.length();
	int cen = 0;
	vector p;
	p.push_back(0);

	for (r = 0; r < len; r++) {
		printf("cen = %d, pcen = %d, r = %d\n", cen, p[cen], r);
		p.push_back(0);
		int cenR = cen + p[cen];
		if (r <= cenR) {
			l = 2 * cen - r;
			p[r] = p[l] < cenR - r ? p[l] : cenR - r;
		} else {
			p[r] = 0;
		}

		// 这里是关键,如果从r为中心长度为0开始向两边扩展,则复杂度为O(n^2), 这里从p[r]开始扩展
		while((r-p[r]-1>=0) && (r+p[r]+1= cenR) {
			cen = r;
		}
	}

	cout << s << endl;
	int i;
	for (i = 0; i < len; i++) {
		printf("%d ", p[i]);
	} printf("\n");
}

int main() {
	solve("11111111111");
	return 0;
}



应用:







你可能感兴趣的:(常用算法,回文串,Manacher,On)