from bisect import * N, M, B = map(int, input().split()) i = 0 S = set() A = pow(M, N, B) D = dict() while i * i < 10 ** 9: S.add(i * i) D[i * i] = i i += 1 L = sorted(list(S)) ind = bisect_right(L, A) ans = [D[L[ind-1]]] A -= L[ind-1] ind = bisect_right(L, A) ans.append(D[L[ind-1]]) A -= L[ind-1] LL = [] for a in L: if a <= A: LL.append(a) N = len(LL) for a in range(N): for b in range(a, N): for c in range(b, N): for d in range(c, N): if LL[a] + LL[b] + LL[c] + LL[d] == A: ans.append(D[LL[a]]) ans.append(D[LL[b]]) ans.append(D[LL[c]]) ans.append(D[LL[d]]) print("YES") print(*ans)