結果
| 問題 | 
                            No.1510 Simple Integral
                             | 
                    
| コンテスト | |
| ユーザー | 
                             magsta
                         | 
                    
| 提出日時 | 2021-03-28 19:20:02 | 
| 言語 | C++14  (gcc 13.3.0 + boost 1.87.0)  | 
                    
| 結果 | 
                             
                                WA
                                 
                             
                            
                         | 
                    
| 実行時間 | - | 
| コード長 | 965 bytes | 
| コンパイル時間 | 734 ms | 
| コンパイル使用メモリ | 84,740 KB | 
| 実行使用メモリ | 6,948 KB | 
| 最終ジャッジ日時 | 2024-09-14 22:57:32 | 
| 合計ジャッジ時間 | 2,038 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge4 / judge6 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | WA * 3 | 
| other | WA * 43 | 
ソースコード
#include <iostream>
#include <string>
#include <algorithm>
#include <utility>
#include <cmath>
#include <cstdlib>
#include <map>
#include <set>
#include <queue>
#include <vector>
using namespace std;
#define MOD 1000000007
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; cin >> N;
	vector<long long> A(N);
	for (int i = 0; i < N; i++) cin >> A[i];
	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 * modinv(2, MOD) % MOD << endl;
}
            
            
            
        
            
magsta