численная точность в рекурсивных функциях

Question

Я озадачен, оценивая следующую функцию, она производит число до F (0,8, 172, 1), но когда я увеличиваю 172 до 173, результат становится бесконечным. Я подозреваю, что существует проблема с числовой точностью?

double F(double d, int c, int t) { 
    // base cases 
    if ((c==1 && t==1) || (c==0 && t==0)) 
     return 1.; 
    if (c==0 || t==0) 
     return 0.; 
    if (t>c) 
     return 0.; 
    return F(d,c-1,t-1) + (c-1 - t*d)*F(d,c-1,t); 
} 

Answer 1

Я не знаю, что ваша функция должна делать, но учитывая аргументы: F(0.8, 172, 1) возвращаемое значение 4.41861e+306, который чуть меньше максимального значения а double может представлять:

// 1.79769e+308 
std::cout << std::numeric_limits<double>::max() << std::endl; 

Когда 172 заменяется на 173, возвращаемое значение превышает максимальное значение a double может представлять и становится положительной бесконечностью. Это можно сделать, изменив тип возврата F на long double, который приводит к значению 7.56466e+308