洛谷P4173 残缺的字符串（FFT）

传送门

话说为什么字符串会和卷积扯上关系呢……到底得脑洞大到什么程度才能想到这种东西啊……大佬太珂怕了……

因为通配符的关系，自动机已经废了

那么换种方式考虑，如果两个字符串每一位对应的编码都相等，那么这两个字符串相等

编码相等就代表$\sum_{i=1}^na[i]-b[i]=0$

然而这是不对的，有可能前面少一点，后面多一点，最好加起来还是$0$

那就平方一下$\sum_{i=1}^n(a[i]-b[i]=0)^2=0$，那就大丈夫了

于是我们得到了比一位一位匹配更麻烦的方法

看到平方……把它展开一下试试……结果……$\sum_{i=1}^na[i]^2+b[i]^2-2a[i]b[i]$

我们考虑把$b$给倒过来……这就是三个卷积啊！直接用FFT啊！加起来之后如果为$0$说明匹配上了

然而……通配符怎么办……如果有一位有通配符，那么这一位代表的编码相减肯定是$0$了

那么再改一下$\sum_{i=1}^na[i]*b[i]*(a[i]-b[i])=0$，如果一个位置是通配符，就把$a[i]/b[i]$设为$0$

然后就没有问题了……求一下卷积……如果第$i$位为$0$那么它就是一个能匹配上的字符串的结尾

 //minamoto
 #include<iostream>
 #include<cstdio>
 #include<cmath>
 using namespace std;
 const int N=;const double Pi=acos(-1.0);
 struct complex{
     double x,y;
     complex(double xx=,double yy=){x=xx,y=yy;}
     complex operator +(complex b){return complex(x+b.x,y+b.y);}
     complex operator -(complex b){return complex(x-b.x,y-b.y);}
     complex operator *(complex b){return complex(x*b.x-y*b.y,x*b.y+y*b.x);}
 }A[N],B[N],C[N];
 int n,m,l,r[N],limit=;double a[N],b[N];
 void FFT(complex *A,int type){
     for(int i=;i<limit;++i)
     if(i<r[i]) swap(A[i],A[r[i]]);
     for(int mid=;mid<limit;mid<<=){
         complex Wn(cos(Pi/mid),type*sin(Pi/mid));
         for(int R=mid<<,j=;j<limit;j+=R){
             complex w(,);
             for(int k=;k<mid;++k,w=w*Wn){
                 complex x=A[j+k],y=w*A[j+mid+k];
                 A[j+k]=x+y,A[j+mid+k]=x-y;
             }
         }
     }
     if(type==-)
     for(int i=;i<limit;++i) A[i]=(int)(A[i].x/limit+0.5);
 }
 char s1[],s2[];int l1,l2,ans[N],res=;
 int main(){
 //    freopen("testdata.in","r",stdin);
     scanf("%d%d",&l2,&l1);
     scanf("%s%s",s2,s1);
     m=l1+l2;
     while(limit<=m) limit<<=,++l;
     for(int i=;i<limit;++i)
     r[i]=(r[i>>]>>)|((i&)<<(l-));
     for(int i=;i<l1;++i)
     a[i]=s1[i]=='*'?:s1[i]-'a'+;
     for(int i=;i<l2;++i)
     b[l2-i-]=s2[i]=='*'?:s2[i]-'a'+;
     for(int i=;i<l1;++i) A[i].x=a[i]*a[i]*a[i];
     for(int i=;i<l2;++i) B[i].x=b[i];
     FFT(A,),FFT(B,);
     for(int i=;i<limit;++i) C[i]=A[i]*B[i];//a^3*b
     for(int i=;i<limit;++i) A[i].x=a[i]*a[i],A[i].y=;
     for(int i=;i<limit;++i) B[i].x=b[i]*b[i],B[i].y=;
     FFT(A,),FFT(B,);
     complex tmp(,);
     for(int i=;i<limit;++i) C[i]=C[i]-A[i]*B[i]*tmp;//-2ab*a*b
     for(int i=;i<limit;++i) A[i].x=a[i],A[i].y=;
     for(int i=;i<limit;++i) B[i].x=b[i]*b[i]*b[i],B[i].y=;
     FFT(A,),FFT(B,);
     for(int i=;i<limit;++i) C[i]=C[i]+A[i]*B[i];//b^3*a
     FFT(C,-);
     for(int i=l2-;i<l1;++i)
     if(C[i].x==0.0) ans[++res]=i-l2+;
     printf("%d\n",res);
     for(int i=;i<=res;++i) printf("%d ",ans[i]);
     return ;
 }