2021.06.04麦森数

2021.06.04麦森数

题目描述

形如$2^P-1$的素数称为麦森数，这时P一定也是个素数。但反过来不一定，即如果P是个素数，$2^P-1$不一定也是素数。到1998年底，人们已找到了37个麦森数。最大的一个是P=3021377，它有909526位。麦森数有许多重要应用，它与完全数密切相关。

任务：从文件中输入P（1000<P<3100000），计算$2^P-1$的位数和最后500位数字（用十进制高精度数表示）

输入格式

文件中只包含一个整数P（1000<P<3100000）

输出格式

第一行：十进制高精度数$2^P-1$的位数。

第2-11行：十进制高精度数$2^P-1$的最后500位数字。（每行输出50位，共输出10行，不足500位时高位补0）

不必验证$2^P-1$与P是否为素数。

样例输入

1279

样例输出

386
00000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000
00000000000000104079321946643990819252403273640855
38615262247266704805319112350403608059673360298012
23944173232418484242161395428100779138356624832346
49081399066056773207629241295093892203457731833496
61583550472959420547689811211693677147548478866962
50138443826029173234888531116082853841658502825560
46662248318909188018470682222031405210266984354887
32958028878050869736186900714720710555703168729087

思路

1. 如果直接用BigInteger的pow，会导致内存超限。但是题目中只需要最后500个数字，所以可以采用modPow。
2. 采用公式推导求出大数的位数。
3. 问，-1和的位数是否不同呢？
   答：位数是相同的。因为对于，末尾数字不可能是0，只可能是2，4，6，8。

代码

	void test() {
		Scanner cin = new Scanner(System.in);
		int p = cin.nextInt();
		int k = (int)(p*Math.log10(2))+1; //求位数
		System.out.println(k);
		//求500数字
		BigInteger t = new BigInteger("2");
		//模pow，余10^500，求次方
		t = t.modPow(new BigInteger(""+p),new BigInteger("10").pow(500)); 
		t = t.subtract(new BigInteger("1"));
		char[] chs = t.toString().toCharArray();
		int cnt = 0;
		for(int i = 500; i > 0; i--) { //输出后500个数字，不足补0
			System.out.print(chs.length-i>=0?chs[chs.length-i]:"0");
			if(++cnt%50==0) System.out.println();
		}
	}
	String n0(int t) {
		StringBuffer sb = new StringBuffer();
		for(int i = 0; i < t; i++)  sb.append("0");
		return sb.toString();
	}