P4781-[模板]拉格朗日插值
生活随笔
收集整理的這篇文章主要介紹了
P4781-[模板]拉格朗日插值
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
正題
題目鏈接:https://www.luogu.com.cn/problem/P4781
題目大意
給出nnn個點確定一個多項式y=f(x)y=f(x)y=f(x),求f(k)f(k)f(k)的值
解題思路
拉格朗日插值:
f(k)=∑i=1nyi∏i≠jk?xjxi?xjf(k)=\sum_{i=1}^ny_i\prod_{i\neq j}\frac{k-x_j}{x_i-x_j}f(k)=i=1∑n?yi?i?=j∏?xi??xj?k?xj??
時間復雜度O(n2)O(n^2)O(n2)
codecodecode
#include<cstdio> #include<cstring> #include<algorithm> #define ll long long using namespace std; const ll N=2e3+10,XJQ=998244353; ll n,k,x[N],y[N],ans; ll power(ll x,ll b){x%=XJQ;ll ans=1;while(b){if(b&1)ans=ans*x%XJQ;x=x*x%XJQ;b>>=1;}return ans; } int main() {scanf("%lld%lld",&n,&k);for(ll i=1;i<=n;i++)scanf("%lld%lld",&x[i],&y[i]);for(ll i=1;i<=n;i++){for(ll j=1;j<=n;j++)if(i!=j)(y[i]*=(k-x[j]+XJQ)%XJQ*power(x[i]-x[j]+XJQ,XJQ-2)%XJQ)%=XJQ;(ans+=y[i])%=XJQ; }printf("%lld",ans); }總結
以上是生活随笔為你收集整理的P4781-[模板]拉格朗日插值的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 游戏多开的电脑配置高吗(游戏多开的电脑配
- 下一篇: jzoj4049-排序【搜索】