<output id="qn6qe"></output>

    1. <output id="qn6qe"><tt id="qn6qe"></tt></output>
    2. <strike id="qn6qe"></strike>

      亚洲 日本 欧洲 欧美 视频,日韩中文字幕有码av,一本一道av中文字幕无码,国产线播放免费人成视频播放,人妻少妇偷人无码视频,日夜啪啪一区二区三区,国产尤物精品自在拍视频首页,久热这里只有精品12

      多項(xiàng)式算法學(xué)習(xí)筆記

      來學(xué)習(xí)一個(gè)多項(xiàng)式全家桶。

      基本算法

      FFT

      先咕著。

      NTT

      先咕著。

      數(shù)學(xué)

      多項(xiàng)式求導(dǎo)

      這其實(shí)是高中數(shù)學(xué)內(nèi)容。

      對于多項(xiàng)式 \(F(x)=u(x)+v(x)\),它的導(dǎo)數(shù)\(F'(x)=u'(x)+v'(x)\)

      再結(jié)合\(u(x)=x^n, u'(x)=nx^{n-1}\),就可以愉快地線性時(shí)間求出導(dǎo)數(shù)。

      多項(xiàng)式積分

      直接根據(jù)公式

      \[\int F(x)dx = (\sum^{n}_{i=0}\frac{f_i}{i+1}x^{i+1})+c \]

      那么右邊就是

      \[f_0x+\frac{f_1}{2}x^2+\frac{f_2}{3}x^3+\frac{f_3}{4}x^4+...+\frac{f_n}{n+1}x^{n+1}+c \]

      原函數(shù)

      \[F(x)=f_0+f_1x+f_2x^2+f_3x^3+...+f_nx^n \]

      \(c\)咋整啊,直接等 \(0\) 吧。

      那么就可以線性做完了。

      多項(xiàng)式求逆

      對于函數(shù) \(F(x)\),求一個(gè)多項(xiàng)式 \(G\)(x),使得在每一項(xiàng)系數(shù)模 \(x^n\) 時(shí),有\(F(x)*G(x) \equiv 1 \pmod {x^n}\)

      如果 \(F\) 只有一項(xiàng),那么就變成了單項(xiàng)式求逆元。

      如果有 \(n\) 項(xiàng)呢?

      假設(shè)我們先求出來了模 \(x^{\lceil \frac{n}{2} \rceil}\) 的逆為 \(G'\),則:

      \[F*G' \equiv 1 \pmod {x^{\lceil \frac{n}{2} \rceil}} \]

      \[F*G \equiv 1 \pmod {x^{\lceil \frac{n}{2} \rceil}} \]

      上減下,除掉 \(F\)

      \[G'-G \equiv 0 \pmod {x^{\lceil \frac{n}{2} \rceil}} \]

      平方

      \[(G'-G)^2 \equiv 0 \pmod {x^n} \]

      拆開

      \[G'^2-2GG'+G^2 \equiv 0 \pmod {x^n} \]

      乘上 \(F\)

      \[FG'^2-2G'+G \equiv 0\pmod{x^n} \]

      移項(xiàng)

      \[G \equiv 2G'-FG'^2 \pmod{x^n} \]

      遞歸算就行了。

      多項(xiàng)式對數(shù)函數(shù)

      多項(xiàng)式 \(F\),求 \(G \equiv ln(F) \pmod {x^n}\)

      設(shè) \(A(x)=ln(x)\),有

      \[G(x)=A(F(x)),A'(x)=\frac{1}{x} \]

      求導(dǎo)

      \[G'(x)=A'(F(x))F'(x)=\frac{F'(x)}{F(x)} \]

      求導(dǎo)求逆乘一起就行了。

      泰勒展開

      求一個(gè)函數(shù) \(f(x)\) 在某一點(diǎn)的值,可以構(gòu)造一個(gè)函數(shù) \(g(x)\),使得

      \[f(x) \approx g(x) = g(0) + \frac{f^1(0)}{1!}x + \frac{f^2(0)}{2!}x^2 + ... + \frac{f^n(0)}{n!}x^n \]

      其中 \(f^n(0)\) 表示對原函數(shù)的圖像上 \(0\) 這個(gè)點(diǎn)進(jìn)行 \(n\) 階求導(dǎo)。

      牛頓迭代

      求函數(shù) \(G(x)\) 的零點(diǎn),即求滿足 \(G(F(z)) \equiv 0 \pmod {z^n}\) 的多項(xiàng)式 \(F(z)\)

      \(n=1\) 時(shí),提前求出\(G(F(z)) \equiv 0 \pmod {z^n}\)

      現(xiàn)在假設(shè)求出了

      \[G(F_0(z)) \equiv 0 \pmod{z^{ \lceil \frac{n}{2} \rceil}} \]

      進(jìn)行泰勒展開

      \[G(F(z))=G(F_0(z))+\frac{G'(F_0(z))}{1!}(F(z)-F_0(z))+\frac{G''(F_0(z))}{2!}(F(z)-F_0(z))^2+... \]

      因?yàn)?\(F(z)\)\(F_0(z)\) 的最后 \(\lceil \frac{n}{2} \rceil\) 項(xiàng)相同,所以 \((F(z)-F_0(z))^2\) 的最低非 \(0\) 項(xiàng)次數(shù)大于 $2 \lceil \frac{n}{2} \rceil $,所以

      \[G(F(z)) \equiv G(F_0(z))+G'(F_0(z))+G'(F_0(z))(F(z)-F_0(z)) \pmod{z^n} \]

      \(G(F(z)) \equiv 0 \pmod{z^n}\),得到

      \[F(z) \equiv F_0(z) - \frac{G(F_0(z))}{G'(F_0(z))} \pmod{z^n} \]

      迭代就可以了。

      多項(xiàng)式指數(shù)函數(shù)

      \(F(x)\equiv e^{A(x)} \pmod{x^n}\)

      取對數(shù)并移項(xiàng)

      \[\ln F(x) - A(x)=0 \]

      設(shè) \(G(F(x)) = \ln F(x) - A(x)\),則求 \(G(F(x))=0\)

      求導(dǎo)

      \[G'(F(x)) = \frac{1}{F(x)} \]

      假設(shè)已經(jīng)求出了 \(F_0(x) \equiv e^{A(x)} \pmod{x^{\lceil\frac{n}{2}\rceil}}\),代入牛頓迭代公式

      \[F(x) \equiv F_0(x) - \frac{G(F_0(x))}{G'(F_0(x))} \pmod{x^n} \\ F(x) \equiv F_0(x)(1 - \ln F_0(x) + A(x)) \pmod{x^n} \]

      \(A(0)=0\),所以 \(F(x)\) 的常數(shù)項(xiàng)為 \(1\)

      也可以遞歸求解了。

      posted @ 2020-02-24 08:28  SamBillon  閱讀(267)  評論(0)    收藏  舉報(bào)
      主站蜘蛛池模板: 中文字幕不卡在线播放| 国产成人一区二区不卡| 亚洲中少妇久久中文字幕| 亚洲综合一区二区精品导航| 久久婷婷成人综合色综合| 色欲AV无码一区二区人妻| 国产极品精品自在线不卡| 麻豆精品久久久久久久99蜜桃| 久久国产成人av蜜臀| 99久久亚洲综合精品成人网| 亚洲精品无码你懂的网站| 国产v综合v亚洲欧美大天堂| 少妇高潮惨叫喷水在线观看| 激情一区二区三区成人文| 依依成人精品视频在线观看| 国产精品国产三级国快看| 亚洲国产欧美在线人成| 在线亚洲午夜理论av大片| 熟女一区| 人妻中文字幕精品一页| 国产网友愉拍精品视频手机| 国精品午夜福利视频| 欧美在线观看www| 久久精品女人的天堂av| 亚洲精品一区二区三区蜜臀| 成年视频人免费网站动漫在线| 99中文字幕国产精品| 线观看的国产成人av天堂| 国产亚洲精品久久久久5区| 色五月丁香五月综合五月| 国产丝袜在线精品丝袜不卡| 国产成人av电影在线观看第一页 | 国产日韩精品一区二区三区在线| 国产高颜值不卡一区二区| 国产精品中文字幕观看| 熟妇的奶头又大又长奶水视频 | 国产日韩欧美| 亚洲狼人久久伊人久久伊| 亚洲真人无码永久在线| 上司的丰满人妻中文字幕| 人妻少妇|