Math 三模乘法

Math 三模乘法,math,multiplication,modulus,integer-overflow,Math,Multiplication,Modulus,Integer Overflow,我计算的金额如下: (a[x]+ma[x-1]+2ma[x-2]+3m*a[x-3]+..%MOD(MOD=1e9+7) 为此,我使用这个循环 long long mulmod(long long a,long long b,long long c) { if (a == 0 || b == 0) return 0; if (a == 1) return b; if (b == 1) return a; long

我计算的金额如下:

(a[x]+ma[x-1]+2ma[x-2]+3m*a[x-3]+..%MOD(MOD=1e9+7)

为此,我使用这个循环

long long mulmod(long long a,long long  b,long long c)
{
    if (a == 0 || b == 0)
        return 0;
    if (a == 1)
        return b;
    if (b == 1)
        return a;
    long long a2 = mulmod(a, b / 2, c);
    if ((b & 1) == 0) 
    {
        return (a2 + a2) % c;
    } 
    else 
    {
        return ((a % c) + (a2 + a2)) % c;
    }
}

res=a[x]%MOD;
for(i=x-1;i>=1;i--)
    res=(res%MOD+mulmod(mulmod(x-i,m,MOD),a[i],MOD))%MOD;
然而,这仍然给了我溢出错误。基本错误,我猜是在(abc)%MOD中


谢谢。

您需要将以下模块算术标识合并到程序中以避免溢出:

(A+B+…)mod C=(A mod C+B mod C+…mod C)mod C


(A*B*…)mod C=(A mod C*B mod C*…mod C)mod C
您需要将以下模块化算术标识合并到程序中以避免溢出:

(A+B+…)mod C=(A mod C+B mod C+…mod C)mod C


(A*B*…)mod C=(A mod C*B mod C*…mod C)mod C

你到底想找什么
1e9+7
看起来像是一个竞赛问题(也许是Euler项目),而这些问题通常都有一种鬼鬼祟祟的方式回到答案中。不过,简单的答案是使用您的语言版本的BigInteger。Long-Long-int不可能适合3个非常大的数字的乘法顺序。所以我需要创建一个计算它的函数。如果任何参数为负,mulmod都不能正常工作。
res=a[x]%MOD
。你确定所有的a[i]都>=0你到底想找什么
1e9+7
看起来像是一个竞赛问题(也许是Euler项目),而这些问题通常都有一种鬼鬼祟祟的方式回到答案中。不过,简单的答案是使用您的语言版本的BigInteger。Long-Long-int不可能适合3个非常大的数字的乘法顺序。所以我需要创建一个计算它的函数。如果任何参数为负,mulmod都不能正常工作。
res=a[x]%MOD
。你确定(i=x-1;i>=1;i--)的所有a[i]都>=0
res=(res%MOD+((x-i)%MOD)*(m%MOD)*(a[i]%MOD))%MOD我已将代码更改为此,但错误仍然存在。这些行:
返回(a2+a2)%c
返回((a%c)+(a2+a2))%c也有溢出的可能性。如果我直接使用上面使用的乘法属性并用这一行替换函数:
res=(res%MOD+((x-I)%MOD)*(m%MOD)*(a[I]%MOD))%MOD)%MOD这不应该工作吗?你的
(x-i)%MOD)
术语不应该是
(a[x-i])%MOD)
用于(i=x-1;i>=1;i--)
res=(res%MOD+((x-i)%MOD)*(m%MOD a[i]%MOD我已将代码更改为此,但错误仍然存在。这些行:
返回(a2+a2)%c
返回((a%c)+(a2+a2))%c也有溢出的可能性。如果我直接使用上面使用的乘法属性并用这一行替换函数:
res=(res%MOD+((x-I)%MOD)*(m%MOD)*(a[I]%MOD))%MOD)%MOD这不应该工作吗?你的
(x-i)%MOD)
术语不应该是
(a[x-i])%MOD)