我的二进制搜索应用程序不在这里

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

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

p=N*N^cN；

其中p是在数千个细胞中测量的生长。

N是经过的天数。

c是与实验毒性水平相关的常数。

当这些细胞达到特定生长时，你的生物学伙伴需要从这些细胞中取出一些组织。他们要求你编写一个程序，告诉他们发生这种情况的确切时间，考虑到毒性水平和所需的生长。

输入

第一行是T（1≤T≤40000），测试用例的数量，然后是T

每个测试用例都是一行，其中有2个整数（Pc），用空格分隔。

p（1≤p≤10¹⁵）

c（1≤c≤5）

输出

对于每个测试用例，您必须以十进制格式输出预期时间。

我所做的是使用二进制搜索来查找天数，如下所示：

#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

附言：我没有张贴图书馆和所有为了避免混乱。此外，一旦得到正确的值，我会将其设置为小数点后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;
        }
    }