USACO以前的1.3.3 [calfflac] manacher算法

最新推荐文章于 2021-07-08 23:06:10 发布

CMTM4

最新推荐文章于 2021-07-08 23:06:10 发布

阅读量536

点赞数

分类专栏：算法研究 USACO 文章标签：算法 USACO manacher算法

本文链接：https://blog.csdn.net/baidu_23081367/article/details/42426239

版权

USACO 同时被 2 个专栏收录

71 篇文章 0 订阅

订阅专栏

算法研究

3 篇文章 0 订阅

订阅专栏

/*
TASK:calfflac
LANG:C++
*/

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;

const int max_n = 200000; 
char origin_string[max_n];//原串
int reflect[max_n * 2];//子串的第i个字母，在原串的位置
int origin_string_length(0), son_string_length(0);//原串长度，子串长度
char son_string[max_n * 2];//子串
int son_string_extend_length[max_n * 2];//子串第i个位置，可以左右扩展的最大长度

bool IF_LETTER(char &T)
{
	if ('a' <= T && T <= 'z')	return true;
	if ('A' <= T && T <= 'Z')	return true;
	return false;
}

bool IF_SAME(char A, char B)
{
	if ('a' <= A && A <= 'z')	A += 'A' - 'a';
	if ('a' <= B && B <= 'z')	B += 'A' - 'a';
	if (A == B)	return true;
	return false;
}

void manacher()
{
	int mid = 0, max_length = 0; //起点位置  最长扩展长度
	for (int i = 1; i != son_string_length - 1; ++ i)
	{
		int len = 0;//从i个地方能扩展的最大长度，默认为0
		if (i < mid + max_length)	len = min(son_string_extend_length[mid - (i - mid)], mid + max_length - i);
		//如果已经在最长的回文内，  那么就是回文串长度边界的位置， 或者就是之前求过的值
		while (IF_SAME(son_string[i + len + 1], son_string[i - len - 1]))	++ len;
		//更新回文
		if (i + len > mid + max_length)	
		{
			mid = i;
			max_length = len;	
			//更新回文串最长位置
		}
		son_string_extend_length[i] = len;//记录
	}
	//下面就是处理答案，因为开头是-开头，所以如果回文串涉及开头的话，那么会导致串减少。所以把所有以#结尾的串的扩展长度-1
	for (int i = 1; i != son_string_length - 1; ++ i)	if (son_string[i + son_string_extend_length[i]] == '#') -- son_string_extend_length[i];
	int ans = 0, P; 
	for (int i = 1; i != son_string_length - 1; ++ i)
		if (son_string_extend_length[i] > ans)
		{
			ans = son_string_extend_length[i];
			P = i;	
		}
	int begin = reflect[P - ans];	
	int end = reflect[P + ans];
	cout<<ans + 1<<endl;
	for (int i = begin; i != end + 1; ++ i)	printf("%c",origin_string[i]);
	cout<<endl;
}

int main()
{
	freopen("calfflac.in", "r", stdin);
	freopen("calfflac.out", "w", stdout);
	son_string[son_string_length++] = '+';
	char ch;
	while ((ch = getchar()) != EOF)	
	{
		origin_string[origin_string_length++] = ch;	
		if (IF_LETTER(ch))	
		{	
			son_string[son_string_length++] = ch;			
			reflect[son_string_length - 1] = origin_string_length - 1;
			son_string[son_string_length++] = '#';
		}
	}
	son_string[son_string_length] = '-';
	manacher();
	return 0;
}

manacher算法，可以在O(n)的时间内解决最长回文子串问题。

这个题显然是最长回文子串问题，只不过要把里面的串抽出来，但是归根结底，还是最长回文子串问题。

下面介绍manacher算法：

在求最长回文串的时候，我们做了很多无用运算，如果能有效利用这些运算，就可以大大的降低时间复杂度。

举一个例子

0123456789ABCDEFG //这样表示第串第i个位置，第一个就是0号位置，第10个就是A位置……

aba??+??aba //这就是串

↑ //b的最长子串就是3。从1位置开始，可以向左或者向右扩展1个单位长度

aba??+??aba

↑ //+这个位置，也就是5位置，可以向左向右扩展5个长度。

aba??+??aba

↑ //假如我们要求这个b字母，也就是9位置，以他为中心的回文子串长度？