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=1i−1r[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=1i−1r[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;
}