拉格朗日插值算法

首先介绍下作用

在平面上给定 n 个点,求一条多项式图像穿过所有的点

(每个点横坐标不同)

由浅入深,首先想一个悬浮点 1 ,它在除 1 以外给定的点的横坐标上函数值都是 0 。已知它的横坐标 \(x_1\) , 怎么用函数图像表示

我们一拍脑袋发现,可以这么表示

\[f1(x) = \prod_{i = 2}^n {x - x_i} \]

这一坨是不是在给定点的横坐标中,只当横坐标是点 \(x_1\) 的时候不为 0

然而实际上这个东西的纵坐标不太可控,所以我们把这玩意再除一下,让横坐标是 \(x_1\) 时纵坐标为 1, 方便后续操作

\[f1(x) = \prod_{i = 2}^n \frac{x - x_i}{x_1 - x_i} \]

下一个目标其实也简单,就是把这个点的纵坐标改变,比如我们现在已知这个点 1 的纵坐标为 \(y_1\) ,式子就变成

\[f1(x) = y_1 \times \prod_{i = 2}^n \frac{x - x_i}{x_1 - x_i} \]

更具有普遍性,一个点 i , 横坐标 \(x_i\) , 纵坐标 \(y_i\),以下这个函数的图像一定经过它

\[f(x) = y_i \times \prod_{j \ne i}^{n} \frac{x - x_j}{x_i - x_j} \]

那么我们可以列举一大堆函数,这些函数都是只在单个点的横坐标上取值为纵坐标,其他给定点上函数值都是 0

比如 刚才的 \(f1(x)\) ,给定点中只在第一个点上取得值是 \(y_i\) ,其他都是 0

那我们是不是就可以用一个求和的形式了捏 ? 这样他就会在不同给定点的坐标自动给出值

\[f(k) = \sum_{i = 1}^{n} y_i \times \prod_{j \ne i}^{n} \frac{k - x_j}{x_i - x_j} \]

这就是经过 n 个点的函数。说到这里,其实板子已经能做了

洛谷P4781 拉格朗日插值

放下代码罢,马蜂诡奇,勿喷拜谢

#include <bits/stdc++.h>#define mod 998244353
#define ri register int
#define N 2007
#define int long longint x[N] , y[N];
namespace P4781 {int n , k , sum = 0;inline int Pow(int a , int b) {int ans = 1;while(b) {if(b & 1) {ans = ans * a % mod;}a = a * a % mod;b >>= 1;}return ans % mod;}void solve() {sum = 0;scanf("%lld%lld",&n,&k);for(ri i = 1; i <= n; ++i) {scanf("%lld%lld", &x[i] , &y[i]);}for(ri i = 1; i <= n; ++i) {int res = 1;for(ri j = 1; j <= n; ++j) {if(i == j) {continue;}res = (res * (k - x[j] + mod) % mod * Pow(x[i] - x[j] , mod - 2) % mod + mod) % mod;}res = res * y[i] % mod;sum = (sum + res) % mod;}printf("%lld",sum % mod);}
}signed main() {P4781 :: solve();return 0;
}

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

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

相关文章

本地部署Gemma3模型

本地部署Gemma3模型 1. 拉取并运行 Ollama 容器 # 拉取 Ollama 镜像 docker pull ollama/ollama# 运行 Ollama 容器 docker run -d --gpus all -v ollama:/root/.ollama -p 11434:11434 --name ollama ollama/ollama2. 进入容器并部署 Gemma3 模型(默认4B版 ) # 进入 Ollama 容…

Anaconda中启动Jupyter lab的方法

Anaconda中启动Jupyter lab的方法: 1、在Anaconda Prompt中,使用命令来启动:jupyter lab 2、在Anaconda界面中点击Jupyter lab下方的launch按钮启动,如下图所示:

React+Next.js+MaterialUI+Toolpad技术栈学习——安装

今天跟大家分享一个React+Next.js+MaterialUI技术栈的前端框架Toolpad。相关资源MaterialUI Toolpad框架效果安装运行安装命令npx create-toolpad-app@latest your-app cd your-app npm run dev文件结构 无身份认证 ├── app │ ├── (dashboard) │ │ ├── layou…

测试驱动开发(TDD)浅析

测试驱动开发(TDD:Test Driven Development)是敏捷开发中的一项核心实践,推崇通过测试来驱动整个开发的进行。TDD有别于传统“先编码,后测试”的开发过程,而是要求在编写业务代码之前,先编写测试用例。TDD的概念大致在上世纪90年代随着极限编程(XP:Extreme Programmin…

(18).命令模式

命令模式 命令模式的核心思想是将请求封装为个对象,将其作为命令发起者和接收者的中介,而抽象出来的命令对象又使得能够对一系列请求进行操作,如对请求进行排队,记录请求日志以及支持可撤销的操作等。命令模式参与者:◇命令的执行者(接收者Receiver):它单纯的只具体实现了功…

学嵌入式C语言,看这一篇就够了(5)

C语言的运算符 学习编程语言,应该遵循“字-->词-->句-->段--->章”,对于一条有意义的语句而言,是离不开标点符号的运算符指明要进行的运算和操作,操作数是指运算符的操作对象,根据运算符操作数的数目不同,C语言标准把运算符分为三种:单目运算符(一元运算符…

20242313 2024-2025-2 《Python程序设计》实验一报告

20242313 2024-2025-2 《Python程序设计》实验一报告 课程:《Python程序设计》 班级:2423 姓名:曾海鹏 学号:20242313 实验教师:王志强 实验日期:2025年3月16日 必修/选修:公选课 1.实验内容 1.熟悉Python开发环境; 2.练习Python运行、调试技能;(编写书中的程序,并…

nn.Embedding()函数详解

nn.Embedding()函数详解 nn.Embedding()函数:随机初始化词向量,词向量在正态分布N(0,1)中随机取值 输入: torch.nn.Embedding(num_embeddings, embedding_dim, padding_idx=None, max_norm=None, norm_type=2.0, scale_grad_by_freq=False, sparse=False, _weight=None) num…

htb Authority

端口扫描 nmap -sC -sV -p- -Pn -T4 10.10.11.222 Starting Nmap 7.92 ( https://nmap.org ) at 2024-10-04 19:42 CST Nmap scan report for 10.10.11.222 (10.10.11.222) Host is up (0.40s latency). Not shown: 65506 closed tcp ports (reset) PORT STATE SERVICE …

蓝桥杯14届省B

蓝桥杯14届省赛B组A:int a[105]; int day[]={0,31,28,31,30,31,30,31,31,30,31,30,31};//记录每个月有多少天 set<int> st;//记录不重复的日期void check(int mm,int dd){if (mm>12||mm<1||dd<1||dd>day[mm]) return;else st.insert(mm*100+dd);//st存日期 …

docker 安装 oracle database 问题记录

pre本地docker (WSL)安装运行 Oracle1. 镜像处理参考链接:https://www.cnblogs.com/wuchangsoft/p/18344847 oracle 镜像获取:https://container-registry.oracle.com/ords/f?p=113:10:::::: (Oracle官网,由于部分问题导致直接pull无法拉取) 阿里云,参考链接里有个个人19…