[Luogu4173]

题解

\(1.\)定义匹配函数

\(2.\)定义完全匹配函数

\(3.\)快速计算每一位的完全匹配函数值

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#define debug(...) fprintf(stderr,__VA_ARGS__)
#define Debug(x) cout<<#x<<"="<<x<<endl
using namespace std;
typedef long long LL;
const int INF=1e9+7;
inline LL read(){
register LL x=0,f=1;register char c=getchar();
while(c<48||c>57){if(c=='-')f=-1;c=getchar();}
while(c>=48&&c<=57)x=(x<<3)+(x<<1)+(c&15),c=getchar();
return f*x;
}
inline char readc(){
register char c=getchar();
while(c==' '||c=='\n'||c=='\t') c=getchar();
return c;
} const int MAXN=2e6+5;
const double Pi=acos(-1); struct cmpx{
double x,y;
inline cmpx(){}
inline cmpx(double _x,double _y){x=_x,y=_y;}
inline friend cmpx operator + (cmpx a,cmpx b){return cmpx(a.x+b.x,a.y+b.y);}
inline friend cmpx operator - (cmpx a,cmpx b){return cmpx(a.x-b.x,a.y-b.y);}
inline friend cmpx operator * (cmpx a,cmpx b){return cmpx(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);}
inline friend cmpx operator * (cmpx a,double b){return cmpx(a.x*b,a.y*b);}
}A[MAXN],B[MAXN],C[MAXN]; char s[MAXN],t[MAXN];
int ans[MAXN],a[MAXN],b[MAXN];
int n,m; namespace F_F_T{
int rev[MAXN],limit,l;
inline void init(int n){
for(limit=1,l=0;limit<=n;limit<<=1) l++;
for(int i=0;i<limit;i++)
rev[i]=(rev[i>>1]>>1)|((i&1)<<(l-1));
}
inline void FFT(cmpx *A,int type){
for(int i=0;i<limit;i++)
if(i<rev[i]) swap(A[i],A[rev[i]]);
for(int len=1;len<limit;len<<=1){
cmpx Wn=(cmpx){cos(Pi/len),type*sin(Pi/len)};
for(int i=0;i<limit;i+=(len<<1)){
cmpx w=(cmpx){1,0};
for(int j=0;j<len;j++,w=w*Wn){
cmpx x=A[i+j],y=w*A[i+len+j];
A[i+j]=x+y;
A[i+len+j]=x-y;
}
}
}
if(type==-1){
for(int i=0;i<limit;i++) A[i].x/=limit;
}
}
}using namespace F_F_T; int main(){
n=read(),m=read();
init(n+m-2);
for(int i=n-1;i>=0;i--){
char c=readc();
a[i]=(c!='*')?c-'a'+1:0;
}
for(int i=0;i<m;i++){
char c=readc();
b[i]=(c!='*')?c-'a'+1:0;
} for(int i=0;i<n;i++) A[i]=cmpx(a[i]*a[i]*a[i],0);
for(int i=0;i<m;i++) B[i]=cmpx(b[i],0);
FFT(A,1);FFT(B,1);
for(int i=0;i<limit;i++) C[i]=A[i]*B[i]; for(int i=0;i<limit;i++) A[i]=B[i]=cmpx(0,0);
for(int i=0;i<n;i++) A[i]=cmpx(a[i],0);
for(int i=0;i<m;i++) B[i]=cmpx(b[i]*b[i]*b[i],0);
FFT(A,1);FFT(B,1);
for(int i=0;i<limit;i++) C[i]=C[i]+A[i]*B[i]; for(int i=0;i<limit;i++) A[i]=B[i]=cmpx(0,0);
for(int i=0;i<n;i++) A[i]=cmpx(a[i]*a[i],0);
for(int i=0;i<m;i++) B[i]=cmpx(b[i]*b[i],0);
FFT(A,1);FFT(B,1);
for(int i=0;i<limit;i++) C[i]=C[i]-A[i]*B[i]*2.0; FFT(C,-1); for(int i=n-1;i<m;i++){
if(fabs(C[i].x)<0.5) ans[++ans[0]]=i-n+2;
//printf("%.3lf\n",C[i].x);
}
printf("%d\n",ans[0]);
for(int i=1;i<=ans[0];i++)
printf("%d ",ans[i]);
}

P4173 残缺的字符串(FFT)的更多相关文章

  1. Luogu P4173 残缺的字符串-FFT在字符串匹配中的应用

    P4173 残缺的字符串 FFT在字符串匹配中的应用. 能解决大概这种问题: 给定长度为\(m\)的A串,长度为\(n\)的B串.问A串在B串中的匹配数 我们设一个函数(下标从\(0\)开始) \(C ...

  2. P4173 残缺的字符串(FFT字符串匹配)

    P4173 残缺的字符串(FFT字符串匹配) P4173 解题思路: 经典套路将模式串翻转,将*设为0,设以目标串的x位置匹配结束的匹配函数为\(P(x)=\sum^{m-1}_{i=0}[A(m-1 ...

  3. 洛谷 P4173 残缺的字符串 (FFT)

    题目链接:P4173 残缺的字符串 题意 给定长度为 \(m\) 的模式串和长度为 \(n\) 的目标串,两个串都带有通配符,求所有匹配的位置. 思路 FFT 带有通配符的字符串匹配问题. 设模式串为 ...

  4. P4173 残缺的字符串 fft

    题意:给你两个字符串,问你第一个在第二个中出现过多少次,并输出位置,匹配时是模糊匹配*可和任意一个字符匹配 题解:fft加速字符串匹配; 假设上面的串是s,s长度为m,下面的串是p,p长度为n,先考虑 ...

  5. luoguP4173 残缺的字符串 FFT

    luoguP4173 残缺的字符串 FFT 链接 luogu 思路 和昨天做的题几乎一样. 匹配等价于(其实我更喜欢fft从0开始) \(\sum\limits_{i=0}^{m-1}(S[i+j]- ...

  6. BZOJ 4259: 残缺的字符串 [FFT]

    4259: 残缺的字符串 题意:s,t,星号任意字符,匹配方案数 和上题一样 多乘上一个\(a_{j+i}\)就行了 #include <iostream> #include <cs ...

  7. 【BZOJ4259】残缺的字符串 FFT

    [BZOJ4259]残缺的字符串 Description 很久很久以前,在你刚刚学习字符串匹配的时候,有两个仅包含小写字母的字符串A和B,其中A串长度为m,B串长度为n.可当你现在再次碰到这两个串时, ...

  8. 洛谷P4173 残缺的字符串(FFT)

    传送门 话说为什么字符串会和卷积扯上关系呢……到底得脑洞大到什么程度才能想到这种东西啊……大佬太珂怕了…… 因为通配符的关系,自动机已经废了 那么换种方式考虑,如果两个字符串每一位对应的编码都相等,那 ...

  9. [Luogu P4173]残缺的字符串 ( 数论 FFT)

    题面 传送门:洛咕 Solution 这题我写得脑壳疼,我好菜啊 好吧,我们来说正题. 这题.....emmmmmmm 显然KMP类的字符串神仙算法在这里没法用了. 那咋搞啊(或者说这题和数学有半毛钱 ...

随机推荐

  1. 用JS 写一个简单的程序,切换七彩盒子背景

    <!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <title> ...

  2. redis的特性

  3. ionic 页面乱码

    是文本编辑器的问题.文本编辑器保存时默认保存成Encoding:ANSI编码格式,保存成utf-8就好了

  4. iOS  地图(自定义地位图标)

    #import <UIKit/UIKit.h> @interface AppDelegate : UIResponder <UIApplicationDelegate> @pr ...

  5. SRA数据转成fastq

    Downloading and installing the SRA Toolkit step1: 下载并安装SRAtoolkit    (Download the Toolkit from the ...

  6. loj10099 矿场搭建

    传送门 分析 我们发现可以将这张图转换为一个联通块来处理.我们求出所有的割点.在求完之后我们我们对于每一个点双连通分量如果它没有割点相连则需要布置两个出口,因为可能有一个出口正好被割掉.而如果有一个割 ...

  7. 构建一个在线ASCII视频流服务

    构建一个在线ASCII视频流服务 2018-03-26  正常的文章 1685 什么是ASCII视频流服务? 其实这个名字是咱胡乱起的,具体叫啥我也不清楚,但效果如下: 大家可以在自己的命令行里试下, ...

  8. Inheritance with EF Code First: Part 3 – Table per Concrete Type (TPC)

    Inheritance with EF Code First: Part 3 – Table per Concrete Type (TPC) This is the third (and last) ...

  9. 《Effective Java》第10章 发并

    第66条:同步访问共享的可变数据 Java语言规范保证读或者写一个变量是原子的(atomic ) ,除非这个变量的类型为long或者double. [java中long和double类型操作的非原子性 ...

  10. [译]Javascript中的do-while循环

    本文翻译youtube上的up主kudvenkat的javascript tutorial播放单 源地址在此: https://www.youtube.com/watch?v=PMsVM7rjupU& ...