拉格朗日插值法

本文介绍了如何使用拉格朗日插值法来构造一个通过特定点集的多项式。通过定义特殊的基函数,确保了每个基函数只在一个指定的点取非零值,从而实现了对任意一组数据点的精确插值。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

介绍

已知 n n n个点 ( x i , y i ) (x_i,y_i) (xi,yi)时,可以唯一地确定一个次数不超过 n n n的多项式 y = f ( x ) y=f(x) y=f(x)。但要怎么求呢?

在其他数学家绞尽脑汁的时候,拉格朗日想到:既然这是唯一的,为什么我不能设计 n n n个函数,使 p i ( x ) p_i(x) pi(x)只在 x = x i x=x_i x=xi时为 y i y_i yi x = x j x=x_j x=xj时值为 0 ( j ≠ i ) 0(j\neq i) 0(j=i),再将其求和呢?

说干就干,令
p i ( x ) = y i ∏ j = 1 , j ≠ i n x − x j x i − x j p_i(x)=y_i\prod\limits_{j=1,j\neq i}^n \dfrac{x-x_j}{x_i-x_j} pi(x)=yij=1,j=inxixjxxj

那么当 x = x i x=x_i x=xi时, p i ( x ) = y i p_i(x)=y_i pi(x)=yi;当 x = x j ( j ≠ i ) x=x_j(j\neq i) x=xj(j=i)时, p i ( x ) = 0 p_i(x)=0 pi(x)=0

那么,很显然地, f ( x ) f(x) f(x)即为所有的 p i ( x ) p_i(x) pi(x)之和。

f ( x ) = ∑ i = 1 n p i ( x ) = ∑ i = 1 n y i ∏ j = 1 , j ≠ i n x − x j x i − x j f(x)=\sum\limits_{i=1}^np_i(x)=\sum\limits_{i=1}^ny_i\prod\limits_{j=1,j\neq i}^n\dfrac{x-x_j}{x_i-x_j} f(x)=i=1npi(x)=i=1nyij=1,j=inxixjxxj

此时对于所有的 i ( 1 ≤ i ≤ n ) , f ( x i ) = y i i(1\leq i\leq n),f(x_i)=y_i i(1in),f(xi)=yi。又因为这个多项式是唯一的,所以 f ( x ) f(x) f(x)即为所求。


例题

P4781 拉格朗日插值

code

#include<bits/stdc++.h>
using namespace std;
int n;
long long k,ans=0,p,q,x[2005],y[2005];
long long mod=998244353;
long long mi(long long t,long long v){
	if(v==0) return 1;
	long long re=mi(t,v/2);
	re=re*re%mod;
	if(v&1) re=re*t%mod;
	return re;
}
int main()
{
	scanf("%d%lld",&n,&k);
	for(int i=1;i<=n;i++){
		scanf("%lld%lld",&x[i],&y[i]);
	}
	for(int i=1;i<=n;i++){
		p=q=1;
		for(int j=1;j<=n;j++){
			if(i==j) continue;
			p=p*(k-x[j])%mod;
			q=q*(x[i]-x[j])%mod;
		}
		ans=(ans+p*mi(q,mod-2)%mod*y[i]%mod+mod)%mod;
	}
	printf("%lld",ans);
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值