結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff #

#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;
}
0