100837D
囤了一个星期。。今天看了下vj上 sysuteam7 三年半之前的代码。。
深刻地认识到了自己智商不足的问题。
先求出来每个点对中心的偏移量。
确实是乱序的,但是我们可以极角排序,这样一定是一个
循环移位匹配到另一个循环移位。
所以枚举n次就可以,check也十分简单
#include <bits/stdc++.h>
using namespace std;
typedef double db;
const db eps = 1e-;
const db pi = acos(-);
int sign(db k){
if(k>eps)return ;else if(k<-eps)return -;return ;
}
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)==&&cmp(y,k1.y)==;}
bool operator<(const point k1)const {
int a = cmp(x,k1.x);
if(a==-)return ;else if(a==) return ;
else return cmp(y,k1.y)==-;
}
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)==-||(sign(y)==&&sign(x)==-);}
};
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))>)||(k1.getP()==k2.getP()&&sign(cross(k1,k2))==&&abs(k1.x)<abs(k2.x));
}
int n;
point a[],b[];
point o1,o2;
int main(){
//freopen("d.in","r",stdin);
//freopen("d.out","w",stdout);
scanf("%d",&n);
if(n==) return *printf("%.11f\n",0.0);
for(int i=;i<n;i++){
scanf("%lf%lf",&a[i].x,&a[i].y);
o1=o1+a[i];
}
for(int i=;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=;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=;i<n;i++){//[0,n)~[i,i+n)
bool f=;
for(int j=;j<n;j++){
if(cmp(a[j].abs(),b[(i+j)%n].abs())!=){
f=;break;
}
}
if(f){
ans=min(ans,acos(dot(a[],b[i])/a[].abs()/b[i].abs()));
}
}
printf("%.11f\n",ans);
}
/**
先求出来每个点对中心的偏移量。
确实是乱序的,但是我们可以极角排序,这样一定是一个
循环移位匹配到另一个循环移位。
所以枚举n次就可以,check也十分简单
*/
100837D的更多相关文章
随机推荐
- mui选择器的坑
mui框架最近比较火,因为在移动端的页面展示效果太好了,web页面相当于APP的效果.连二年级的小明同学都知道了..你别说你不知道哦 但是这毕竟是一个不成熟的框架,维护和解决方案都跟不上,因此新手入坑 ...
- mysql用户管理与备份
用户管理 我们知道在Mysql中root用户是最高权限的用户,其他用户的创建和权限授予都是通过root用户来操作的 查看用户 在root用户界面下 select user,host,password ...
- Javaweb学习笔记——(二十七)——————泛型、泛型的通配符、反射泛型信息、反射注解、注解
泛型 1.泛型类:具有一个或多个类型变量的类,称之为泛型类 class A<T>{ } 2.在创建泛型实例时,需要为其类型变量赋值 A<String> a = new ...
- SpringBoot之解决云服务器VPS在所处云端集群的内网不能解析域名的问题:java.net.UnknownHostException:abc.cn: Temporary failure in name resolution
一.起因与原因分析过程 前端小伙伴儿告诉我,说服务器崩了. 请求数据接口,接口有响应,但报的json提示指向:数据库异常错误. 遂登陆云主机查看日志,核心记录显示如下: 2018-11-09 22:1 ...
- sessionStorage:写入记事本功能[内容写入sessionStorage中,读取,删除]
知识点: 1.设置sessionStorage----setItem:sessionStorage.setItem(key,data); 存储数据使用key是唯一,不可重复,每触发都生成:如用一个固定 ...
- 第二章 Android系统与嵌入式开发
第二章 Android系统与嵌入式开发 第二章首先要先了解Android和嵌入式Lnux系统有什么区别和联系,嵌入式Linux系统是在嵌入式设备中运行Linux系统:Android系统是在嵌入式设备中 ...
- win10设置vscode的终端为管理员权限
右击vscode 点击属性选择兼容性,勾选 “以管理员身份运行此程序” 确定即可!
- tmux 后台训练
参考链接:https://blog.csdn.net/u014381600/article/details/54588531
- Java开发之@PostConstruct和@PreDestroy注解
从Java EE5规范开始,Servlet增加了两个影响Servlet生命周期的注解(Annotation):@PostConstruct和@PreConstruct.这两个注解被用来修饰一个非静态的 ...
- (七)Knockout 创建自定义绑定
创建自定义绑定 你可以创建自己的自定义绑定 – 没有必要非要使用内嵌的绑定(像click,value等).你可以你封装复杂的逻辑或行为,自定义很容易使用和重用的绑定.例如,你可以在form表单里自定义 ...