結果
| 問題 |
No.1938 Lagrange Sum
|
| コンテスト | |
| ユーザー |
noya2
|
| 提出日時 | 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;
}
noya2