結果

問題 No.1938 Lagrange Sum
ユーザー noya2noya2
提出日時 2021-08-10 18:53:13
言語 C++17
(gcc 11.2.0 + boost 1.78.0)
結果
AC  
実行時間 589 ms / 3,000 ms
コード長 994 bytes
コンパイル時間 3,874 ms
使用メモリ 9,772 KB
最終ジャッジ日時 2023-02-21 11:56:48
合計ジャッジ時間 13,690 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge13
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 589 ms
7,564 KB
testcase_01 AC 585 ms
7,596 KB
testcase_02 AC 584 ms
7,596 KB
testcase_03 AC 2 ms
7,584 KB
testcase_04 AC 1 ms
7,720 KB
testcase_05 AC 368 ms
7,592 KB
testcase_06 AC 405 ms
7,684 KB
testcase_07 AC 2 ms
7,612 KB
testcase_08 AC 82 ms
7,668 KB
testcase_09 AC 26 ms
9,684 KB
testcase_10 AC 538 ms
7,576 KB
testcase_11 AC 91 ms
9,772 KB
testcase_12 AC 543 ms
7,596 KB
testcase_13 AC 552 ms
7,612 KB
testcase_14 AC 190 ms
9,692 KB
testcase_15 AC 16 ms
7,596 KB
testcase_16 AC 507 ms
7,660 KB
testcase_17 AC 179 ms
9,752 KB
testcase_18 AC 182 ms
7,652 KB
testcase_19 AC 522 ms
7,592 KB
testcase_20 AC 522 ms
9,700 KB
testcase_21 AC 8 ms
7,568 KB
testcase_22 AC 2 ms
9,700 KB
testcase_23 AC 7 ms
9,692 KB
testcase_24 AC 510 ms
7,648 KB
testcase_25 AC 2 ms
7,656 KB
testcase_26 AC 1 ms
7,564 KB
testcase_27 AC 1 ms
7,604 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