【bzoj4259】 残缺的字符串 FFT

又是一道FFT套路题

思路可以参考bzoj4503，题解

我们对串S和串T中出现的*处全部赋值为0。

反正最终的差异度式子大概就是

$C[i]=\sum_{j=0}^{|T|-1}S[i+j]T[j](S[i+j]-T[j])^2$

然后和上一题一样的展开方式，将T串reverse一下做FFT再统计下即可。

然后这题卡常，FFT的长度是100W，所以用NTT会被卡常（我就T了）

然后就没了

 #include<bits/stdc++.h>
 #define PI acos(-1)
 #define zero(x) (fabs(x)<0.5)
 #define M 1<<20
 using namespace std;
 struct cp{
     double r,i;
     cp(){r=i=;} cp(int x){r=x; i=;}
     cp(double rr,double ii){r=rr;i=ii;}
     friend cp operator +(cp a,cp b){return cp(a.r+b.r,a.i+b.i);}
     friend cp operator -(cp a,cp b){return cp(a.r-b.r,a.i-b.i);}
     friend cp operator *(cp a,cp b){return cp(a.r*b.r-a.i*b.i,a.r*b.i+a.i*b.r);}
 };

 cp a[M]={},aa[M]={},aaa[M]={};
 cp b[M]={},bb[M]={},bbb[M]={};
 int n; cp ans[M]={};

 void change(cp a[],int n){
     for(int i=,j=;i<n-;i++){
         if(i<j) swap(a[i],a[j]);
         int k=n>>;
         while(j>=k) j-=k,k>>=;
         j+=k;
     }
 }
 void FFT(cp a[],int n,int on){
     change(a,n);
     for(int h=;h<=n;h<<=){
         cp wn=cp(cos(*PI/h),on*sin(*PI/h));
         for(int j=;j<n;j+=h){
             cp w=cp(,);
             for(int k=j;k<j+(h>>);k++){
                 cp u=a[k],t=w*a[k+(h>>)];
                 a[k]=u+t;
                 a[k+(h>>)]=u-t;
                 w=w*wn;
             }
         }
     }
     if(on==-) for(int i=;i<n;i++) a[i].r/=n;
 }

 char s[M]={},c[M]={};
 int lens,lenc,len=;
 int main(){
     scanf("%d%d",&lens,&lenc);
     scanf("%s%s",c,s);
     lens=strlen(s); lenc=strlen(c);
     while(len<lens+lenc) len<<=;
     reverse(c,c+lenc);
     for(int i=;i<lens;i++) a[i]=(s[i]=='*'?:s[i]-'a'+),aa[i]=a[i]*a[i],aaa[i]=aa[i]*a[i];
     for(int i=;i<lenc;i++) b[i]=(c[i]=='*'?:c[i]-'a'+),bb[i]=b[i]*b[i],bbb[i]=bb[i]*b[i];
     FFT(a,len,); FFT(aa,len,); FFT(aaa,len,);
     FFT(b,len,); FFT(bb,len,); FFT(bbb,len,);
     for(int i=;i<len;i++) ans[i]=aaa[i]*b[i]-*aa[i]*bb[i]+a[i]*bbb[i];
     FFT(ans,len,-);
     int sum=;
     for(int i=lenc-;i<lens;i++)
     if(zero(ans[i].r)) sum++;
     cout<<sum<<endl;
     for(int i=lenc-;i<lens;i++)
     if(zero(ans[i].r)) printf("%d ",i-lenc+);
 }