BZOJ4503 两个串

18 篇文章 0 订阅
4 篇文章 0 订阅
本文深入探讨了一种利用FFT快速傅里叶变换进行字符串匹配的方法,通过将字符串转换为复数序列,采用FFT进行卷积计算,实现了高效地查找模式串在主串中的所有出现位置。文章详细介绍了算法原理,包括如何计算匹配函数f(x),并通过实例演示了如何使用FFT加速计算过程。
摘要由CSDN通过智能技术生成

原题链接:https://www.lydsy.com/JudgeOnline/problem.php?id=4503

两个串

Description

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

Input

两行两个字符串,分别代表S和T

Output

第一行一个正整数k,表示T在S中出现了几次。

接下来k行正整数,分别代表T每次在S中出现的开始位置。按照从小到大的顺序输出,S下标从0开始。

Sample Input

bbabaababaaaaabaaaaaaaabaaabbbabaaabbabaabbbbabbbbbbabbaabbbababababbbbbbaaabaaabbbbbaabbbaabbbbabab

a?aba?abba

Sample Output

0

HINT

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

题解

n = ∣ S ∣ , m = ∣ T ∣ n=|S|,m=|T| n=S,m=T,定义函数 f ( x ) = ∑ i = 0 m − 1 ( S x + i − T i ) 2 × T i f(x)=\sum_{i=0}^{m-1}(S_{x+i}-T_i)^2\times T_i f(x)=i=0m1(Sx+iTi)2×Ti,其中小写字母的值直接用 A C S I I ACSII ACSII码,’ ? ? ? ’ 的值为 0 0 0,这样如果 f ( x ) = 0 f(x)=0 f(x)=0,那么说明可以从 x x x开始匹配。

考虑如何计算 f ( x ) f(x) f(x),将 ( S x + i − T i ) 2 × T i (S_{x+i}-T_i)^2\times T_i (Sx+iTi)2×Ti展开,得到 S x + i 2 T i − 2 S x + i T i 2 + T i 3 S^2_{x+i}T_i-2S_{x+i}T^2_i+T_i^3 Sx+i2Ti2Sx+iTi2+Ti3,前两个为乘积形式, ∑ \sum 起来是平行的:

1.png

这样我们没法快速计算,所以我们把 T T T翻转过来:

2.png

运算变成了卷积形式,我们就可以用 F F T \mathcal{FFT} FFT快速计算了。

代码
#include<bits/stdc++.h>
#define db double
using namespace std;
const int M=(1<<18)+5;
const db pi=acos(-1.0);
struct cpx{db x,y;}s[M],t[M],s2[M],t2[M];
cpx operator +(cpx a,cpx b){return (cpx){a.x+b.x,a.y+b.y};}
cpx operator -(cpx a,cpx b){return (cpx){a.x-b.x,a.y-b.y};}
cpx operator *(cpx a,cpx b){return (cpx){a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x};}
int n,m,len,mx,rev[M],ans[M],tot,t3;
char ch[M];
void fft(cpx *f,int typ)
{
	cpx wn,w,x,y;int i,j,k,mid;
	for(i=0;i<mx;++i)if(i<rev[i])swap(f[i],f[rev[i]]);
	for(mid=1;mid<mx;mid<<=1)for(j=0,wn=(cpx){cos(pi/mid),typ*sin(pi/mid)};j<mx;j+=mid<<1)
	for(k=0,w=(cpx){1,0};k<mid;++k,w=w*wn)x=f[j+k],y=w*f[j+mid+k],f[j+k]=x+y,f[j+mid+k]=x-y;
}
void in()
{
	scanf("%s",ch);for(int i=n=strlen(ch)-1;i>=0;--i)s[i].x=ch[i]-'a'+1;
	scanf("%s",ch);for(int i=m=strlen(ch)-1;i>=0;--i)t[i].x=ch[i]=='?'?0:ch[i]-'a'+1;
}
void ac()
{
	for(mx=1;mx<=n+m;mx<<=1,++len);for(int i=0;i<mx;++i)rev[i]=rev[i>>1]>>1|((i&1)<<(len-1));
	reverse(t,t+m+1);
	for(int i=0;i<mx;++i)s2[i]=s[i]*s[i],t2[i]=(cpx){2,0}*t[i]*t[i],t3+=(t[i]*t[i]*t[i]).x;
	fft(s2,1),fft(t2,1),fft(s,1),fft(t,1);
	for(int i=0;i<mx;++i)s2[i]=s2[i]*t[i]-s[i]*t2[i];
	fft(s2,-1);
	for(int i=0;i<=n-m;++i)if(int(s2[m+i].x/mx+t3+0.5)==0)ans[++tot]=i;
	printf("%d\n",tot);for(int i=1;i<=tot;++i)printf("%d\n",ans[i]);
}
int main(){in();ac();}
BZOJ 2908 题目是一个数据下载任务。这个任务要求下载指定的数据文件,并统计文件小于等于给定整数的数字个数。 为了完成这个任务,首先需要选择一个合适的网址来下载文件。我们可以使用一个网络爬虫库,如Python的Requests库,来帮助我们完成文件下载的操作。 首先,我们需要使用Requests库的get()方法来访问目标网址,并将目标文件下载到我们的本地计算机。可以使用以下代码实现文件下载: ```python import requests url = '目标文件的网址' response = requests.get(url) with open('本地保存文件的路径', 'wb') as file: file.write(response.content) ``` 下载完成后,我们可以使用Python内置的open()函数打开已下载的文件,并按行读取文件内容。可以使用以下代码实现文件内容读取: ```python count = 0 with open('本地保存文件的路径', 'r') as file: for line in file: # 在这里实现对每一行数据的判断 # 如果小于等于给定整数,count 加 1 # 否则,不进行任何操作 ``` 在每一行的处理过程,我们可以使用split()方法将一行数据分割成多个字符串,并使用int()函数将其转换为整数。然后,我们可以将该整数与给定整数进行比较,以判断是否小于等于给定整数。 最后,我们可以将统计结果打印出来,以满足题目的要求。 综上所述,以上是关于解决 BZOJ 2908 数据下载任务的简要步骤和代码实现。 希望对您有所帮助。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

ShadyPi

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值