————FFT
BAJim_H
比孤独更可悲的事情,就是根本不知道自己很孤独,或者分明很孤独,却把自己都骗得相信自己不孤独。
展开
-
多项式系列模板
FFTstruct Z{ double x,y; Z(double _x=0,double _y=0):x(_x),y(_y){};};Z operator +(Z x,Z y) {return Z(x.x+y.x,x.y+y.y);}Z operator -(Z x,Z y) {return Z(x.x-y.x,x.y-y.y);}Z operator *(Z ...原创 2018-04-23 15:07:58 · 718 阅读 · 0 评论 -
【杂题】【CometOJ Contest #5】E:迫真大游戏【概率】【排列组合】【多项式】
Description有一个n个点的环,有一个指针会从1号点开始向后扫描,每次扫描有p的概率删除当前点询问每个点最后一个被删除的概率。答案对998244353取模n<=200000Solution直接计算要考虑前后两部分比较麻烦,我们不妨先考虑1号点如何计算。记fnf_nfn表示n个点的环1号点最后一个被删除的概率。枚举一轮中删去了几个fn=∑i=0n−1pi(1−p)n...原创 2019-06-15 21:13:00 · 275 阅读 · 0 评论 -
【PKUSC2019】线弦图【计数】【树形DP】【分治FFT】
Description定义线图为把无向图的边变成点,新图中点与点之间右边当且仅当它们对应的边在原图中有公共点,这样得到的图。定义弦图为不存在一个长度大于3的纯环,纯环的定义是在环上任取两个不相邻的点,它们之间都没有边,也就是不存在没有弦的环的无向图。现在给出一棵n个点的树,你可以在上面添加任意多条边(不能重边),要求得到的图的线图是弦图,求加边的方案数。n<=200000Solut...原创 2019-05-29 20:08:12 · 482 阅读 · 0 评论 -
[JZOJ6088] [BZOJ5376] [loj #2463]【2018集训队互测Day 1】完美的旅行【线性递推】【多项式】【FWT】
DescriptionSolution我们考虑将问题一步步拆解第一步求出FS,iF_{S,i}FS,i表示一次旅行按位与的值为S,走了i步的方案数。第二步答案是FS,iF_{S,i}FS,i的二维重复卷积,记答案为SS,iS_{S,i}SS,i,那么FS,i×ST,jF_{S,i}\times S_{T,j}FS,i×ST,j能够贡献到SS&T,i+jS_{S\...原创 2019-03-28 22:00:49 · 277 阅读 · 0 评论 -
[JZOJ6065]【NOI2019模拟2019.3.18】One?One!【数学】【多项式】
DescriptionSolution这道题虽然相对套路,但还是蛮考人的我们要求的实际上就是∑d⌊n111...1⌋(d个1)\sum\limits_{d} {\lfloor{n\over 111...1}\rfloor}(d个1)d∑⌊111...1n⌋(d个1)有一个套路就是上下同时乘9令N=9nN=9nN=9n那就变成了∑d≥2⌊N10d−1⌋\sum\limits_{d...原创 2019-03-18 22:41:04 · 437 阅读 · 0 评论 -
[JZOJ6054]【NOI2019模拟2019.3.12】Z的礼物【计数】【斯特林反演】【分治FFT】
DescriptionSolution设pi=∑j=1ibjp_i=\sum\limits_{j=1}^{i}b_jpi=j=1∑ibj由题意可以列出式子ai=∑j=1iS(i,j)pja_i=\sum\limits_{j=1}^{i}S(i,j)p_jai=j=1∑iS(i,j)pj其中S(i,j)为第二类斯特林数根据斯特林反演的式子F(n)=∑i=0nS(i,j)...原创 2019-03-13 21:53:16 · 273 阅读 · 0 评论 -
【杂题】[LibreOJ #2409]【THUPC 2017】小 L 的计算题 / Sum 【数学】【多项式】
Description给出n个数a1...ana_1...a_na1...an求它们的111~nnn次方和n≤200000,a≤109n\leq 200000,a\leq 10^9n≤200000,a≤10920组数据,n总和不超过400000Solution此题的关键是一个叫牛顿恒等式的东西。考虑多项式F(x)=anxn+an−1xn−1+⋯+a1x+a0F(x)=a_nx^n...原创 2019-03-01 22:26:13 · 430 阅读 · 0 评论 -
[JZOJ5998]【WC2019模拟2019.1.14】操作
DescriptionSolutioin推一波式子,容易将一个操作写成一次函数的形式。现在只需要对于每个操作计算所有可能在它后面的操作对它的影响,也就是斜率乘积乘上这种情况出现的概率。考虑分治NTT但此时我们发现需要对于每一个操作求和,且自己不能选。有一个很妙的做法,分治NTT的时候对于两边分别维护全选上的和缺一个的和左边缺一个 * 右边全选 + 右边缺一个 * 左边全选就是整个...原创 2019-01-21 09:21:02 · 349 阅读 · 1 评论 -
【学习小记】常系数齐次线性递推
问题引入给出数列ggg,满足当n&amp;amp;amp;gt;mn&amp;amp;amp;gt;mn&amp;amp;gt;m时gn=∑i=1mgn−i×aig_n=\sum\limits_{i=1}^{m}g_{n-i}\times a_ign=i=1∑mgn−i×ai当n&amp;amp;amp;lt;=mn&amp;amp;amp;lt;=mn&amp;amp;lt;=m时,gn=c原创 2019-01-03 22:29:19 · 1488 阅读 · 0 评论 -
[JZOJ5529] 朱老师的难题 【多项式】【生成函数】
Description有一个n个数的严格递增序列a定义一组数为序列a的一个非空子集,也就是说这里面的元素不能重复,也没有顺序之分。一组数SSS的权值为(−1)∣S∣+1∏i∈Sai(-1)^{|S|+1}\prod\limits_{i\in S} a_i(−1)∣S∣+1i∈S∏ai你可以选择若干组数(组与组之间有顺序,一个元素可以同时出现在很多组)作为一个选择方案,这个方案的权值是所...原创 2018-12-21 22:25:15 · 538 阅读 · 0 评论 -
多项式多点求值
Preface首先我们需要的前置知识有:多项式取模余式定理多项式取模可以看这里余式定理的内容是这样的:对于多项式F(x),其在x=x0处的点值等于F(x) mod (x-x0)因为除式是个一次式,那么余式一定是一个常数Procedure考虑分治假设我们要求多项式F(x)F(x)F(x)在(x0…xn-1)这n个点的点值设Pl,r(x)=∏i=lr(x−xi),Fl,r...原创 2018-12-21 12:28:11 · 1745 阅读 · 0 评论 -
多项式求逆与多项式除法/取模
多项式求逆多项式求逆是多项式模块中的一个重要操作(“操作”这个词看出如今多项式题是多么的工业化,犹如毒瘤8操作LCT),在做生成函数/多项式除法、多项式取模/多项式多点求值中均有应用对于一个n次多项式F(x)F(x)F(x),我们希望求出一个m-1次多项式G(x)G(x)G(x),满足F(x)G(x)≡1(modxm)F(x)G(x)\equiv 1\pmod {x^m}F(x)G(x)≡1(...原创 2018-10-17 08:04:33 · 6164 阅读 · 0 评论 -
[JZOJ3303][BZOJ3456] 城市规划【多项式】【生成函数】
Description刚刚解决完电力网络的问题, 阿狸又被领导的任务给难住了.刚才说过, 阿狸的国家有n 个城市, 现在国家需要在某些城市对之间建立一些贸易路线, 使得整个国家的任意两个城市都直接或间接的连通.为了省钱, 每两个城市之间最多只能有一条直接的贸易路径. 对于两个建立路线的方案, 如果存在一个城市对, 在两个方案中是否建立路线不一样, 那么这两个方案就是不同的, 否则就是相同的...原创 2018-10-18 12:28:53 · 357 阅读 · 0 评论 -
【杂题】【树形DP】【NTT】[未知来源] 树 【无实现】
Description 给定一棵有n个节点的树,每个点有一个点权a[i],每条边有一个出现的概率p(u,v),那么随机出现若干个联通块,给定一个整数k,对于每个联通块的特征值记为(∑a[i])k(∑a[i])k(\sum a[i])^k,对于整棵树的特征值我们记为联通块的特征值的和,求整棵树的特征值的期望。 答案对998244353取模 1&lt;=n,k&lt;=1000...原创 2018-09-04 17:32:08 · 285 阅读 · 0 评论 -
[JZOJ4203] Tree
Description 注意:本题空间限制为7168KBSolution我们按照每层操作涉及的区间长度来考虑,以包含1的区间为例,[1,1],[1,2],[1,4],[1,8]….. 可以发现这很像FFT中的蝴蝶变换 前半段转移是or,后半段转移是and于是只需要开一个2n2n2^n的数组,类似FFT的过程写出来,最后变换到顶层每个位置就是所对应的i 直接累加答...原创 2018-06-19 16:16:13 · 204 阅读 · 0 评论 -
[JZOJ5666]【GDOI2018Day2模拟4.18】法力风暴(分治NTT 模板)
Description 有2≤n≤105,0≤Ai,k≤1092≤n≤105,0≤Ai,k≤1092 ≤ n ≤ 10^5 ,0 ≤ A_i ,k ≤ 10^9 Solution注意到一次操作打出的伤害就是原来A的乘积减去操作后A的乘积那么题目转化为求原来A的乘积减去最终A的乘积的期望设a[i]a[i]a[i]最终被减去了b[i]b[i]b[i]那么最终期望为 E(...原创 2018-04-23 16:18:20 · 422 阅读 · 0 评论 -
【杂题】[LibreOJ #6608] 无意识的石子堆【容斥原理】【FFT】
DescriptionSolution943718401=225*2^22+1显然每行必须有两个,我们不妨枚举有k列有2个石子,那么有2(n-k)列有1个石子。Ans=∑k=0n(mk)(m−k2(n−k))SkAns=\sum\limits_{k=0}^{n}{m\choose k}{m-k\choose 2(n-k)}S_kAns=k=0∑n(km)(2(n−k)m−k)Sk...原创 2019-06-10 21:16:26 · 395 阅读 · 1 评论