poj 1673 EXOCENTER OF A TRIANGLE

垂心。

/*
 * Author: stormdpzh
 * Created Time: 2012/7/19 15:12:17
 * File Name: poj_1673.cpp
 */
#include <iostream>
#include <cstdio>
#include <sstream>
#include <cstring>
#include <string>
#include <cmath>
#include <vector>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <algorithm>
#include <functional>

#define sz(v) ((int)(v).size())
#define rep(i, n) for(int i = 0; i < n; i++)
#define repf(i, a, b) for(int i = a; i <= b; i++)
#define repd(i, a, b) for(int i = a; i >= b; i--)
#define out(n) printf("%d\n", n)
#define mset(a, b) memset(a, b, sizeof(a))
#define wh(n) while(1 == scanf("%d", &n))
#define whz(n) while(1 == scanf("%d", &n) && n != 0)
#define lint long long

using namespace std;

const int MaxN = 1 << 30;

struct Point 
{
    double x, y;
    Point() {}
    Point(double _x, double _y) : x(_x), y(_y) {}
};

Point get_intersection(const Point &p1, const Point &p2, const Point &p3, const Point &p4)
{
    Point tmp = p1;
    double t = ((p1.x - p3.x) * (p3.y - p4.y) - (p1.y - p3.y) * (p3.x - p4.x)) / 
                ((p1.x - p2.x) * (p3.y - p4.y) - (p1.y - p2.y) * (p3.x - p4.x));
    tmp.x += (p2.x - p1.x) * t;
    tmp.y += (p2.y - p1.y) * t;
    return tmp;
}

double x[5], y[5];
Point p1, p2;

void gao()
{
    double xa = x[0] - y[2] + y[1];
    double ya = y[0] + x[2] - x[1];
    double xb = x[1] - y[2] + y[0];
    double yb = y[1] + x[2] - x[0];
    Point res = get_intersection(Point(x[0], y[0]), Point(xa, ya), Point(x[1], y[1]), Point(xb, yb));
    printf("%.4lf %.4lf\n", res.x, res.y);
}

int main()
{
    int t;
    scanf("%d", &t);
    while(t--) {
        rep(i, 3) scanf("%lf%lf", &x[i], &y[i]);
        gao();
    }
    return 0;
}