#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;

int gcd(int u, int v) {
	if(u==v) return u;
	if(u==0) return v;
	if(v==0) return u;
	if(~u & 1) {
		if(v & 1) return gcd(u>>1, v);
		else return gcd(u>>1, v>>1)<<1;
	}
	if(~v & 1) return gcd(u, v>>1);
	if(u>v) return gcd((u-v)>>1, v);
	return gcd((v-u)>>1, u);
}

int main() {
	int a, b, x, y, g, Bx, By, B;
	scanf("%d%d", &a, &b);
	for(int i=0; i<a; i++) {
		for(int j=0; j<a; j++) {
			if(pow(i, 2) + pow(j, 2)==pow(a, 2)) {
				x=i, y=j;
				g=gcd(x, y);
				Bx=-y/g;
				By=x/g;
				B=sqrt(pow(Bx, 2) + pow(By, 2));
				if(b%B==0 && sqrt(pow(By*b/B-y, 2))!=0 && sqrt(pow(Bx*b/B-x, 2))!=0) {
					printf("YES\n");
					printf("0 0\n");
					printf("%d %d\n", x, y);
					printf("%d %d\n", Bx*b/B, By*b/B);
					return 0;
				}
			}
		}
	}
	printf("NO");
	return 0;
}