結果
問題 | No.1938 Lagrange Sum |
ユーザー |
![]() |
提出日時 | 2021-08-10 18:53:13 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 559 ms / 3,000 ms |
コード長 | 994 bytes |
コンパイル時間 | 3,794 ms |
コンパイル使用メモリ | 253,572 KB |
最終ジャッジ日時 | 2025-01-23 17:24:21 |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 25 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/all> #define rep(i,n) for (int i = 0; i < int(n); ++i) #define repb(i,n) for (int i = int(n)-1; i >= 0; --i) using namespace std; using namespace atcoder; using ll = long long; using mint = modint998244353; int main(){ int n; cin >> n; ll X; cin >> X; mint z = X; vector<mint> x(n); vector<mint> y(n); rep(i,n){ ll a, b; cin >> a >> b; x[i] = a; y[i] = b; } vector<mint> lrui(n+1,1); vector<mint> rrui(n+1,1); rep(i,n) lrui[i+1] = lrui[i] * (z - x[i]); repb(i,n) rrui[i] = rrui[i+1] * (z - x[i]); vector<mint> dxs(n,1); rep(i,n) rep(j,n){ if (i == j) continue; dxs[i] *= x[i] - x[j]; } rep(i,n) dxs[i] = mint(1) / dxs[i]; mint a = 0, b = 0, c = 0; rep(i,n) a += y[i] * lrui[i] * rrui[i+1] * dxs[i]; a *= n; rep(i,n) b += y[i] * dxs[i]; rep(i,n) c += lrui[i] * rrui[i+1]; mint ans = a - b * c; cout << ans.val() << endl; }