#include <iostream>
#include <cmath>
#include <numbers>

using namespace std;

int sq(int a) {
    return a * a;
}

double sq(double a) {
    return a * a;
}


double getDistance(pair<double, double> p1, pair<double, double> p2) {
    auto [x1, y1] = p1;
    auto [x2, y2] = p2;
    return sqrt(sq(x1-x2) + sq(y1-y2));
}

bool lineOnCircle(int x1, int y1, int x2, int y2, int xc, int yc, int r) {
    int a = sq(x2-x1) + sq(y2-y1);
    int b = 2 * (((x2 - x1) * (x1 - xc)) + ((y2 - y1) * (y1 - yc)));
    int c = sq(x1-xc) + sq(y1-yc) - sq(r);
    int disc = sq(b) - 4 * a * c;
    return (disc >= 0);
}

const int BS_ITERS = 10000;

int main() {
    int T;
    cin >> T;
    while (T--) {
        int xa, ya, xb, yb, xc, yc, r;
        cin >> xa >> ya >> xb >> yb >> xc >> yc >> r;
        pair<int, int> a = {xa, ya};
        pair<int, int> b = {xb, yb};
        double ans;
        double theta = M_PI/2;
        pair<double, double> currPos = {xc, yc+r};
        for (int _ = 0; _ < BS_ITERS; _++) {
            auto [currX, currY] = currPos;

            double relX = currX-xc;
            double relY = currY-yc;

            pair<double, double> ccw = {xc + (relX*cos(theta)-relY*sin(theta)),
                                        yc + (relX*sin(theta)+relY*cos(theta))};

            pair<double, double> cw = {xc + (relX*cos(-theta)-relY*sin(-theta)),
                                        yc + (relX*sin(-theta)+relY*cos(-theta))};

            double ccwD = getDistance(a, ccw) + getDistance(b, ccw);
            double cwD = getDistance(a, cw) + getDistance(b, cw);

            if (ccwD < cwD) {
                currPos = ccw;
            } else {
                currPos = cw;
            }
            theta /= 2;
            ans = getDistance(a, currPos) + getDistance(b, currPos);
        }
        cout << ans << "\n";
    }


    return 0;
}