首页 > 其他分享 >拉格朗日插值

拉格朗日插值

时间:2024-09-20 19:45:18浏览次数:8  
标签:拉格朗 frac 插值 多项式 sum prod sim neq

  • 应用范围:求一个 \(n\) 次多项式过 \((x_1,y_1)\sim (x_n,y_n)\)

  • 构造思想:

    • 设 \(f_i(x)\) 使得 对于 \(x_i\neq x_j\),\(f(x_j)=0\),且 \(f(x_i)=1\),注意并不是对全体 \(R\) 满足。
    • 由上 \(F(x)=\sum y_if_i(x)\) 即为所求。
    • 构造方法:\(f_i(x)=\frac{\prod_{j=1,j\neq i}^n(x-x_j)}{\prod_{j=1,j\neq i}^n(x_i-x_j)}\)

也就是:

\[\begin{aligned} F(x)=\sum_{i=1}^ny_i\prod_{j\neq i}\frac{(x-x_j)}{(x_i-x_j)} \end{aligned} \]

计算时可以直接带值 \(O(n^2)\) 计算,也可以利用背包的操作先求出 \(\prod_{i=1}^n(x-x_i)\) 的具体系数 \(g_0\sim g_n\),然后除掉一个 \((x-x_i)\) 后的 \(g'\),其实相当于 \(-x_ig'_k+g'_{k-1}\)。根据这个式子可以直接确定 \(g'_0\),进而逐步确定 \(g’_{1\sim n}\)

在实际应用中 往往维护 \(x_i\) 对应的 \(y_i\),也就是点值,最后一步再插回原多项式

我们可以在 \(O(n^2)\) 的复杂度内还原原多项式。

例题 CF995F,code

标签:拉格朗,frac,插值,多项式,sum,prod,sim,neq
From: https://www.cnblogs.com/spdarkle/p/18423178

相关文章

  • 牛顿插值法-C++【可直接复制粘贴/欢迎评论点赞】
    牛顿插值法是一种基于给定数据点集构造插值多项式的方法,用于近似未知函数的值。该方法通过构造差商表并利用该表逐步构建插值多项式。相较于拉格朗日插值法,牛顿插值法的一个显著优势是,当需要增加插值点时,只需重附上一项即可,无需重新计算所有插值点的值。基本概念牛顿插值法的......
  • vue优点/插值表达式/强制绑定
    1.Vue.js的优点体积小:压缩后只有33k;更高的运行效率:基于虚拟DOM,一种可以预先通过JavaScript进行各种计算,把最终的DOM操作计算出来并优化的技术,由于这种DOM操作属于预处理操作,并没有真实的操作DOM,所以叫做虚拟DOM;双向数据绑定:让开发者不用再去操作DOM对象,把更多的精力投入到业务......
  • Lagrange 插值
    给定\(n\)个横坐标不同的点,求过这\(n\)个点的\(n-1\)次多项式。算法引入这可以直接用高斯消元做,但是时间复杂度\(\mathcalO(n^3)\)不可接受,我们需要优化。我们令\((x_1,y_1),(x_2,y_2),\dots,(x_t,y_t)\)为这些点。考虑构造一个函数\(\ell_j(x)\)满足\[\ell_......
  • 如何通过插值法计算cell delay?
            我们知道,celldelay是根据inputtransition和outputload计算得到的。如图所示,为X8驱动的buffer的timing查找表。由于buffer是正单边类型cell,那么当一个1->0 翻转的信号经过buffer时,计算timingdelay应该去查找cell_fall这个表格。假设inputtransition为0.......
  • Numba最近邻插值(CPU+ GPU + Z轴切块 + XYZ轴切块 + 多线程)
    文章目录最近邻插值(加速方法)(1)scipy.ndimage.zoom(2)Numba-CPU加速(3)Numba-GPU加速(4)Numba-CPU加速(Z轴切块)(5)Numba-CPU加速(XYZ轴切块)(6)Numba-CPU加速(XYZ轴切块)+多线程输入数据插值倍数时耗scipy.ndimage.zoom(1024,1024,512)4172.16sNumba-CPU(1024,1024,512)456.58sN......
  • 拉格朗日插值
    CSR:又拉又插的东西(又垃圾,又傻叉)JCY:什么你拉插了一晚上?(我学习拉插学了一晚上)什么是拉插给定一些点值,是否可以求出一个函数,使得函数图像穿过这些点,并求出给定的\(x\)所对应的\(y\)。初步思路前置我们想到两个点肯定可以确定一条直线,而三个点肯定可以确定一条抛物线,所以我们......
  • matlab中的插值与拟合(代码)
    目录1.对均匀数据的插值与拟合2.对散点数据的拟合(如ANSYSfluent导出的节点数据)1.对均匀数据的插值与拟合interp1:一维插值。这是最常用的插值函数之一,用于对一维数据进行插值。它可以执行线性插值、最近邻插值、样条插值等多种类型的插值。%已知数据点x=1:5;y......
  • 拉格朗日插值优化 DP 做题笔记
    本来想在洛谷题单里找斜率优化DP的,然后发现了一个拉格朗日插值优化DP的题单,就点进去尝试了一下。题单。于是先看了雨兔的题解,学了CF995F的做法,然后A了这个题。雨兔题解的链接和我的代码见CF上的提交记录。现在正在做后面的题。P3643[APIO2016]划艇\(a_i,b_i......
  • [转]插值-样条插值 - 努力的孔子 - 博客园
    百度百科定义插值:在离散数据的基础上插补连续函数,使得这条连续曲线经过全部离散点,同时也可以估计出函数在其他点的近似值。样条插值:一种以可变样条来作出一条经过一系列点的光滑曲线的数学方法。插值样条是由一些多项式组成的,每一个多项式都是由相邻的两个数据点决定的,这样,任......
  • 有符号浮点运算的基本步骤:以双线性插值为例
    参考:韩彬的图像处理书、无双软件学院方法。步骤一:无损定点化浮点数在硬件计算中首先需要做的便是定点化,一般是左移一定位宽,可以是2048或4096;这个过程要注意保障无损;步骤二:运算和位宽匹配;要确定所有参与计算的数小数位位宽是匹配的,否则无法进行任何层次的计算;需要特别注意很......