Mersenne素数

Mersenne素数

如果一个素数可以写成2^p-1的形式，那么该素数就称为Mersenne素数，其中p是一个正整数。编写程序找出p<=100的所有Mersenne素数，然后显示如下输出。图片中的输出数据并不完整（数据太大，程序要运行几个小时，就不弄全图啦）。

下面是代码：

package java实验1019;
import java.math.BigInteger;
public class java_1019 {

	/**
	 * @param args
	 */
	final static int s=2;//定义常量，便于算式中a^p-1的a的修改
	public static void main(String[] args) {
		// TODO Auto-generated method stub
		BigInteger a=BigInteger.valueOf(s);//定义一个数值为2的BigInteger对象
		BigInteger b=BigInteger.valueOf(s);
		BigInteger c=new BigInteger("1");
		System.out.println("p\t\t2^p-1");
		for(int i=1;i<=100;)
		{
			if(i==1)//由于1和2均不是素数，所以从3开始；p=2时，a=3
			{
			a=a.multiply(b).subtract(c);//a=a*2-1;
			}
			else
			{
				a=a.add(c).multiply(b).subtract(c);//a=(a+1)*2-1
			}
			i++;
			if(素数(a)==true)
				{
				System.out.println(i+"\t\t"+a);
				}
		}

	}
	private static BigInteger BigInteger(String string) {
		// TODO Auto-generated method stub
		return null;
	}
	public static boolean 素数(BigInteger a)
	{
		int j=2;
		BigInteger i=BigInteger.valueOf(j);
		BigInteger s=new BigInteger("0");
		BigInteger b=new BigInteger("1");
		while(i.compareTo(a)<0)//从2开始计算，如果a的数值能整除i的数值，就不是素数
		{
			if(a.remainder(i).compareTo(s)==0)//BigInteger.remainder()计算余数
			{
				return false;
			}
			i=i.add(b);//i=i+1；
		}
		return true;
		
	}
}