---多项式
BAJim_H
比孤独更可悲的事情,就是根本不知道自己很孤独,或者分明很孤独,却把自己都骗得相信自己不孤独。
展开
-
拉格朗日插值法
Text对于一个k次多项式函数f(x)=∑i=0kaixif(x)=\sum\limits_{i=0}^{k}a_ix^i有两种表示其的方法。 可以用传统的每一项的系数来表示。同时由于只要有K+1个(x,fx)的点对,就能唯一确定这个多项式,因此把这种表示法叫点值表示法。系数表示法转换成点值表示法叫求值,反过来就叫插值。当然可以高斯消元,但是复杂度比较大。主要讨论拉格朗日插值法。假设已经求得了K+原创 2017-10-06 22:03:50 · 966 阅读 · 0 评论 -
[JZOJ6044]【NOI2019模拟2019.3.7】联邦【状压DP】【子集卷积】【FWT】
DescriptionSolution城市规划想必都做过,求无向连通图的个数这题的思路是类似的设G[S],F[S]G[S],F[S]G[S],F[S]分别表示选取了S中的点所有图的c的乘积的和,以及连通图的c的乘积的和。G[S]G[S]G[S]可以简单O(2n∗n2)O(2^n*n^2)O(2n∗n2)预处理对于F,考虑总的减去不连通的,枚举编号最小/最大的点所在的连通块,令其为i...原创 2019-03-08 21:35:45 · 359 阅读 · 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 · 279 阅读 · 0 评论 -
[JZOJ6080]【GDOI2019模拟2019.3.23】IOer【生成函数】【数学】
Description有m天时间,每天可以完成任意数量的事件,第i天可以完成的事件有(ui+v)种(可以同一种完成多次)问m天总共完成n个事件的方案数。注意每天完成的时间是有顺序的,也就是说同一天内1,2和2,1是两种不同的方案。T组数据。n≤1018,m≤200000,u,v≤109,T≤5n\leq 10^{18},m\leq 200000,u,v\leq 10^9,T\leq 5n...原创 2019-03-24 22:33:37 · 717 阅读 · 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 · 440 阅读 · 0 评论 -
[JZOJ6082]【GDOI2019模拟2019.3.25】染色问题【多项式】【DP】
Description考虑一个长度为n的序列,你需要用颜色1~m依次对序列染色每次染色只能选择一段区间来染,要求最后整个序列没有未染色的元素。问最终的颜色序列的种类数,答案对998244353取模。Solution我们考虑这样一个DP设F[i][j]F[i][j]F[i][j]表示当前染了前i种颜色,已经覆盖了jjj个元素。有F[i][j]=F[i−1][j]+∑k=0j−1F[i−...原创 2019-03-26 08:11:00 · 523 阅读 · 0 评论 -
[JZOJ6090]【GDOI2019模拟2019.3.27】圆【数学】
Description有nnn场比赛,第i场比赛的前aia_iai名晋级有q次询问,每次询问一个区间[l,r][l,r][l,r],和一个xxx,表示只参加[l,r][l,r][l,r]的比赛,每场的排名会独立的在[1..x][1..x][1..x]中随机。回答至少有一场晋级的概率。没有模数答案与标准答案相差10−610^{-6}10−6以内即为正确。Soluton很有意思的一题。...原创 2019-03-28 15:10:39 · 185 阅读 · 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 评论 -
【WC2019】数树【计数】【DP】【多项式】
Description此题含有三个子问题问题1:给出n个点的两棵树,记m为同时在两棵树中的边的个数,求ymy^mym问题2:给出n个点的一棵树,另外一棵树任意生成,求所有方案总的ymy^mym的和问题3:两棵树均任意生成,求所有方案总的ymy^mym的和Solution留坑待填Code#include <bits/stdc++.h>#define fo(i,a,...原创 2019-04-20 21:48:44 · 260 阅读 · 0 评论 -
[LibreOJ 3120]【CTS2019】珍珠 【生成函数】【计数】
DescriptionSolution有一个直观的思路是考虑每种颜色个数的奇偶性,奇数个数的颜色不能超过n−2mn-2mn−2m因此若n−2m≥Dn-2m\geq Dn−2m≥D则答案一定是DnD^nDn否则由于每种颜色其实没有区别,我们考虑一种颜色为奇数和为偶数的指数型生成函数奇数是ex−e−x2e^x-e^{-x}\over 22ex−e−x,偶数是ex+e−x2e^x+e^{...原创 2019-05-21 19:59:09 · 2440 阅读 · 2 评论 -
【PKUSC2019】线弦图【计数】【树形DP】【分治FFT】
Description定义线图为把无向图的边变成点,新图中点与点之间右边当且仅当它们对应的边在原图中有公共点,这样得到的图。定义弦图为不存在一个长度大于3的纯环,纯环的定义是在环上任取两个不相邻的点,它们之间都没有边,也就是不存在没有弦的环的无向图。现在给出一棵n个点的树,你可以在上面添加任意多条边(不能重边),要求得到的图的线图是弦图,求加边的方案数。n<=200000Solut...原创 2019-05-29 20:08:12 · 485 阅读 · 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 · 280 阅读 · 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 · 432 阅读 · 0 评论 -
[JZOJ5998]【WC2019模拟2019.1.14】操作
DescriptionSolutioin推一波式子,容易将一个操作写成一次函数的形式。现在只需要对于每个操作计算所有可能在它后面的操作对它的影响,也就是斜率乘积乘上这种情况出现的概率。考虑分治NTT但此时我们发现需要对于每一个操作求和,且自己不能选。有一个很妙的做法,分治NTT的时候对于两边分别维护全选上的和缺一个的和左边缺一个 * 右边全选 + 右边缺一个 * 左边全选就是整个...原创 2019-01-21 09:21:02 · 352 阅读 · 1 评论 -
[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 · 427 阅读 · 0 评论 -
多项式系列模板
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 · 719 阅读 · 0 评论 -
[JZOJ3303][BZOJ3456] 城市规划【多项式】【生成函数】
Description刚刚解决完电力网络的问题, 阿狸又被领导的任务给难住了.刚才说过, 阿狸的国家有n 个城市, 现在国家需要在某些城市对之间建立一些贸易路线, 使得整个国家的任意两个城市都直接或间接的连通.为了省钱, 每两个城市之间最多只能有一条直接的贸易路径. 对于两个建立路线的方案, 如果存在一个城市对, 在两个方案中是否建立路线不一样, 那么这两个方案就是不同的, 否则就是相同的...原创 2018-10-18 12:28:53 · 358 阅读 · 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 · 6193 阅读 · 0 评论 -
快速沃尔什变换(FWT)
Preface我们知道,对于下面的式子f(n)=∑i⊕j=ng(i)h(j)f(n)=\sum\limits_{i⊕j=n}g(i)h(j)f(n)=i⊕j=n∑g(i)h(j)当⊕⊕⊕为+运算时,可以用FFT优化当⊕⊕⊕为位运算(or,and,xor)时是否也有一种优化呢答案是肯定的。FWT首先,对于多项式,我们定义A∗BA*BA∗B为某一种位运算卷积显然它满足交换律、结合律...原创 2018-10-24 11:01:37 · 1159 阅读 · 0 评论 -
【杂题】[LibreOJ 2541] 【PKUWC2018】猎人杀【生成函数】【概率与期望】
Description猎人杀是一款风靡一时的游戏“狼人杀”的民间版本,他的规则是这样的:一开始有 n个猎人,第 i 个猎人有仇恨度 wi。每个猎人只有一个固定的技能:死亡后必须开一枪,且被射中的人也会死亡。然而向谁开枪也是有讲究的,假设当前还活着的猎人有[i1...im][i_1...i_m][i1...im],那么有wik∑j=1mwijw_{i_k}\over \sum\limits...原创 2018-11-08 15:33:47 · 362 阅读 · 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 · 1751 阅读 · 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 · 542 阅读 · 0 评论 -
【学习小记】常系数齐次线性递推
问题引入给出数列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 · 1489 阅读 · 0 评论 -
【2018.12.28模拟赛】再一次张开翅膀 【多项式】【生成函数】【线性递推】(无实现)
Descritption给出n,mn,mn,m,求在[1..n][1..n][1..n]中选出m个互不相同的数,它们能作为一个凸多边形的边长的方案数。n≤109,3≤m≤200n\leq 10^9,3\leq m\leq 200n≤109,3≤m≤200Solution我们发现,若干个数能构成凸多边形的条件就是最大值小于其他数之和。考虑补集转化,计算和不超过最大值的方案数,再用总数(nm...原创 2019-01-04 22:22:40 · 276 阅读 · 0 评论 -
[JZOJ5520] Every one will meet some difficult [JZOJ5985] 多项式【TCO2013 3A】【组合数学】【多项式】【生成函数】
Description求∑i=1mai≤S\sum\limits_{i=1}^{m}a_i\leq Si=1∑mai≤S且∀i,ai&gt;0\forall i,a_i&gt;0∀i,ai>0∀i≤n,ai≤T\forall i\leq n,a_i\leq T∀i≤n,ai≤T的a的解的组数n≤m≤1e9,T≤105,n∗T≤S≤1018,m−n≤1000...原创 2019-01-05 16:19:46 · 508 阅读 · 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 · 398 阅读 · 1 评论