import math def f(N): # Nを4平方和で表す O(N^1.5) for a in range(N+1): if 4*a*a>N: break for b in range(a,N+1): if a*a+3*b*b>N: break for c in range(b,N+1): if a*a+b*b+2*c*c>N: break dd=N-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()) print("YES") X=pow(M,N,B) a=math.isqrt(X) X-=a*a b=math.isqrt(X) X-=b*b print(*f(X),a,b)