Vijos 1237 隐形的翅膀

题目

背景

小杉终于进入了天堂。他看到每个人都带着一双隐形翅膀，他也想要。
（小杉是怎么看到的？……）

描述

天使告诉小杉，每只翅膀都有长度，两只翅膀的长度之比越接近黄金分割比例，就越完美。
现在天使给了小杉N只翅膀，小杉想挑出一对最完美的。

格式

输入格式

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

输出格式

对每组测试数据输出两个整数，表示小杉挑选出来的一对翅膀。
注意，比较短的在前，如果有多对翅膀的完美程度一样，请输出最小的一对。

样例1

样例输入1

4
2 3 4 6

样例输出1

2
3

限制

每个测试点1s

提示

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

题解

离散（排序）+二分
注意：C++中实数运算要求表达式中至少有一个变量是实数，得出结果才能是实数

代码

#include <algorithm>
#include <stdio.h>
#include <cmath>

using namespace std;

int a[30005];
double const gold=0.6180339887498949;
int g,h;

void best(int i,int j)
{
    double c=a[i],d=g;
    if (fabs(c/a[j]-gold)<fabs(d/h-gold)) 
    {
        g=a[i];h=a[j];
    }
}

int main()
{
    int n;
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
       scanf("%d",&a[i]);
    sort(a+1,a+n+1);
    g=a[1];h=a[2];
    for (int i=1;i<=n;i++)
    {
        int l=1,r=i-1,mid;
        double c;
        do
        {
            mid=(l+r)/2;
            c=a[mid];
            if (c/a[i]<gold) l=mid; else r=mid;
        }   
        while (l<r-1);
        best(l,i);best(r,i);
        l=i+1;r=n;
        do
        {
            mid=(l+r)/2;
            c=a[mid];
            if (a[i]/c<gold) r=mid; else l=mid;
        }
        while (l<r-1);

        best(i,l);best(i,r);
    }
    printf("%d\n%d",g,h);
    return 0;
}