FFT做字符串匹配即可,详见代码
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#define debug(...) fprintf(stderr,__VA_ARGS__)
#define Debug(x) cout<<#x<<"="<<x<<endl
using namespace std;
typedef long long LL;
const int INF=1e9+7;
inline LL read(){
register LL x=0,f=1;register char c=getchar();
while(c<48||c>57){if(c=='-')f=-1;c=getchar();}
while(c>=48&&c<=57)x=(x<<3)+(x<<1)+(c&15),c=getchar();
return f*x;
}
inline char readc(){
register char c=getchar();
while(c==' '||c=='\n'||c=='\t') c=getchar();
return c;
}
const int MAXN=6e5+5;
const double Pi=acos(-1);
struct cmpx{
double x,y;
inline cmpx(){}
inline cmpx(double _x,double _y){x=_x,y=_y;}
inline friend cmpx operator + (cmpx a,cmpx b){return cmpx(a.x+b.x,a.y+b.y);}
inline friend cmpx operator - (cmpx a,cmpx b){return cmpx(a.x-b.x,a.y-b.y);}
inline friend cmpx operator * (cmpx a,cmpx b){return cmpx(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);}
}A[MAXN],B[MAXN];
namespace F_F_T{
int rev[MAXN],limit=1,l;
inline void init(int n){
for(limit=1;limit<=n;limit<<=1) l++;
for(int i=0;i<limit;i++)
rev[i]=(rev[i>>1]>>1)|((i&1)<<(l-1));
}
inline void FFT(cmpx *A,int type){
for(int i=0;i<limit;i++)
if(i<rev[i]) swap(A[i],A[rev[i]]);
for(int len=1;len<limit;len<<=1){
cmpx Wn=(cmpx){cos(Pi/len),type*sin(Pi/len)};
for(int i=0;i<limit;i+=(len<<1)){
cmpx w=(cmpx){1,0};
for(int j=0;j<len;j++,w=w*Wn){
cmpx x=A[i+j],y=w*A[i+len+j];
A[i+j]=x+y;
A[i+len+j]=x-y;
}
}
}
if(type==-1){
for(int i=0;i<limit;i++) A[i].x/=limit;
}
}
}using namespace F_F_T;
char s[MAXN],t[MAXN],id[4]={'A','G','C','T'};
int sum[4][MAXN],ans[MAXN];
int n,m,K,Ans;
int main(){
n=read(),m=read(),K=read();
for(int i=0;i<n;i++){
s[i]=readc();
for(int k=0;k<4;k++)
if(s[i]==id[k]) sum[k][i+1]++;//+1是为了实现前缀和中的初始值0
}
for(int i=0;i<m;i++) t[i]=readc();
for(int k=0;k<4;k++){
for(int i=1;i<=n;i++) sum[k][i]+=sum[k][i-1];
}
init(n+m-2);//最大下标的和
for(int k=0;k<4;k++){
for(int i=0;i<limit;i++) A[i]=B[i]=cmpx(0,0);
for(int i=0;i<n;i++)
if(sum[k][min(i+K+1,n)]-sum[k][max(i-K,0)])//前缀和
A[i].x=1;
for(int i=0;i<m;i++)
if(t[i]==id[k]) B[m-i-1].x=1;//翻转一个
FFT(A,1);FFT(B,1);
for(int i=0;i<limit;i++) A[i]=A[i]*B[i];
FFT(A,-1);
for(int i=m-1;i<n;i++) ans[i-m+1]+=(int)(A[i].x+0.5);
}
for(int i=0;i<n;i++)
if(ans[i]==m) Ans++;
printf("%d\n",Ans);
}