Codeforces17EPalisection

manacher+差分

我们定义l[i]为以i开头的回文串个数,r[i]为以i结尾的回文串个数,于是我们发现,对于一个点i, ∑ k = 1 i − 1 r [ k ] ∗ l [ i ] \sum_{k=1}^{i-1}r[k]*l[i] k=1i1r[k]l[i]就是i之前不与以i开头的回文串相交的个数,于是我们得到所有不相交的回文串个数为 ∑ i = 1 n ( l [ i ] ∗ ∑ k = 1 i − 1 r [ k ] ) \sum_{i=1}^n(l[i]*\sum_{k=1}^{i-1}r[k]) i=1n(l[i]k=1i1r[k]),然后我们就可以先做一遍马拉车,得到每个点的回文半径,然后差分,对于每个i,对于l数组在他的i-半径+1处+1,在i+1处-1,因为i之前的半径个点都能当做一个起点,对于r数组在i处+1,在i+半径处-1,这些是可以结尾的位置,然后暴力统计一遍就好了

代码

//By AcerMo
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
const int mod=51123987;
const int M=4000005;
using namespace std;
int n,f[M],l[M],r[M];
char a[M],b[M];
int main()
{
    scanf("%d%s",&n,a),b[0]=1;
    for(int i=0;i<n;i++) b[i*2]='@',b[i*2+1]=a[i];
    b[n*2]='@';b[n*2+1]='%';
    int mx=0,p=0,len=strlen(b)-1;
    for (int i=0;i<len;i++)
    {
        if (mx>i) f[i]=min(mx-i,f[2*p-i]);
        else f[i]=1;
        while (b[i-f[i]]==b[i+f[i]]) f[i]++;
        if (f[i]+i>mx) mx=i+f[i],p=i;
    }
    int ans=0,que=0;
    for(int i=0;i<=len;i++)
    {
        l[i-f[i]+1]++,l[i+1]--;
        r[i]++,r[i+f[i]]--;
        ans=(ans+f[i]/2)%mod;
    }
    ans=1ll*ans*(ans-1)/2%mod;
    for(int i=0;i<len;i++)
    {
        l[i]+=l[i-1],r[i]+=r[i-1];
        if(i%2) ans=(ans-1ll*que*l[i]%mod)%mod,que=(que+r[i])%mod;
    }
    cout<<(ans+mod)%mod;
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值