C++打印前5个完全数程序异常输出2096128，求排查

嘿，我来帮你排查这个问题！首先明确一下：前5个完全数应该是6、28、496、8128、33550336，所以2096128混进来肯定是代码逻辑出了问题。结合你给出的错误数值，我几乎可以确定问题出在梅森素数的判断环节，下面给你拆解细节：

为什么会出现2096128？ #

完全数和梅森素数强相关：如果2^p - 1是梅森素数，那么2^(p-1)*(2^p -1)就是一个完全数。前几个符合条件的p是2、3、5、7、13，对应生成前5个正确的完全数。

而2096128正好是p=11时的计算结果：2^(10)*(2^11 -1) = 1024*2047 = 2096128。但问题在于，2^11 -1 = 2047是合数（2047=23*89），根本不是梅森素数！所以你的代码肯定是错误地把2047判定成了素数，才生成了这个非完全数。

常见的错误点排查 #

你可以对照自己的代码检查以下几个地方：

• 素数判断函数的逻辑漏洞：比如你的素数判断是否只检查到了sqrt(num)之前的数？如果循环只跑到num/2或者中途提前终止，就可能漏掉23、89这样的因子。举个反例：如果你的素数判断函数是这样的，就会误判2047为素数：

  bool isPrime(int num) {
      if (num <=1) return false;
      for(int i=2; i<num/2; i++){ // 这里应该改成i*i <= num，否则会漏掉大因子
          if(num%i ==0) return false;
      }
      return true;
  }
  

• 整数溢出问题：如果你的代码用int存储梅森数，当p增大时可能会溢出，但2047还在32位int的范围内，所以这个可能性较低，但如果是计算过程中溢出导致数值错乱，也可能出现类似问题，建议用long long存储大数。
• 真因子和的计算错误：如果你的代码不是用梅森素数生成完全数，而是直接计算每个数的真因子和，那可能是在求和时把自身算进去了，或者漏掉了某些因子。比如循环到n而不是n-1，或者没有处理平方数的因子重复问题。

修复示例 #

给你一个简单的正确实现参考，重点是素数判断的正确性：

#include <iostream>
using namespace std;

bool isPrime(long long num) {
    if (num <= 1) return false;
    if (num == 2) return true;
    if (num % 2 == 0) return false;
    // 只检查到sqrt(num)，且跳过偶数
    for (long long i = 3; i * i <= num; i += 2) {
        if (num % i == 0) return false;
    }
    return true;
}

void findFirst5PerfectNumbers() {
    int found = 0;
    int p = 2;
    while (found < 5) {
        long long mersenne = (1LL << p) - 1; // 用1LL避免移位溢出
        if (isPrime(mersenne)) {
            long long perfect = (1LL << (p-1)) * mersenne;
            cout << "第" << found+1 << "个完全数：" << perfect << endl;
            found++;
        }
        p++;
    }
}

int main() {
    findFirst5PerfectNumbers();
    return 0;
}

这个代码会正确输出前5个完全数，不会出现2096128。

内容的提问来源于stack exchange，提问作者Catz