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;
		
	}
}