結果
問題 | No.2847 Birthday Attack |
ユーザー | 👑 rin204 |
提出日時 | 2024-08-23 22:39:45 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,144 bytes |
コンパイル時間 | 365 ms |
コンパイル使用メモリ | 82,092 KB |
実行使用メモリ | 77,348 KB |
最終ジャッジ日時 | 2024-08-23 22:40:14 |
合計ジャッジ時間 | 26,442 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 736 ms
61,104 KB |
testcase_01 | AC | 706 ms
60,920 KB |
testcase_02 | AC | 726 ms
60,300 KB |
testcase_03 | AC | 1,051 ms
76,540 KB |
testcase_04 | WA | - |
testcase_05 | WA | - |
testcase_06 | WA | - |
testcase_07 | WA | - |
testcase_08 | WA | - |
testcase_09 | WA | - |
testcase_10 | WA | - |
testcase_11 | WA | - |
testcase_12 | WA | - |
testcase_13 | WA | - |
testcase_14 | WA | - |
testcase_15 | WA | - |
testcase_16 | WA | - |
testcase_17 | WA | - |
ソースコード
from math import gcd x, y, MOD = map(int, input().split()) ans = 0 ox = x // 2 ex = x - ox oy = y // 2 ey = y - oy ans += ox * (ox - 1) // 2 * y ans += ex * (ex - 1) // 2 * y ans += oy * (oy - 1) // 2 * x ans += ey * (ey - 1) // 2 * x ans %= MOD for n in range(1, 5000000): for m in range(n + 1, 5000000): if (m - n) % 2 == 0 or gcd(m, n) != 1: continue a = m * m - n * n b = 2 * m * n c = m * m + n * n if b > max(x, y): break a *= 2 b *= 2 k = min(x // a, y // b) ans += x * y * k ans %= MOD ans -= b * x * (k * (k + 1) // 2) % MOD ans %= MOD ans -= a * y * (k * (k + 1) // 2) % MOD ans %= MOD ans += a * b * (k * (k + 1) * (2 * k + 1) // 6) % MOD ans %= MOD k = min(x // b, y // a) ans += x * y * k ans %= MOD ans -= a * x * (k * (k + 1) // 2) % MOD ans %= MOD ans -= b * y * (k * (k + 1) // 2) % MOD ans %= MOD ans += a * b * (k * (k + 1) * (2 * k + 1) // 6) % MOD ans %= MOD ans *= 2 print(ans % MOD)