def norm(a, b, c):
    return (a**2+b**2+c**2)**0.5


eps = 10 ** -9
def solve(a, b, c):#a,b,c <= 1/2
    ans = []

    
    d = min(a, b, c)
    dd = d / (1-d)
    ans.append((1, 1, 1, dd * norm(1-a,1-b,1-c)))
    a,b,c = a - (1-a) * dd,b - (1-b) * dd,c - (1-c) * dd    
    
    dd = a / (1 - a)
    ans.append((1, 0, 0, dd * norm(1-a, b, c)))
    a, b, c = a - (1 - a) * dd, b + b * dd, c + c * dd
    #print(a, b, c)
    
    dd = b / (1 - b)
    ans.append((0, 1, 0, dd * norm(a, 1-b, c)))
    a, b, c = a + a * dd, b - (1 - b) * dd, c + c * dd
    
    
    dd = c / (1 - c)
    ans.append((0, 0, 1, dd * norm(a, b, 1-c)))
    a, b, c = a + a * dd, b + b * dd, c - (1 - c) * dd
    
    assert max(abs(a), abs(b), abs(c)) <= 10 ** -8
    return ans[::-1]

def check(ans, A,B,C):
    x,y,z = 0, 0, 0
    for r,g,b, t in ans:
        v = norm(r-x,g-y,b-z)
        x = x + (r - x) * t / v
        y = y + (g - y) * t / v
        z = z + (b - z) * t / v
        # print(x, y, z)
    # print(max(abs(x-a), abs(b-y), abs(c-z)))
    return max(abs(x-A), abs(B-y), abs(C-z)) <= 10 ** -6

m = 255

for _ in range(int(input())):
    a,b,c = map(int, input().split())
    A,B,C = a/m, b/m, c/m
    # print(A, B, C)

    flag = [0,0,0]
    ans = []
    if m - a < a:
        flag[0] = 1
        a = m - a
        ans.append((*flag, 1))
    if m - b < b:
        flag[1] = 1
        b = m - b
        ans.append((*flag, 1))
    if m - c < c:
        flag[2] = 1
        c = m - c
        ans.append((*flag, 1))
    # print(flag)
    ans += [(r^flag[0], flag[1]^g, b^flag[2], t) for r,g,b,t in solve(a/m, b/m, c/m)]
    print(len(ans))
    for r,g,b, t in ans:
        # print(r,g,b,t)
        print(r * m, g * m, b * m, t * m)
    assert check(ans, A,B,C)