結果

問題 No.1510 Simple Integral
ユーザー magstamagsta
提出日時 2021-04-06 21:06:41
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 1,248 bytes
コンパイル時間 6,973 ms
コンパイル使用メモリ 226,096 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-09-14 22:58:25
合計ジャッジ時間 8,302 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample RE * 3
other RE * 43
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <string>
#include <algorithm>
#include <utility>
#include <cmath>
#include <cstdlib>
#include <map>
#include <set>
#include <queue>
#include <vector>
#include "testlib.h"
using namespace std;
#define MOD 998244353

long long modinv(long long a, long long mod) {
    long long b = mod, u = 1, v = 0;
    while (b) {
        long long t = a / b;
        a -= t * b; swap(a, b);
        u -= t * v; swap(u, v);
    }
    u %= mod;
    if (u < 0) u += mod;
    return u;
}

int main() {
    long long N = inf.readInt(1, 1000);
    inf.readEoln();
    vector<long long> A(N);
    A[0] = inf.readInt(1, 1000000);
    for (int i = 1; i < N; i++) {
        inf.readSpace();
        A[i] = inf.readInt(1, 1000000);
    }
    inf.readEoln();
    inf.readEof();
    for (int i = 0; i < N - 1; i++) {
        if (A[i] >= A[i + 1]) cout << -1 << endl;
    }
    long long sum = 0;
    for (long long i = 0; i < N; i++) {
        long long mul = modinv(A[i], MOD);
        for (long long j = 0; j < N; j++) {
            if (i == j) continue;
            mul = mul * modinv(abs(A[j] * A[j] - A[i] * A[i]) % MOD, MOD) % MOD;
        }
        if (i % 2 == 1) mul = MOD - mul;
        sum += mul;
    }
    cout << sum % MOD << endl;
}
0