bzoj4403 两个串
Description
Input
Output
两个串a,b相等(b中有通配符)当且仅当Σ(a[i]-b[i])2b[i]=0,其中a[i],b[i]为对应字符的对应编号,且通配符对应0
翻转S串并把T串用0补至与S等长后上式可化为卷积形式用fft计算
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
int N,P;
const double pi=3.14159265358979323846;
int r[];
struct cplx{
double a,b;
inline cplx(double r=,double i=){a=r;b=i;}
inline cplx operator+(cplx x){return cplx(a+x.a,b+x.b);}
inline cplx operator-(cplx x){return cplx(a-x.a,b-x.b);}
inline cplx operator*(cplx x){return cplx(a*x.a-b*x.b,a*x.b+b*x.a);}
}a[],b[];
void dft(cplx*a,int t){
for(int i=;i<N;i++)if(i>r[i])std::swap(a[i],a[r[i]]);
for(int i=;i<N;i<<=){
cplx w(cos(pi/i),t*sin(pi/i));
for(int j=;j<N;j+=i<<){
cplx e(),*b=a+j,*c=b+i;
for(int k=;k<i;k++,e=e*w){
cplx x=b[k],y=c[k]*e;
b[k]=x+y,c[k]=x-y;
}
}
}
}
char s1[],s2[];
int v1[],v2[];
int h1[],h2[];
int as[],ap=;
int l1,l2,h3=;
int main(){
input:{
scanf("%s%s",s1,s2);
l1=strlen(s1)-;
l2=strlen(s2)-;
for(int i=;i<=l1;i++)v1[i]=s1[l1-i]-'a'+;
for(int i=;i<=l2;i++)v2[i]=s2[i]!='?'?s2[i]-'a'+:;
for(int i=;i<=l2;i++)h3+=v2[i]*v2[i]*v2[i];
}
init:{
for(N=,P=;N<=l1+;N<<=,++P);
N<<=,++P;
for(int i=;i<N;i++)r[i]=r[i>>]>>|(i&)<<P;
}
calc:{
for(int i=;i<=l1;i++)a[i]=cplx(v1[i]*v1[i]);
for(int i=;i<=l1;i++)b[i]=cplx(v2[i]);
for(int i=l1+;i<N;i++)a[i]=b[i]=cplx();
dft(a,),dft(b,);
for(int i=;i<N;i++)a[i]=a[i]*b[i];
dft(a,-);
for(int i=;i<=l1;i++)h1[i]=(int)(a[i].a/N+0.5); for(int i=;i<=l1;i++)a[i]=cplx(v1[i]);
for(int i=;i<=l1;i++)b[i]=cplx(v2[i]*v2[i]);
for(int i=l1+;i<N;i++)a[i]=b[i]=cplx(,);
dft(a,),dft(b,);
for(int i=;i<N;i++)a[i]=a[i]*b[i];
dft(a,-);
for(int i=;i<=l1;i++)h2[i]=(int)(a[i].a/N+0.5);
}
output:{
for(int i=l1;i>=l2;i--)if(h1[i]-*h2[i]+h3==)as[ap++]=l1-i;
printf("%d\n",ap);
for(int i=;i<ap;i++)printf("%d\n",as[i]);
}
return ;
}
bzoj4403 两个串的更多相关文章
- bzoj 4503 两个串
Description 兔子们在玩两个串的游戏.给定两个字符串S和T,兔子们想知道T在S中出现了几次, 分别在哪些位置出现.注意T中可能有“?”字符,这个字符可以匹配任何字符. Input 两行两个字 ...
- BZOJ4503: 两个串
Description 兔子们在玩两个串的游戏.给定两个字符串S和T,兔子们想知道T在S中出现了几次, 分别在哪些位置出现.注意T中可能有“?”字符,这个字符可以匹配任何字符. Input 两行两个字 ...
- SPOJ 1811 Longest Common Substring (后缀自动机第一题,求两个串的最长公共子串)
题目大意: 给出两个长度小于等于25W的字符串,求它们的最长公共子串. 题目链接:http://www.spoj.com/problems/LCS/ 算法讨论: 二分+哈希, 后缀数组, 后缀自动机. ...
- BZOJ 4503: 两个串 [FFT]
4503: 两个串 题意:兔子们在玩两个串的游戏.给定两个只含小写字母的字符串S和T,兔子们想知道T在S中出现了几次, 分别在哪些位置出现.注意T中可能有"?"字符,这个字符可以匹 ...
- 【BZOJ4503】两个串(FFT)
[BZOJ4503]两个串(FFT) 题面 给定串\(S\),以及带通配符的串\(T\),询问\(T\)在\(S\)中出现了几次.并且输出对应的位置. \(|S|,|T|<=10^5\),字符集 ...
- bzoj4503: 两个串 bitset
目录 题目链接 题解 代码 题目链接 bzoj4503: 两个串 题解 暴一发bitset f[i][j] 表示 S[1..i] 是否有个后缀能匹配 T[1..j] 那么假设 S[i+1] 能匹配 T ...
- URAL 1517 Freedom of Choice (后缀数组 输出两个串最长公共子串)
版权声明:本文为博主原创文章,未经博主同意不得转载. https://blog.csdn.net/whyorwhnt/article/details/34075603 题意:给出两个串的长度(一样长) ...
- 【BZOJ 4503】4503: 两个串 (FFT)
4503: 两个串 Time Limit: 20 Sec Memory Limit: 256 MBSubmit: 497 Solved: 226 Description 兔子们在玩两个串的游戏.给 ...
- SPOJ 1811 Longest Common Substring(求两个串的最长公共子串 || 或者n个串)
http://www.spoj.com/problems/LCS/ 题目:求两个串的最长公共子串 参考:https://www.cnblogs.com/autoint/p/10345276.html: ...
随机推荐
- awk根据指定的字符串分割字符串
以从字符串"hello-kitty-red-for-you"中获取-for前面的内容为例: echo "hello-kitty-red-for-you" |aw ...
- WebApi Session支持
代码: WebApiConfig using System; using System.Collections.Generic; using System.Linq; using System.Net ...
- JMeter的下载以及安装使用
下载 https://jmeter.apache.org/download_jmeter.cgi 安装 无须安装,解压之后即可使用. 解压到C:\Program Files\apache-jmeter ...
- Exception.StackTrace
Exception中的StackTrace属性 执行堆栈跟踪在给定时刻正在执行的所有方法. 对方法调用的跟踪称为堆栈跟踪. 堆栈跟踪列表提供了一种循着调用堆叠跟踪到方法中异常发生处行号的手段.Stac ...
- 【max_result_window大小】 Result window is too large的问题
方法一: 如果需要搜索分页,可以通过from size组合来进行.from表示从第几行开始,size表示查询多少条文档.from默认为0,size默认为10, 如果搜索size大于10000,需要设置 ...
- JS浏览器BOM
浏览器对象模型 (BOM) BOM的核心是window,而window对象又具有双重角色,它既是通过js访问浏览器窗口的一个接口,又是一个Global(全局)对象.这意味着在网页中定义的任何对象,变 ...
- 解决Github Desktop Repo publish 失败问题
参考: src refspec master does not match any when pushing commits in git Github上传项目步骤和常见问题 Trying to gi ...
- 原生DOM操作vs框架虚拟DOM比较
1. 原生 DOM 操作 vs. 通过框架封装操作. 这是一个性能 vs. 可维护性的取舍.框架的意义在于为你掩盖底层的 DOM 操作,让你用更声明式的方式来描述你的目的,从而让你的代码更容易维护.没 ...
- data.table 中的动态作用域
data.table 中最常用的语法就是 data[i, j, by],其中 i.j 和 by 都是在动态作用域中被计算的.换句话说,我们不仅可以直接使用列,也可以提前定义诸如 .N ..I 和 .S ...
- vue双向数据绑定最最最最最简单直观的例子
vue双向数据绑定最最最最最简单直观的例子 一.总结 一句话总结:双向绑定既不仅model可以影响view的数据,view也可以影响model的数据 view model 数据 1.vue双向数据绑定 ...