結果

問題 No.2189 六平方和
ユーザー FromBooskaFromBooska
提出日時 2023-06-02 18:06:33
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 1,548 bytes
コンパイル時間 304 ms
コンパイル使用メモリ 87,288 KB
実行使用メモリ 409,676 KB
最終ジャッジ日時 2023-08-28 02:24:08
合計ジャッジ時間 6,354 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 68 ms
71,508 KB
testcase_01 AC 67 ms
71,472 KB
testcase_02 AC 69 ms
71,332 KB
testcase_03 AC 69 ms
71,280 KB
testcase_04 AC 73 ms
70,960 KB
testcase_05 AC 70 ms
71,176 KB
testcase_06 AC 68 ms
71,368 KB
testcase_07 AC 68 ms
71,208 KB
testcase_08 AC 68 ms
71,568 KB
testcase_09 AC 67 ms
71,496 KB
testcase_10 AC 68 ms
71,284 KB
testcase_11 AC 75 ms
79,336 KB
testcase_12 AC 68 ms
71,396 KB
testcase_13 AC 67 ms
71,096 KB
testcase_14 AC 69 ms
71,536 KB
testcase_15 AC 67 ms
71,428 KB
testcase_16 AC 68 ms
71,368 KB
testcase_17 AC 68 ms
71,468 KB
testcase_18 AC 77 ms
77,592 KB
testcase_19 AC 78 ms
77,668 KB
testcase_20 TLE -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

# 前半の作戦はこれまでと同じ
# まずM**N (mod B) = Rをpowで計算する
# xi**2 (mod B)は常に0 for xi=B, or 1 for xi=1にできる
# x1からR以下一番近いところにx1**2 (mod B)がなる値にx1を固定
# 同様にx2を決める
# testerの4平方定理defを後半4つに使う

import math
def f(L):
    # Lを4平方和で表す O(L^1.5)
    for a in range(L+1):
        if 4*a*a > L: 
            break
        for b in range(a, L+1):
            if a*a + 3*b*b > L: 
                break
            for c in range(b, L+1):
                if a*a + b*b + 2*c*c > L: 
                    break
                dd = L - a*a - b*b - c*c
                d = math.isqrt(dd)
                if d*d == dd: 
                    return (a, b, c, d)

N, M, B = map(int, input().split())

R = pow(M, N, B)
R_remainder = R
INF = 10**20
X = []
for i in range(2):
    #print('i', i, 'R_remainder', R_remainder, X)
    if R_remainder == 0:
        X.append(B)
        continue
    
    diff = INF
    visited = set()
    for j in range(1, B+1):
        calc = (j**2)%B
        if calc not in visited:
            visited.add(calc)
            if calc <= R_remainder and R_remainder-calc < diff:
                diff = R_remainder-calc
                num = j
        elif calc in visited:
            break
    X.append(num)
    R_remainder -= (num**2)%B

if R_remainder == 0:
    for i in range(4):
        X.append(B)
    print('YES')
    print(*X)
else: 
    X2 = f(R_remainder)
    X.extend(X2)
    print('YES')
    print(*X)
0