CF528D Solution
题解
差分求出\([i-k,i+k]\)区间(后文称其为 容错范围 )中每种字母出现的次数。设bitset数组\(pos[i][j]\)表示\(s\)串第\(j\)位的容错范围内是/否(\(1/0\))有第\(i\)种字母。对于\(t_i\),设\(tmp=pos[t_i]>>i\),\(tmp_j\)表示\(s_j\)为区间起点时是/否(\(1/0\))满足存在\(t_i\)在容错范围内,将\(pos[t_i]>>i\)并在一起即为答案。
代码
#include<bits/stdc++.h>
using namespace std;
const int N=2e5+10;
char s[N],t[N];
int cha[N][5],mp[200];
bitset<N> pos[5],ans;
inline int read()
{
int s=0,w=1; char ch=getchar();
while(ch<'0' || ch>'9') {if(ch=='-') w=-1; ch=getchar();}
while(ch>='0' && ch<='9') s=s*10+ch-'0',ch=getchar();
return s*w;
}
int main()
{
int n=read(),m=read(),k=read(),res=0;
scanf("%s%s",s+1,t+1);
mp['A']=1,mp['T']=2,mp['G']=3,mp['C']=4;
for(int i=1;i<=n;i++)
cha[max(1,i-k)][mp[s[i]]]++,cha[min(n+1,i+k+1)][mp[s[i]]]--;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=4;j++)
{
if(i) cha[i][j]+=cha[i-1][j];
if(cha[i][j]) pos[j][i]=1;
else pos[j][i]=0;
}
}
for(int i=0;i<=n;i++) ans[i]=1;
for(int i=1;i<=m;i++) ans&=(pos[mp[t[i]]]>>i);
for(int i=0;i<n;i++)
if(ans[i]==1) res++;
printf("%d",res);
return 0;
}