UVa 11178 Morley's Theorem (向量旋转)

http://uva.onlinejudge.org/index.php?option=onlinejudge&page=show_problem&problem=2119

/*0.025s*/#include<cstdio>#include<cmath>struct P{double x, y;P(double x = 0.0, double y = 0.0): x(x), y(y) {}void read() {scanf("%lf%lf", &x, &y);}void output() {printf("%f %f", x, y);}};typedef P Vector;Vector operator + (const Vector &A, const Vector &B) {return Vector(A.x + B.x, A.y + B.y);}Vector operator - (const P &A, const P &B) {return Vector(A.x - B.x, A.y - B.y);}Vector operator * (const Vector &A, double p) {return Vector(A.x * p, A.y * p);}Vector Rotate(const Vector &A, double rad) {return Vector(A.x * cos(rad) - A.y * sin(rad), A.x * sin(rad) + A.y * cos(rad));}inline double Dot(const Vector &A, const Vector &B) {return A.x * B.x + A.y * B.y;}inline double Cross(const Vector &A, const Vector &B) {return A.x * B.y - A.y * B.x;}inline double Length(const Vector &A) {return hypot(A.x, A.y);}inline double Angle(const Vector &A, const Vector &B) {return acos(Dot(A, B) / Length(A) / Length(B));}inline P GetLineIntersection(const P &p1, const Vector &s1, const P &p2, const Vector &s2){return p1 + s1 * (Cross(s2, p1 - p2) / Cross(s1, s2));}P getP(P A, P B, P C){Vector v1 = Rotate(C - B, Angle(A - B, C - B) / 3);Vector v2 = Rotate(B - C, -Angle(A - C, B - C) / 3); /// 负数表示顺时针旋转return GetLineIntersection(B, v1, C, v2);}int main(){int T;P A, B, C;scanf("%d", &T);while (T--){A.read(), B.read(), C.read();getP(A, B, C).output(), putchar(32);getP(B, C, A).output(), putchar(32);getP(C, A, B).output(), putchar(10);}return 0;}