题面
求环中的长度为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 回文检测