投资管理有限公司注册要求,石家庄网络seo推广,福州网站制作推广,做素材类的网站赚钱吗传送门 文章目录题意#xff1a;思路#xff1a;题意#xff1a;
给你两个串a,ba,ba,b长度分别为n,mn,mn,m#xff0c;你需要输出m1m1m1个数#xff0c;第iii个数表示当允许有i−1i-1i−1个数可以不匹配时aaa中长度为mmm的子串与bbb匹配的数量#xff0c;匹配的意思就是…传送门
文章目录题意思路题意
给你两个串a,ba,ba,b长度分别为n,mn,mn,m你需要输出m1m1m1个数第iii个数表示当允许有i−1i-1i−1个数可以不匹配时aaa中长度为mmm的子串与bbb匹配的数量匹配的意思就是可以有i−1i-1i−1个位置不同其他位置相同。 n≤2e5,a∈(0,1,...,9,∗)n\le2e5,a\in {(0,1,...,9,*)}n≤2e5,a∈(0,1,...,9,∗)其中∗*∗代表通配符即与任何其他字符相同。
思路
不多bb直接上思路 由于字符集很小考虑每个字符对答案的影响。 考虑aibja_ib_jaibj那么在aaa中以im−jim-jim−j的位置为结尾的长度为mmm的子串中他的贡献是111由于其满足im−jkim-jkim−jk所以容易想到用FFTFFTFFT快速计算。为了处理m−jm-jm−j考虑将其reversereversereverse一下这样其实就变成了ijkijkijk直接FFTFFTFFT卷一下就可以啦。 假设我们已经卷出来fff数组了该怎么用它呢 先考虑没有通配符的情况我们可以发现m−f[i]m-f[i]m−f[i]就是他不能匹配的数量所以给ans[m−f[i]]ans[m-f[i]]ans[m−f[i]]让后求一个前缀和即可。 对于通配符我们考虑容斥的来求即acntbcnt−fcnta_{cnt}b_{cnt}-f_{cnt}acntbcnt−fcntfcntf_{cnt}fcnt即对∗*∗卷一下得到的答案。 有一个小细节可以优化就是我们发现最终卷完只有指数为nnn以内的项有用所以在卷的时候卷到nnn即可。
// Problem: C - Forgiving Matching
// Contest: Virtual Judge - 2021多校第三场补题
// URL: https://vjudge.net/contest/449636#problem/C
// Memory Limit: 524 MB
// Time Limit: 12000 ms
//
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize(Ofast,no-stack-protector,unroll-loops,fast-math)
//#pragma GCC target(sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tunenative)
#pragma GCC optimize(2)
#includecstdio
#includeiostream
#includestring
#includecstring
#includemap
#includecmath
#includecctype
#includevector
#includeset
#includequeue
#includealgorithm
#includesstream
#includectime
#includecstdlib
#includerandom
#includecassert
#define X first
#define Y second
#define L (u1)
#define R (u1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].ltr[u].r)1)
#define Len(u) (tr[u].r-tr[u].l1)
#define random(a,b) ((a)rand()%((b)-(a)1))
#define db puts(---)
using namespace std;//void rd_cre() { freopen(d://dp//data.txt,w,stdout); srand(time(NULL)); }
//void rd_ac() { freopen(d://dp//data.txt,r,stdin); freopen(d://dp//AC.txt,w,stdout); }
//void rd_wa() { freopen(d://dp//data.txt,r,stdin); freopen(d://dp//WA.txt,w,stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pairint,int PII;const int N1000010,mod1e97,INF0x3f3f3f3f;
const double eps1e-6,PIacos(-1);int n,m;
int rev[N];
int bit,limit;
int f[N],pre[N],ans[N];
char s1[N],s2[N];
double p1[N],p2[N];struct Complex {double x,y;Complex operator (const Complex t) const { return {xt.x,yt.y}; }Complex operator - (const Complex t) const { return {x-t.x,y-t.y}; }Complex operator * (const Complex t) const { return {x*t.x-y*t.y,x*t.yy*t.x}; }
}a[N],b[N];void fft(Complex a[],int inv) {for(int i0;ilimit;i) if(irev[i]) swap(a[i],a[rev[i]]);for(int mid1;midlimit;mid1) {Complex w1Complex({p1[mid],inv*p2[mid]});for(int i0;ilimit;imid*2) {Complex wkComplex({1,0});for(int j0;jmid;j,wkwk*w1) {Complex xa[ij],ywk*a[ijmid];a[ij]xy; a[ijmid]x-y;}}}
}int main()
{
// ios::sync_with_stdio(false);
// cin.tie(0);int _; scanf(%d,_);while(_--) {int cnt0;scanf(%d%d%s%s,n,m,s1,s2);pre[0]0; ans[m]0; ans[m1]0;for(int i0;in;i) {if(i0) pre[i]pre[i-1];pre[i]s1[i]*;if(im) cnts2[i]*;ans[i]0; f[i]0;}reverse(s2,s2m);bit0;while((1bit)n) bit;limit1bit;for(int i0;ilimit;i) rev[i](rev[i1]1)|((i1)(bit-1));for(int mid1;midlimit;mid1) p1[mid]cos(PI/mid),p2[mid]sin(PI/mid);for(int op0;op10;op) {char nowop0;if(op10) now*;for(int i0;ilimit;i) a[i]b[i]{0,0};for(int i0;in;i) a[i]{s1[i]now? 1.0:0,0};for(int i0;im;i) b[i]{s2[i]now? 1.0:0,0};fft(a,1); fft(b,1);for(int i0;ilimit;i) a[i]a[i]*b[i];fft(a,-1);for(int im-1;in;i) {int val(int)(a[i].x/limit0.5);if(op10) f[i]val;else f[i]-val;}}for(int im-1;in;i) {int allpre[i];if(im-1) all-pre[i-m];allcntf[i];ans[m-all];}for(int i1;im;i) ans[i]ans[i-1];for(int i0;im;i) printf(%d\n,ans[i]);}return 0;
}
/**/
