B00007 快速模幂运算的两个C语言程序

最新推荐文章于 2023-09-05 02:56:29 发布

weixin_33895604

最新推荐文章于 2023-09-05 02:56:29 发布

阅读量69

点赞数

CC 4.0 BY-SA版权

文章标签： c/c++

原文链接：http://www.cnblogs.com/tigerisland/p/7564856.html

本文提供了两种实现64位整数模乘运算的方法。一种通过位操作实现，另一种利用浮点运算来减少溢出风险。这两种方法适用于需要进行大整数运算但不超过64位整数范围的应用场景。

代码来自维基百科的Modular arithmetic。

这两段代码都不是大整数计算的程序，是2进制64整数的计算程序，数据不能大于2进制63位。两段代码分别如下：

uint64_t mul_mod(uint64_t a, uint64_t b, uint64_t m)
{
   uint64_t d = 0, mp2 = m >> 1;
   int i;
   if (a >= m) a %= m;
   if (b >= m) b %= m;
   for (i = 0; i < 64; ++i)
   {
       d = (d > mp2) ? (d << 1) - m : d << 1;
       if (a & 0x8000000000000000ULL)
           d += b;
       if (d > m) d -= m;
       a <<= 1;
   }
   return d;
}

uint64_t mul_mod(uint64_t a, uint64_t b, uint64_t m)
{
   long double x;
   uint64_t c;
   int64_t r;
   if (a >= m) a %= m;
   if (b >= m) b %= m;
   x = a;
   c = x * b / m;
   r = (int64_t)(a * b - c * m) % (int64_t)m;
   return r < 0 ? r + m : r;
}

转载于:https://www.cnblogs.com/tigerisland/p/7564856.html