题目传送门（内部题86）
输入格式
第一行$n$
接下来$n$行，每行四个浮点数，分别表示变换前的坐标和变换后的坐标
输出格式
第一行浮点数$\theta$以弧度制表示
第二行浮点数$scale$
第三行两个浮点数$d_x,d_y$
我们将用$SPJ$以$10^{-3}$的绝对误差来判断变换结果是否正确，建议输出$10$位小数以上。
样例
样例输入1：
5
0 0 -1 1
0 1 -2 1
1 0 -1 2
1 1 0 0
2 1 1 0
样例输出1：
1.5707963268
1
-1 1
样例输入2：
5
0 0 1 1
0 1 1 2
1 0 2 1
1 1 0 0
2 1 1 0
样例输出2：
0
1
1 1
数据范围与提示
对于$30\%$的数据，保证存在可行的答案$\theta=0,scale=1$
对于另$30\%$的数据，$n\leqslant 500$
对于所有数据$n\leqslant 100,000$，坐标属于$[-10^9,10^9]$
题解
正解是随机化额……
不过还是有正确性的。
如果我们知道两组$(x_i,y_i)$，那么就可以解出$\theta,scale,d_x,d_y$了，然后将解出来的这四个值带入检验即可。
因为至少有一半是正确的，所以每次可以解出的概率就是$\frac{1}{4}$，解不出来的概率就是$\frac{3}{4}$，那么解$k$次，解不出来的概率就是$\frac{1}{4}^k$，$k=50$时概率已经$<10^{-5}$了。
对于官方题解里说的高斯消元，因为就两组，没有必要。
实在不行多交几遍（其实我就脸黑了……）
时间复杂度：$\Theta($玄学$)$。
期望得分：$100$分。
实际得分：$0\sim 100$分（脸黑别怪我）。
代码时刻
#include<bits/stdc++.h>
using namespace std;
const double eps=1e-5;
struct rec{double x0,y0,x2,y2;}e[100001];
pair<double,double> d[100001];
int n;
double dx,dy,Cos,Sin,scale,theta;
void solve(int x,int y)
{
double d1=e[x].x0-e[y].x0;
double d2=e[x].y0-e[y].y0;
double d3=e[x].x2-e[y].x2;
double d4=e[x].y2-e[y].y2;
double COS=(d1*d3+d2*d4)/(d1*d1+d2*d2);
double SIN=(d4-d2*COS)/d1;
theta=atan(SIN/COS);
dx=e[x].x2-e[x].x0*COS+e[x].y0*SIN;
dy=e[x].y2-e[x].x0*SIN-e[x].y0*COS;
Cos=cos(theta);
Sin=sin(theta);
scale=SIN/Sin;
if(scale<0.0)
{
theta+=M_PI;
Cos=cos(theta);
Sin=sin(theta);
scale=SIN/Sin;
}
}
bool judge()
{
int sum=0;
for(int i=1;i<=n;i++)
{
double nowx=e[i].x0*Cos-e[i].y0*Sin;
double nowy=e[i].x0*Sin+e[i].y0*Cos;
nowx=nowx*scale+dx;nowy=nowy*scale+dy;
if(fabs(nowx-e[i].x2)<eps&&abs(nowy-e[i].y2)<eps)sum++;
}
return (2*sum>=n);
}
int main()
{
srand(time(NULL));
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
scanf("%lf%lf%lf%lf",&e[i].x0,&e[i].y0,&e[i].x2,&e[i].y2);
d[i]=make_pair(e[i].x2-e[i].x0,e[i].y2-e[i].y0);
}
sort(d+1,d+n+1);
int sum=0;
for(int i=1;i<=n;i++)
{
if(fabs(d[i].first-d[i-1].first)>eps||fabs(d[i].second-d[i-1].second)>eps)sum=0;
sum++;
if(2*sum>=n)
{
puts("0\n1");
printf("%.10lf %.10lf",d[i].first,d[i].second);
return 0;
}
}
while(1)
{
int x,y;
do
{
x=rand()%n+1;
y=rand()%n+1;
}while(x==y);
solve(x,y);
if(judge())
{
printf("%.10lf\n%.10lf\n%.10lf %.10lf",theta,scale,dx,dy);
return 0;
}
}
return 0;
}
rp++