梅森素数与四素数：算法实现与优化策略

1. 梅森素数与四素数：算法实现与优化策略

在计算机科学和数学的交叉领域，素数研究一直是个经典课题。今天我将分享两个有趣的素数相关问题：梅森素数的查找和四素数的统计。这两个问题看似简单，但涉及到了算法选择、性能优化等核心编程技巧。

2. 梅森素数：概念与实现

2.1 梅森素数的数学定义

梅森数是指形如2^p-1的一类数，其中指数p是素数。当这个梅森数本身也是素数时，我们称它为梅森素数。例如：

• M(2) = 2²-1 = 3（是素数）
• M(3) = 2³-1 = 7（是素数）
• M(5) = 2⁵-1 = 31（是素数）

2.2 算法实现解析

题目要求输入一个长整型数n，输出不大于n的所有梅森素数。参考代码采用了以下策略：

cpp复制#include<bits/stdc++.h>
using namespace std;

// 判断是否为素数的函数
int fun(long n){
    if(n <= 1) return 0;
    for(int i = 2; i * i <= n; i ++){
        if(n % i == 0){
            return 0;
        }
    }
    return 1;
}

int main(){
    long n;
    cin>>n;
    for(int i = 2; i < n; i ++){
        long long p = pow(2, i) - 1; 
        if(p > n){
            break;
        }
        if(fun(i) && fun(p)) {
            cout<<"M("<<i<<")="<<p<<endl;
        }
    }
    return 0;
}

2.2.1 关键点解析

1. 素数判断函数fun()：采用试除法，检查从2到√n的所有整数是否能整除n
2. 主循环逻辑：
   • 遍历可能的指数i（从2开始）
   • 计算梅森数p = 2^i - 1
   • 检查i和p是否都是素数
   • 当p超过n时提前终止循环

注意：这里使用pow(2