設(shè)a[i]=bool(s[i]=='a')，b[i]=bool(s[i]=='b')，考慮a和a、b和b的卷積，由于卷積是對稱的，就可以統(tǒng)計出不連續(xù)回文子串個數(shù)了。可能說得比較簡略。再用manacher算出連續(xù)回文子串個數(shù)并減去。

#include<bits/stdc++.h> using namespace std; const int p=1e9+7; const int N=1<<18; typedef double flo; const flo pi=acos((flo)-1); struct vec{flo x,y;}; vec operator+(vec a,vec b){return{a.x+b.x,a.y+b.y};} vec operator-(vec a,vec b){return{a.x-b.x,a.y-b.y};} vec operator*(vec a,vec b){return{a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x};} void fft(vec*c,int n,int d){static int f[N];f[1]=n>>1;for(int i=2;i<n;++i){f[i]=f[i>>1]>>1|f[i&1];if(i>f[i])swap(c[i],c[f[i]]);}for(int i=1;i<n;i*=2){vec w={cos(pi/i*d),sin(pi/i*d)};for(int j=0;j<n;j+=i*2){vec s={1};for(int k=j;k<j+i;++k){vec v=s*c[k+i];c[k+i]=c[k]-v,c[k]=c[k]+v,s=s*w;}}}if(!~d)for(int i=0;i<n;++i)c[i].x/=n; } char z[N]; vec a[N],b[N]; int c[N],f[N]; int main(){scanf("%s",z);int n=strlen(z);int m=4<<__lg(n);for(int i=0;i<n;++i){a[i].x=98-z[i];b[i].x=z[i]-97;}fft(a,m,1);fft(b,m,1);for(int i=0;i<m;++i)a[i]=a[i]*a[i]+b[i]*b[i];fft(a,m,-1);for(int i=1;i<=n;++i)c[i]=(c[i-1]*2+1)%p;int s=0;for(int i=0;i<m;++i){int j=a[i].x+.5;(s+=c[(j+1)/2])%=p;}for(int i=n-1;~i;--i){z[i*2+2]=z[i];z[i*2+3]='#';}z[1]='#';z[0]='^';for(int i=1,j=0;z[i];++i){int k=j+f[j];f[i]=i<k?min(f[j*2-i],k-i):1;if(i+f[i]>=k){j=i;while(z[i-f[i]]==z[i+f[i]])++f[i];}(s+=p-f[i]/2)%=p;}printf("%d\n",s); }

轉(zhuǎn)載于:https://www.cnblogs.com/f321dd/p/5533930.html

與50位技術(shù)專家面對面20年技術(shù)見證，附贈技術(shù)全景圖

總結(jié)

以上是生活随笔為你收集整理的BZOJ3160: 万径人踪灭的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。