#include<bits/stdc++.h>
using namespace std;
#define db double
const int O=255;
const db eps=1e-9;
bool ap(db x,db y){return abs(x-y)<eps;}
db p2(db x){return x*x;}
int ans;
int X[12],Y[12],Z[12];
db D[12];
void xp(int x,int y,int z,db d){
    ans++,X[ans]=x,Y[ans]=y,Z[ans]=z,D[ans]=d;
}
void sol(){
    ans=0;
    db x,y,z;
    scanf("%lf%lf%lf",&x,&y,&z);
    x/=O,y/=O,z/=O;
    if(ap(x,1)&&ap(y,1)&&ap(z,1))xp(1,1,1,1000);
    else{
        db D=sqrt(p2(1-x)+p2(1-y)+p2(1-z));
        db te=1e9+7;
        if(!ap(x,1))te=min(te,x/(1-x));if(!ap(y,1))te=min(te,y/(1-y));if(!ap(z,1))te=min(te,z/(1-z));
        xp(1,1,1,D*te);
        x-=te*(1-x),y-=te*(1-y),z-=te*(1-z);
        int w=-1;
        if(ap(x,0))w=0;
        else if(ap(y,0))w=1,swap(x,y);
        else w=2,swap(x,z);
        if(ap(y,1)&&ap(z,1))xp(0,1,1,1000);
        else{
            D=sqrt(p2(1-y)+p2(1-z)),te=1e9+7;
            if(!ap(y,1))te=min(te,y/(1-y));if(!ap(z,1))te=min(te,z/(1-z));
            xp(0,1,1,D*te),y-=te*(1-y),z-=te*(1-z);
            if(ap(y,0))xp(0,0,1,z);
            else xp(0,1,0,y);
        }
        for(int i=2;i<=ans;i++){
            if(w==1)swap(X[i],Y[i]);
            if(w==2)swap(X[i],Z[i]);
        }
    }
    printf("%d\n",ans);
    for(int i=ans;i;i--)printf("%d %d %d %.10lf\n",X[i]*O,Y[i]*O,Z[i]*O,D[i]*O);
}
int main(){
    int T;scanf("%d",&T);while(T--)sol();
    return 0;
}