#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