【BZOJ】3621: 我想那还真是令人高兴啊

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

题意：给两个三角形，问A能否通过旋转伸缩到B。

#include <bits/stdc++.h>
using namespace std;
struct cp {
	static const double eps=1e-4;
	double x, y;
	cp(double _x=0, double _y=0) : x(_x), y(_y) {}
	void scan() { scanf("%lf%lf", &x, &y); }
	void P() const { printf("%.9f %.9f
", x, y); }
	cp operator - (const cp &a) { return cp(x-a.x, y-a.y); }
	cp operator * (const cp &a) { return cp(x*a.x-y*a.y, x*a.y+y*a.x); }
	cp operator / (const cp &a) { static double c; c=a.x*a.x+a.y*a.y; return cp((x*a.x+y*a.y)/c, (y*a.x-x*a.y)/c); }
	bool operator == (const cp &a) { return abs(x-a.x)<eps && abs(y-a.y)<eps; }
}a[2][3], ans;
int vis[3], b[3];
bool check() {
	static cp T, P;
	T=(a[1][b[0]]-a[1][b[1]])/(a[0][0]-a[0][1]);
	P=(a[1][b[0]]-a[0][0]*T)/(cp(1)-T);
	if((a[0][2]-P)*T==a[1][b[2]]-P) {
		ans=P;
		return 1;
	}
	return 0;
}
bool dfs(int x) {
	if(x==3) return check();
	for(int i=0; i<3; ++i) if(!vis[i]) {
		vis[i]=1; b[x]=i;
		if(dfs(x+1)) return 1;
		vis[i]=0;
	}
	return 0;
}
int main() {
	int T;
	scanf("%d", &T);
	while(T--) {
		for(int i=0; i<2; ++i)
			for(int j=0; j<3; ++j)
				a[i][j].scan();
		memset(vis, 0, sizeof vis);
		if(!dfs(0)) puts("-1");
		else ans.P();
	}
	return 0;
}

　　

复数乘法的几何意义就是极角相加，长度相乘。因此当做复平面来做就行辣= =