题目大意就是字符串匹配,不过有一个门限k而已

之前有提到过fft做字符串匹配,这里和之前那种有些许不同

因为只有A,C,G,T四种字符,所以就考虑构造4个01序列

例如,模板串a关于'A'的01序列中,1代表这个位置可以匹配,而0则代表不能匹配。

这样构造出4个序列后,再对匹配串b做同样的处理

下面用a['A']代表a关于'A'的01序列,b同理

然后可以知道a['A'][i]&b['A'][i]如果是1则代表可以匹配,如果是0则代表不能匹配。

那么在位置i两个串能否匹配就可以写做

for(x = i ~ i+lenb) ans += a['A'][i]&b['A'][i]

如果ans恰好等于b中'A'的个数,那么就说明'A'匹配成功,接下来做'C','G','T'的情况

如果都可以匹配成功,那么这个位置就可以匹配了

如何用fft做这个呢?实际上也很简单

把b串颠倒一下就变成了a['A'][i]&b['A'][lenb-i]

就是一个卷积形式,于是就可以fft了

(测试感觉stl里的complex比较慢,非递归比递归快很多)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <bitset>
#include <cmath>
#include <complex>
using namespace std;
const double pi = acos(-);
const int maxn = ;
int A[][maxn], B[][maxn], ANS[maxn];
struct cp {
double a,b;
cp() { a = b = ; }
cp(double _a, double _b):a(_a), b(_b) {}
cp operator+(const cp&y)const { return cp(a+y.a, b+y.b); }
cp operator-(const cp&y)const { return cp(a-y.a, b-y.b); }
cp operator*(const cp&y)const { return cp(a*y.a-b*y.b,a*y.b+b*y.a); }
cp operator!()const { return cp(a,-b); }
};
int T[];
int la, lb, k;
char str[maxn];
class FFT {
int n, L, R[maxn];
cp a[maxn], b[maxn];
void DFT(cp *a,int n,int f) {
for(int i = ; i < n; i++) if(i < R[i]) swap(a[i], a[R[i]]);
for(int i = ; i < n; i <<= ) {
cp t, x, wn(cos(pi/i), sin(pi*f/i));
for(int j = ; j < n; j += (i<<)) {
cp w(, );
for(int k = ; k < i; k++,w = w*wn) {
x = a[j+k];
t = w*a[j+k+i];
a[j+k] = x+t;
a[j+k+i] = x-t;
}
}
}
}
public:
int c[maxn];
FFT() { memset(R, , sizeof(R)); }
void init(int *A, int *B, int n1, int n2) {
memset(a, , sizeof(a));
memset(b, , sizeof(b));
for(int i = ; i <= n1; i++) a[i] = cp(A[i], );
for(int i = ; i <= n2; i++) b[i] = cp(B[i], );
n2 += n1;
for(n = , L = ; n <= n2; n <<= ) L++;
for(int i = ; i < n; i++) R[i] = (R[i>>]>>)|((i&)<<(L-));
}
void calc() {
DFT(a, n, );
DFT(b, n, );
for(int i = ; i <= n; i++) a[i] = a[i]*b[i];
DFT(a, n, -);
for(int i = ; i <= n; i++) c[i] = (a[i].a/n + 0.5);
}
} fft; int main() {
cin>>la>>lb>>k;
T['A'] = ; T['C'] = ; T['G'] = ; T['T'] = ;
scanf("%s", str);
for(int i = ; i < la; i++) {
int ty = T[str[i]];
A[ty][i] = ;
for(int j = i+; j <= min(i+k, la-); j++) {
if(ty == T[str[j]]) break;
A[ty][j] = ;
}
for(int j = i-; j >= max(, i-k); j--) {
if(A[ty][j]) break;
A[ty][j] = ;
}
}
scanf("%s", str);
for(int i = ; i < lb; i++) {
int ty = T[str[i]];
B[ty][lb-i-] = ;
}
for(int i = lb-; i <= la+lb-; i++) ANS[i] = ;
for(int i = ; i < ; i++) {
fft.init(A[i], B[i], la, lb);
fft.calc();
int t = ;
for(int j = ; j < lb; j++) if(B[i][j]) t++;
for(int j = lb-; j <= la+lb-; j++)
ANS[j] &= (fft.c[j] == t);
}
int ans = ;
for(int i = lb-; i <= la+lb-; i++) ans += ANS[i];
cout<<ans<<endl;
}

51nod 1565模糊搜索(FFT)的更多相关文章

  1. 51nod 1565 模糊搜索 FFT

    这...好强啊\(QwQ\) 思路:卷积?\(FFT\)? 提交:\(5\)次 错因:一开始的预处理写错了(竟然只错了最后几个大点)闹得我以为\(FFT\)写挂了\(QwQ\) 题解: 对四种字符分开 ...

  2. 51NOD 1565:模糊搜索——题解

    http://www.51nod.com/onlineJudge/questionCode.html#problemId=1565&noticeId=445588 有两个基因串S和T,他们只包 ...

  3. 【51nod】1565 模糊搜索

    题解 这个字符集很小,我们可以把每个字符拿出来做一次匹配,把第一个字符串处理每个出现过的该字符处理成一个区间加,即最后变成第一个字符串的该位置能够匹配某字符 例如对于样例 10 4 1 AGCAATT ...

  4. 51nod 算法马拉松 34 Problem D 区间求和2 (FFT加速卷积)

    题目链接  51nod 算法马拉松 34  Problem D 在这个题中$2$这个质数比较特殊,所以我们先特判$2$的情况,然后仅考虑大于等于$3$的奇数即可. 首先考虑任意一个点对$(i, j)$ ...

  5. FFT/NTT [51Nod 1028] 大数乘法 V2

    题目链接:51Nod 传送门 没压位,效率会低一点 1.FFT #include <cstdio> #include <cstring> #include <algori ...

  6. 51NOD 1258 序列求和 V4 [任意模数fft 多项式求逆元 伯努利数]

    1258 序列求和 V4 题意:求\(S_m(n) = \sum_{i=1}^n i^m \mod 10^9+7\),多组数据,\(T \le 500, n \le 10^{18}, k \le 50 ...

  7. 51nod 1172 Partial Sums V2 卡精度的任意模数FFT

    卡精度的任意模数fft模板题……这道题随便写个表就能看出规律来(或者说考虑一下实际意义),反正拿到这题之后,很快就会发现他是任意模数fft模板题.然后我就去网上抄了一下板子……我打的是最土的任意模数f ...

  8. 51nod 1028 大数乘法 V2 【FFT模板题】

    题目链接 模板题.. #include<bits/stdc++.h> using namespace std; typedef int LL; typedef double db; nam ...

  9. 51Nod 快速傅里叶变换题集选刷

    打开51Nod全部问题页面,在右边题目分类中找到快速傅里叶变换,然后按分值排序,就是本文的题目顺序. 1.大数乘法问题 这个……板子就算了吧. 2.美妙的序列问题 长度为n的排列,且满足从中间任意位置 ...

随机推荐

  1. python的字典数据类型及常用操作

    字典的定义与特性 字典是Python语言中唯一的映射类型. 定义:{key1: value1, key2: value2} 1.键与值用冒号“:”分开: 2.项与项用逗号“,”分开: 特性: 1.ke ...

  2. Apache Maven(三):POM

    什么是 POM? POM (Project Object Model) 项目对象模型.它是一个XML文件,其中包含有关Maven用于构建项目的项目和配置细节的信息.它包含大多数项目的默认值.例如,构建 ...

  3. Apache httpd Server 配置正向代理

    背景 代理(Proxy),位于客户端与实际服务端之间,当客户端需要请求服务端内容时,先向代理发起请求,代理将请求转发到实际的服务器,再原路返回.也可以在代理服务器设置缓存,将实际服务器上不常变化的内容 ...

  4. typecho博客组插件:openSug.js百度搜索框下拉提示免费代码

      Typecho候选搜索增强插件:安装openSug插件即可获得带有“搜索框提示”功能的搜索框,让Typecho搜索更便捷! 支持百度.谷歌.雅虎.Yandex.360好搜.UC神马.酷狗.优酷.淘 ...

  5. Laravel5.x 封装的上传图片类

    图片缩放需要用conposer安装 ImageManagerStatic类 可参考下面的地址安装: https://www.jb51.net/article/128159.htm 控制器里: 控制器里 ...

  6. lnmp+phpmyadmin+laravel 环境配置

    腾讯云 Ubuntu16.04 添加用户 useradd 与 adduser Ubuntu下useradd不会在/home下自动创建与用户名同名的用户目录,而且不会自动选择shell版本,也没有设置密 ...

  7. Python学习:for 循环 与 range()函数

    for 循环   For … in 语句是另一种循环语句,其特点是会在一系列对象上进行迭代(Iterates),即它会遍历序列中的每一个项目 注意:  1.else 部分是可选的.当循环中包含它时,它 ...

  8. 转译符,re模块,random模块

    一, 转译符 1.python 中的转译符 正则表达式中的内容在Python中就是字符串 ' \n ' : \ 转移符赋予了这个n一个特殊意义,表示一个换行符 ' \ \ n' :  \ \  表示取 ...

  9. ionic 打包apk Failure [INSTALL_FAILED_USER_RESTRICTED: Install canceled by user]

    错误日志如下: Built the following apk(s): /Users/hongye0/Documents/project/haitoujiaApp/platforms/android/ ...

  10. Spring事务:一种编程式事务,三种声明式事务

    事务隔离级别 隔离级别是指若干个并发的事务之间的隔离程度.TransactionDefinition 接口中定义了五个表示隔离级别的常量: TransactionDefinition.ISOLATIO ...