#include <iostream>
#include <complex>
#include <cmath>
#include <cstdio>

using namespace std;

const double eps = 1e-12;

int main() {
    int n;
    cin >> n;
    for (int i = 0; i < n; i++) {
        complex<double> p1, p2, p3, m1, m2, d1, d2;
        cin >> p1.real() >> p1.imag() >> p2.real() >> p2.imag() >> p3.real() >> p3.imag();
        m1 = (p1 + p2) / 2.0;
        m2 = (p2 + p3) / 2.0;
        d1 = polar(1.0, arg(p2 - p1) + M_PI_2);
        d2 = polar(1.0, arg(p3 - p2) + M_PI_2);
        
        double t;
        t = (d1.real() * m1.imag() + d1.imag() * m2.real() - m1.real() * d1.imag() - d1.real() * m2.imag()) / (d1.real() * d2.imag() - d2.real() * d1.imag());
        
        double x = m2.real() + t * d2.real();
        double y = m2.imag() + t * d2.imag();
        double r = abs(complex<double>(x,y) - p1);
        
        printf("%.3lf %.3lf %.3lf\n", x, y, r);
    }
    return 0;
}