結果
問題 | No.1844 Divisors Sum Sum |
ユーザー | lloyz |
提出日時 | 2022-02-18 23:31:31 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 922 ms / 3,000 ms |
コード長 | 758 bytes |
コンパイル時間 | 371 ms |
コンパイル使用メモリ | 86,940 KB |
実行使用メモリ | 78,904 KB |
最終ジャッジ日時 | 2023-09-11 20:10:38 |
合計ジャッジ時間 | 16,430 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 858 ms
78,904 KB |
testcase_01 | AC | 912 ms
78,436 KB |
testcase_02 | AC | 922 ms
78,564 KB |
testcase_03 | AC | 916 ms
78,268 KB |
testcase_04 | AC | 922 ms
78,180 KB |
testcase_05 | AC | 908 ms
78,432 KB |
testcase_06 | AC | 910 ms
78,520 KB |
testcase_07 | AC | 908 ms
78,456 KB |
testcase_08 | AC | 919 ms
78,552 KB |
testcase_09 | AC | 910 ms
78,660 KB |
testcase_10 | AC | 841 ms
78,648 KB |
testcase_11 | AC | 266 ms
78,628 KB |
testcase_12 | AC | 666 ms
78,608 KB |
testcase_13 | AC | 185 ms
78,328 KB |
testcase_14 | AC | 367 ms
78,380 KB |
testcase_15 | AC | 69 ms
71,272 KB |
testcase_16 | AC | 70 ms
71,268 KB |
testcase_17 | AC | 71 ms
71,208 KB |
testcase_18 | AC | 70 ms
71,224 KB |
testcase_19 | AC | 71 ms
70,840 KB |
testcase_20 | AC | 70 ms
71,168 KB |
testcase_21 | AC | 71 ms
71,032 KB |
testcase_22 | AC | 70 ms
71,220 KB |
testcase_23 | AC | 69 ms
71,216 KB |
testcase_24 | AC | 71 ms
71,196 KB |
testcase_25 | AC | 71 ms
70,836 KB |
testcase_26 | AC | 71 ms
71,012 KB |
testcase_27 | AC | 71 ms
71,276 KB |
testcase_28 | AC | 70 ms
71,176 KB |
testcase_29 | AC | 71 ms
71,272 KB |
testcase_30 | AC | 70 ms
71,036 KB |
testcase_31 | AC | 70 ms
71,168 KB |
testcase_32 | AC | 72 ms
71,008 KB |
testcase_33 | AC | 70 ms
71,124 KB |
testcase_34 | AC | 69 ms
71,268 KB |
testcase_35 | AC | 71 ms
71,276 KB |
testcase_36 | AC | 71 ms
70,840 KB |
testcase_37 | AC | 71 ms
71,132 KB |
ソースコード
def xgcd(a, b): x0, y0, x1, y1 = 1, 0, 0, 1 while b != 0: q, a, b = a // b, b, a % b x0, x1 = x1, x0 - q * x1 y0, y1 = y1, y0 - q * y1 return a, x0, y0 def modinv(a, m): g, x, y = xgcd(a, m) if g != 1: raise Exception('modular inverse does not exist') else: return x % m n = int(input()) mod = 10**9 + 7 ans = 1 for _ in range(n): p, e = map(int, input().split()) temp1 = (e + 1) * (pow(p, e + 1, mod) - 1) % mod temp1 *= modinv(p - 1, mod) temp1 %= mod temp2 = 1 - ((e + 1) * pow(p, e, mod) % mod) + (e * pow(p, e + 1, mod) % mod) temp2 %= mod temp2 *= modinv((p - 1)**2, mod) * p % mod temp2 %= mod ans *= (temp1 - temp2) % mod ans %= mod print(ans)