隐形的翅膀

**

隐形的翅膀

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

这道理如果采用枚举的方法的话，写两重循环显然超时 O(n^2)
我们先用进行快速排序简单的二分就可以完成了，比较的时候采用绝对值的方法来进行比较
但要注意的是，不一定会正好把这个区间二分，说以要左右比较，但我们先把比较前的数先存起来
举个例子：
排序之后的数据为
2 3 4 6 8
下标mid=2
2/4=0.5<gold
mid=0+(1-0)/2
2/2=1 更加离谱，所以2的partner就是4啦
以此类推下去

#include<iostream>
#include<algorithm>
#include<cmath>
const double N=0.6180339887498949;
using namespace std;
int a[30010];

int main()
{
    int n,i,l,r,mid,ans,z,y;
    double minn=100000.0,ha;
    cin>>n;
    for(i=1;i<=n;i++)
     cin>>a[i];
    sort(a+1,a+n+1);
    for(i=1;i<n;i++)
     {
        l=i+1;
        r=n;
        while(l<=r)
        {
            mid=(l+r)/2;
            ha=double(a[i])/(double(a[mid]));
            if(fabs(ha-N)<fabs(minn-N))//先来比较一下再来更新最合适的数
            {
             minn=ha;
             z=a[i];
             y=a[mid];
            }
            if(ha<N)
             r=mid-1;
            else
             l=mid+1;
        }
     }
     cout<<z<<endl<<y;
     return 0;
}

|# 状态 耗时 内存占用
#1 Accepted 2ms 216.0 KiB
#2 Accepted 1ms 232.0 KiB
#3 Accepted 16ms 420.0 KiB
#4 Accepted 15ms 420.0 KiB
#5 Accepted 16ms 424.0 KiB
#6 Accepted 16ms 420.0 KiB
#7 Accepted 15ms 420.0 KiB
#8 Accepted 15ms 416.0 KiB
#9 Accepted 16ms 408.0 KiB
#10 Accepted 15ms 424.0 KiB

学会程序和算法，走遍天下都不怕

大理古城