c - 递归调和函数返回 NaN

Question

我编写了以下示例代码来查找 N. (1+1/2+1/3+...1/N) 的谐波值。阅读用 BOLD 编写的代码中的注释，帮助我找出为什么会发生这种情况。

#include <stdio.h>

float harmonic(float n, float har) {

    if(n==0) {
        return 0;
    }

    if(n==1) {
        printf("%f\n", har+1.0f);***/* This prints value 1.5000*/***
        return har+1.0f;
    }else{
        harmonic(n-1, (har+(1/n)));
    } 
} 

int main() 
{ 
    printf("%f\n", harmonic(2, 0.0f)); **/* But this prints value nan(Not a  Number)*/**  
    return 0; 
}

谢谢，纳迦

Answer 1

我想你想做：

return harmonic(n-1, (har+(1/n)));

Answer 2

我的第一个想法是您几乎不应该将浮点数与简单的相等性进行比较，因此“if(n==0)”应该是“if(n<=EPSILON)”，而“if(n==1)”应该是“if( n<= 1.0f + EPSILON)" 其中 EPSILON 是一个小的正分数，可能是 1.0e-5。取决于您可以依赖多少精度。

但后来我意识到n应该是一个int。在除法之前将其转换为浮点数。由于与“n”的比较意味着你冒着无限递归的风险。

考虑使用双精度而不是浮点数。

Matthew Flaschen 的回答说明了您收到 NaN 消息的真正原因。原始代码不会从“else”返回任何内容，因此调用者可能正在从堆栈中读取垃圾。因此，NaN。