【计算思维作业】A.隐形的翅膀

问题

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

输入数据

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

输出数据

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

样例输入

4
2 3 4 6

样例输出

2
3

样例说明

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

思路

1.思路，TLE

最朴素的思路，两层for循环，得到所有的比率情况，然后从中筛选出最接近黄金分割比的。用数学语言描述就是

但是根据题目所给的N的范围，最大可到，平方后已经到了的量级，有超时的风险，提交到oj上也确实超时了。所以这个思路肯定不是最简单的。

2.思路，AC

换个角度思考，如果我们求最接近黄金比的比率，是不是可以反过来，用（小翅膀长/黄金比）来求出理想大翅膀长，然后去用二分查找法去搜索和理想大翅膀长最接近的翅膀长度。

这样的话

        1.输入保存数组是

        2.sort快排是

        3.二分查找

三项相加，整体的时间复杂度就是

代码

1.自己写的TLE

错误的代码也贴出来供对比

#include<iostream>
#include<vector>
#include<cmath>
#include<algorithm>
#define goldratio 0.6180339887498949
using namespace std;

//tle 自己写的
int main()
{
    int n,biggertmp,smallertmp,biggerwings,smallerwings;
    double ratio,tmpratio;
    ratio=100;
    cin>>n;
    vector<int> vec(n,0);
    for(int i=0;i<n;i++)
        cin>>vec[i];
    for(int i=0;i<n-1;i++)    //两层for循环，超时了
    {
        for(int j=i+1;j<n;j++)
        {
            smallertmp=vec[i];
            biggertmp=vec[j];
            tmpratio=smallertmp*1.0/biggertmp;
            if(fabs(tmpratio-goldratio)<fabs(ratio-goldratio)){
                ratio=tmpratio;
                biggerwings=biggertmp;
                smallerwings=smallertmp;
            }
        }
    }
    cout<<smallerwings<<endl;
    cout<<biggerwings<<endl;
    return 0;
}

2.自己写的AC

#include<iostream>
#include<vector>
#include<cmath>
#include<algorithm>
#define goldratio 0.6180339887498949
using namespace std;
int main()
{
    int n,right,left,mid,biggerwings,tmpbigger,smallerwings,tmpsmaller;
    double tmpratio,ratio,bestwings;
    ratio=2;
    cin>>n;
    vector<int> vec(n,0);
    for(int i=0;i<n;i++)
        cin>>vec[i];
    sort(vec.begin(),vec.end());
    for(int i=0;i<n-1;i++)
    {
        left=i+1;
        right=n-1;
        bestwings=vec[i]/goldratio;
        while(right-left>1)
        {
            mid=(left+right)/2;
            if(vec[mid]>bestwings) right=mid;
            else if(vec[mid]<bestwings) left=mid;
        }
        if(fabs(vec[left]-bestwings)<fabs(vec[right]-bestwings))
            tmpbigger=vec[left];
        else tmpbigger=vec[right];
        tmpsmaller=vec[i];
        tmpratio=tmpsmaller*1.0/tmpbigger;
        if(fabs(tmpratio-goldratio)<fabs(ratio-goldratio))
        {
            biggerwings=tmpbigger;
            smallerwings=tmpsmaller;
            ratio=tmpratio;
        }
    }
    cout<<smallerwings<<endl;
    cout<<biggerwings<<endl;
    return 0;
}