Math 三模乘法
我计算的金额如下: (a[x]+ma[x-1]+2ma[x-2]+3m*a[x-3]+..%MOD(MOD=1e9+7) 为此,我使用这个循环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
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)
?