結果

問題 No.1938 Lagrange Sum
ユーザー pockynypockyny
提出日時 2022-05-13 23:26:39
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 558 ms / 3,000 ms
コード長 1,251 bytes
コンパイル時間 907 ms
コンパイル使用メモリ 74,128 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-22 03:46:25
合計ジャッジ時間 8,819 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 557 ms
6,812 KB
testcase_01 AC 558 ms
6,940 KB
testcase_02 AC 556 ms
6,940 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 2 ms
6,944 KB
testcase_05 AC 348 ms
6,944 KB
testcase_06 AC 389 ms
6,940 KB
testcase_07 AC 2 ms
6,940 KB
testcase_08 AC 79 ms
6,940 KB
testcase_09 AC 26 ms
6,940 KB
testcase_10 AC 516 ms
6,944 KB
testcase_11 AC 88 ms
6,940 KB
testcase_12 AC 520 ms
6,940 KB
testcase_13 AC 526 ms
6,944 KB
testcase_14 AC 182 ms
6,944 KB
testcase_15 AC 16 ms
6,940 KB
testcase_16 AC 492 ms
6,940 KB
testcase_17 AC 169 ms
6,940 KB
testcase_18 AC 171 ms
6,940 KB
testcase_19 AC 492 ms
6,944 KB
testcase_20 AC 501 ms
6,940 KB
testcase_21 AC 7 ms
6,944 KB
testcase_22 AC 1 ms
6,944 KB
testcase_23 AC 6 ms
6,940 KB
testcase_24 AC 484 ms
6,944 KB
testcase_25 AC 2 ms
6,940 KB
testcase_26 AC 2 ms
6,940 KB
testcase_27 AC 2 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <atcoder/modint>

using namespace std;
using namespace atcoder;
using mint = modint998244353;
mint x[10010],y[10010],z[10010],sumL[10010],sumR[10010],dif[10010];
int main(){
    int i,j,n; cin >> n;
    int a; cin >> a; mint X = a;
    for(i=0;i<n;i++){
        int b,c; cin >> b >> c; x[i] = b; y[i] = c;
    }
    for(i=0;i<n;i++){
        dif[i] = 1;
        for(j=0;j<n;j++){
            if(j==i) continue;
            dif[i] *= (x[i] - x[j]);
        }
    }
    sumL[0] = 1;
    for(i=1;i<=n;i++) sumL[i] = sumL[i - 1]*(X - x[i - 1]);
    sumR[n] = 1;
    for(i=n - 1;i>=0;i--) sumR[i] = sumR[i + 1]*(X - x[i]);
    mint ans = 0;
    for(i=0;i<n;i++){
        ans += y[i]*sumL[i]*sumR[i + 1]/dif[i];
        //cout << ans.val() << endl;
    }
    ans *= (mint)(n - 1);
    //cout << ans.val() << endl;
    for(i=0;i<n;i++) z[i] = y[i]/dif[i];
    mint sum = 0;
    for(i=0;i<n;i++) sum += z[i];
    for(i=0;i<n;i++){
        ans -= (sum - z[i])*sumL[i]*sumR[i + 1];
    }
    cout << ans.val() << "\n";
    /*for(i=0;i<n;i++) cout << dif[i].val() << " ";
    cout << endl;
    for(i=0;i<=n;i++) cout << sumL[i].val() << " ";
    cout << endl;
    for(i=0;i<=n;i++) cout << sumR[i].val() << " ";
    cout << endl;*/
}
0