快速乘法模2 ^ 16 1

前端之家收集整理的这篇文章主要介绍了快速乘法模2 ^ 16 1前端之家小编觉得挺不错的,现在分享给大家,也给大家做个参考。
IDEA密码使用乘法模2 ^ 16 1.是否有算法在没有通用模运算符的情况下执行此操作(仅模2 ^ 16(截断))?在IDEA的上下文中,零被解释为2 ^ 16(这意味着零不是我们乘法的参数,它不能是结果,因此我们可以保存一位并将值存储为2 ^ 16作为位模式0000000000000000).我想知道如何在不使用标准模运算符的情况下有效地实现它(或者是否有可能).

解决方法

您可以利用(N-1)%N == -1这一事实.

因此,(65536 * a)%65537 == -a%65537.
此外,-a%65537 == -a 1(mod 65536),当0被解释为65536时

  1. uint16_t fastmod65537(uint16_t a,uint16_t b)
  2. {
  3. uint32_t c;
  4. uint16_t hi,lo;
  5. if (a == 0)
  6. return -b + 1;
  7. if (b == 0)
  8. return -a + 1;
  9.  
  10. c = (uint32_t)a * (uint32_t)b;
  11. hi = c >> 16;
  12. lo = c;
  13. if (lo > hi)
  14. return lo-hi;
  15. return lo-hi+1;
  16. }

这里唯一的问题是如果hi == lo,结果将是0.幸运的是,测试套件确认,它实际上不能……

  1. int main()
  2. {
  3. uint64_t a,b;
  4. for (a = 1; a <= 65536; a++)
  5. for (b = 1; b <= 65536; b++)
  6. {
  7. uint64_t c = a*b;
  8. uint32_t d = (c % 65537) & 65535;
  9. uint32_t e = m(a & 65535,b & 65535);
  10. if (d != e)
  11. printf("a * b % 65537 != m(%d,%d) real=%d m()=%d\n",(uint32_t)a,(uint32_t)b,d,e);
  12. }
  13. }

输出:无

猜你在找的C&C++相关文章