結果

問題 No.2847 Birthday Attack
ユーザー 👑 seekworserseekworser
提出日時 2024-07-09 23:01:42
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 524 ms / 3,000 ms
コード長 742 bytes
コンパイル時間 1,228 ms
コンパイル使用メモリ 82,484 KB
実行使用メモリ 76,712 KB
最終ジャッジ日時 2024-07-09 23:01:51
合計ジャッジ時間 7,316 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 41 ms
54,180 KB
testcase_01 AC 40 ms
53,008 KB
testcase_02 AC 44 ms
53,052 KB
testcase_03 AC 165 ms
76,080 KB
testcase_04 AC 40 ms
52,572 KB
testcase_05 AC 505 ms
76,580 KB
testcase_06 AC 292 ms
76,360 KB
testcase_07 AC 475 ms
75,976 KB
testcase_08 AC 201 ms
76,108 KB
testcase_09 AC 333 ms
76,172 KB
testcase_10 AC 264 ms
75,992 KB
testcase_11 AC 290 ms
76,712 KB
testcase_12 AC 456 ms
76,120 KB
testcase_13 AC 524 ms
76,208 KB
testcase_14 AC 376 ms
76,388 KB
testcase_15 AC 409 ms
76,356 KB
testcase_16 AC 326 ms
76,184 KB
testcase_17 AC 423 ms
76,484 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from math import gcd


x,y,md = map(int, input().split())

def calc_total(x) -> int:
    xk = x // 2
    xtotal = (xk + 1) * xk
    return (xk * x - xtotal) % md
ans = (calc_total(x) * 2 * y + calc_total(y) * 2 * x) % md
mx = max(x, y) // 2 + 1
for m in range(1, 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: ans = (ans + (x - 2*a) * (y - 2*b) * 4) % md
            if b*2 < x and a*2 < y: ans = (ans + (x - 2*b) * (y - 2*a) * 4) % md
            mul += 1
print(ans % md)
0