传送门

两种做法,一种是依次考虑每种字符,然后如果某个位置是该字符或者是\(*\)对应的值就是1,否则是0,然后把第一个串倒过来,fft卷积起来,最后看对应位置的值是否为m

然而上面那个做法在字符集大小过大的时候会GG,所以考虑一次性匹配,如果不考虑通配符\(*\),设(开头)位置i的匹配函数为\(f(i)=\sum_{j=1}^{m}(a_j-b_{i+m-j})^2\),\(f(i)=0\)说明i可以匹配

因为通配符可以匹配任何字符,也就是会让某组匹配一定是0,那么\(f(i)\)可以改成\(f(i)=\sum_{j=1}^{m}(a_j-b_{i+m-j})^2a_jb_{i+m-j}\),化简得到\(f(i)=\sum_{j=1}^{m}{a_j}^3b_{i+m-j}+a_j{b_{i+m-j}}^3-2{a_j}^2{b_{i+m-j}}^2\),然后三遍fft就好了

然后我ntt被fft吊着打qwq

// luogu-judger-enable-o2
#include<bits/stdc++.h>
#define LL long long
#define db double
#define il inline
#define re register using namespace std;
const int N=(1<<20)+10,mod=998244353;
il int rd()
{
int x=0,w=1;char ch=0;
while(ch<'0'||ch>'9') {if(ch=='-') w=-1;ch=getchar();}
while(ch>='0'&&ch<='9') {x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
return x*w;
}
int rdr[N];
il int fpow(int a,int b){int an=1;while(b){if(b&1) an=1ll*an*a%mod;a=1ll*a*a%mod,b>>=1;} return an;}
void ntt(int *a,int n,bool op)
{
int l=0,x,y;
while((1<<l)<n) ++l;
for(int i=0;i<n;++i) rdr[i]=(rdr[i>>1]>>1)|((i&1)<<(l-1));
for(int i=0;i<n;++i) if(i<rdr[i]) swap(a[i],a[rdr[i]]);
for(int i=1;i<n;i<<=1)
{
int ww=fpow(op?3:332748118,(mod-1)/(i<<1));
for(int j=0;j<n;j+=i<<1)
for(int k=0,w=1;k<i;++k,w=1ll*w*ww%mod)
x=a[j+k],y=1ll*a[j+k+i]*w%mod,a[j+k]=(x+y)%mod,a[j+k+i]=(x-y+mod)%mod;
}
if(!op) for(int i=0,w=fpow(n,mod-2);i<n;++i) a[i]=1ll*a[i]*w%mod;
}
int n,m,hh[N],a[N],b[N],an[N],ta;
char cc[N],ss[N]; int main()
{
n=rd(),m=rd();
int l=1;
while(l<=n+m) l<<=1;
scanf("%s%s",cc,ss);
for(int i=0;i<n;++i)
{
int x=cc[n-i-1]=='*'?0:cc[n-i-1]-96;
a[i]=x*x*x;
}
for(int i=0;i<m;++i)
{
int x=ss[i]=='*'?0:ss[i]-96;
b[i]=x;
}
ntt(a,l,1),ntt(b,l,1);
for(int i=0;i<l;++i) a[i]=1ll*a[i]*b[i]%mod;
ntt(a,l,0);
for(int i=1;i<=m;++i) hh[i]=(hh[i]+a[n+i-2])%mod;
for(int i=0;i<l;++i) a[i]=b[i]=0;
for(int i=0;i<n;++i)
{
int x=cc[n-i-1]=='*'?0:cc[n-i-1]-96;
a[i]=x;
}
for(int i=0;i<m;++i)
{
int x=ss[i]=='*'?0:ss[i]-96;
b[i]=x*x*x;
}
ntt(a,l,1),ntt(b,l,1);
for(int i=0;i<l;++i) a[i]=1ll*a[i]*b[i]%mod;
ntt(a,l,0);
for(int i=1;i<=m;++i) hh[i]=(hh[i]+a[n+i-2])%mod;
for(int i=0;i<l;++i) a[i]=b[i]=0;
for(int i=0;i<n;++i)
{
int x=cc[n-i-1]=='*'?0:cc[n-i-1]-96;
a[i]=x*x;
}
for(int i=0;i<m;++i)
{
int x=ss[i]=='*'?0:ss[i]-96;
b[i]=x*x;
}
ntt(a,l,1),ntt(b,l,1);
for(int i=0;i<l;++i) a[i]=1ll*a[i]*b[i]%mod;
ntt(a,l,0);
for(int i=1;i<=m;++i) hh[i]=((hh[i]-a[n+i-2]-a[n+i-2])%mod+mod)%mod;
for(int i=0;i<l;++i) a[i]=b[i]=0;
for(int i=1;i<=m-n+1;++i)
if(!hh[i]) an[++ta]=i;
printf("%d\n",ta);
for(int i=1;i<=ta;++i) printf("%d ",an[i]);
return 0;
}

luogu P4173 残缺的字符串的更多相关文章

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

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

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

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

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

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

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

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

  5. P4173 残缺的字符串 fft

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

  6. 洛谷 P4173 残缺的字符串

    (不知道xjb KMP可不可以做的说) (假设下标都以0开头) 对于有一定偏移量的序列的 对应位置 匹配或者数值计算的题,这里是有一种套路的,就是把其中一个序列翻转过来,然后卷积一下,所得到的新序列C ...

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

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

  8. 洛谷P4173 残缺的字符串

    题目大意: 两个带通配符的字符串\(a,b\),求\(a\)在\(b\)中出现的位置 字符串长度\(\le 300000\) 考虑魔改一发\(kmp\),发现魔改不出来 于是考虑上网搜题解 然后考虑\ ...

  9. P4173 残缺的字符串

    题目链接 题意分析 啥 ? ? ? \(FFT\)做字符串匹配 可是就是这样 我们定义匹配函数 我们定义\(A\)是匹配串 \(B\)是被匹配串 我们当前到达\(B\)串的\(x\)位置 \[P(x) ...

随机推荐

  1. Injection with CDI (Part I)

    官方参考:http://docs.jboss.org/weld/reference/latest/en-US/html/index.html https://antoniogoncalves.org/ ...

  2. 系统API执行没效果,可以检查一下是否与 360安全卫士 有关?!

    今天在写一个工具软件,使用到一个系统API: mouse_event() 代码如下: mouse_event(MOUSEEVENTF_ABSOLUTE | MOUSEEVENTF_MOVE,  x2, ...

  3. bigdata learning unit two--Spark environment setting

    1.下载 Spark安装之前的准备 文件的解压与改名 tar -zxvf spark-2.2.0-bin-hadoop2.7.tgz rm -rf spark-2.2.0-bin-hadoop2.7. ...

  4. HTML学习笔记Day15

    一.CSS3渐变 (一).CSS3渐变(gradient)可以让你在两个或多个指定的颜色之间显示平稳的过度:渐变效果比使用图片在放大时看起来效果更好,因为渐变(gradient)是由浏览器生成的 1. ...

  5. (转载)python: getopt的使用;

    注: 该文转载于https://blog.csdn.net/tianzhu123/article/details/7655499python中 getopt 模块, 该模块是专门用来处理命令行参数的 ...

  6. day-04(jquery)

    回顾: js: 组成部分: ECMAScript BOM DOM 变量声明 var 变量名=初始化值; 数据类型: 原始类型 Undefined Null String Number Boolean: ...

  7. 列举一些 MacBook Pro 必需的外设和应用程序推荐

    来源:知乎 文章收录于:风云社区SCOEE,提供上千款mac软件下载 基于从事Apps设计或开发者,使用 MacBook Pro,以下罗列一些必需的外设和应用程序推荐. Retina 256GB SS ...

  8. 干货分享:互联网运营 学习SEO从零开始 SEO深度解析学习笔记

    最近在自学SEO,互联网运营,把做的笔记干货分享给大家啊! 希望能帮到大家,如有好的建议可以关注我[磨人的小妖精]或留言,大家一起探讨. 之前还写过一篇文章互联网运营+SEO:推荐必看的5本书籍,学习 ...

  9. bzoj1233 单调队列优化dp

    https://www.lydsy.com/JudgeOnline/problem.php?id=1233 数据结构优化dp的代码总是那么抽象 题意:奶牛们讨厌黑暗. 为了调整牛棚顶的电灯的亮度,Be ...

  10. chrome浏览器安装不上的惨痛经历

    项目在赶进度的时候,手贱把chrome的一些文件删除了,整个浏览器都没法打开 决定重装下,但是连卸载的功能都打不开了 上网重新下载了个安装包,发现安装包都打不来 很绝望,查了很多资料 很多人说要删除注 ...