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

using namespace std;

typedef complex<double> Point;

bool eqv(double d1, double d2) {
    return abs(d1-d2) < 1e-5;
}

int main() {
    int N;
    double x1, x2, x3, y1, y2, y3;

    while(true) {
        cin >> N;
        if(!N) break;

        cin >> x1 >> y1 >> x2 >> y2 >> x3 >> y3;
        Point p1(x1, y1), p2(x2, y2), p3(x3, y3);

        double R = abs(abs(p3-p2) / sin(arg(p2-p1) - arg(p3-p1))) / 2;
        Point pc = (p1+p2) / 2.;
        Point ph1 = (pc-p1) * Point(0, 1);
        double L = sqrt(R*R - norm(pc-p1));
        Point center = pc + ph1 * L / abs(ph1);
        if(!eqv(abs(center-p1), abs(center-p2)) ||
           !eqv(abs(center-p1), abs(center-p3)) ||
           !eqv(abs(center-p3), abs(center-p2))) center = pc - ph1*L/abs(ph1);


        printf("%.3f %.3f %.3f\n", center.real(), center.imag(), R);
    }
    return 0;
}