結果
問題 | No.886 Direct |
ユーザー | maspy |
提出日時 | 2020-03-28 10:32:19 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 297 ms / 4,000 ms |
コード長 | 1,241 bytes |
コンパイル時間 | 1,237 ms |
コンパイル使用メモリ | 87,048 KB |
実行使用メモリ | 81,720 KB |
最終ジャッジ日時 | 2023-08-30 17:29:30 |
合計ジャッジ時間 | 7,873 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge15 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 100 ms
72,296 KB |
testcase_01 | AC | 101 ms
72,384 KB |
testcase_02 | AC | 101 ms
72,604 KB |
testcase_03 | AC | 124 ms
77,496 KB |
testcase_04 | AC | 98 ms
72,352 KB |
testcase_05 | AC | 100 ms
72,432 KB |
testcase_06 | AC | 99 ms
72,348 KB |
testcase_07 | AC | 100 ms
72,168 KB |
testcase_08 | AC | 100 ms
72,584 KB |
testcase_09 | AC | 101 ms
72,420 KB |
testcase_10 | AC | 100 ms
72,224 KB |
testcase_11 | AC | 101 ms
72,420 KB |
testcase_12 | AC | 104 ms
72,560 KB |
testcase_13 | AC | 102 ms
72,448 KB |
testcase_14 | AC | 101 ms
72,528 KB |
testcase_15 | AC | 99 ms
72,412 KB |
testcase_16 | AC | 99 ms
72,528 KB |
testcase_17 | AC | 129 ms
78,664 KB |
testcase_18 | AC | 135 ms
78,312 KB |
testcase_19 | AC | 138 ms
78,060 KB |
testcase_20 | AC | 130 ms
78,288 KB |
testcase_21 | AC | 122 ms
77,732 KB |
testcase_22 | AC | 134 ms
78,580 KB |
testcase_23 | AC | 213 ms
79,112 KB |
testcase_24 | AC | 231 ms
79,620 KB |
testcase_25 | AC | 200 ms
79,272 KB |
testcase_26 | AC | 208 ms
79,140 KB |
testcase_27 | AC | 297 ms
81,080 KB |
testcase_28 | AC | 281 ms
80,016 KB |
testcase_29 | AC | 98 ms
72,368 KB |
testcase_30 | AC | 105 ms
72,356 KB |
testcase_31 | AC | 254 ms
80,876 KB |
testcase_32 | AC | 266 ms
81,616 KB |
testcase_33 | AC | 293 ms
81,720 KB |
testcase_34 | AC | 267 ms
81,440 KB |
testcase_35 | AC | 257 ms
80,812 KB |
ソースコード
#!/usr/bin/ python3.8 import sys read = sys.stdin.buffer.read readline = sys.stdin.buffer.readline readlines = sys.stdin.buffer.readlines from functools import lru_cache H, W = map(int, readline().split()) if H < W: H, W = W, H def S1(n): return n * (n + 1) // 2 def S2(n): return n * (n + 1) * (2 * n + 1) // 6 @lru_cache(None) def f(H, W): """ 1 <= x <= H, 1 <= y <= W, gcd(x,y) = 1 に対して、sum 1, sum x, sum y, sum xyをかえす""" SH = H * (H + 1) // 2 SW = W * (W + 1) // 2 sum_1 = H * W sum_x = SH * W sum_y = SW * H sum_xy = SH * SW d = 2 while True: Hd = H // d Wd = W // d if not Wd: break next_d = min(H // Hd, W // Wd) + 1 s1, sx, sy, sxy = f(Hd, Wd) sum_1 -= s1 * (next_d - d) sum_x -= sx * (S1(next_d - 1) - S1(d - 1)) sum_y -= sy * (S1(next_d - 1) - S1(d - 1)) sum_xy -= sxy * (S2(next_d - 1) - S2(d - 1)) d = next_d return sum_1, sum_x, sum_y, sum_xy sum_1, sum_x, sum_y, sum_xy = f(H, W) answer = sum_1 * H * W - W * sum_x - H * sum_y + sum_xy answer *= 2 # まよこ、まうえ answer += (W - 1) * H answer += (H - 1) * W MOD = 10 ** 9 + 7 print(answer % MOD)