CodeForces - 17E Palisection【Manacher】【前缀和】【差分】

最新推荐文章于 2023-07-01 17:05:35 发布

invoKer-

最新推荐文章于 2023-07-01 17:05:35 发布

阅读量183

点赞数

CC 4.0 BY-SA版权

分类专栏： codeforces 差分 Manacher

本文链接：https://blog.csdn.net/inv00ker/article/details/108943762

codeforces 同时被 3 个专栏收录

145 篇文章

订阅专栏

差分

3 篇文章

订阅专栏

Manacher

1 篇文章

订阅专栏

这篇博客主要介绍了如何运用Manacher算法高效地找出字符串中的最长回文子串，并给出了C++实现的详细代码。算法核心是利用回文串的对称性质减少重复计算，提高求解效率。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目链接：https://codeforces.com/contest/17/problem/E

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
static const int MAXN=2e6+10;
static const int MAXM=4e6+10;
static const int MOD=51123987;
char s[MAXN];
char s_new[MAXM];
int p[MAXM];
int h[MAXM],t[MAXM];
int n;
int init()
{
    int len=strlen(s);
    s_new[0]='$'; s_new[1]='#';
    int j=2;
    for (int i=0;i<len;i++)
    {
        s_new[j++]=s[i]; s_new[j++]='#';
    }
    s_new[j++]= '^'; s_new[j]='\0';    
    return j;
}
void manacher()
{
    int len=init();
    int max_len=-1;
    int id,mx=0;
    for (int i=1;i<len;i++)
    {
        if(i<mx) p[i]=min(p[2*id-i],mx-i);
        else p[i]=1;
        while(s_new[i-p[i]]==s_new[i+p[i]]) p[i]++;
        if(mx<i+p[i])
        {
            id=i; mx=i+p[i];
        }
    }
}
int main()
{
    scanf("%d%s",&n,s);
    manacher();
    long long cnt=0;
    for(int i=1;i<=2*n+1;i++)
    {
        h[i-p[i]+1]++; h[i+1]--;
        t[i]++; t[i+p[i]]--; 
        cnt=(cnt+p[i]/2)%MOD;
    }
    for(int i=1;i<=2*n+1;i++) h[i]+=h[i-1],t[i]+=t[i-1];
    cnt=cnt*(cnt-1)/2%MOD;
    long long tott=0;
    for(int i=2;i<=2*n-2;i+=2)
    {
        tott=(tott+t[i])%MOD;
        cnt=(cnt-tott*h[i+2]%MOD)%MOD;
    }
    printf("%lld\n",(cnt+MOD)%MOD);
    return 0;
}