Question

我用JavaScript编写了这段代码，但是对于大输入，它返回了错误的数字它应该用模数（mo）计算指数（ex）幂的基数。
我用C编写了相同的代码并且正在工作。请有人告诉我有什么问题。
尝试以模10 ^ 9 +7来测试费马定理。

function powFun(base, ex, mo) {
    var r;
    if(ex === 0) 
        return 1;
    else if(ex % 2 === 0) {
        r = powFun(base, ex/2, mo) % mo ;
        return (r * r) % mo;
    }else 
        return (base * powFun(base, ex - 1, mo)) % mo;
}

Answer 1

此行发生溢出：

return (r * r) % mo;

有关实施this question的一些简单算法，请参阅a*a mod n without overflow的答案。我已将one of the answers改编为以下Javascript：

＆＃13;

function addmod(x, y, n)
{
    // Precondition: x<n, y<n
    // If it will overflow, use alternative calculation
    if (x + y <= x) x = x - (n - y) % n;
    else x = (x + y) % n;
    return x;
}

function sqrmod(a, n)
{
    var b;
    var sum = 0;

    // Make sure original number is less than n
    a = a % n;

    // Use double and add algorithm to calculate a*a mod n
    for (b = a; b != 0; b >>= 1) {
        if (b & 1) {
            sum = addmod(sum, a, n);
        }
        a = addmod(a, a, n);
    }
    return sum;
}

function powFun(base, ex, mo) {
    var r;
    if(ex === 0) 
        return 1;
    else if(ex % 2 === 0) {
        r = powFun(base, ex/2, mo) % mo ;
        // return (r * r) % mo;
        return sqrmod(r, mo);
    }else 
        return (base * powFun(base, ex - 1, mo)) % mo;
}

result = powFun(4, 1000000006, 1000000007);

alert(result);

＆＃13;

要支持更大的数字，请使用支持大整数的库。

Javascript自制动力与模数

1 个答案: