CF17E Palisection（manacher)

题意

给出一个长度为N的字符串S，问S中有多少个回文子串对(i,j)使得i,j在S中的位置相交？
（N<=2*10⁶)

题解

 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #include<cmath>
 #include<algorithm>
 using namespace std;
 const long long N=;
 const long long mod=;
 char s[N*],str[N*];
 long long p[N*],n,m,f[N*],g[N*],sum,ans,tot,num;
 void init(){
     str[]=str[]='#';
     for(long long i=;i<=n;i++){
         str[i*]=s[i];
         str[i*+]='#';
     }
     m=n*+;
 }
 void manacher(){
     long long mx=,id;
     for(long long i=;i<=m;i++){
         if(mx>i)p[i]=min(p[id-(i-id)],p[id]+id-i);
         else p[i]=;
         while(str[i-p[i]]==str[i+p[i]])p[i]++;
         if(i+p[i]->mx)mx=i+p[i]-,id=i;
     }
 }
 int main(){
     scanf("%lld",&n);
     scanf("%s",s+);
     init();
     manacher();
     for(long long i=;i<=m;i++){
         f[i]++;
         f[i+p[i]-+]--;
         g[i]++;
         g[i-p[i]+-]--;
     }
     for(long long i=;i<=m;i++){
         f[i]+=f[i-];
         f[i]%=mod;
     }
     for(long long i=m;i>=;i--){
         g[i]+=g[i+];
         g[i]%=mod;
     }
     for(long long i=;i<=n;i++){
         ans+=g[i*]*sum;
         ans%=mod;
         sum+=f[i*];
         sum%=mod;
         num+=f[i*];
         num%=mod;
     }
     printf("%lld",((num-)*num/%mod-ans+mod)%mod);
     return ;
 }