#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;
}
				I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y3N0ZGlvPgojaW5jbHVkZSA8Y21hdGg+CnVzaW5nIG5hbWVzcGFjZSBzdGQ7CgppbnQgZ2NkKGludCB1LCBpbnQgdikgewoJaWYodT09dikgcmV0dXJuIHU7CglpZih1PT0wKSByZXR1cm4gdjsKCWlmKHY9PTApIHJldHVybiB1OwoJaWYofnUgJiAxKSB7CgkJaWYodiAmIDEpIHJldHVybiBnY2QodT4+MSwgdik7CgkJZWxzZSByZXR1cm4gZ2NkKHU+PjEsIHY+PjEpPDwxOwoJfQoJaWYofnYgJiAxKSByZXR1cm4gZ2NkKHUsIHY+PjEpOwoJaWYodT52KSByZXR1cm4gZ2NkKCh1LXYpPj4xLCB2KTsKCXJldHVybiBnY2QoKHYtdSk+PjEsIHUpOwp9CgppbnQgbWFpbigpIHsKCWludCBhLCBiLCB4LCB5LCBnLCBCeCwgQnksIEI7CglzY2FuZigiJWQlZCIsICZhLCAmYik7Cglmb3IoaW50IGk9MDsgaTxhOyBpKyspIHsKCQlmb3IoaW50IGo9MDsgajxhOyBqKyspIHsKCQkJaWYocG93KGksIDIpICsgcG93KGosIDIpPT1wb3coYSwgMikpIHsKCQkJCXg9aSwgeT1qOwoJCQkJZz1nY2QoeCwgeSk7CgkJCQlCeD0teS9nOwoJCQkJQnk9eC9nOwoJCQkJQj1zcXJ0KHBvdyhCeCwgMikgKyBwb3coQnksIDIpKTsKCQkJCWlmKGIlQj09MCAmJiBzcXJ0KHBvdyhCeSpiL0IteSwgMikpIT0wICYmIHNxcnQocG93KEJ4KmIvQi14LCAyKSkhPTApIHsKCQkJCQlwcmludGYoIllFU1xuIik7CgkJCQkJcHJpbnRmKCIwIDBcbiIpOwoJCQkJCXByaW50ZigiJWQgJWRcbiIsIHgsIHkpOwoJCQkJCXByaW50ZigiJWQgJWRcbiIsIEJ4KmIvQiwgQnkqYi9CKTsKCQkJCQlyZXR1cm4gMDsKCQkJCX0KCQkJfQoJCX0KCX0KCXByaW50ZigiTk8iKTsKCXJldHVybiAwOwp9