結果

問題 No.1938 Lagrange Sum
ユーザー pockynypockyny
提出日時 2022-05-13 23:26:39
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 576 ms / 3,000 ms
コード長 1,251 bytes
コンパイル時間 626 ms
コンパイル使用メモリ 73,896 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-09-29 09:09:30
合計ジャッジ時間 9,363 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge14
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 575 ms
4,380 KB
testcase_01 AC 573 ms
4,376 KB
testcase_02 AC 576 ms
4,376 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 1 ms
4,376 KB
testcase_05 AC 362 ms
4,380 KB
testcase_06 AC 397 ms
4,376 KB
testcase_07 AC 2 ms
4,376 KB
testcase_08 AC 81 ms
4,376 KB
testcase_09 AC 26 ms
4,380 KB
testcase_10 AC 530 ms
4,376 KB
testcase_11 AC 90 ms
4,376 KB
testcase_12 AC 535 ms
4,376 KB
testcase_13 AC 541 ms
4,376 KB
testcase_14 AC 187 ms
4,380 KB
testcase_15 AC 16 ms
4,380 KB
testcase_16 AC 498 ms
4,376 KB
testcase_17 AC 177 ms
4,376 KB
testcase_18 AC 180 ms
4,376 KB
testcase_19 AC 514 ms
4,376 KB
testcase_20 AC 513 ms
4,376 KB
testcase_21 AC 8 ms
4,380 KB
testcase_22 AC 2 ms
4,380 KB
testcase_23 AC 6 ms
4,376 KB
testcase_24 AC 501 ms
4,380 KB
testcase_25 AC 2 ms
4,380 KB
testcase_26 AC 2 ms
4,380 KB
testcase_27 AC 2 ms
4,380 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