遞歸函數中的數值精度

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