[复数 解析几何 || 随机化] Ural 1046 Geometrical Dreams & POJ 2600 Geometrical dreams

这个可以先设出a1 然后 一圈递推出an+1 an+1==a1 解这个方程就好了

这里用复数来实现 因为复数乘法是天然的旋转 

复数旋转就是类似这样的东西

以下来自 POJ Discuss

point getnp(point ori, point cen, double angle)  //向量cen->ori绕cen旋转逆时针旋转angle弧度
{
	//复数平面向量旋转
	point v;
	v.x = ori.x - cen.x;
	v.y = ori.y - cen.y;
	point res;
	res.x = v.x*cos(angle) - v.y*sin(angle) + cen.x;
	res.y = v.x*sin(angle) + v.y*cos(angle) + cen.y;
	return res;
}

唯一坑的就是poj上精度 输出要保留零位

#include<cstdio>  
#include<cstdlib>
#include<algorithm>
#include<complex>
#include<cmath>
using namespace std;  
  
typedef complex<double> E;

const int N=55;
const double PI=acos(-1.0);

int n;  
E M[N],ans[N];  
double a[N];
  
int main(){
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  scanf("%d",&n);  
  double ix,iy;
  for (int i=1;i<=n;i++)
    scanf("%lf%lf",&ix,&iy),M[i]=E(ix,iy);
  for (int i=1;i<=n;i++)
    scanf("%lf",&a[i]),a[i]=a[i]/180.0*PI; 
  E x(1,0),y(0,0),tmp;
  for (int i=1;i<=n;i++){
    tmp=E(cos(a[i]),sin(a[i]));  
    x*=tmp;  
    y-=M[i],y*=tmp,y+=M[i];  
  }
  ans[1]=(-y)/(x-E(1,0));  
  x=E(1,0),y=E(0,0);  
  for (int i=1;i<=n;i++){
    tmp=E(cos(a[i]),sin(a[i]));  
    x*=tmp;  
    y-=M[i],y*=tmp,y+=M[i];  
    ans[i+1]=x*ans[1]+y;  
  }
  for (int i=1;i<=n;i++)
    printf("%.2lf %.2lf\n",ans[i].real(),ans[i].imag());  
  return 0;  
}  

论文：刘家骅《浅谈随机化在信息学竞赛中的应用》 描述了一种随机化的做法 类似于爬山

让我们从另外一个角度思考问题，从题目的条件易知，只要确定了一个顶点的坐标，多边形的其他顶点的坐标就能够通过简单计算得到，那么问题就转化为确定多边形的一个顶点的坐标。

如何确定一个顶点的坐标呢，枚举和二分等常用算法都无法为我们解决问题，于是，我们想到了随机化。

我们开始时将第一个点放在原点，通过计算能够得出第N+1个顶点的坐标，如果第N+1个顶点和第1个顶点重合，这个多边形就是所求，但是显然这样的可能性非常渺茫，于是我们需要调整第一点的位置。显然，第一个点距离第N+1个点的距离越小，其位置越接近其实际位置。我们每次可以在暂时确定第一个点的位置附近随机一个点，判断第一个点放在这里这个位置时与原来相比第一个点与第N+1个点的距离是否比原来更小，如果是，则将第一个点的位置暂时定在这个位置，然后继续上述操作，直至第一个点与第N+1个点重合，我们就确定了多边形的顶点的坐标了。