当p是素数并且Mp也是素数时,称Mp为默森尼数,
例如,p=5,Mp=2^5-1=31,5和31都是素数,因此31是默森尼数;
而Mp=2^11-1=2047,11是素数,2047不是素数,因此2047不是默森尼数
#include <stdio.h>
#include <math.h>
int JudgePrime(int p) //判断一个数是否为素数,是则返回0,不是返回-1
{
int flag = 0;
int i = 0;
for (i = 2; i < p; i++)
{
if (p % i == 0)
{
flag = 1;
}
}
if (flag == 0)
{
return 0;
}
return -1;
}
int MersenneNum(long *Arr) //求Mersenne数,放入数组中
{
long int Mp = 0;
int p = 0;
int count = 0;
for (p = 2; p <= 10000; p++)
{
if (JudgePrime(p) == 0) //判断p是不是素数
{
Mp = pow(2, p) - 1;
if (JudgePrime(Mp) == 0)
{
Arr[count++] = Mp; //判断Mp是不是素数,是则放入数组中
if (count >= 7)
{
break; //算出7个则退出
}
}
}
}
return count;
}
void PrintArr(long *Arr, int count) //打印结果
{
int i = 0;
for (i = 0; i < count; i++)
{
printf ("%ld ", Arr[i]);
}
printf ("\n");
}
int main()
{
long int Arr[100] = {0}; //存放求得的Mersenne number
int count = 0;
count = MersenneNum(Arr);
PrintArr(Arr, count);
return 0;
}