Educational Codeforces Round 1 C. Nearest vectors

寻找最近向量

最新推荐文章于 2017-07-20 08:38:23 发布

原创最新推荐文章于 2017-07-20 08:38:23 发布 · 411 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#codeforces

---------------------计算几何同时被 2 个专栏收录

6 篇文章

订阅专栏

极角排序

1 篇文章

订阅专栏

部署运行你感兴趣的模型镜像

C. Nearest vectors

题意:输入1e5范围内的点，每个点表示原点到该点的向量，并且向量的编号就是输入的序号(从1开始)；

原本一道很水的题，就不贴了，但是这道题的坑点真是2333，坑精度。。使用double被被人hack了几次，不贴出来难解心头之恨啊。。。
思路:使用atan2(y,x)来得到与原点的夹角(范围为-PI~PI,即y > 0时，表示逆时针否则为负的顺时针)，long double 存储角度。。。还要注意就是最小的角度和最大的角度之间要形成一个“环”;

62ms
#include<bits/stdc++.h>
using namespace std;
#define rep(i,n) for(int i = 1;i <= (n);i++)
typedef pair<long double,int> PII;
#define A first
#define B second
const int MAXN = 100100;
const long double PI = acos(-1.0)*2;
PII p[MAXN];
int main()
{
    int n,i,x,y;
    cin>>n;
    rep(i,n){
        scanf("%d%d",&x,&y);
        long double angle = atan2(y,x);
        if(angle < 0) angle += PI;
        p[i] = PII{angle,i};
    }
    sort(p+1,p+1+n);
    p[n+1] = {p[1].A+PI,p[1].B};
    int l,r;long double mn = PI;
    rep(i,n){
        if(p[i+1].A - p[i].A < mn){
            mn = p[i+1].A - p[i].A;
            l = p[i].B,r = p[i+1].B;
        }
    }
    cout<<l<<" "<<r;
}