BZOJ 3621: 我想那还真是令人高兴啊 计算几何 复数

https://www.lydsy.com/JudgeOnline/problem.php?id=3621

给定两个三角形，其中一个可以通过以某点为中心旋转并放缩的方式得到另一个，求这个中心

https://blog.csdn.net/PoPoQQQ/article/details/44346555

↑题解见这个博客。

get到了复数的特殊用法，原来旋转缩放都可以用复数实现的。因为复数相乘是极角相加、长度相乘，所以把点放到复平面上之后可以用一个复数点T代表缩放和旋转的系数，一个复数点P为中心点，列方程验证就行了。

我太智障了，最开始写的代码的验证方式特别复杂精度还有问题（我为什么要算出来三个中心点然后把三个中心点都比较一遍啊有毒吧），看了dalao的代码发现自己真实智障完全不用比较那么多，算出来两个比一下就可以了。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<complex>
using namespace std;
#define LL long long
const double eps=0.0001;
typedef complex< double >cd;
cd a[2][3];
inline double mabs(double x){return x<0?(-x):x;}
int main(){
    //freopen("a.in","r",stdin);
    int T;scanf("%d",&T);
    while(T-->0){
        double x,y;
        for(int i=0;i<2;++i){
            for(int j=0;j<3;j++){
                scanf("%lf%lf",&x,&y);
                a[i][j]=cd(x,y);
            }
        }
        cd ans=0;bool f=0;
        for(int i=0;i<3;i++){
            for(int j=0;j<3;j++){
                if(i==j)continue;
                cd z=(a[1][i]-a[1][j])/(a[0][0]-a[0][1]);
                int w=3-i-j;
                cd po=(a[0][0]*z-a[1][i])/(z-1.0);
                cd zz=(a[0][2]-po)*z-(a[1][w]-po);
                if(mabs(zz.real())<eps&&mabs(zz.imag())<eps){ans=po;f=1;break;}
            }
            if(f)break;
        }
        printf("%.4f %.4f\n",ans.real(),ans.imag());
    }
    return 0;
}