ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

拉格朗日插值模板题 luoguP4871

2020-04-26 10:05:34  阅读:255  来源: 互联网

标签:const 插值 res int 1ll luoguP4871 include 模板 mod


学习博客

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<map>
#include<queue>
#include<vector>
#include<string>
#include<fstream>

using namespace std;
typedef long long ll;
const int N = 4e6 + 105;
const int mod = 998244353;
const double Pi = acos(- 1.0);
const int INF = 0x3f3f3f3f;

int n, k;
int x[N], y[N];

int sub(const int& a, const int& b){
	int t = a - b;
	return t < 0 ? t + mod : t;
}


int inv(int x){
	int res = 1;
	int p = mod - 2;
	for(; p; p >>= 1, x = 1ll * x * x % mod)
		if(p & 1) res = 1ll * res * x % mod;
	return res;
}


int Lagrange(int n, int x[], int y[], int k){
	int res = 0;
	for(int i = 1; i <= n; ++ i){
		int s1 = 1, s2 = 1;	//s1:分子,s2:分母
		for(int j = 1; j <= n; ++ j){
			if(i != j){
				s1 = 1ll * s1 * sub(k, x[j]) % mod;
				s2 = 1ll * s2 * sub(x[i], x[j]) % mod;
			}
		}
		res = (res + 1ll * y[i] * s1 % mod * inv(s2) % mod) % mod; 
	}
	return res;
}

int main()
{
	scanf("%d%d",&n,&k);
	for(int i = 1; i <= n; ++ i) scanf("%d%d",&x[i],&y[i]);
	printf("%d\n", Lagrange(n, x, y, k));
	return 0;
}

标签:const,插值,res,int,1ll,luoguP4871,include,模板,mod
来源: https://www.cnblogs.com/A-sc/p/12777358.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有