# この作戦でどうだ # まず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, x3、ーーーと決めていく # どんどん差は小さくなり、差を埋められればYes # WA出た、かなり通ったけど。 # ラグランジュの四平方定理により、全ての自然数が高々四個の平方数の和で表される # つまりこの問題にNOはないのだと思う # 最初の2つのxを求めたらあとは全探索にしてはどうだろう 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: for x3 in range(0, R_remainder+1): for x4 in range(0, R_remainder+1): for x5 in range(0, R_remainder+1): for x6 in range(0, R_remainder+1): if (x3**2 + x4**2 + x5**2 + x6**2)%B == R_remainder: if x3 == 0: x3 = B if x4 == 0: x4 = B if x5 == 0: x5 = B if x6 == 0: x6 = B X.extend([x3, x4, x5, x6]) print('YES') print(*X) exit()