计数(5):多项式相关

news/2025/9/30 16:20:45/文章来源:https://www.cnblogs.com/thirty-two/p/19085590

计数中我目前只见过两类与多项式相关的问题,一类是采用 gf 来刻画,一类是用拉插来优化掉 dp 的很大一个维度,还有一类就是神秘多项式问题,但是考虑 fft 已经被 cnoi 除名了所以我就不学了(

这部分写的非常浅显,我本身的多项式水平是负数,如果大家有题目推荐欢迎狠狠地踹主包的说。


拉格朗日插值法

想到一个很有意思的事情,我第一次见到拉格朗日插值法被名字吓跑了,后来有一天看到老高中数学课本,看到 CRT 居然有用拉格朗日插值法的理解,然后就去写板子了。

基本形式

给定 \(n\) 个点 \((x_1, y_1)\sim (x_n, y_n)\),求一个 $(n - 1) $ 次经过所有点的多项式中 \(x = k\) 的函数值。

无脑的做法是高斯消元求解系数。

但是注意到我们并不在乎系数只在乎求值。因此有:

\[f(k) = \sum\limits_{i = 1}^n y_i\prod_{j \not= i}\dfrac{k - x_j}{x_i - x_j} \]

感性认知一下:我们给每个 \(y\) 配上一个系数,对于 \(x_i\) 那么 \(y_i\) 系数为 \(1\),否则为 \(0\)。因此将 \(x_i\) 带入进来的时候我们希望首先 \(j\not= i\),其次总是有分子分母相等,同样的如果 \(x_i\) 与这一项不相等那么就等于 \(0\)。得到系数 \(\prod\limits_{j \not = i}\dfrac{k - x_j}{x_i - x_j}\)


拓展

值域连续时

\(x_i\) 的值域连续时,以 \(x_i\)\(1\sim n\) 为例,此时:

\[f(k) = \sum\limits_{i = 1}^n y_i\prod\limits_{j \not= i}\dfrac{(k - j)}{i - j} \\ = \sum\limits_{i = 1}^n y_i \dfrac{\frac{(k - 1)^{\underline{n}}}{k - i}}{(i-1)!(-1)^{\underline{n - i}}} \]

必须注意的是,\(k\le n\) 时,\((k - 1)^{\underline{n}}\) 可能等于 \(0\) 要特判!

预处理一些阶乘即可。

重心拉格朗日插值

有的时候我们要对单点 \(k\) 求取值,但是点会不断增加。

变换式子:

\[w = \prod\limits_{j = 1}^n (k - x_j)\\f(k) = w\sum\limits_{i = 1}^n \dfrac{1}{k - x_i}\prod\limits_{j \not =i}^n\dfrac{y_i}{x_j - x_i} \]

那么令 \(t_i = \prod\limits_{j \not= i}\dfrac{y_i}{x_j - x_i}\),新加入一个点只需要对所有 \(t_i\) 进行更新,计算新的 \(t_n\)\(w\),最后计算贡献。加入一个点的代价是 \(O(n)\) 的。


例题

AT_abc208_f

关键是要注意到组合数也可以看作是多项式,后面忘了

题解 https://www.luogu.com.cn/article/95bxguv3

P8290

我被卡常了。

首先考虑暴力,最大值与最小值的差很难做,考虑钦定范围在 \([L, R]\) 内其中 \(R - L = k\)。这样会算重,常见的套路是钦定端点被选择,比如这里我们再减去 \([L+1, R]\) 的情况即可。那么每个点能选择的就是 \([l_i, r_i]\cap [L, R]\) 这个区间。第一问直接 dp,第二问要做一个换根,对于一个区间可以 \(O(n)\) 求。最终时间复杂度是 \(O(nV)\) 的。现在考虑优化。首先如果我们将 \(l_i - 1, l_i, r_i, r_i+1\) 这些可能的端点离散化,那么两个端点之间的情况显然是相同的,要么然完全包含,要么然是与 \(l_i, r_i\) 形成一次函数,总共乘起来是 \((n+1)\) 次的函数,最后再做前缀和就是 \((n+2)\) 次函数,要枚举前 \((n+3)\) 项算前缀和。拉差即可,时间复杂度 \(O(n^3)\)

P3270

当时应该是省选之前和 yzy 一起看的这题,我俩当时怎么都不会啊,暴力都不会。

为了方便下问认为 \(n \gets n - 1, R_i = n - R_i + 1\)

注意到直接约束“恰好 \(K\) 个人全部科目小于等于 B 的”很难做,所以考虑容斥,钦定有 \(x\) 个人被吊打了就很容易了。这是一个明显的二项式反演的形式,\(g(x)\) 为钦定 \(x\) 个人被碾压了,那么就有:

\[g(x) = \binom{n}{x}\prod\limits_{i = 1}^m(\sum\limits_{t = 1}^{U_i}\binom{n - x}{R_i - x}t^{R_i}(U_i - t)^{n - R_i})\\ f(K) = \sum\limits_{i = K}^n (-1)^{i - K}\binom{i}{K}g(i) \]

瓶颈在于快速计算 \(\sum\limits_{t=1}^{U_i}t^{R_i}(U_i - t)^{n - R_i}\),拉插即可。


生成函数相关

形式幂级数也是幂级哈!

利用形式幂级数分析问题,而不是使用多项式科技

P5333

这也太难了

不会,我学。

首先考虑刻画哈密尔顿回路,通过手玩不难发现,本质上是将这些树分别剖成若干条不交链,哈密尔顿回路可以看作这些链的拼接,要求:

  1. 相邻的两条链来自不同树。
  2. 第一条链必须来自第一棵树,并且必须包含 \(1\) 点。
  3. 最后一条链必须不来自第一棵树。

假设已经钦定了链的剖分方式,如果一条链大小大于 \(1\) 那么产生 \(2\) 方案,否则只产生 \(1\) 方案。设 \(f[i, j]\) 代表第 \(i\) 棵树剖分成 \(j\) 条链的权和,若钦定了 \(j_1, j_2, \dots, j_m\),那么总权和就是 \(\prod_{i = 1}^n f[i, j_i]\)。那么只需要计算链的排列方式。

首先考虑第一条约束。

考虑容斥。独立考虑一棵树的影响,对于树 \(i\),剖分成 \(j\) 条链,钦定其中有 \(k\) 个位置最后相邻,容斥系数则为 \((-1)^k\),排列一共应当由 \(j!\binom{j - 1}{k}\) 中选择方法。于是我们得到了新的 \(x^{j - k}\) 条链。最后我们是按照这个 \((j - k)\) 之和相关计算答案的,不难想到构造出 egf:

\[G_{i} = \sum\limits_{j = 0}f[i, j]j!\sum\limits_{k = 0}^j\binom{j - 1}{k}(-1)^k\dfrac{x^{j - k}}{(j - k)!} \]

如果没有别的约束,将 \(G_i\) 全部相乘计算每一项系数即可。

对于第二条约束,强制第一条链,注意,链内部只要将 \(j!\) 改成 \((j - 1!)!\) 代表第一条链固定,即可,合并链的形式幂级数则要去掉这条链。即 \(G_1\) 变成:

\[G_1 = \sum\limits_{j = 1}f[1, j](j - 1)!\sum\limits_{k = 0}^j\binom{j - 1}{k}(-1)^k\dfrac{x^{j - k - 1}}{(j - k - 1)!} \]

即可。

对于第三条约束,我们正难则反,减去不合法,即对于第一棵树强制第一条链和最后一条链。注意,最后一条链是任选的,所以 \((j - 1)!\) 仍然不变,代表钦定了第一条,剩下任意选,最后还剩下一个不准动。唯一的修改也是在最后的形式幂级数中:

\[G_1 = \sum\limits_{j = 1}f[1, j](j - 1)!\sum\limits_{k = 0}^j\binom{j - 1}{k}(-1)^k\dfrac{x^{j - k - 2}}{(j - k - 2)!} \]

计算答案全部卷起来就行了。时间复杂度 \(O(\sum\limits k^2)\)


接下来三题都是形如:\(K\) 次操作取反之类的问题,前两个问题我们都采用生成函数刻画,但是最后一个问题并不采用这个方法而是直接考虑集合幂级数卷积。一方面是最后一道题用生成函数推导有点太困难了我还不会 QAQ,一方面也是提个醒不要被生成函数卷积给局限住啦。

zr3272

\(n\) 个点的树,每条边拥有权 \(p_i\) 与存在状况 \(c_i\)。无形的大手子进行了 \(K\) 次操作,每次有 \(\dfrac{p_i}{\sum_{j = 1}^{n-1}p_j}\) 取反 \(i\) 这条边的存在状况。求最终树上联通块大小乘积的期望。

猫生首次见到这种问题,猫很震撼!令 \(S = \sum\limits p\),接下来我们默认 \(S^{K}\) 最后除以。

首先考虑已知每条边存在情况计算概率。因为每个边有可能被取反多次,最后操作数量和为 \(K\),并不好直接刻画,所以我们采用 egf,这样卷积的过程天然为我们做好了这个背包。

对于第 \(i\) 条边,它取反情况记作 \(G_{i, 1}\),未被取反的情况记作 \(G_{i, 1}\),那么有:

\[G_1 = \dfrac{e^{px} - e^{-px}}{2}\\G_0 = \dfrac{e^{px} +e^{-px}}{2} \]

盘面的 egf 为所有 \(G\) 的乘积 \(F\),那么最终概率为 \(K![x^K]F\)。令 \(F = \sum\limits_{i\in [-S, S]} c_ie^{ix}\),那么有:

\[K![x^K]F = \sum\limits_{i\in [-S, S]} c_ii^K \]

维护 \(c_i\) 即可。

回到题目当中来,树形 dp 合并子树的时候考虑 \(u-v\) 选择哪一个 gf,重点在于我们无法轻易的维护“联通块乘积”这一个信息。使用组合意义,联通块乘积看作联通块内部选一个点方案数。这个可以增量计算,设 \(f[u, i, 0/1]\) 代表 \(u\) 所在联通块选/没选,当前维护的系数为 \(c_i\) 的系数只和。直接转移即可。

AT_abc231_g

离开事件发生独立性谈论期望的可乘性是无稽之谈。用生成函数刻画。

建立一个 gf,\([x^k]G_i(x)\) 代表 \(i\) 位置选中了 \(k\) 次的期望,也就是 \((a_i+k)\dfrac{1}{n^k}\)。外面套上 egf 计算可重集合,最后相乘卷积起来即可。

具体而言对于每个 \(i\) 的 ogf 为 \(G_i(x) = \sum\limits_{m = 0}{(a_i+{m})\dfrac{(\frac{x}{n})^m}{m!}}\)。最后计算 \(n![x^k]\prod_{i = 1}^nG_i(x)\) 即可。\(k\) 很大所以没有什么直接卷积的余地。考虑求出封闭形式后直接计算。

\[G_i(x)=a_i\sum\limits_{m = 0}\dfrac{(\frac{x}{n})^m}{m!} + \sum\limits_{m=1}\dfrac{(\frac{x}{n})^m}{(m-1)!}\\=a_ie^{\frac{x}{n}}+\dfrac{x}{n}e^{\frac{x}{n}} \]

注意到 \(\frac{x}{n}e^{\frac{x}{n}}\) 看上去很难维护,但是对于每个 \(i\) 都有这一项,所以可以直接枚举乘积中选了多少项这个。考虑计算答案的贡献形式,答案形式为 \(t(\frac{x}{n}e^{\frac{x}{n}})^me^{\frac{(n-m)}{n}x}\),展开计算一下:

\[t(\frac{x}{n}e^{\frac{x}{n}})^me^{(n-m)x} = t\dfrac{x^m}{n^m}e^{x}\\ = \dfrac{t}{n^m}\sum\limits_{i=0}\dfrac{x^{m+i}}{i!} \]

故第 \(k\)\(\dfrac{tk!}{n^m(k-m)!}\)。对于相同 \(m\) 记录 \(t\) 和即可。直接 dp 即可。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/923059.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

最新WTAPI开发微信机器人教程说明

最新WTAPI开发微信机器人教程说明 微信二次开发机器人接口开发:框架提供了开放接口,可以通过这些接口对微信进行操作。如接收用户消息、发送消息、操作朋友圈等。 自然语言处理:机器人需要能够理解自然语言,以便能…

线性DP - 学习笔记

动态规划:用空间代替重复计算。 有些递归在展开计算时,总是重复调用一个子问题的解,这种重复调用的递归变成动态规划很有收益。 如果每次展开都是不同的解,或者重复调用的现象很少,那么没有改动态规划的必要。 任…

2025 年最新制氮机厂家权威推荐排行榜:聚焦行业优质厂商综合实力,助力企业精准选购优质设备制氮机产生氮气/氮气纯化/设备改造/维修/保养/半导体用制氮机厂家推荐

在工业生产持续升级的当下,制氮机作为保障多行业生产安全、提升生产效率的关键设备,市场需求日益旺盛。但当前制氮机市场品牌繁杂,产品质量参差不齐,部分厂商技术落后导致设备能耗高、纯度波动大,还有些厂商售后体…

南京的网站建设高端网站的建设

一、效果图 1.点击显示菜单栏&#xff0c;点击x号关闭&#xff1b; 2.点击一级菜单&#xff0c;展开显示二级&#xff0c;并且加号变为减号&#xff1b; 3.点击其他一级导航&#xff0c;自动收起展开的导航。 二、代码实现 <!DOCTYPE html> <html><head>&…

现在asp做网站ai效果图网站

0x01应用场景 获取JavaCard卡内剩余空间&#xff0c;一方面是在评估一张卡的时候需要用到&#xff0c;另一方面是在应用个人化或者运行时需要用到。 例如&#xff1a;应用提供商为了保证自己的应用在卡内运行期间能够不受空间影响&#xff0c;一般会在个人化(安装应用)的时候先…

2025 年氨分解设备厂家最新推荐榜单:含半导体 / 冶金 / 稀土行业专用设备及品牌综合实力排名

随着双碳目标深入推进,氨分解制氢技术在氢能产业中的应用愈发广泛,2025 年全球市场规模预计达 28.7 亿美元。然而,行业发展面临诸多痛点:催化剂性能与成本难以平衡,铁基催化剂转化率低、钌基催化剂成本高;设备能…

idea打包推送maven仓库及同时推送到不同的maven仓库,本地和云上的腾讯云

idea打包推送maven仓库及同时推送到不同的maven仓库,本地和云上的腾讯云idea打包推送maven仓库及同时推送到不同的maven仓库,本地和云上的腾讯云 idea客户端打包推送jar报错如下:Execution failed for task :scoote…

2025 阳台柜厂家最新推荐榜:企业资质 / 材质 / 服务全景解析,选对品牌少走弯路杭州阳台柜/浙江阳台柜厂家推荐

阳台已从传统功能区升级为 “家庭第二客厅”,但消费者选购阳台柜时仍面临多重困扰:部分产品在高湿高晒环境下易发霉变形,寿命不足 3 年;多数品牌缺乏全链条服务,设计与落地严重脱节;中间商溢价使性价比大打折扣,…

广州个人网站搭建天津招标信息网

前面两篇文章已经整理了CompletableFuture大部分的特性&#xff0c;本文会整理完CompletableFuture余下的特性&#xff0c;以及将它跟RxJava进行比较。 3.6 Either Either 表示的是两个CompletableFuture&#xff0c;当其中任意一个CompletableFuture计算完成的时候就会执行。 …

2025 年除湿机厂家最新权威推荐排行榜:实力厂家技术口碑评测及场景适配选购指南吊顶/泳池/车库/防爆/调温/新风除湿机厂家推荐

当前除湿机市场需求持续攀升,工业生产对湿度精度要求愈发严苛,民用场景也对环境舒适度有了更高追求,但行业乱象却让用户选购难度大增。一方面,大量产品缺乏核心技术,除湿效率低、能耗高,难以适配电子、制药、家用…

2025 年液氨蒸发器厂家联系方式,众众电热:多领域加热设备供应与定制化解决方案提供商

在政策支持与产业升级的双重驱动下,液氨蒸发器等工业加热设备市场持续扩容,2024 年上半年市场规模已突破 52.7 亿元,预计 2025 年末将增至 68.5 亿元。化工合成、制冷储能等领域的需求增长,以及智能化、低能耗的技…

debian11不进入桌面环境打开chromiu

debian11不进入桌面环境打开chromiu在Debian非桌面环境中,使用Wayland(Weston)打开Chromium并打开特定网页,可以按照以下步骤进行:安装Weston和Chromium启动Weston合成器在Weston环境中启动Chromium 1,安装chrom…

【Batch】批量修改文件后缀

批量修改文件后缀 @echo off for %%F in (*.jpg) do (ren "%%F" "%%~nF.png" )-------------------------------------------------------------------------------------如果万事开头难 那请结局…

【solace】基于docker部署solace环境

一、solace https://www.solace.dev/ solace是一款高性能的消息中间件二、部署脚本docker run -d \-p 8080:8080 \-p 55555:55555 \-p:8008:8008 \-p:1883:1883 \-p:8000:8000 \-p:5672:5672 \-p:9000:9000 \-p:2222:2…

2025 年阿里巴巴开通实力商家店铺开通代运营 / 阿里巴巴新手开店全托管代运营公司推荐:南鑫粤网络 —— 全域运营解决方案与实战服务优势解析

2025 年的阿里巴巴 1688 平台已成为超 2.5 亿活跃采购商与千万级供应商汇聚的核心 B2B 阵地,实力商家与超级工厂模式成为店铺基础建设的关键方向。但随着平台商家突破 2000 万,竞争从 “入驻” 转向 “运营”,企业普…

Vue-element-admin开发指南 - 教程

Vue-element-admin开发指南 - 教程pre { white-space: pre !important; word-wrap: normal !important; overflow-x: auto !important; display: block !important; font-family: "Consolas", "Monaco&…

2025 年国内工作服厂家最新推荐排行榜:聚焦工艺设计与服务,精选权威榜单助企业采购冬季/春季/工人/车间/防静电/餐饮/劳保工作服厂家推荐

当前工作服市场产品质量参差不齐,部分产品面料不达标、耐用性差,且多数供应商缺乏个性化定制能力,难以满足不同行业企业的特殊需求,如医疗行业的抗菌需求、户外作业的防水保暖需求等。同时,企业采购时还常面临交付…

邯郸做网站的公司wordpress表情按钮

高维高斯分布基础 多位高斯分布的几何理解 多维高斯分布表达式为&#xff1a; p(x∣μ,Σ)1(2π)p/2∣Σ∣1/2e−12(x−μ)TΣ−1(x−μ)p(x|\mu,\Sigma)\frac{1}{(2\pi)^{p/2}|\Sigma|^{1/2}}e^{-\frac{1}{2}(x-\mu)^{T}\Sigma^{-1}(x-\mu)} p(x∣μ,Σ)(2π)p/2∣Σ∣1/21​…

ClickHouse 窗口函数使用详解(一) - 若

ROW_NUMBER() OVER (PARTITION BY department ORDER BY salary DESC) 详解 一句话概括 "在每个部门内部,按工资从高到低给员工编号" 逐词解析 1. ROW_NUMBER()作用:生成顺序编号特点:从1开始,连续不重复…

学做网站容易吗源代码下载

一、文件读取 Client向NameNode发起RPC请求&#xff0c;来确定请求文件block所在的位置&#xff1b;NameNode会视情况返回文件的部分或者全部block列表&#xff0c;对于每个block&#xff0c;NameNode 都会返回含有该 block 副本的 DataNode 地址&#xff1b; 这些返回的 DN 地…