博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
[学习笔记] 形式幂级数
阅读量:6092 次
发布时间:2019-06-20

本文共 1313 字,大约阅读时间需要 4 分钟。

很早就开始学了吧 一直没有写学习笔记..

candy? 的博客写的很不错啊

多项式一系列操作复杂度\(T(n) = T(n/2) + \mathcal O(nlogn) = \mathcal O(nlogn)\)

多项式求逆

倍增的思想

已知\(A(x)B_0(x) \equiv 1 \pmod {x^n}\)

\(B(x)\)满足\(A(x)B(x) \equiv 1 \pmod {x^{2n}}\)

多项式有逆元的前提是常数项有逆元

两式相减 然后平方展开 左右同时乘上\(A(x)\)可以得到

\[B(x) \equiv B_0(x) (2 - B_0(x)A(x))\pmod {x^{2n}}\]

多项式开方

同样倍增

已知\(B_0^{2}(x) \equiv A(x) \pmod {x^n}\)

\(B(x)\)满足\(B^2(x) \equiv A(x) \pmod {x^{2n}}\)

\[ (B_0^{2}(x) - A(x)) ^ 2 \equiv 0 \pmod {x^{2n}}\]

\[ (B_0^{2}(x) + A(x)) ^ 2 \equiv 4 B_0^{2}(x) A(x) \pmod {x^{2n}}\]

\[ (\frac {B_0^{2}(x) + A(x)}{2B_0(x)}) ^ 2 \equiv A(x)\pmod {x^{2n}}\]

同样需要多项式求逆元

多项式求ln

\[B(x) = ln(A(x))\]

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

求一下导积分积回去

牛顿迭代

已知\(g(x)\)\(f(x)\) 满足 \(g(f(x)) = 0\)

已知前\(n\)\(f(x) \equiv f_0(x) \pmod {x^n}\)

\(g(f_0(x))\)泰勒展开 .. 这里感觉好强

得到

\[f(x) = f_0(x) - \frac {g(f_0(x))}{g'(f_0(x))}\]

多项式求exp

定义\(f(x) = e^{A(x)}\)

以及 \(g(f(x)) = ln(f(x)) - A(x) = 0\)

代入牛顿迭代式子中去

\[f(x) = f_0(x) - \frac {ln(f_0(x)) - A(x)}{1/(f_0(x))}\]

整理得

\[f(x) = f_0(x) (1 - ln(f_0(x)) + A(x))\]

bzoj3456

定义\(G(x)\) 为无向(不一定连通)图数目 有

\[G(x) = \sum_{n \ge 0} 2^{C_n^2} \frac {x^n} {n!}\]

根据指数生成函数 设\(F(x)\)为无向连通图数目 有

\[G(x) = \sum_{n \ge 0} \frac {F(x) ^ n} {n!}\]

因为无向图是由若干无向连通图组成的集合

于是\(G=e^{F}\) \(F = ln(G)\)

求一个ln就好了

转载于:https://www.cnblogs.com/foreverpiano/p/8831342.html

你可能感兴趣的文章
Tyvj 1728 普通平衡树
查看>>
[Usaco2015 dec]Max Flow
查看>>
javascript性能优化
查看>>
多路归并排序之败者树
查看>>
java连接MySql数据库
查看>>
转:Vue keep-alive实践总结
查看>>
android studio修改新项目package名称
查看>>
深入python的set和dict
查看>>
C++ 11 lambda
查看>>
Hadoop2.5.0 搭建实录
查看>>
实验吧 recursive write up
查看>>
High-speed Charting Control--MFC绘制图表(折线图、饼图、柱形图)控件
查看>>
go test命令參数问题
查看>>
linux 搜索文本
查看>>
超实用Mac软件分享(二)
查看>>
Android JSON数据解析
查看>>
DEV实现日期时间效果
查看>>
java注解【转】
查看>>
Oracle表分区
查看>>
centos 下安装g++
查看>>