結果

問題 No.2847 Birthday Attack
ユーザー 👑 seekworserseekworser
提出日時 2024-07-09 23:12:29
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 548 ms / 3,000 ms
コード長 870 bytes
コンパイル時間 261 ms
コンパイル使用メモリ 82,148 KB
実行使用メモリ 76,672 KB
最終ジャッジ日時 2024-07-09 23:12:42
合計ジャッジ時間 6,752 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 42 ms
52,096 KB
testcase_01 AC 41 ms
52,224 KB
testcase_02 AC 42 ms
52,096 KB
testcase_03 AC 173 ms
75,904 KB
testcase_04 AC 41 ms
52,480 KB
testcase_05 AC 548 ms
76,180 KB
testcase_06 AC 305 ms
76,032 KB
testcase_07 AC 507 ms
76,672 KB
testcase_08 AC 197 ms
76,416 KB
testcase_09 AC 347 ms
76,288 KB
testcase_10 AC 268 ms
76,416 KB
testcase_11 AC 304 ms
76,460 KB
testcase_12 AC 473 ms
76,416 KB
testcase_13 AC 537 ms
76,672 KB
testcase_14 AC 411 ms
76,160 KB
testcase_15 AC 423 ms
76,416 KB
testcase_16 AC 339 ms
76,288 KB
testcase_17 AC 442 ms
76,672 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from math import gcd


x,y,md = map(int, input().split())
assert 1 <= x <= 4_000_000 and 1 <= y <= 4_000_000 and 100_000_000 <= md <= 1_000_000_000

ans = 0
mx = max(x, y) // 2 + 1
for m in range(0, mx+1):
    for n in range(m+1, mx+1):
        if (m % 2 == n % 2) or gcd(m, n) != 1: continue
        ai = n*n - m*m
        bi = 2*m*n
        if (ai > mx) or (bi > mx): break
        mul = 1
        while (mul * max(ai, bi) <= mx):
            a = ai*mul
            b = bi*mul
            if a*2 < x and b*2 < y:
                if (b == 0): ans = (ans + (x - 2*a) * (y - 2*b) * 2) % md
                else: ans = (ans + (x - 2*a) * (y - 2*b) * 4) % md
            if b*2 < x and a*2 < y:
                if (b == 0): ans = (ans + (x - 2*b) * (y - 2*a) * 2) % md
                else: ans = (ans + (x - 2*b) * (y - 2*a) * 4) % md
            mul += 1
print(ans % md)
0