BZOJ 4259 [FFT]

最新推荐文章于 2024-07-31 00:30:00 发布
最新推荐文章于 2024-07-31 00:30:00 发布
阅读量734 收藏 1
点赞数 2
分类专栏： DFT 文章标签： fft
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Vectorxj/article/details/70214774
版权

Description

很久很久以前，在你刚刚学习字符串匹配的时候，有两个仅包含小写字母的字符串A和B，其中A串长度为m，B串长度为n。可当你现在再次碰到这两个串时，这两个串已经老化了，每个串都有不同程度的残缺。你想对这两个串重新进行匹配，其中A为模板串，那么现在问题来了，请回答，对于B的每一个位置i，从这个位置开始连续m个字符形成的子串是否可能与A串完全匹配?

Solution

将星号看成0；匹配就等价于 ∑aibi(ai−bi)2==0 翻转后FFT即可。
为什么我的FFT这么sb卡常卡了好久，最后用最原始的FFT过掉了，难道这个有问题。。。

（这是过的）

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;

inline char get(void) {
    static char buf[100000], *S = buf, *T = buf;
    if (S == T) {
        T = (S = buf) + fread(buf, 1, 100000, stdin);
        if (S == T) return EOF;
    }
    return *S++;
}
inline void read(int &x) {
    static char c; x = 0;
    for (c = get(); c < '0' || c > '9'; c = get());
    for (; c >= '0' && c <= '9'; c = get()) x = x * 10 + c - '0';
}
inline void read(char *s) {
    static int len; len = 0; static char c;
    for (c = get(); (c < 'a' || c > 'z') && c != '*'; c = get());
    for (; (c >= 'a' && c <= 'z') || c == '*'; c = get()) s[len++] = c;
}

const int N = 1203030;
const double PI = acos(-1);

struct C {
    double r, i;
    inline C(double R = 0, double I = 0):r(R), i(I) {}
    friend C operator +(const C &b, const C &a) {
        return C(b.r + a.r, b.i + a.i);
    }
    friend C operator -(const C &b, const C &a) {
        return C(b.r - a.r, b.i - a.i);
    }
    friend C operator *(const C &b, const C &a) {
        return C(b.r * a.r - b.i * a.i, b.r * a.i + b.i * a.r);
    }
    friend C operator /(const C &b, const double &a) {
        return C(b.r / a, b.i / a);
    }
};
int n1, n2, L, num, m, n, ans;
int a[N], b[N];
C A[N], B[N], F[N];
int R[N], Ans[N];
char s[N];

inline void FFT(C *a, int n, int f) {
    static C x, y;
    for (int i = 0; i < n; i++)
        if (R[i] > i) swap(a[R[i]], a[i]);
    for (int i = 1; i < n; i <<= 1) {
        C wn(cos(PI / i), f * sin(PI / i));
        for (int j = 0; j < n; j += (i << 1)) {
            C w(1, 0);
            for (int k = 0; k < i; k++) {
                x = a[j + k]; y = a[j + k + i] * w;
                a[j + k] = x + y; a[j + k + i] = x - y; w = w * wn;
            }
        }
    }
    if (f == -1) for (int i = 0; i < n; i++) a[i] = a[i] / n;
}

int main(void) {
    freopen("1.in", "r", stdin);
    freopen("1.out", "w", stdout);
    read(n1); read(n2);
    read(s); n1--;
    for (int i = 0; i <= n1; i++)
        a[i] = (s[i] == '*' ? 0 : s[i] - 'a' + 1);
    read(s); n2--;
    for (int i = 0; i <= n2; i++)
        b[i] = (s[i] == '*' ? 0 : s[i] - 'a' + 1);
    n = n2 * 2; n += 2; reverse(a, a + n1 + 1);
    for (m = 1; m <= n; m <<= 1) L++; L--;
    for (int i = 1; i < m; i++) R[i] = (R[i >> 1] >> 1) | ((i & 1) << L);
    for (int i = 0; i < m; i++) {
        A[i] = a[i] * a[i] * a[i]; B[i] = b[i];
    }
    FFT(A, m, 1); FFT(B, m, 1);
    for (int i = 0; i < m; i++) F[i] = A[i] * B[i];
    for (int i = 0; i < m; i++) {
        A[i] = a[i]; B[i] = b[i] * b[i] * b[i];
    }
    FFT(A, m, 1); FFT(B, m, 1);
    for (int i = 0; i < m; i++) F[i] = F[i] + A[i] * B[i];
    for (int i = 0; i < m; i++) {
        A[i] = a[i] * a[i]; B[i] = b[i] * b[i];
    }
    FFT(A, m, 1); FFT(B, m, 1);
    for (int i = 0; i < m; i++) F[i] = F[i] - A[i] * B[i] * 2;
    FFT(F, m, -1);
    for (int i = n1; i <= n2; i++)
        if ((int)(F[i].r + 0.4) == 0) Ans[++ans] = i;
    printf("%d\n", ans);
    for (int i = 1; i <= ans; i++) printf("%d ", Ans[i] - n1 + 1);
    putchar('\n');
}

（这有很大问题啊）

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;

inline char get(void) {
    static char buf[100000], *S = buf, *T = buf;
    if (S == T) {
        T = (S = buf) + fread(buf, 1, 100000, stdin);
        if (S == T) return EOF;
    }
    return *S++;
}
inline void read(int &x) {
    static char c; x = 0;
    for (c = get(); c < '0' || c > '9'; c = get());
    for (; c >= '0' && c <= '9'; c = get()) x = x * 10 + c - '0';
}
inline void read(char *s) {
    static int len; len = 0; static char c;
    for (c = get(); (c < 'a' || c > 'z') && c != '*'; c = get());
    for (; (c >= 'a' && c <= 'z') || c == '*'; c = get()) s[len++] = c;
}

const int N = 1203030;
const double PI = acos(-1);

struct C {
    double r, i;
    inline C(double R = 0, double I = 0):r(R), i(I) {}
    friend C operator +(const C &b, const C &a) {
        return C(b.r + a.r, b.i + a.i);
    }
    friend C operator -(const C &b, const C &a) {
        return C(b.r - a.r, b.i - a.i);
    }
    friend C operator *(const C &b, const C &a) {
        return C(b.r * a.r - b.i * a.i, b.r * a.i + b.i * a.r);
    }
    friend C operator /(const C &b, const double &a) {
        return C(b.r / a, b.i / a);
    }
};
int n1, n2, L, num, m, n, ans;
int a[N], b[N];
C A[N], B[N], F[N];
int R[N], Ans[N];
char s[N];
C g;
C w[2][N];

void Prep(int n) {
    num = n; g = C(cos(PI * 2 / num), sin(PI * 2 / num));
    w[0][0] = w[1][0] = C(1, 0);
    for (int i = 1; i < num; i++) w[1][i] = w[1][i - 1] * g;
    for (int i = 1; i < num; i++) w[0][i] = w[1][num - i];
}
inline void FFT(C *a, int n, int f) {
    static C x, y;
    for (int i = 0; i < n; i++)
        if (R[i] > i) swap(a[R[i]], a[i]);
    for (int i = 1; i < n; i <<= 1)
        for (int j = 0; j < n; j += (i << 1))
            for (int k = 0; k < i; k++) {
                x = a[j + k]; y = a[j + k + i] * w[f][num / (i << 1) * k];
                a[j + k] = x + y; a[j + k + i] = x - y;
            }
    if (!f) for (int i = 0; i < n; i++) a[i] = a[i] / n;
}

int main(void) {
    freopen("1.in", "r", stdin);
    freopen("1.out", "w", stdout);
    read(n1); read(n2);
    read(s); n1--;
    for (int i = 0; i <= n1; i++)
        a[i] = (s[i] == '*' ? 0 : s[i] - 'a' + 1);
    read(s); n2--;
    for (int i = 0; i <= n2; i++)
        b[i] = (s[i] == '*' ? 0 : s[i] - 'a' + 1);
    n = n2 * 2; n += 2; reverse(a, a + n1 + 1);
    for (m = 1; m <= n; m <<= 1) L++; L--;
    Prep(m);
    for (int i = 1; i < m; i++) R[i] = (R[i >> 1] >> 1) | ((i & 1) << L);
    for (int i = 0; i < m; i++) {
        A[i] = a[i] * a[i] * a[i]; B[i] = b[i];
    }
    FFT(A, m, 1); FFT(B, m, 1);
    for (int i = 0; i < m; i++) F[i] = A[i] * B[i];
    for (int i = 0; i < m; i++) {
        A[i] = a[i]; B[i] = b[i] * b[i] * b[i];
    }
    FFT(A, m, 1); FFT(B, m, 1);
    for (int i = 0; i < m; i++) F[i] = F[i] + A[i] * B[i];
    for (int i = 0; i < m; i++) {
        A[i] = a[i] * a[i]; B[i] = b[i] * b[i];
    }
    FFT(A, m, 1); FFT(B, m, 1);
    for (int i = 0; i < m; i++) F[i] = F[i] - A[i] * B[i] * 2;
    FFT(F, m, 0);
    for (int i = n1; i <= n2; i++)
        if ((int)(F[i].r + 0.4) == 0) Ans[++ans] = i;
    printf("%d\n", ans);
    for (int i = 1; i <= ans; i++) printf("%d ", Ans[i] - n1 + 1);
    putchar('\n');
}
Vectorxj
关注
专栏目录
洛谷P4173 残缺的字符串(FFT匹配字符串)
qq_21829533的博客
03-13 319
首先考虑完整两串的匹配（我知道可以kmp，但是模糊匹配的思路需要从这里引出 假设短串a长为lena长串b长为lenb 我们知道两串如果完整匹配，他们每一个位置都要相同，即a[i]−b[y+i]=0a[i]-b[y+i]=0a[i]−b[y+i]=0 所以说a要完整匹配要满足： ∑i=0lenaa[i]−b[y+i]=0\sum_{i=0}^{lena} a[i]-b[y+i]=0i=0∑lena​...
【模板】通配符单模式串匹配 -FFT bzoj4259: 残缺的字符串
远行客
09-20 577
传送门:bzoj4259 题解 带通配符的字符串匹配无法有效地用kmpkmpkmp处理，这时大常数的FFTFFTFFT派上了用场。 这题已经升级为一种套路/模板了，暂且引用ebola’s题解的称呼：带通配符的单模式串匹配。 设模式串为TTT，文本串为SSS。“*”对应0，a,b,...,za,b,...,za,b,...,z分别对应1,2,...,261,2,...,261,2,...,26。 ...
P4173 残缺的字符串（FFT）
weixin_30799995的博客
03-10 87
[Luogu4173] 题解 \(1.\)定义匹配函数 \(2.\)定义完全匹配函数 \(3.\)快速计算每一位的完全匹配函数值 #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #include<cmath> #define debug(....
数字信号处理之 快速傅里叶变换（FFT）
最新发布
mohen_777的博客
07-31 1131
快速傅里叶变换（FFT）并不是与DFT不同的另一种变换，而是为了减少计算次数的一种快速有效的算法。将基-2和基-4分解揉合在一起，提出一种分裂基FFT算法。以此二分法继续下分，可得到一个按级分解的“按时间抽取法”。基-2算法、基-4算法、基-8算法、混合基算法。着重讨论算法原理，并提供了实现算法的详细细节。为2的整数幂的FFT，也称基-2FFT。，上述的FFT算法皆可直接用来运算IDFT。直接算法对FFT算法的计算时间之比近似为。DIT算法的输入是反序，输出是正序的。，则有两个快速傅里叶变换的方法。
luogu P4173 残缺的字符串
zsyz_ZZY的博客
12-17 209
背景： FFT的字符串假题。 题目传送门： https://www.luogu.org/problemnew/show/P4173 题意： 定义∗*∗表示一个任意字符，求一个字符串从哪里开始能在模式串完全匹配。 思路： 考虑转化为数字问题。 另000表示∗*∗，111表示aaa，222表示bbb… 考虑一一对应。 那么假设两个位置匹配，就会有(ai−bi)2aibi=0(a_i-b_i...
洛谷P4173：残缺的字符串（FFT、通配符匹配）
BUG_Creater_jie的博客
01-08 360
利用FFT解决字符串匹配问题
FFT字符串匹配（解决通配符问题）
lifehappy
08-16 752
FFT字符串匹配 定义字符串下标从000，开始，有文本串AAA长度为nnn，模式串BBB长度为mmm，我们可以考虑一个函数f(x,y)=A(x)−B(y)f(x, y) = A(x) - B(y)f(x,y)=A(x)−B(y)。 我们设F(x)(x≥m−1)=∑i=0m−1f(x−m+1+i,i)F(x)(x \ge m - 1) = \sum\limits_{i = 0} ^{m - 1} f(x - m + 1 + i, i)F(x)(x≥m−1)=i=0∑m−1​f(x−m+1+i,i)，由定义显然
[BZOJ4259] 残缺的字符串 （FFT）
AZUI _ oi memories
02-21 2301
题意：定义*号可匹配任意字符。给出A串，B串，求A串在B中完全匹配的所有位置。 将*号视作0，则两个等长的串可匹配当且仅当Σ(a[i]-b[i])^2*a[i]*b[i]==0。将A串和B串最左边对齐，每次上式都要重算一次，不科学。所以讲A串先反转过来然后补上*号，就是卷积啦。 原来的万径人踪灭一题当中用FFT来求了回文串，这里又能玩字符串匹配，真是太6了。 #include #includ
bzoj FFT 的模版
12-20
bzoj FFT 的模版
BZOJ2179: FFT快速傅立叶 FFT实现高精度乘法
05-14
这道题目是要求用 FFT 实现高精度乘法。我们可以先将两个高精度数转换为多项式，然后对这两个多项式进行相乘，最后将结果转换为高精度数输出即可。 具体来说，将两个高精度数 $a$ 和 $b$ 分别转换为多项式 $A(x)$ ...
FFT简介（洛谷P1919/3803/4245、BZOJ2179）
forezxl的博客
07-08 577
学习自FFT详解。 很久前就想学，然而一直不能理解，这两天稍微懂了一些。 含义及用途 FFT（Fast Fourier Transformation）是离散傅氏变换（DFT）的快速算法。即为快速傅氏变换。它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的，把多项式乘法的复杂度从O(n2)O(n2)O(n^2)降到了O(nlogn)O(nlog⁡n)O(n\log...
[Luogu P4173] [BZOJ 4259] 残缺的字符串
LPA的博客
09-29 251
洛谷传送门 BZOJ传送门 题目描述 很久很久以前，在你刚刚学习字符串匹配的时候，有两个仅包含小写字母的字符串AAA和BBB，其中AAA串长度为mmm，BBB串长度为nnn。可当你现在再次碰到这两个串时，这两个串已经老化了，每个串都有不同程度的残缺。 你想对这两个串重新进行匹配，其中AAA为模板串，那么现在问题来了，请回答，对于BBB的每一个位置iii，从这个位置开始连续mmm个字符形成的子串是否...
洛谷 P4173 残缺的字符串 (FFT)
G19850917的博客
09-23 232
题目链接：P4173 残缺的字符串 题意 给定长度为 \(m\) 的模式串和长度为 \(n\) 的目标串，两个串都带有通配符，求所有匹配的位置。 思路 FFT 带有通配符的字符串匹配问题。 设模式串为 \(p\)，目标串为 \(t\)，将两个串的内容都根据字母先后顺序映射到 \(1\) 到 \(26\)。 如果不带有通配符，那么 \(t\) 以第 \(k\) 位结束的长度为 \(|p|\...
洛谷 P4173 残缺的字符串 fft
liangzihao1的博客
08-24 273
题目描述 很久很久以前，在你刚刚学习字符串匹配的时候，有两个仅包含小写字母的字符串 AAA 和 BBB ，其中 AAA 串长度为 mmm ， BBB 串长度为 nnn 。可当你现在再次碰到这两个串时，这两个串已经老化了，每个串都有不同程度的残缺。 你想对这两个串重新进行匹配，其中 AAA 为模板串，那么现在问题来了，请回答，对于 BBB 的每一个位置 iii ，从这个位置开始连续 mmm 个字...
洛谷4173（fft带通配符字符串匹配）
llmxby
11-23 660
题目链接：https://www.luogu.org/problemnew/show/P4173 思路：就是FFT的匹配，具体可以去洛谷看题解（PS：第二个大佬的题解写的挺好的emmm） #include &lt;cstdio&gt; #include &lt;cstdlib&gt; #include &lt;cstring&gt; #include &lt;bitset&gt; #incl...
BZOJ 4503: 两个串 FFT 通配符匹配
BlackJack
08-06 1062
4503: 两个串 Time Limit: 20 Sec  Memory Limit: 256 MB Submit: 630  Solved: 277 [Submit][Status][Discuss] Description 兔子们在玩两个串的游戏。给定两个字符串S和T，兔子们想知道T在S中出现了几次， 分别在哪些位置出现。注意T中可能有“?”字符，这个字符可以匹配任何字符。
[openjudge]带通配符的字符串匹配（dp）
zyf2000
10-23 783
题目描述传送门题解感觉这道题刚开始令我迷惑的是一个*可能匹配多个字符，这样的话看起来没法搞啊。 实际上不用想那么多。只需要考虑从哪里转移就可以了。 f(i,j)表示第一个串的前i个字符能否匹配第二个串的前j个字符。那么可以进行分类讨论： ①如果s1[i]=’?’或者s1[i]=s2[j]的话，f(i,j)可以从f(i-1,j-1)转移过来。 ②如果s1[i]=∗*的话，又有几种情况。假设*代
BZOJ3513[MUTC2013]idiots——FFT+生成函数
weixin_33862993的博客
02-18 128
题目描述 给定n个长度分别为a_i的木棒，问随机选择3个木棒能够拼成三角形的概率。 输入 第一行T(T<=100)，表示数据组数。 接下来若干行描述T组数据，每组数据第一行是n，接下来一行有n个数表示a_i。 3≤N≤10^5,1≤a_i≤10^5 输出 T行，每行一个整数，四舍五入保留7位小数。 样例输入 2 4 1 3 3 4 4 2 3...
FFT
weixin_33929309的博客
10-25 85
1 问题描述多项式的系数表示法：设关于x的n-1次多项式，次数为i的系数为ai。 \[A\left( x \right) = \sum\limits_{j = 0}^{n - 1} {{a_j}{x^j}} \] 多项式的点值表示法：用n个点可以确定一条曲线， \[\left( {{x_0},{y_0}} \right)\left( {{x_1},{y_1}} \right) \ldots \le...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值