c++ 梅森數(shù)源碼示例解析
題目:
要求:按從小到大的順序輸出所有不超過2……n−1的梅森數(shù),每行一個。如果完全沒有,則輸出“None”。
別人的例子
#include <math.h>
int main() {
int n = 0, m = 0, e = 0,h=0;
int i = 0;
scanf("%d", &n);
int a = (int)pow(2, n) - 1;//最大數(shù)
for (i = 2; i < a; i++) {//數(shù)字挨個增加
m = 0;
for (e = 2; e <= sqrt(i); e++) {//判斷是不是為素數(shù)
if (i % e == 0) {
m++;
break;
}
}
if (m == 0) {
for (e = 1; e < n; e++) {
if (2.0 == pow(i + 1, 1.0 / (1.0 * e))) {//判斷是否滿足另一個條件
printf("%d\n", i);
h++;
break;
}
}
}
}
if (h == 0) {
printf("None");
}
return 0;
}
發(fā)現(xiàn)報錯: 分析錯誤: 我感覺可能是求梅森數(shù)的思路不對。
1.for循環(huán)里的i是否需要在循環(huán)外聲明。
2.判斷i是否為素數(shù)的方式可以改為使用線性篩法。
3.在判斷梅森數(shù)時,使用pow函數(shù)求2的冪次方會影響精度,建議使用位運算進行優(yōu)化。
4.如果完全沒有梅森數(shù),應該在循環(huán)外輸出“None”。
我的代碼
#include <stdio.h>
#include <math.h>
int main() {
int n = 0, p = 0, flag = 0;
scanf("%d", &n);
for (p = 2; p <n; p++) {
int i = 2, m = 2; //m為當前計算出的2^p-1的值
//線性篩法求2^p-1是否為素數(shù)
while (i <= sqrt(m)) {
if (m % i == 0) {
break; //不是素數(shù),跳出循環(huán)
}
i++;
}
//如果i>sqrt(m)說明m是素數(shù)
if (i > sqrt(m)) {
printf("%d\n", m);
flag = 1; //標記是否有梅森數(shù)
}
m = (m << 1) + 1; //計算下一個2^p-1的值
}
if (!flag) {
printf("None");
}
return 0;
}
顯示部分正確修改過的代碼
#include <iostream>
#include <cmath>
using namespace std;
int main()
{
int n;
cin >> n;
bool flag = false; // 記錄是否有梅森數(shù)
for (int i = 2; i < n; i++)
{
bool isPrime = true; // 記錄i是否為素數(shù)
for (int j = 2; j <= sqrt(i); j++)
{
if (i % j == 0)
{
isPrime = false;
break;
}
}
if (isPrime) // 如果i為素數(shù)
{
int m = i; // 記錄i的值
int sum = 0; // 記錄i的位數(shù)和
while (m > 0)
{
sum += m % 10;
m /= 10;
}
if (sum == i) // 如果i的位數(shù)和等于i本身,則為梅森數(shù)
{
cout << i << endl;
flag = true; // 有梅森數(shù)
}
}
}
if (!flag) // 如果沒有梅森數(shù)
{
cout << "None" << endl;
}
return 0;
}以上就是c++ 梅森數(shù)源碼示例解析的詳細內(nèi)容,更多關(guān)于c++ 梅森數(shù)源的資料請關(guān)注腳本之家其它相關(guān)文章!
相關(guān)文章
C++設(shè)計與實現(xiàn)ORM系統(tǒng)實例詳解
這篇文章主要為大家介紹了C++設(shè)計與實現(xiàn)ORM系統(tǒng)實例詳解,有需要的朋友可以借鑒參考下,希望能夠有所幫助,祝大家多多進步,早日升職加薪2022-09-09
C++ 中靜態(tài)成員函數(shù)與非靜態(tài)成員函數(shù)的區(qū)別
這篇文章主要介紹了C++ 中靜態(tài)成員函數(shù)與非靜態(tài)成員函數(shù)的區(qū)別的相關(guān)資料,需要的朋友可以參考下2017-05-05

