我的二分搜索应用程序不在这里吗？

Question

所以，我正在解决这个问题：http://www.spoj.com/problems/immersed/

一个奇妙的发现即将在生物学领域发生，你是进行研究的团队的一员。这项研究是关于测量细胞在无氧和有毒物质存在的环境中的生长。研究小组得出了一个精确的假设，分析的数据告诉他们：生长、天数和毒性；由以下公式关联：

P=n*n^cn；

第一行是T(1≤T≤400000)，测试用例的数量，然后是T个测试用例。

每个测试用例都是一条用空格分隔的2个整数(P，c)的线。

P(1≤P≤10¹⁵)

#define eps 1e-7
const double cont = 14.0;
double p,c;
double binary (double start, double end);
double func (double n);
int main (void)
{
    int t;
    cin>>t;
    while (t != 0)
    {
        cin>>p>>c;
        double mid,ans,start=0,end=cont;
        ans = binary(start,end);
        cout.precision(6);
        cout<<fixed<<ans<<"\n";         
        t--;
    }
    return 0;
}

double func (double n)
{
    double ret = n*pow(n,c*n);
    return ret;
}

double binary (double start, double end)
{
    if (start <= end)
    {
        double mid = start + (end-start)/2.0;
        if (func(mid)-p <= eps)
            return mid;
        else if (func(mid)-p > eps)
            return binary(start,mid);
        else
            return binary(mid,end);
    }
}

Input:
3
1 1
3 4
100 1
Output:
1.000000
1.207384
3.086308

My output (for the above input)

0.875
0.875
1.75

PS：我没有张贴库和所有，以避免杂乱。此外，我将设置为小数点后6位，一旦我得到正确的值。我只想知道，是我的逻辑不正确还是我的二分搜索被错误地实现了？

编辑：我提交的新代码

double binary (double start, double end)
    {
        if (start <= end)
        {
            double mid = start + (end-start)/2.0;
            double delta = func(mid) - p;
            if (delta < -1*eps)
                return binary(mid,end);
            else if (delta > eps)
                return binary(start,mid);
            else
                return mid;
        }
    }

Answer 1

您正在以不合理的顺序进行测试：

if (func(mid)-p <= eps)
    return mid;
else if (func(mid)-p > eps)
    return binary(start,mid);
else
    return binary(mid,end);

试试看

if (func(mid)-p < -eps)
    return binary(mid,end);
else if (func(mid)-p > eps)
    return binary(start,mid);
else
    return mid;

我也不确定这两个递归调用。我保留了你的逻辑，哪种情况下叫哪种（因为我可能误解了你的公式），但他们回头看我。

double binary (double start, double end)
{
    for (;;)
    {
        double mid = (start + end) * .5;
        double delta=func(mid)-p;

        if ( delta < -eps)
            start = mid;
        else if (delta > eps)
            end = mid;
        else
            return mid;
    }
}