大数场景下Mod()返回0导致斐波那契数判断异常求助

问题根源 #

你遇到的是双精度浮点数的精度限制问题。Matlab默认使用64位双精度浮点数，它只能精确表示小于等于2^53（约9.007e15）的整数。当数值超过这个范围时，整数无法被精确存储——比如10^20+1在双精度中会被存储为和10^20完全相同的近似值，导致mod(10^20+1,2)返回0，而非预期的1。

你的原代码中，对大数计算sqrt(5*a.^2±4)时，浮点数的精度丢失会让结果出现偏差，进而导致mod(...,1)~=0的判断完全失效，错误地将大数判定为斐波那契数。

解决思路 #

核心是用整数类型替代浮点数运算，避免精度丢失：

• 将输入转换为64位无符号整数（uint64），确保所有运算在整数域内进行。
• 改用整数平方根判断5n²±4是否为完全平方数，避免浮点数sqrt的近似误差。

修改后的代码 #

function y = fib_count(x)
    % 转换为64位无符号整数，规避浮点数精度问题
    a = unique(uint64(x));
    is_fib = false(size(a));
    
    for i = 1:length(a)
        n = a(i);
        val1 = 5 * n * n + 4;
        val2 = 5 * n * n - 4;
        
        % 使用整数平方根函数isqrt（Matlab R2020b及以上支持），返回不超过目标值的最大整数
        sqrt_val1 = isqrt(val1);
        sqrt_val2 = isqrt(val2);
        
        % 判断是否为完全平方数
        if sqrt_val1 * sqrt_val1 == val1 || sqrt_val2 * sqrt_val2 == val2
            is_fib(i) = true;
        end
    end
    
    y = sum(is_fib);
end

注意事项 #

• 如果你的Matlab版本低于R2020b，没有isqrt函数，可以用uint64(floor(sqrt(double(val))))临时替代，但要注意：当val超过2^53时，double(val)仍会有精度丢失。这种情况下可以自己实现整数平方根的二分法，确保全程用整数运算。
• 如果输入包含负数或0，可根据Problem 44481的题目要求添加前置过滤逻辑，因为斐波那契数通常定义为正整数序列。

内容的提问来源于stack exchange，提问作者Trà Đức Thịnh