Morley’s Theorem - UVa 11178 求直线交点

Output

Foreachlineofinputyoushouldproduceonelineofoutput.Thislinecontainssixfloatingpointnumbers

separatedbyasinglespace.Thesesixfloating-point
actuallymeansthattheCartesiancoordinatesofD,EandFare

respectively.Errorslessthan

will
beaccepted.

SampleInputOutputforSampleInput

2

112212

0010005050

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
usingnamespacestd;
structPoint
{
doublex,y;
Point(doublex=0,doubley=0):x(x),y(y){}
};
typedefPointVector;
Vectoroperator+(VectorA,VectorB){returnVector(A.x+B.x,A.y+B.y);}
Vectoroperator-(VectorA,VectorB){returnVector(A.x-B.x,A.y-B.y);}
Vectoroperator*(VectorA,doublep){returnVector(A.x*p,A.y*p);}
Vectoroperator/(VectorA,doublep){returnVector(A.x/p,A.y/p);}
doubleDot(VectorA,VectorB){returnA.x*B.x+A.y*B.y;}
doubleLength(VectorA){returnsqrt(Dot(A,A));}
doubleAngle(VectorA,VectorB){returnacos(Dot(A,B)/Length(A)/Length(B));}
doubleCross(VectorA,VectorB){returnA.x*B.y-A.y*B.x;}
doubleArea(PointA,PointB,PointC){returnCross(B-A,C-A);}
VectorRotate(VectorA,doublerad)
{returnVector(A.x*cos(rad)-A.y*sin(rad),A.x*sin(rad)+A.y*cos(rad));}
PointGetLineInterSection(PointP,Vectorv,PointQ,Vectorw)
{
Vectoru=P-Q;
doublet=Cross(w,u)/Cross(v,w);
returnP+v*t;
}
Pointsolve(PointA,PointB,PointC)
{
Vectorv1=C-B;
doublea1=Angle(A-B,v1);
v1=Rotate(v1,a1/3);

Vectorv2=B-C;
doublea2=Angle(A-C,v2);
v2=Rotate(v2,-a2/3);
returnGetLineInterSection(B,v1,C,v2);
}
intmain()
{
intT,t,i,j,k;
PointA,B,C,D,E,F;
scanf("%d",&T);
for(t=1;t<=T;t++)
{
scanf("%lf%lf%lf%lf%lf%lf",&A.x,&A.y,&B.x,&B.y,&C.x,&C.y);
D=solve(A,B,C);
E=solve(B,C,A);
F=solve(C,A,B);
printf("%lf%lf%lf%lf%lf%lf\n",D.x,D.y,E.x,E.y,F.x,F.y);
}
}