BZOJ 4259 殘缺的字串 FFT

2022-05-08 12:48:10 字數 1063 閱讀 4855

很久很久以前,在你剛剛學習字串匹配的時候,有兩個僅包含小寫字母的字串a和b,其中a串長度為m,b串長度為n。可當你現在再次碰到這兩個串時,這兩個串已經老化了,每個串都有不同程度的殘缺。

你想對這兩個串重新進行匹配,其中a為模板串,那麼現在問題來了,請回答,對於b的每乙個位置i,從這個位置開始連續m個字元形成的子串是否可能與a串完全匹配?

第一行包含兩個正整數m,n(1<=m<=n<=300000),分別表示a串和b串的長度。

第二行為乙個長度為m的字串a。

第三行為乙個長度為n的字串b。

兩個串均僅由小寫字母和*號組成,其中*號表示相應位置已經殘缺。

第一行包含乙個整數k,表示b串中可以完全匹配a串的位置個數。

若k>0,則第二行輸出k個正整數,從小到大依次輸出每個可以匹配的開頭位置(下標從1開始)。

題解:構造多項式 $c$,設 $c_=\sum_^(a-b_)^2\times a_\times b$.

不難發現,如果 $c_=0$,則 $b_$ 開始的一段長度為 $m$ 的字串會與 $a$ 串對應完全相等.

將式子展開,得 $c_=\sum_^b_^+a_^b_-2a_^b_^$.

將 $a$ 下標反轉,來三次 fft,最終加和即可. 

#include#define maxn 1200000       

using namespace std;

void setio(string s)

#define pi 3.1415926535898

struct cpx

cpx operator+(const cpx b)

cpx operator-(const cpx b)

cpx operator*(const cpx b)

}a1[maxn],a2[maxn],a3[maxn],b1[maxn],b2[maxn],b3[maxn];

void fft(cpx *a,int n,int flag)

for(int mid=1;midq;

int main()

return 0;

}

BZOJ4259 殘缺的字串

給出兩個字串,乙個模式串,乙個匹配串,問匹配串中哪些位置可以匹配上模式串,其中 可以作任意字元。這個可以轉化為多項式,我們可以把 看作0,其他字母看作各個數字,然後發現如果兩個字串相同,當且僅當 i 0n 1 a i b i 2 a i b i 0 sum 0 i 0n 1 a i b i 2 a ...

BZOJ4259 殘缺的字串

其實大部分字串的題都可以用多項式來想,包括這道題。於是,我們可以嘗試去構造兩個多項式,使其乘後的係數為0即可。相等為0,那麼我們可以用減法表示。可 可以匹配所有的符號,我們又該咋辦?那不就相當於乘個0嘛。於是,我們得到了下式 然後就可以卷積計算了。include include include in...

bzoj 4259 殘缺的字串

這題好神啊,居然是fft,表示一直在往資料結構上想。把 當成0,那麼兩個串可以匹配當且僅當 sum a i b i 2 times a i times b i 0 我們可以把平方拆開,然後就變成了幾個乘積相加的形式,那就大力翻轉乙個串然後跑fft。因為最開始mle了所以複製貼上了好多東西。1 inc...