100837D

囤了一个星期。。今天看了下vj上 sysuteam7 三年半之前的代码。。

深刻地认识到了自己智商不足的问题。

先求出来每个点对中心的偏移量。
确实是乱序的，但是我们可以极角排序，这样一定是一个
循环移位匹配到另一个循环移位。
所以枚举n次就可以，check也十分简单

#include <bits/stdc++.h>
using namespace std;
typedef double db;
const db eps = 1e-6;
const db pi = acos(-1);
int sign(db k){
    if(k>eps)return 1;else if(k<-eps)return -1;return 0;
}
int cmp(db k1,db k2){ return sign(k1-k2);}
struct point{
    db x,y;
    point operator+(const point &k1)const { return point{x+k1.x,y+k1.y};}
    point operator-(const point &k1)const { return point{x-k1.x,y-k1.y};}
    point operator *(db k1)const{ return point{x*k1,y*k1};}
    point operator/(db k1)const { return point{x/k1,y/k1};}
    int operator == (const point &k1) const{return cmp(x,k1.x)==0&&cmp(y,k1.y)==0;}
    bool operator<(const point k1)const {
        int a = cmp(x,k1.x);
        if(a==-1)return 1;else if(a==1) return 0;
        else return cmp(y,k1.y)==-1;
    }
    db abs(){ return sqrt(x*x+y*y);}
    db abs2(){return x*x+y*y;}
    db dis(point k1){ return ((*this)-k1).abs();}
    point unit(){db w = abs(); return point{x/w,y/w};}
    point turn90(){return point{-y,x};}
    db getP()const{return sign(y)==-1||(sign(y)==0&&sign(x)==-1);}
};
db cross(point k1,point k2){return k1.x*k2.y-k1.y*k2.x;}
db dot(point k1,point k2){return k1.x*k2.x+k1.y*k2.y;}
int compareangle (point k1,point k2){//[0,2pi]
    return k1.getP()<k2.getP()||(k1.getP()==k2.getP()&&sign(cross(k1,k2))>0)||(k1.getP()==k2.getP()&&sign(cross(k1,k2))==0&&abs(k1.x)<abs(k2.x));
}
int n;
point a[1005],b[1005];
point o1,o2;
int main(){
    //freopen("d.in","r",stdin);
    //freopen("d.out","w",stdout);
    scanf("%d",&n);
    if(n==1) return 0*printf("%.11f\n",0.0);
    for(int i=0;i<n;i++){
        scanf("%lf%lf",&a[i].x,&a[i].y);
        o1=o1+a[i];
    }
    for(int i=0;i<n;i++){
        scanf("%lf%lf",&b[i].x,&b[i].y);
        o2=o2+b[i];
    }
    o1=o1/n,o2=o2/n;
    for(int i=0;i<n;i++){
        a[i]=a[i]-o1,b[i]=b[i]-o2;
    }
    sort(a,a+n,compareangle);
    sort(b,b+n,compareangle);
//    printf("%.11f\n",a[0].abs());
//    printf("%.11f\n",b[2].abs());
    db ans = pi;
    for(int i=0;i<n;i++){//[0,n)~[i,i+n)
        bool f=1;
        for(int j=0;j<n;j++){
            if(cmp(a[j].abs(),b[(i+j)%n].abs())!=0){
                f=0;break;
            }
        }
        if(f){
            ans=min(ans,acos(dot(a[0],b[i])/a[0].abs()/b[i].abs()));
        }
    }
    printf("%.11f\n",ans);
}
/**
 先求出来每个点对中心的偏移量。
 确实是乱序的，但是我们可以极角排序，这样一定是一个
 循环移位匹配到另一个循环移位。
 所以枚举n次就可以，check也十分简单
 */