B - Fuzzy Search （FFT）

题目链接：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+;
 struct complex
 {
     double r,i;
     complex(double _r = ,double _i = )
     {
         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 = , j = len/; i < len-; i++)
     {
         if(i < j)
             swap(y[i],y[j]);
         k = len/;
         while( j >= k)
         {
             j -= k;
             k /= ;
         }
         if(j < k)
             j += k;
     }
 }
 void fft(complex y[],int len,int on)
 {
     change(y,len);
     for(int h = ; h <= len; h <<= )
     {
         complex wn(cos(-on**PI/h),sin(-on**PI/h));
         for(int j = ; j < len; j += h)
         {
             complex w(,);
             for(int k = j; k < j+h/; k++)
             {
                 complex u = y[k];
                 complex t = w*y[k+h/];
                 y[k] = u+t;
                 y[k+h/] = u-t;
                 w = w*wn;
             }
         }
     }
     if(on == -)
         for(int i = ; i < len; i++)
             y[i].r /= len;
 }
 char str1[maxn],str2[maxn];
 complex x1[maxn],x2[maxn];
 int len1,len2,len=;
 int vis[maxn][],id[maxn],cnt[];
 ll ans[maxn];
 int main()
 {
     int n,m,k;
     scanf("%d %d %d",&n,&m,&k);
     scanf("%s",str1+);
     scanf("%s",str2+);
     len1=strlen(str1+);
     len2=strlen(str2+);
     while(len<len1*||len<len2*)
         len<<=;
    //     cout<<len<<endl;
     int l=,r=;
     id['A']=,id['G']=,id['T']=,id['C']=;
     for(int i=; 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=; j<=; j++)
             if(cnt[j])
                 vis[i][j]=;
     }
     for(int i=; i<=; i++)
     {
         for(int j=; j<len; j++)
         {
             x1[j]=complex(,);
             x2[j]=complex(,);
         }
         for(int j=; j<=n; j++)
         {
             if(vis[j][i])
                 x1[j-]=complex(,);
         }
         for(int j=; j<=m; j++)
         {
             if(id[(int)str2[j]]==i)
                 x2[m-j]=complex(,);
         }
         fft(x1,len,);
         fft(x2,len,);
         for(int j=; j<len; j++)
         {
             x1[j]=x1[j]*x2[j];
         }
         fft(x1,len,-);
         for(int j=; j<len; j++)
         {
             ans[j]+=(ll)(x1[j].r+0.5);
         }
     }
     int num=;
     for(int i=; i<len; i++)
     {
     //    cout<<i<<" "<<ans[i]<<endl;
         if(ans[i]==m)
             num++;
     }
     printf("%d\n",num);
     return ;
 }