morley's theorem uva11178

1、问题描述
2、只用利用直线交点求解即可，注意三分线可以用Rotate函数更新

#include <iostream>#include<cmath>using namespace std;/* run this program using the console pauser or add your own getch, system("pause") or input loop */struct point{    double x,y;    point(double x=0,double y=0):x(x),y(y){}};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);}bool operator <(const point&a,const point&b){    return a.x<b.x||(a.x==b.x&&a.y<b.y);}const double eps=1e-10;int dcmp(double x){    if(fabs(x)<eps)return 0;    else return x<0?-1:1;}bool operator ==(const point&a,const point&b){    return dcmp(a.x-b.x)==0&&dcmp(a.y-b.y)==0;}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 Angle(Vector a,Vector b){    return acos(Dot(a,b)/Length(a)/Length(b));}double Cross(Vector a,Vector b){    return a.x*b.y-a.y-b.x;}double Area2(point a,point b,point c){    return Cross(b-a,c-a);}Vector Rotate(Vector a,double 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;}double distancetoline(point p,point a,point b){    Vector v1=b-a,v2=p-a;    return fabs(Cross(v1,v2))/Length(v1);}point getd(point a,point b,point c){    Vector v1=c-b;    double a1=Angle(a-b,v1);    v1=Rotate(v1,a1/3);    Vector v2=b-c;    double a2=Angle(a-c,v2);    v2=Rotate(v2,-a2/3);    return Getlineintersection(b,v1,c,v2);}point read_point(){    int x,y;    cin>>x>>y;    point a(x,y);    return a;}int main(int argc, char** argv) {    int t;    point a,b,c,d,e,f;    cin>>t;    while(t--){        a=read_point();        b=read_point();        c=read_point();        d=getd(a,b,c);        e=getd(b,c,a);        f=getd(c,a,b);        cout<<endl;    }    return 0;}