ll@n,@m,@b,x[5],y;
m=powmod(m,n,b);
rrep(y,b){
	ll m2=(m-y*y)%%b;
	rep(i,5){
		x[i]=Isqrt(m2);
		m2-=x[i]*x[i];
	}
	if(m2==0){
		break;
	}
}
wt("YES");
wt(x(5),y);