CF528D. Fuzzy Search

題意：DNA序列，在母串s中匹配模式串t，對(duì)于s中每個(gè)位置i，只要s[i-k]到s[i+k]中有c就認(rèn)為匹配了c。求有多少個(gè)位置匹配了t

---

預(yù)處理\(f[i][j]\)表示位置i可以匹配字符j
分別考慮每一個(gè)字符c，對(duì)s的每個(gè)位置i求出用\(s[i,i+m-1]\)匹配t，這個(gè)字符匹配了幾次
用\(a_i=[s的位置i匹配c]，\ b_i=[t_i=c]\)
那么c的匹配次數(shù)就是\(c_j=\sum\limits_{i=0}^{m-1}a_{j+i}b_i\)，位置i匹配了t當(dāng)且僅當(dāng)四種字符的匹配次數(shù)和等于t的長(zhǎng)度m

這時(shí)候就可以考慮bitset暴力過了

一個(gè)常用技巧是，反轉(zhuǎn)模式串(或母串)，然后就成了卷積的形式：
\[ c_j=\sum\limits_{i=0}^{m-1}a_{j+i}b_{m-1-i}=D_{m+j-1} \]
這樣計(jì)算是沒有問題的，因?yàn)閎只有\([0,m-1]\)有值其他地方為0

注意處理每個(gè)字符前memset a和b！！！！！

#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #include <cmath> using namespace std; typedef long long ll; const int N=(1<<20)+5, INF=1e9; const double PI=acos(-1); inline int read(){char c=getchar();int x=0,f=1;while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}return x*f; }struct meow{double x, y;meow(double a=0, double b=0):x(a), y(b){} }; meow operator +(meow a, meow b) {return meow(a.x+b.x, a.y+b.y);} meow operator -(meow a, meow b) {return meow(a.x-b.x, a.y-b.y);} meow operator *(meow a, meow b) {return meow(a.x*b.x-a.y*b.y, a.x*b.y+a.y*b.x);} meow conj(meow a) {return meow(a.x, -a.y);} typedef meow cd;namespace FFT{int n, rev[N];void ini(int lim) {n=1; int k=0;while(n<lim) n<<=1, k++;for(int i=0; i<n; i++) rev[i] = (rev[i>>1]>>1) | ((i&1)<<(k-1));}void dft(cd *a, int flag) {for(int i=0; i<n; i++) if(i<rev[i]) swap(a[i], a[rev[i]]);for(int l=2; l<=n; l<<=1) {int m=l>>1; cd wn = meow(cos(2*PI/l), flag*sin(2*PI/l));for(cd *p=a; p!=a+n; p+=l) {cd w(1, 0);for(int k=0; k<m; k++) {cd t = w*p[k+m];p[k+m] = p[k] - t;p[k] = p[k] + t;w=w*wn;}}}if(flag==-1) for(int i=0; i<n; i++) a[i].x/=n;}void mul(cd *a, cd *b) {dft(a, 1); dft(b, 1);for(int i=0; i<n; i++) a[i]=a[i]*b[i];dft(a, -1);} }using FFT::mul; using FFT::ini;int n, m, k, lim, f[N][5], cnt[5], id[300]; cd a[N], b[N], c[N]; char s[N], t[N]; int ans[N]; void solve(int now) { memset(a, 0, sizeof(a)); memset(b, 0, sizeof(b));for(int i=0; i<n; i++) a[i].x = f[i][now];for(int i=0; i<m; i++) b[m-1-i].x = id[(int)t[i]]==now; mul(a, b);for(int i=0; i<n; i++) ans[i] += int(a[m-1+i].x+0.5); } int main() {freopen("in","r",stdin);n=read(); m=read(); k=read(); lim=n+m-1; ini(lim);scanf("%s%s",s,t);id['A']=0; id['T']=1; id['C']=2; id['G']=3;int l=0, r=0; cnt[ id[(int)s[0]] ]++;for(int i=0; i<n; i++) {while(l<i-k) cnt[ id[(int)s[l++]] ]--;while(r<n-1 && r<i+k) cnt[ id[(int)s[++r]] ]++;for(int j=0; j<4; j++) if(cnt[j]) f[i][j]=1;}for(int i=0; i<4; i++) solve(i);int sum=0;for(int i=0; i<n; i++) if(ans[i]==m) sum++;printf("%d",sum); }

總結(jié)

以上是生活随笔為你收集整理的CF528D. Fuzzy Search [FFT]的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。