前幾天在網上看到有人曬一道面試題，寫一個函數求平方根，該函數帶2個參數，第一個參數是目標數字，第二個參數是精度。即：

double sqrt(double target, double g);
a = sqrt(t, g);
// 要求 |a^2 - t| < g

原文看到題目，就沒有接著往下看，想了一會兒，就有思路了，覺得拿來做面試題不錯。可是沒想到面了好幾個人，都沒有人能給出稍微好看點的代碼。難道大家工作久了之後，寫代碼反而生疏了嗎？一直手癢，今晚就花了幾分鐘實現了一下，感覺還是比較簡單的。實現完成之後，感覺不錯，給出的答案都滿足要求了，程序應該是對了。但是後來仔細一思考，又打了些程序處理過程的trace出來，發現不對了，有bug. 什麽bug呢？收斂比較慢，和原來二分查找的設想不一樣。雖然也能給出結果，但是收斂慢了，程序的實現和自己的思路不一致。想了一會兒，給出了正確的代碼。該sqrt函數實現大約24行左右，如下：

#include <iostream>
using namespace std;
 
 
#define abs(X) ((X) > 0 ? (X) : -(X))
 
double sqrt(double target, double g)
{
    if (target < 0 || g < 0) return -1;
    
    double result = target;
    double small  = target > 1 ? 0 : target;
    double big    = target > 1 ? target : 1;
    
    while
 
 (true) {
        double diff = result * result - target;
        
        if (abs(diff) < g) {
            break;
        }
        else if (diff > 0) {
            big = result;
        }
        else if (diff < 0) {
            small = result;
        }
        result = (big + small) / 2;
    }
    return
 
 result;
}
 
int main()
{
    cout << "sqrt(10, 0.1) = " << sqrt(10, 0.1) << endl;
    cout << "sqrt(25, 0.001) = " << sqrt(25,0.001) << endl;
    
    cout << "sqrt(0.09, 0.001) = " << sqrt(0.09, 0.001) << endl;
    cout << "sqrt(0.4, 0.01) = " << sqrt(0.4, 0.01) << endl;
    
    return 0;
}

求平方根的面試題