uva 11178 Morley's Theorem 三角形内角三等分线交点

#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
using namespace std;

#define all(x) (x).begin(), (x).end()
#define for0(a, n) for (int (a) = 0; (a) < (n); (a)++)
#define for1(a, n) for (int (a) = 1; (a) <= (n); (a)++)
typedef long long ll;
typedef pair<int, int> pii;
const int INF =0x3f3f3f3f;
//const int maxn=    ;
const double eps=1e-10;
struct Point
{
double x,y;
Point(double x=0,double y=0):x(x),y(y) {};
void show()
{
cout<<x<<" "<<y<<endl;
}
}A,B,C,D,E,F;

typedef Point Vector;
Vector operator +(Vector A,Vector B) {return Vector(A.x+B.x,A.y+B.y); }
Vector operator -(Vector A,Vector B) {return Vector(A.x-B.x,A.y-B.y); }
Vector operator *(Vector A,double p) {return Vector(A.x*p,A.y*p); }
Vector operator /(Vector A,double p) {return Vector(A.x/p,A.y/p); }
Vector operator -(Vector A)  {return  Vector(-A.x,-A.y);}
double Dot(Vector A,Vector B)//点乘
{
return A.x*B.x+A.y*B.y;
}

double Length(Vector A)
{
return sqrt(Dot(A,A));
}

double Cross(Vector A,Vector B)//叉乘
{
return A.x*B.y-A.y*B.x;
}

double angle(Vector A,Vector B)
{
return acos(Dot(A,B)/Length(A)/Length(B));
}
Vector Rotate(Vector A,double rad)//向量A 逆时针旋转rad度
{
return Vector(A.x*cos(rad)-A.y*sin(rad),A.x*sin(rad)+A.y*cos(rad)     );
}
Point GetLineIntersection(Point P,Vector v,Point Q,Vector w)
{
Vector u=P-Q;
double t=Cross(w,u)/Cross(v,w);
return P+v*t;

}

int main()
{
int T;scanf("%d",&T);
while(T--)
{
scanf("%lf%lf%lf%lf%lf%lf",&A.x,&A.y,&B.x,&B.y,&C.x,&C.y);
Vector AB=B-A,AC=C-A,BC=C-B;
Vector BA=-AB,CA=-AC,CB=-BC ;
double aA=angle(AB,AC);
double aB=angle(BA,BC);
double aC=angle(CA,CB);
Vector BF=Rotate(BA,-aB/3 );//旋转
Vector AF=Rotate(AB,aA/3  );
F=GetLineIntersection(B,BF,A,AF);
Vector CE=Rotate(CA,aC/3);
Vector AE=Rotate(AC,-aA/3);
E=GetLineIntersection(C,CE,A,AE);
Vector BD=Rotate(BC,aB/3);
Vector CD=Rotate(CB,-aC/3);
D=GetLineIntersection(B,BD,C,CD);
printf("%.6f %.6f %.6f %.6f %.6f %.6f\n",D.x,D.y,E.x,E.y,F.x,F.y);
}

return 0;
}