結果
問題 | No.2847 Birthday Attack |
ユーザー | shobonvip |
提出日時 | 2024-08-23 23:27:17 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 477 ms / 3,000 ms |
コード長 | 771 bytes |
コンパイル時間 | 443 ms |
コンパイル使用メモリ | 82,360 KB |
実行使用メモリ | 76,876 KB |
最終ジャッジ日時 | 2024-08-23 23:27:24 |
合計ジャッジ時間 | 6,677 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 68 ms
60,320 KB |
testcase_01 | AC | 66 ms
60,148 KB |
testcase_02 | AC | 66 ms
59,780 KB |
testcase_03 | AC | 235 ms
76,036 KB |
testcase_04 | AC | 66 ms
59,800 KB |
testcase_05 | AC | 477 ms
76,772 KB |
testcase_06 | AC | 329 ms
76,564 KB |
testcase_07 | AC | 458 ms
76,596 KB |
testcase_08 | AC | 262 ms
76,832 KB |
testcase_09 | AC | 338 ms
76,812 KB |
testcase_10 | AC | 287 ms
76,804 KB |
testcase_11 | AC | 283 ms
76,500 KB |
testcase_12 | AC | 415 ms
76,772 KB |
testcase_13 | AC | 437 ms
76,704 KB |
testcase_14 | AC | 418 ms
76,876 KB |
testcase_15 | AC | 341 ms
76,464 KB |
testcase_16 | AC | 323 ms
76,688 KB |
testcase_17 | AC | 329 ms
76,692 KB |
ソースコード
import math x, y, mod = map(int,input().split()) ans = 0 for p in range(1, 10 ** 7 + 1): for q in range(0, p): if p ** 2 + q ** 2 > max(x, y): break if p % 2 == q % 2: continue if math.gcd(p, q) != 1: continue a = p ** 2 - q ** 2 b = 2 * p * q c = p ** 2 + q ** 2 k = 1 while True: if x - a * k * 2 <= 0: break if y - b * k * 2 <= 0: break ans += (x - 2 * a * k) * (y - 2 * b * k) % mod * ((a > 0) + 1) % mod * ((b > 0) + 1) % mod ans %= mod k += 1 a = 2 * p * q b = p ** 2 - q ** 2 c = p ** 2 + q ** 2 k = 1 while True: if x - a * k * 2 <= 0: break if y - b * k * 2 <= 0: break ans += (x - 2 * a * k) * (y - 2 * b * k) % mod * ((a > 0) + 1) % mod * ((b > 0) + 1) % mod ans %= mod k += 1 print(ans % mod)