Vijos隐形的翅膀

Vijos隐形的翅膀

背景

小杉终于进入了天堂。他看到每个人都带着一双隐形翅膀，他也想要。

（小杉是怎么看到的？……）

描述

天使告诉小杉，每只翅膀都有长度，两只翅膀的长度之比越接近黄金分割比例，就越完美。

现在天使给了小杉N只翅膀，小杉想挑出一对最完美的。

格式

输入格式

每组测试数据的 第一行有一个数N(2<=N<=30000) 第二行有N个不超过1e5的正整数，表示N只翅膀的长度。

20%的数据N<=100

输出格式

对每组测试数据输出两个整数，表示小杉挑选出来的一对翅膀。

注意，比较短的在前，如果有多对翅膀的完美程度一样，请输出最小的一对。

样例1

样例输入1

4
2 3 4 6

Copy

样例输出1

2
3

Copy

限制

每个测试点1s

提示

你可以认为黄金分割比就是0.6180339887498949

分析

我们可以将所有数据从小到大快排，然后对于数组中的每个数a[i]，去除以黄金分割比可以得到一个数，再在数组中找和所得数相近的数，用找到的数除以a[i],最后找和黄金分割比最相近的一组数。

代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath> using namespace std;

const double xyz = 0.6180339887498949;
int n;
double a[30005];
double ans = 9999999999.0;
int x, y;

int search(double e, int p)
{
    int l = p + 1, r = n;
    while (l <= r)
    {
        int mid = (l + r) / 2;
        if (a[mid] <= e)
            l = mid + 1;
        else
            r = mid - 1;
    }
    double w = fabs(a[l] - e);
    double b = fabs(a[l - 1] - e);
    double c = fabs(a[l + 1] - e);
    if (w < b && w < c)
        return l;
    else if (b < w && b < c)
        return l - 1;
    else
        return l + 1;
}

int main()
{
    cin >> n;
    for (int i = 1; i <= n; i++)
        cin >> a[i];
    sort(a + 1, a + n + 1);
    for (int i = 1; i <= n; i++)
    {
        double t = (a[i] / (xyz));
        int k = search(t, i);
        double temp = fabs(a[i] / a[k] - xyz);
        if (temp < ans)
        {
            ans = temp;
            x = a[i];
            y = a[k];
        }
    }
    cout << x << endl
         << y << endl;
    return 0;
}