首页 > 其他分享 >一些多项式常用操作的公式推导

一些多项式常用操作的公式推导

时间:2024-01-18 21:00:26浏览次数:30  
标签:推导 pmod dfrac times 公式 多项式 aligned equiv

怕我以后忘记了看不懂自己的板子(((

牛顿迭代

用于求解函数零点的近似值。

设函数 \(f(x)\) 的零点近似值为 \(x_0\),过点 \((x_0,f(x_0))\) 作 \(f(x)\) 的切线,切线与 \(x\) 轴交点的横坐标即为新的近似值。

切线解析式为 \(y=f'(x_0)(x-x_0)+f(x_0)\),当 \(y=0\) 时 \(x=x_0-\dfrac{f(x_0)}{f'(x_0)}\)。

同理,若要求满足 \(G(F(x))=0\) 的多项式 \(F(x)\),则 \(F(x)=F_0(x)-\dfrac{G(F_0(x))}{G'(F_0(x))}\),每迭代一次精度翻倍。


P4238 多项式求逆

\[\begin{aligned} A(x)\times F(x)&\equiv 1\pmod{x^n}\\ \\ \dfrac{1}{F(x)}&\equiv A(x)\pmod {x^n}\\ \\ \dfrac{1}{F(x)}-A(x)&\equiv 0\pmod {x^n}\\ \end{aligned}\]

令 \(G(F(x))=\dfrac{1}{F(x)}-A(x)\),求 \(\dfrac{1}{A(x)}\) 即为求 \(G(F(x))\) 的零点。

代入牛顿迭代的式子:

\[F(x)=F_0(x)-\dfrac{\dfrac{1}{F(x)}-A(x)}{G'(F_0(x))} \]

将 \(G(F(x))\) 中的 \(F(x)\) 看作变量,\(A(x)\) 看作常数,可得 \(G'F(x)=-\dfrac{1}{F(x)^2}\),代入得

\[\begin{aligned} F(x)&=F_0(x)+(\dfrac{1}{F_0(x)}-A(x))\times F_0(x)^2\\ \\ F(x)&=2F_0(x)-A(x)\times F_0(x)^2\\ \end{aligned}\]

初始时 \(F_0(x)=\dfrac{1}{A(0)}\)。


P4512 多项式除法

对于任意 \(n\) 次多项式 \(A(x)\),设 \(A_r(x)=\displaystyle\sum_{i=0}^n[x^{n-i}]A(x)x^i\)。

\(\displaystyle A_r(x)=\sum_{i=0}^n[x^{n-i}]A(x)(\dfrac{1}{x})^{n-i}x^n=x^nA(\dfrac{1}{x})\)

将 \(\dfrac{1}{x}\) 代入题目中式子:

\[\begin{aligned} F(\dfrac{1}{x})&=Q(\dfrac{1}{x})\times G(\dfrac{1}{x})+R(\dfrac{1}{x})\\ \\ x^n\times F(\dfrac{1}{x})&=x^n\times Q(\dfrac{1}{x})\times G(\dfrac{1}{x})+x^n\times R(\dfrac{1}{x})\\ \\ x^nF(\dfrac{1}{x})&=x^{n-m} Q(\dfrac{1}{x})\times x^mG(\dfrac{1}{x})+x^{n-m+1}\times x^{m-1}R(\dfrac{1}{x})\\ \\ F_r(x)&= Q_r(x)\times G_r(x)+x^{n-m+1}R_r(x)\\ \\ F_r(x)&\equiv Q_r(x)\times G_r(x)\pmod {x^{n-m+1}}\\ \\ Q_r(x)&\equiv F_r(x)\times G_r^{-1}(x)\pmod {x^{n-m+1}}\\ \\ \end{aligned}\]

\(Q_r(x)\) 翻转即为 \(Q(x)\),\(R(x)=F(x)-Q(x)\times G(x)\)。


P4725 多项式 ln

\[\begin{aligned} F(x)&\equiv\ln A(x)\\ \\ F'(x)&\equiv\dfrac{1}{A(x)}\times A'(x)\\ \\ F(x)&\equiv\int\dfrac{1}{A(x)}\times A'(x)\text{ d}x\\ \\ \end{aligned}\]


P4726 多项式 exp

\[F(x)\equiv e^{A(x)}\pmod{x^n} \]

两边求对数后移项得

\[\ln F(x)-A(x)\equiv 0\pmod{x^n} \]

令 \(G(F(x))=\ln F(x)-A(x)\),代入牛顿迭代的式子:

\[F(x)=F_0(x)-\dfrac{\ln F(x)-A(x)}{G'(F(x))} \]

\(G'F(x)=F^{-1}(x)\),代入得

\[F(x)=F_0(x)(1-\ln F(x)+A(x)) \]

初始时 \(F_0(x)=e^0=1\)。

标签:推导,pmod,dfrac,times,公式,多项式,aligned,equiv
From: https://www.cnblogs.com/ConvergentSeries/p/17973399

相关文章

  • 多项式相关
    板子:1constintN=4e5+5;23structNTT{4constintmod=998244353,G=3,invG=332748118;5intrev[N],res[N],Inv[N],Res[N],eres[N],Eres[N];67voidpreinv(intn){8Inv[0]=Inv[1]=1;9......
  • 多项式计数杂谈
    多项式计数杂谈-command_block的博客-洛谷博客command_block的博客导航切换首页文章command_block的博客多项式计数杂谈postedon2020-02-1118:16:01|under算法|Ox00.前面的话&前置芝士本文Markd......
  • 多项式求值软件下载Polynomial evaluation software mus 2025 download
    本软件是Windows下64位软件。本软件能计算如a0+a1*x+a2*x^2+......+an*x^n的式子的对b1的求值结果。具体的方法就是在多多项式系数区输入a0到an的值,然后点击计算多项式的结果即可在结果栏算出结果。最大项数为1000项。多项式系数输入时1项1行,从上到下是a0到an,中间不能空行。T......
  • 如何使用Java在Excel中添加动态数组公式?
    本文由葡萄城技术团队发布。转载请注明出处:葡萄城官网,葡萄城为开发者提供专业的开发工具、解决方案和服务,赋能开发者。前言动态数组公式是Excel引入的一项重要功能,它将Excel分为两种风格:Excel365和传统Excel(2019或更早版本)。动态数组功能允许用户从单个单元格中的公式......
  • 切比雪夫多项式
    切比雪夫多项式通常我们使用切比雪夫多项式时都在范围[-1,1]之间。定义切比雪夫多项式在[-1,1]上的定义是:\(T_n(x)=cos(narccos(x)),-1\leqx\leq1\),其中,T_n(x)是阶数为n的切比雪夫多项式。性质\(T_n(x)\)是n阶多项式。\(T_n(x)\)的奇偶性和n的奇偶性一致。\(T_n(x)\)在区......
  • Python列表差异值统计:集合操作、列表推导式、对称差集详解
     在Python中,统计两个列表的差异值有多种方法,其中包括使用集合操作、列表推导式等。下面我将通过实例详细讲解几种常见的方法,并提供相应的实例源代码。方法一:使用集合操作list1=[1,2,3,4,5]list2=[3,4,5,6,7]#找到在list1中而不在list2中的元素difference1......
  • 多项式基础
    FFT/NTT板子,就不说了。分治FFT给定序列\(g_{1\dotsn-1}\),求序列\(f_{0\dotsn-1}\)。其中\(f_i=\sum_{j=1}^if_{i-j}g_j\),边界为\(f_0=1\)。Solution1朴素地做是\(\mathcal{O}(n^2)\)的。观察到后面的项依赖前面的项。于是我们可以考虑用cdq分治实现。具体......
  • 多项式定积分计算软件2025 64位WIN版下载Polynomial definite integral calculation s
    本软件功能强大,价格实惠,欢迎试用本软件的WIN64位版本。本软件能计算如a0+a1*x+a2*x^2+......+an*a^n的式子的对b1和b2的积分的结果。具体的方法就是在多多项式系数区输入a0到an的值,然后点击计算积分结果即可在结果栏算出结果。最大项数为1000项。多项式系数输入时1项1行,从上......
  • python经典有序序列的list列表推导式
    生成一个数据列表#初始化一个列表list_1=[]#使用循环生成一个列表数据forvalinrange(0,20,1):#加入集合list_1.append(val)#打印列表数据print(list_1)#[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19]使用列表推导式生......
  • 2d物理引擎学习 - 基于约束的公式解决接触稳定性问题
    先看下直接用弹性碰撞的公式,会出现的问题:Box落在地面上后,没有停在地面上,而是还在不断的下沉。弹性碰撞公式处理碰撞后弹开没有大问题,但是处理物体碰撞后的接触存在不稳定问题。 如何解决?目前物理引擎最主流的解决方法是:基于约束来组织物理公式,而不是直接套用物理公式。什......