BZOJ3160 万径人踪灭 FFT+manacher

对不连续的对称序列计数显然可以转化为对称序列-连续对称序列

连续对称序列显然可以用manacher算法计算得到：P[i]的和即为连续对称序列的个数

所有对称序列总是关于某条对称轴对称，所以对称点的下标和为定值

另f[i]表示以i为对称轴的点对个数（含自身对称）

显然

其中下标为按manacher规则翻倍之后的下标

而这个布尔表达式我处理不来QAQ。。。

膜了PoPoQQQ题解才知道可以强行定下字母之后转为01乘积（这是一种处理01问题的有趣手段）

由于FFT的性质，每对点都会被统计两次（除了自身对称点），实际f[i]=f[i]+1>>1;

最终，其中-1是为了排除空串

注意数组不要开小了。。一定要开够2的XXX次幂才可以。。

#include
#include
#include
#include
using namespace std;
typedef long long LL;
const int Mn=300005;
const LL Mod=1000000007;
const double Pi=acos(-1);
int len,n,b[Mn],f[Mn],p[Mn];
char s[Mn],a[Mn];
complexA[Mn],ft[Mn];
LL bin[Mn],ans=0;
void FFT(complexA[],int n,int ty){
	 int i,j,k,m;
	 complext0,t1;
	 for(i=0;i>=1);
	 	if(i(cos(Pi/m),sin(Pi/m)*ty);
	 	for(j=1;ji)p[i]=min(mx-i,p[2*k-i]);
	 	while(a[i-p[i]]==a[i+p[i]])p[i]++;
	 	if(i+p[i]>mx){mx=i+p[i];k=i;}
	 	ans=(ans-p[i]/2)%Mod;
	 }
	 ans=(ans+Mod)%Mod;
	 printf("%lld\n",ans);
}
int main(){
    init();
    manacher();
	return 0;
}