日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

P4781-[模板]拉格朗日插值

發布時間:2023/12/3 编程问答 35 豆豆
生活随笔 收集整理的這篇文章主要介紹了 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=1n?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-[模板]拉格朗日插值的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。