C++语言: 高亮代码由发芽网提供

/**
*返回0表示x==0，返回-1表示 x < 0, 返回1表示 x > 0
*accuracy(a-b) == 0, a == b 或者 fabs(a-b) < eps
*accuracy(a-b) != 0, a != b 或者 fabs(a-b) > eps
*accuracy(a-b) <  0, a <  b 或者 a - b < -eps
*accuracy(a-b) >  0, a >  b 或者 a - b > eps
*accuracy(a-b) <= 0, a <= b 或者 a - b < eps
*accuracy(a-b) >= 0, a >= b 或者 a - b > -eps
*/
#define eps 1e-12
int accuracy(double x){ return x < -eps?-1:((x < eps)?0:1);}
int isOk(double x){
    //尝试这个高度
    double total = 0;
    for(int i = 0;i < m;i++){
        for(int j = 0; j < n;j++){
            if(accuracy(x-square[i][j]) > 0){
                total += 100*(x - square[i][j]);
            }
        }
    }
    return accuracy(total-all);
}
double binaryTrial(double left,double right){
    double mid;
    while(accuracy(left-right) < 0){
        mid = (left + right)*1.0/2;
        int t  = isOk(mid);
        if(t <= 0){left = mid;}
        else{right = mid;}
    }
    return mid;
}