题目大意就是字符串匹配,不过有一个门限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. yaml文件 .yml

    YAML文件简介 我们可能在spring配置文件里见到过.yml格式的东东,配置文件不都是.propertie或者.xml文件吗?.yml是什么鬼,今天我带你们来一探究竟. YAML(Yet Anot ...

  2. python网络编程,通过服务名称和会话类型(tcp,udp)获取端口号,简单的异常处理

    作为一个php程序员,同时有对网络方面感兴趣,php就比较蛋疼了,所以就抽了些时间看python 之前学python基础因为工作原因,断断续续的看了个基础,差不多是可以写代码了 最近在看<pyt ...

  3. linux几条基本命令和解释

    pwd 查看当前目录/     根目录ls    查看当前目录所包含文件ls -l    查看当前目录所包含文件的详细信息d rwx rwx r-x 1 root root1  2     3   4 ...

  4. python快速改造:基础知识

    改造"Hacking"并不同于破坏"cracking" python快速改造:基础知识 一行就是一行,不管多少,不用加分号 交互式python解释器可以当作计算 ...

  5. android 界面切换教程

    Android的界面切换其实并没有什么难度,就和java的界面切换一样简单,但是java的使用的是Swing框架进行的切换,而Android使用的是Intent 来切换框架的,这是一个全新的使用方法, ...

  6. UVA - 12230

    #include <bits/stdc++.h> using namespace std; int n; double d; double p,l,v,ret,sum; ; /* 村庄A, ...

  7. 如何打war包

    1. 利用jdk里的工具 例如我们要打包的文件在D:\myHome\dist: 运行 cmd: cd D:\myHome\dist 进入D:\myHome\dist 然后输入 D:\myHome\di ...

  8. Win10启动不了的问题处理记录

    前几天电脑突然出现蓝屏的情况,而且使用Win10自带的修复功能根本没有卵用,修复不了,很郁闷,死活进不了系统了,说什么“INACCESSABE BOOT DEVICE”,好像是引导设备不可用. 到网上 ...

  9. Google序列化库FlatBuffers 1.1发布,及与protobuf的比较

    个人总结: FlatBuffer相对于Protobuffer来讲,优势如下: 1. 由于省去了编解码的过程,所以从速度上快于Protobuffer,个人测试结果100w次编解码,编码上FlatBuff ...

  10. 原生js实现五子棋

    为什突然做这个,因为这是个笔试题,拖了一个月才写(最近终于闲了O(∩_∩)O),废话不多说,说说这个题吧 题目要求 编写一个单机[五子棋]游戏,要求如下: 1.使用原生技术实现,兼容 Chrome 浏 ...