C程序求梅森数

形如2 的n次-1素数称为梅森数(Mersenne Number)。例如2 

2∧2−1=3、2 ∧3 −1=7都是梅森数。1722年,双目失明的瑞士数学大师欧拉证明了2 ∧31

 −1=2147483647是一个素数,堪称当时世界上“已知最大素数”的一个记录。

本题要求编写程序,对任一正整数n(n<20),输出所有不超过2 ∧n−1的梅森数。

输入格式:

输入在一行中给出正整数n(n<20)。

输出格式:

按从小到大的顺序输出所有不超过2 ∧n −1的梅森数,每行一个。如果完全没有,则输出“None”。

代码如下。

#include<stdio.h>

#include<math.h>

int sushu(int n);

int sushu(int n)

{

    int i;

    for(i=2;i<=sqrt(n);i++)

    {

        if(n%i==0)

        {

            return 0;

            break;

        }

    }

    return 1;

}//前为定义函数判断是否为素数,是素数反回1为真,不是则反回0为假。

int main()

{

    int n,s=0,t,ret;

    scanf("%d",&n);

    for(t=2;t<=n;t++)

    {

        ret=pow(2,t);

        if(sushu(ret-1))

        {

            s=1;//设置标志判断是否有梅森数。

            printf("%d\n",ret-1);

        }

    }//循环找出梅森数

    if(s==0)

       printf("None");

    return 0;

P2000梅森数一个特殊的数学概念,指的是形如\(2^p - 1\)的素数,其中\(p\)也是质数。梅森数以其发现者罗伯特·梅森命名,它们的寻找是数论中的一个重要课题。在C语言中,由于梅森数的数量非常稀少,直接计算出所有梅森数并不是一个高效的实践,特别是对于较大的\(p\)值。 如果你想在C语言中处理梅森数,通常会涉及到编写算法来测试这些候选数是否为梅森数。这通常涉及幂运算和模运算的高效利用。以下是一个简单的C程序片段,用于测试一个给定的数是否可能是梅森数: ```c #include <stdio.h> #include <stdbool.h> // 基于Miller-Rabin素数检验的函数,不是绝对精确的判断,但概率高 bool is_probably_mersenne(int n) { if (n <= 1 || n % 2 == 0) return false; int d = n - 1; // d should be a power of 2 while (!(d & 1)) d >>= 1; // d now is the largest odd factor of n-1 for (int i = 0; i < 50; i++) { // 选择足够次数的测试,增加可信度 int a = rand() % (n - 1) + 1; // 选择随机的a作为测试因子 bool witness = pow(a, d) % n == 1; // 初步测试 if (!witness) { for (int r = 1; r < d; r++) { witness = !witness; witness = witness * pow(a, 2 * r) % n; if (witness == n - 1) break; // 如果找到一个r使得a^(2^r) mod n = n-1, 则可能不是梅森数 } } if (witness) return false; // 如果所有测试都未否定,那么n可能是个梅森数 } return true; } int main() { int candidate = 2047; // 比如测试一个数 if (is_probably_mersenne(candidate)) { printf("%d可能是梅森数.\n", candidate); } else { printf("%d不可能是梅森数.\n", candidate); } return 0; } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值