結果

問題 No.1938 Lagrange Sum
ユーザー noya2noya2
提出日時 2021-08-10 18:53:13
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 568 ms / 3,000 ms
コード長 994 bytes
コンパイル時間 4,434 ms
コンパイル使用メモリ 265,072 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-22 00:01:30
合計ジャッジ時間 13,674 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 568 ms
5,248 KB
testcase_01 AC 568 ms
5,376 KB
testcase_02 AC 568 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 1 ms
5,376 KB
testcase_05 AC 358 ms
5,376 KB
testcase_06 AC 388 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 81 ms
5,376 KB
testcase_09 AC 26 ms
5,376 KB
testcase_10 AC 530 ms
5,376 KB
testcase_11 AC 90 ms
5,376 KB
testcase_12 AC 525 ms
5,376 KB
testcase_13 AC 534 ms
5,376 KB
testcase_14 AC 185 ms
5,376 KB
testcase_15 AC 16 ms
5,376 KB
testcase_16 AC 499 ms
5,376 KB
testcase_17 AC 174 ms
5,376 KB
testcase_18 AC 178 ms
5,376 KB
testcase_19 AC 504 ms
5,376 KB
testcase_20 AC 498 ms
5,376 KB
testcase_21 AC 9 ms
5,376 KB
testcase_22 AC 2 ms
5,376 KB
testcase_23 AC 6 ms
5,376 KB
testcase_24 AC 497 ms
5,376 KB
testcase_25 AC 2 ms
5,376 KB
testcase_26 AC 2 ms
5,376 KB
testcase_27 AC 2 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

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