from collections import defaultdict, deque, Counter
import copy
from itertools import combinations, permutations, product, accumulate, groupby, chain
from heapq import heapify, heappop, heappush
import math
import bisect
from pprint import pprint
from random import randint
import sys
# sys.setrecursionlimit(700000)
input = lambda: sys.stdin.readline().rstrip('\n')
inf = float('inf')
mod1 = 10**9+7
mod2 = 998244353
def ceil_div(x, y): return -(-x//y)

#################################################

Q = int(input())
for _ in range(Q):
    d, k = map(int, input().split())
    if d == 0:
        if k == 1:
            print("Yes")
            print(0, 0)
        else:
            print("No")
    else:
        if k <= 4*d:
            print("Yes")
            if d%2 == 0:
                x = ceil_div(k, 12)
                y = (k-1)%12+1
                if y <= 4:
                    print(d//2+2*(x-1), d//2-2*(x-1))
                else:
                    print(d//2+(2*(x-1)+1), d//2-(2*(x-1)+1))
            else:
                if k <= d*4-4:
                    x = ceil_div(k, 8)
                    print(ceil_div(d, 2)+(x-1), ceil_div(d, 2)-(x-1))
                else:
                    print(d, 0)
        else:
            print("No")