![](https://img-blog.csdnimg.cn/20201014180756923.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
多项式
Freopen
这个作者很懒,什么都没留下…
展开
-
模拟赛 我的朋友们(分治NTT)
题意:有一个长度为nnn的序列,序列上每个位置有一个物品,每个物品在每次被询问到的时候有pip_ipi的概率被拒绝,考虑这样一个过程:先取前LLL个物品,然后每次对这LLL个物品都询问,如果有xxx个物品被拒绝,则将前xxx个物品丢弃,然后从未被选取的物品中选取最靠前的xxx个物品重复这个过程,当总物品数<L\lt L<L时停止,问期望询问多少组次。n,L≤1e5n , L \leq 1e5n,L≤1e5, mod 998244353\bmod 998244353mod998244353原创 2020-08-03 19:57:37 · 478 阅读 · 0 评论 -
多项式多点求值的小常数解法
设MULT(f(z),g(z))=∑i=0zi∑j=0([zi+j]f(z))([zj]g(z))MUL^T(f(z),g(z)) = \sum_{i=0} z^i\sum_{j=0} ([z^{i+j}]f(z))([z^j]g(z))MULT(f(z),g(z))=∑i=0zi∑j=0([zi+j]f(z))([zj]g(z))也就是差卷积。可以发现F(x0)=∑i=0nx0i[xi]F(x)=[x0]MULT(F(x),11−x0z)F(x_0) =\sum_{i=0}^n x_0^i[x^原创 2020-07-28 22:47:38 · 463 阅读 · 0 评论 -
MTT
合并FFTFFTFFT:三次翻转一次:共轭单位根两次:代入共轭单位根后的值需要再取共轭三次:A(x)+iB(x)A(x) + i B(x)A(x)+iB(x)其中A(x)=p−q2A(x) = \frac {p - q} 2A(x)=2p−q ,则B(x)=iq−p2B(x) = i\frac {q - p}{2}B(x)=i2q−p。任意模数多项式卷积:#include<bits/stdc++.h>#define maxn 300005#define rep(i,j,k) f原创 2020-07-22 22:38:16 · 456 阅读 · 0 评论 -
多项式题单
LOJ #556. 「Antileaf’s Round」咱们去烧菜吧求混合背包(某一体积的物品可能有无限个也可能有有限个),得到体积和为1...n1...n1...n的方案数。热身题。如果一个体积为vvv的物品有无限个,那么其关于体积的普通生成函数为∑i=0xvi=11−xv\sum_{i=0}x^{vi} = \frac 1{1-x^v}∑i=0xvi=1−xv1如果一个体积为vvv的物品有KKK个,那么其关于体积的普通生成函数为∑i=0Kxvi=1−xv(K+1)1−xv\sum_{i=0}原创 2020-06-30 15:33:00 · 658 阅读 · 1 评论 -
无标号无根树计数
EulerEulerEuler变换ε(F(x))=∑i=1(∑j=0xij)[xi]F(x)\varepsilon(F(x)) = \sum_{i=1} (\sum_{j=0} x^{ij})^{[x^i]F(x)}ε(F(x))=∑i=1(∑j=0xij)[xi]F(x)可以用等比数列求和公式+ln\lnln+exp\expexp求得ε(F(x))=e∑i=1F(xi)i\varepsilon(F(x)) = e^{\sum_{i=1} \frac {F(x^i)}{i}}ε(F(x))=e∑原创 2020-06-01 22:47:10 · 1297 阅读 · 0 评论 -
多项式牛顿迭代
以前以为可以不写,现在发现不得不写数域的变化重要的问题:什么是x?x?x?如果我们的数域是实数,那么xxx是一个自变量,我们可以求dxdf(x)\frac {{\rm d }x}{{\rm d}f(x)}df(x)dx,可以在x=Cx = Cx=C处泰勒展开,可以牛顿迭代求出xxx的值。但是我们现在的数域是形式幂级数,xxx的地位,和实数域中的111一样,是一个常数。求dxdf(x)?\frac {{\rm d }x}{{\rm d}f(x)}?df(x)dx?显然对于常数xxx,它等于xx原创 2020-06-01 14:31:36 · 750 阅读 · 0 评论 -
LOJ #6363. 「地底蔷薇」(拉格朗日反演)
题目首先可以看到然后如果你做过点双连通图计数就会发现这个条件让你可以直接多次拉格朗日反演求出大小∈S\in S∈S的所有点双连通图的方案数,设他的生成函数为C(x)C(x)C(x)。然后和点双联通图计数一样的方法设答案为F(x)F(x)F(x)。那么考虑包含根的点双,把这些点双的边删掉之后,所有的点都可以挂一个我们正在求的图也即F(x)F(x)F(x),然后自由组合。所以可以得到方程:F(x)=xexpC(F(x))F(x) = x \exp C(F(x))F(x)=xexpC(F(x))原创 2020-05-31 15:32:16 · 515 阅读 · 0 评论 -
CF1349F Slime and Sequences (欧拉数,扩展拉格朗日反演)
这篇博客大部分抄自此人好序列为对于一个长度为nnn的序列,其中若有ai=k(k>1)a_i= k(k>1)ai=k(k>1),则必存在aj∈[1,i−1]=k−1a_{j\in[1,i-1]} = k-1aj∈[1,i−1]=k−1求对于所有长度为nnn的好序列,分别求出k=1...nk = 1...nk=1...n的出现次数和。第一步神仙映射:把所有的好序列映射到排列中,发现可以一一对应。对应方法如下:排列到好序列:对于一个排列pip_ipi,如果p1,p2...pi原创 2020-05-27 11:55:34 · 667 阅读 · 3 评论 -
烷烃计数(Burnside引理,分治FFT)
在本学期的期末考试中, 连原电池正负极都分不清的你深感自己要爆零了没错正是在下. 就在这时, 你看到了一道附加题:33.\texttt{33.}33.(本题共计 100\texttt{100}100 分)23333\texttt{23333}23333 烷的同分异构体个数为 _______。(不考虑立体异构)是时候表演真正的计数了! 现在, 请你快速计算化学式为 CnH2n+2C_nH_{2n...原创 2020-03-31 18:57:00 · 2533 阅读 · 2 评论 -
牛客网OJ 54842 序列(分段多项式转移方程)
题目好野。写了一发分段多项式。感觉要想清楚才好写。转移是(定)积分的形式,阶乘最后再除也没问题(要记次数)。转移后的常数项,就是Ri≤xR_i\le xRi≤x的部分的贡献。AC Code\mathrm {AC \ Code}AC Code#include<bits/stdc++.h>#define maxn 405#define P 9982...原创 2020-01-17 22:35:41 · 337 阅读 · 0 评论 -
多项式的各种算法\cdot极
极\cdot多项式假之前的一篇感觉太不系统了。这次将用链套链的方式写好这些博文造好这些轮子。0.前置知识:0.1 极限0.2 生成函数1.FFT\mathrm {FFT}FFT2.NTT\mathrm {NTT}NTT3.MTT暂且跳过。4.多项式求逆5.多项式对数函数6.多项式指数函数7.多项式开根8.多项式快速幂9.多项式除法10....原创 2019-12-22 17:04:08 · 2852 阅读 · 0 评论 -
多项式卷积与DFT与FFT
之所以取这个标题。。。是因为我自己以前只认为FFT只是一种精妙的构造。并没能意识到这只是信息学中庞大多项式理论中的the tip of the iceberg{\rm the\ tip \ of \ the \ iceberg}the tip of the iceberg查漏补缺格物致知。多项式:多项式可以是...原创 2019-05-17 15:06:58 · 1641 阅读 · 0 评论 -
特征多项式及Cayley-Hamilton定理
大佬blog大佬blog大佬blog学数学的时候,可能会接触到一个叫做特征根法,特征根方程的东西,当时不觉明历。实际上这是和线性代数中的特征多项式离不开的。学OI的时候,可能会接触到矩阵快速幂求解常系数齐次线性递推的东西,懂了但是只会当模板用。实际上这也和特征多项式有着紧密的联系。言归正传。特征多项式是指对一个nnn阶(转移)方阵AAA , f(λ)=det(λE−A)=λN+b...原创 2019-05-21 22:07:03 · 4903 阅读 · 2 评论 -
LOJ #2462. 「2018 集训队互测 Day 1」完美的集合(组合数取(奇怪模数的)模 + DFS序DP+容斥)
题目part1\texttt{part1}part1DP求出完美集合的价值之和,这种树上联通块的DP一般可以用DFS序DP,把DFS序拿出来从右到左,如果当前dfs序为i的点选,则可以选择dfs序为i+1的点,否则跳过整颗子树,即f[i]=max(f[i+1],f[i+siz[pos[i]]])f[i] = max(f[i+1] , f[i + siz[pos[i]]])f[i]=max(f[...原创 2019-05-20 22:19:47 · 856 阅读 · 0 评论