結果

問題 No.2189 六平方和
ユーザー FromBooskaFromBooska
提出日時 2023-06-02 18:06:33
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 1,548 bytes
コンパイル時間 317 ms
コンパイル使用メモリ 81,920 KB
実行使用メモリ 516,432 KB
最終ジャッジ日時 2024-12-28 15:43:00
合計ジャッジ時間 11,953 ms
ジャッジサーバーID
(参考情報)
judge1 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 25 TLE * 2 MLE * 1
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#
# M**N (mod B) = Rpow
# xi**2 (mod B)0 for xi=B, or 1 for xi=1
# x1Rx1**2 (mod B)x1
# x2
# tester4def4使
import math
def f(L):
# L4 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)
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0