计算 2 的极大幂

2022-09-03 03:26:02

我用Java制作了一个计算2次幂的程序,但它似乎非常低效。对于较小的功率(例如2 ^ 4000),它可以在不到一秒钟的时间内完成。但是,我正在计算2 ^ 43112609,它比已知最大的素数大一个。超过1200万位数字,需要很长时间才能运行。这是我到目前为止的代码:

import java.io.*;

public class Power
{
 private static byte x = 2;
 private static int y = 43112609;
 private static byte[] a = {x};
 private static byte[] b = {1};
 private static byte[] product;
 private static int size = 2;
 private static int prev = 1;
 private static int count = 0;
 private static int delay = 0;
 public static void main(String[] args) throws IOException
 {
  File f = new File("number.txt");
  FileOutputStream output = new FileOutputStream(f);
  for (int z = 0; z < y; z++)
  {
   product = new byte[size];
   for (int i = 0; i < a.length; i++)
   {
    for (int j = 0; j < b.length; j++)
    {
     product[i+j] += (byte) (a[i] * b[j]);
     checkPlaceValue(i + j);
    }
   }
   b = product;
   for (int i = product.length - 1; i > product.length - 2; i--)
   {
    if (product[i] != 0)
    {
     size++;
     if (delay >= 500) 
     {
      delay = 0;
      System.out.print(".");
     }
     delay++;
    }
   }
  }
  String str = "";
  for (int i = (product[product.length-1] == 0) ? 
   product.length - 2 : product.length - 1; i >= 0; i--)
  {
   System.out.print(product[i]);
   str += product[i];
  }
  output.write(str.getBytes());
  output.flush();
  output.close();
  System.out.println();
 }

 public static void checkPlaceValue(int placeValue)
 {
  if (product[placeValue] > 9)
  {
   byte remainder = (byte) (product[placeValue] / 10);
   product[placeValue] -= 10 * remainder;
   product[placeValue + 1] += remainder;
   checkPlaceValue(placeValue + 1);
  }
 }  
}

这不是为了学校项目或任何东西;只是为了好玩。任何关于如何使这更有效率的帮助将不胜感激!谢谢!

凯尔

附言:我没有提到输出应该是以10为基数,而不是二进制。


答案 1

这里的关键是要注意:

2^2 = 4
2^4 = (2^2)*(2^2)
2^8 = (2^4)*(2^4)
2^16 = (2^8)*(2^8)
2^32 = (2^16)*(2^16)
2^64 = (2^32)*(2^32)
2^128 = (2^64)*(2^64)
... and in total of 25 steps ...
2^33554432 = (2^16777216)*(16777216)

然后从:

2^43112609 = (2^33554432) * (2^9558177)

你可以找到使用相同方法的剩余部分,并且由于,您可以找到使用相同的方法,并且由于,您可以找到使用相同的方法,依此类推......(2^9558177)(2^9558177 = 2^8388608 * 2^1169569)2^1169569(2^1169569 = 2^1048576 * 2^120993)2^120993

编辑:以前本节中有一个错误,现在已修复:

此外,通过注意以下几点进一步简化和优化:

2^43112609 = 2^(0b10100100011101100010100001)
2^43112609 = 
      (2^(1*33554432))
    * (2^(0*16777216))
    * (2^(1*8388608))
    * (2^(0*4194304))
    * (2^(0*2097152))
    * (2^(1*1048576))
    * (2^(0*524288))
    * (2^(0*262144))
    * (2^(0*131072))
    * (2^(1*65536))
    * (2^(1*32768))
    * (2^(1*16384))
    * (2^(0*8192))
    * (2^(1*4096))
    * (2^(1*2048))
    * (2^(0*1024))
    * (2^(0*512))
    * (2^(0*256))
    * (2^(1*128))
    * (2^(0*64))
    * (2^(1*32))
    * (2^(0*16))
    * (2^(0*8))
    * (2^(0*4))
    * (2^(0*2))
    * (2^(1*1))

另请注意,2^(0*n) = 2^0 = 1

使用此算法,您可以计算出 、 、 、 ... 在25个乘法中。然后,您可以转换为其二进制表示形式,并使用少于25个乘法轻松找到。总共,您需要使用少于50个乘法来查找0到67108864之间的任何位置。2^12^22^42^82^162^33554432431126092^431126092^nn


答案 2

以二进制文件显示它既简单又快速 - 尽可能快地写入磁盘!100000...... :D


推荐