「Luogu P4987」回文项链 解题报告

题面

求环中的长度为k(k为奇数)且回文中心不同的回文串个数

思路:

刚学manacher算法,就送上一道模板题,此题注重对manacher算法的理解

Manacher,但是不用插入其他符号,因为k是奇数,中心一定在字符上

不知道Manacher?

洛谷日报上有讲,但是比较难懂,建议上B站更深入了解下( o)

可以先把代码抄下来,然后一边看讲解一边理解代码含义,这样更能理解

反正我第一遍已经看蒙了

补充和纠正:

关于前几篇题解,一些小细节纠正一下

文末的代码更加简洁易懂

1、这是一个环,要做Manacher就应该拆环为链,第一篇题解说:应该重复复制三次(扩展为原来的三倍),但是实际上两次(扩展为原来的两倍)就够了

2、p数组(p[i])表示的是以i为中心的最长回文串的半径(包括i这个字符),原本的Manacher算法因为加入的额外的#号,而是p[i]表示的是就是原本最长回文串的长度+1,因此求Ans的时候应该用p[i]-1,但是现在我们不加#号了,p[i]表示的就是以i为中心的最长回文串的半径,此时长度就应该表示为p[i]*2-1,意思是半径*2-回文中心重复计算的字符

Code:

#include<bits/stdc++.h>
#define N 1000010//N<<1表示长度扩展为原来的两倍
using namespace std;
int n,k,p[N<<1],ans,res;
char s[N<<1];
void manacher()//标准的Manacher算法,就是没有加'#'号,而且起点为1
{
	s[0]='?';
	s[2*n+1]='!';
	int id=0,mx=0;//id表示目前最长回文串的中心,mx则是右边界
	for(int i=1;i<=2*n;i++)
	{
		if(i<mx)
			p[i]=min(p[2*id-i],mx-i);//如果中心在mx范围内,用mx去更新p[i]
		else
			p[i]=1;
		while(s[i-p[i]]==s[i+p[i]])//然后暴力
			p[i]++;
		if(mx<i+p[i])//更新mx和id
		{
			id=i;
			mx=i+p[i];
		}
	}
	return;
}
int main()
{
	int i;
	scanf("%d%d",&n,&k);
	scanf("%s",s+1);
	for(i=1;i<=n;i++)//复制
		s[i+n]=s[i];
	manacher();//跑一遍Manacher
	res=(k+1)>>1;//表示起始的中心i的位置,这个位置是第一个有可能长度为k的回文串的中心(串s[1-k]的中心)
	for(i=res;i<=res+n-1;i++)//对于每一个点遍历一次,所以是res+n-1
		if(p[i]*2-1>=k)//真正的长度只要大于等于k就是可以的(大于k的可以把他砍成k啊~)
			ans++;
	printf("%d",ans);
	return 0;
}

推荐题目:

一道有趣的黄题 P1210 回文检测

原文地址:https://www.cnblogs.com/hovny/p/10172247.html