bzoj4503 两个串

标签: bzoj
10人阅读 评论(0) 收藏 举报
分类:

Description


兔子们在玩两个串的游戏。给定两个字符串S和T,兔子们想知道T在S中出现了几次,
分别在哪些位置出现。注意T中可能有“?”字符,这个字符可以匹配任何字符。

S下标从0开始。
S 长度不超过 10^5, T 长度不会超过 S。 S 中只包含小写字母, T中只包含小写字母和“?”

来自 https://www.lydsy.com/JudgeOnline/problem.php?id=4503

Solution


原来fft还能这么用

我们定义两个长度为n的字符串距离为i=1ns2[i](s1[i]s2[i])2
其中当s2=’?’时为0,这样就能保证这是通配符惹

拆一下可以发现两个项可以FFT,另一个项直接求就好了

Code


#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <complex>

typedef std:: complex <double> com;
typedef double db;

const db pi=3.1415926535897932;
const int N=300005;

char str1[N],str2[N];

int a[N],b[N],rev[N],ans[N];

com c[N],d[N],e[N],f[N];

void FFT(com *a,int len,db f) {
    for (int i=0;i<len;i++) if (i<rev[i]) std:: swap(a[i],a[rev[i]]);
    for (int i=1;i<len;i*=2) {
        com wn(cos(pi/i),f*sin(pi/i));
        for (int j=0;j<len;j+=i*2) {
            com w(1,0);
            for (int k=0;k<i;k++) {
                com u=a[j+k],v=w*a[j+k+i];
                a[j+k]=u+v; a[j+k+i]=u-v;
                w*=wn;
            }
        }
    }
    if (f==-1) for (int i=0;i<len;i++) a[i]/=len;
}

int main(void) {
    scanf("%s",str1+1); scanf("%s",str2+1);
    int len1=strlen(str1+1),len2=strlen(str2+1),tmp=0;
    for (int i=1;i<=len1;i++) a[i]=str1[i]-'a'+1;
    for (int i=1;i<=len2;i++) b[len2-i+1]=(str2[i]=='?')?(0):(str2[i]-'a'+1);

    int len,lg; for (len=1,lg=0;len<=len1*2;len*=2,lg++);
    for (int i=0;i<=len;i++) rev[i]=(rev[i/2]/2)|((i&1)<<(lg-1));

    for (int i=1;i<=len2;i++) tmp+=b[i]*b[i]*b[i];
    for (int i=1;i<=len1;i++) c[i]=a[i]*a[i];
    for (int i=1;i<=len2;i++) d[i]=b[i];

    for (int i=1;i<=len1;i++) e[i]=a[i];
    for (int i=1;i<=len2;i++) f[i]=b[i]*b[i]*2;

    FFT(c,len,1); FFT(d,len,1); FFT(e,len,1); FFT(f,len,1);
    for (int i=0;i<len;i++) c[i]*=d[i],e[i]*=f[i];
    FFT(c,len,-1); FFT(e,len,-1);

    for (int i=0;i<len;i++) ans[i]=(int)(0.1+c[i].real())-(int)(0.1+e[i].real());
    int prt=0;
    for (int i=len2+1;i<=len1+1;i++) if (ans[i]+tmp==0) prt++;
    printf("%d\n", prt);
    for (int i=len2+1;i<=len1+1;i++) if (ans[i]+tmp==0) printf("%d\n", i-len2-1);
    return 0;
}
查看评论

BZOJ4503 两个串

orz xuruifan发discuss问题解…… 题解可以看这题discuss 至于如何用FFT算卷积,我们写一下高精乘,用字母代表每一位,发现乘完了其实就是卷积的形式 #includ...
  • neither_nor
  • neither_nor
  • 2016-04-30 13:49:36
  • 1319

bzoj 4503: 两个串 fft

没想到这种字符串题居然可以用fft。。。        设有两个串S1,S2,定义S1和S2的距离为:Σ(i=1,n) (s1[i]-s2[i])^2,那么两个串相同当且仅当距离=0。那么现在S2存在...
  • lych_cys
  • lych_cys
  • 2016-04-11 19:32:10
  • 1411

Bzoj4503:两个串:FFT,构造

题目链接:4503:两个串 听说这样的通配符匹配问题常用FFT搞一搞?好吧蒟蒻并不会 对于这道题如果没有通配符,我们构造一个函数 这样如果s1和s2相等当且仅当f[x]=0 但是这道题有通配符...
  • qq_34025203
  • qq_34025203
  • 2016-04-19 07:09:44
  • 1118

[BZOJ4503]两个串(FFT)

题目描述传送门题解受到上一道题的影响,这题想了很久… 因为这题的思路和上题真的是相似又互异的 令F(i)F(i)表示把小串结尾怼到大串的第i个字符上是否能匹配,0表示能匹配,非0表示不能匹配 对...
  • Clove_unique
  • Clove_unique
  • 2017-03-09 18:34:36
  • 296

Bzoj4503 两个串 FFT

问题描述给定两个串s1,s2,问你s2在s1中出现的次数和位置,s2中有‘?’字符,可以匹配任何字符。字符串长度为105{10}^{5}。题解我们考虑没有‘?’时,两个长度相等的字符串相等,当且仅当∑...
  • Man_God
  • Man_God
  • 2016-12-18 18:30:47
  • 136

[bzoj4503]两个串

4503: 两个串Time Limit: 20 Sec Memory Limit: 256 MB Submit: 190 Solved: 106 [Submit][Status][Discus...
  • FZHvampire
  • FZHvampire
  • 2016-06-02 15:17:04
  • 514

bzoj4503 两个串(dp+FFT)

Description兔子们在玩两个串的游戏。给定两个字符串S和T,兔子们想知道T在S中出现了几次,分别在哪些位置出现。注意T中可能有“?”字符,这个字符可以匹配任何字符。 Input两行两个字符串...
  • wu_tongtong
  • wu_tongtong
  • 2018-03-01 16:01:02
  • 57

bzoj4503两个串 快速傅里叶变换(FFT)

题目传送门:bzoj4503题目大意:给出两个串S和T,T中有通配符(可以与任意字符匹配),求T在S中匹配的次数及位置。题目一看,怎么做啊,kmp?AC自动机?好像很困难啊。 但是,我们转换一下思路...
  • Nickwzk
  • Nickwzk
  • 2017-02-22 19:40:58
  • 826

BZOJ4503:两个串 (构造+FFT/随机化+NTT)

权限题传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4503 题目大意:给你两个串S和T,询问T在S中出现了几次,以及在哪些位置出现。S...
  • KsCla
  • KsCla
  • 2018-02-17 02:08:52
  • 105

bzoj4503 两个串(FFT好题)

bzoj4503 两个串 原题地址:http://www.lydsy.com/JudgeOnline/problem.php?id=4503 题意: 兔子们在玩两个串的游戏。给定两个字符串S和...
  • Bfk_zr
  • Bfk_zr
  • 2018-02-04 22:28:51
  • 68
    个人资料
    持之以恒
    等级:
    访问量: 12万+
    积分: 8130
    排名: 3102
    乱七八糟
    文章分类
    最新评论