B

题目链接：https://cn.vjudge.net/contest/281959#problem/B

题目大意:给你n,m,k。然后输入两个字符串，n代表第一个字符串s1，m代表第二个字符串s2，然后问你第二个字符串在第一个字符串能匹配的次数（选定第一个字符串的位置之后，任意s2中一个字符串，都能在s1对应的位置左右k个都能找到相同的字符）。

具体思路：和 这一篇的思路基本使用一样的。

FFT（Rock Paper Scissors Gym - 101667H）

AC代码：

#include<iostream>
#include<cstring>
#include<string>
#include<cmath>
#include<algorithm>
#include<stdio.h>
using namespace std;
# define ll long long
const double PI = acos(-1.0);
const int maxn = 8e5+100;
struct complex
{
    double r,i;
    complex(double _r = 0,double _i = 0)
    {
        r = _r;
        i = _i;
    }
    complex operator +(const complex &b)
    {
        return complex(r+b.r,i+b.i);
    }
    complex operator -(const complex &b)
    {
        return complex(r-b.r,i-b.i);
    }
    complex operator *(const complex &b)
    {
        return complex(r*b.r-i*b.i,r*b.i+i*b.r);
    }
};
void change(complex y[],int len)
{
    int i,j,k;
    for(i = 1, j = len/2; i < len-1; i++)
    {
        if(i < j)
            swap(y[i],y[j]);
        k = len/2;
        while( j >= k)
        {
            j -= k;
            k /= 2;
        }
        if(j < k)
            j += k;
    }
}
void fft(complex y[],int len,int on)
{
    change(y,len);
    for(int h = 2; h <= len; h <<= 1)
    {
        complex wn(cos(-on*2*PI/h),sin(-on*2*PI/h));
        for(int j = 0; j < len; j += h)
        {
            complex w(1,0);
            for(int k = j; k < j+h/2; k++)
            {
                complex u = y[k];
                complex t = w*y[k+h/2];
                y[k] = u+t;
                y[k+h/2] = u-t;
                w = w*wn;
            }
        }
    }
    if(on == -1)
        for(int i = 0; i < len; i++)
            y[i].r /= len;
}
char str1[maxn],str2[maxn];
complex x1[maxn],x2[maxn];
int len1,len2,len=1;
int vis[maxn][5],id[maxn],cnt[10];
ll ans[maxn];
int main()
{
    int n,m,k;
    scanf("%d %d %d",&n,&m,&k);
    scanf("%s",str1+1);
    scanf("%s",str2+1);
    len1=strlen(str1+1);
    len2=strlen(str2+1);
    while(len<len1*2||len<len2*2)
        len<<=1;
   //     cout<<len<<endl;
    int l=0,r=0;
    id['A']=1,id['G']=2,id['T']=3,id['C']=4;
    for(int i=1; i<=n; i++)
    {
        while(l<n&&l<i-k)//判断这个区间内的字符。
            cnt[id[(int)str1[l++]]]--;
        while(r<n&&r<i+k)
            cnt[id[(int)str1[++r]]]++;
        for(int j=1; j<=4; j++)
            if(cnt[j])
                vis[i][j]=1;
    }
    for(int i=1; i<=4; i++)
    {
        for(int j=0; j<len; j++)
        {
            x1[j]=complex(0,0);
            x2[j]=complex(0,0);
        }
        for(int j=1; j<=n; j++)
        {
            if(vis[j][i])
                x1[j-1]=complex(1,0);
        }
        for(int j=1; j<=m; j++)
        {
            if(id[(int)str2[j]]==i)
                x2[m-j]=complex(1,0);
        }
        fft(x1,len,1);
        fft(x2,len,1);
        for(int j=0; j<len; j++)
        {
            x1[j]=x1[j]*x2[j];
        }
        fft(x1,len,-1);
        for(int j=0; j<len; j++)
        {
            ans[j]+=(ll)(x1[j].r+0.5);
        }
    }
    int num=0;
    for(int i=0; i<len; i++)
    {
    //    cout<<i<<" "<<ans[i]<<endl;
        if(ans[i]==m)
            num++;
    }
    printf("%d
",num);
    return 0;
}