結果
問題 | No.2847 Birthday Attack |
ユーザー | 👑 rin204 |
提出日時 | 2024-08-23 22:38:06 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,133 bytes |
コンパイル時間 | 327 ms |
コンパイル使用メモリ | 82,420 KB |
実行使用メモリ | 848,976 KB |
最終ジャッジ日時 | 2024-08-23 22:38:16 |
合計ジャッジ時間 | 9,053 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 591 ms
60,604 KB |
testcase_01 | AC | 569 ms
60,828 KB |
testcase_02 | AC | 571 ms
60,476 KB |
testcase_03 | AC | 1,346 ms
279,916 KB |
testcase_04 | WA | - |
testcase_05 | MLE | - |
testcase_06 | -- | - |
testcase_07 | -- | - |
testcase_08 | -- | - |
testcase_09 | -- | - |
testcase_10 | -- | - |
testcase_11 | -- | - |
testcase_12 | -- | - |
testcase_13 | -- | - |
testcase_14 | -- | - |
testcase_15 | -- | - |
testcase_16 | -- | - |
testcase_17 | -- | - |
ソースコード
from math import gcd x, y, MOD = map(int, input().split()) ans = 0 lst = [] for n in range(1, 4000000): for m in range(n + 1, 4000000): 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 lst.append((a, b, c)) 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 a, b, c in lst: 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)