https://www.lydsy.com/JudgeOnline/problem.php?id=3625 愉快地列式子.设\(F[i]\)表示权值为\(i\) 的子树的方案数,\(A[i]\)为\(i\)在不在集合中. \[ F[n]=\sum_{i=0}^n \sum_{j=0}^{n-i}F[i]\cdot F[j]\cdot A[n-i-j] \] 初始状态\(F[0]=1\). 我们把\(F,A\)看成多项式. \[ F(x)-1=F^2(x)\cdot A(x)\\ A(x)\cdo…

Code: #include <bits/stdc++.h> #define N 1000010 #define mod 998244353 #define setIO(s) freopen(s".in","r",stdin) using namespace std; namespace NTT { int a[N],b[N],f[N],g[N]; int qpow(int x,int y) { int tmp = 1; while(y) { if(y&…

(首先要%miskcoo,这位dalao写的博客(这里)实在是太强啦qwq大部分多项式相关的知识都是从这位dalao博客里面学的,下面这篇东西是自己对其博客学习后的一些总结和想法,大部分是按照其博客里面的思路来分析的,并添加了一些自己的理解) 多项式求逆(元) 定义 对于一个多项式\(A(x)\),如果存在一个多项式\(B(x)\),满足\(B(x)\)的次数小于等于\(A(x)\)且\(A(x)B(x)\equiv 1(mod\ x^n)\),那么我们称\(B(x)\)为\(A(x)\)在模\…

出处0.0用到第二类斯特林数的性质,做法好像很多,我打的是直接ntt,由第二类斯特林数的容斥公式可以推出,我们可以对于每一个i,来一次ntt求出他与所有j组成的第二类斯特林数的值,这个时候我们是O(n^2logn)的,还不如暴力,但是我们发现,对于刚刚提到的容斥的式子,将其化为卷积形式后,其一边的每一项对于每一个i都相同,另一边的每一项是对于所有的i形成一个n项的等比数列,这样我们可以把成等比数列的一边求和,用固定的一边去卷他们的和,这时候的答案的每一项就是所有的i的这一项的和,然后我们再O(n…

参考:https://www.cnblogs.com/2016gdgzoi509/p/8999460.html 列出生成函数方程,g(x)是价值x的个数 \[ f(x)=g(x)*f^2(x)+1 \] +1是f[0]=1 根据公式解出 \[ f(x)=\frac{1+(-)\sqrt{1-4*g(x)}}{2*g(x)} \] 舍去+的答案,分式上下同乘\( 1-\sqrt{1-4*g(x)} \) \[ f(x)=\frac{2}{1+\sqrt{1-4*g(x)}} \] 然后套多项式开跟…

http://blog.miskcoo.com/2015/05/polynomial-division 好神啊! 通过翻转多项式消除余数的影响,主要原理是商只与次数不小于m的项有关. #include<iostream> #include<cstdio> #include<cmath> #include<cstdlib> #include<cstring> #include<algorithm> using namespace std…

题意:求多项式的逆 题解:多项式最高次项叫度deg,假设我们对于多项式\(A(x)*B(x)\equiv 1\),已知A,求B 假设度为n-1,\(A(x)*B(x)\equiv 1(mod x^{\lceil \frac{n}{2} \rceil})\),\(A(x)*B'(x)\equiv 1(mod x^{\lceil \frac{n}{2} \rceil})\) 两式相减得\(B(x)-B'(x)\equiv 0(mod x^{\lceil \frac{n}{2} \rceil})\),…

题目链接 设多项式\(f(x)\)在模\(x^n\)下的逆元为\(g(x)\) \[f(x)g(x)\equiv 1\ (mod\ x^n)\] \[f(x)g(x)-1\equiv 0\ (mod\ x^n)\] \[f^2(x)g^2(x)-2f(x)g(x)+1\equiv 0\ (mod\ x^{2n})\] \[2f(x)g(x)-f^2(x)g^2(x)\equiv 1\ (mod\ x^{2n})\] \[2f(x)g(x)-f^2(x)g^2(x)\equiv f(x)g'(x)…

http://172.20.6.3/Problem_Show.asp?id=2043 最开始用了FFT,交上去全tle和wa了(tle的比较多),测了一组数据发现求逆元的过程爆double了(毕竟系数的指数幂也是指数增长的,科学计数法也撑不住). 然后问了出题人,发现出题人忘了给用来%的P(想用钢丝球刷出题人QnQ),所以其实是ntt,改完ntt加了个快读a了. 题解:http://blog.miskcoo.com/2015/05/polynomial-division 依然是推式子,系数反转求…

题目大意:给定长度为$n-1$的数组$g_{[1,n)}$,求$f_{[0,n)}$,要求: $$f_i=\sum_{j=1}^if_{i-j}g_j\\f_0=1$$ 题解:分治$FFT$博客,发现这道题就是求$f*g=f-1$($f-1$就是没有常数项的$f$),改写一下式子:$$f*g\equiv f-1\pmod{x^n}\\f-f*g\equiv1\pmod{x^n}\\f*(1-g)\equiv1\pmod{x^n}\\f\equiv(1-g)^{-1}\pmod{x^n}$$ 卡点…