BZOJ3621: 我想那还真是令人高兴啊 复数运算

题意：给出一对三角形，已知其中一个是由另一个绕某一点旋转放缩而得到，求这个点的坐标。每个测试点有T(<=10)组数据。
一开始以为是计算几何，然而向量的旋转和放缩是两种操作，复数却可以一步完成。
对于复数a和b，a*b的几何意义为a与b长度相乘，极角相加。
那么本题就可以看成从旋转点出发，指向其中一个三角形顶点的三个复数向量乘以同一个复数后变成了从旋转点出发，指向另一个三角形顶点的三个复数向量。设出旋转点坐标和同乘的复数，枚举顶点的六种可能对应关系。列出三个方程，用前两个方程解出解，代到第三个方程中检验是否正确即可。
解方程需要用到复数除法。对于(a+bi)/(c+di)，分子分母同乘(c-di)可进行分母有理化，
整理得结果为((ac+bd)/(cc+dd))+((bc-ad)/(cc+dd))i。
然而复数除法非常掉精度，取eps=1e-5都会WA，要1e-4才能过。。。

#include<cstdio>#include<cmath>using namespace std;#define y1 __y__struct cpx{    double a,b;    cpx(const double &a=0,const double &b=0):a(a),b(b){}    cpx(FILE* f){fscanf(f,"%lf%lf",&a,&b);}    #define c ano.a    #define d ano.b    inline cpx operator + (const cpx &ano) const    {return cpx(a+c,b+d);}    inline cpx operator - (const cpx &ano) const    {return cpx(a-c,b-d);}    inline cpx operator * (const cpx &ano) const    {return cpx(a*c-b*d,a*d+b*c);}    inline cpx operator / (const cpx &ano) const    {        double temp=c*c+d*d;        return cpx((a*c+b*d)/temp,(b*c-a*d)/temp);    }    inline bool operator == (const cpx &ano) const    {return fabs(a-c)<1e-4&&fabs(b-d)<1e-4;}    inline void operator() (FILE* f){fprintf(f,"%lf %lf\n",a,b);}}ans,a1,a2,a3;int t;inline bool check(const cpx &b1,const cpx &b2,const cpx &b3){    cpx __t=(b1-b2)/(a1-a2);    ans=(a1*__t-b1)/(__t-cpx(1.0,0.0));    return (a3-ans)*__t==(b3-ans);}int main(){    scanf("%d",&t);    while(t--)    {        a1=stdin,a2=stdin,a3=stdin;        cpx b1(stdin),b2(stdin),b3(stdin);        if(check(b1,b2,b3)) goto successed_to_match;        if(check(b1,b3,b2)) goto successed_to_match;        if(check(b2,b1,b3)) goto successed_to_match;        if(check(b2,b3,b1)) goto successed_to_match;        if(check(b3,b1,b2)) goto successed_to_match;        if(check(b3,b2,b1)) goto successed_to_match;        successed_to_match:        ans(stdout);    }    return 0;}