【寻找第n个默尼森数】利用python找出第n个默尼森数

经典程序设计问题：找第n个默尼森数。P是素数且M也是素数，并且满足等式M=2**P-1，则称M为默尼森数。例如，P=5，M=2**P-1=31，5和31都是素数，因此31是默尼森数。（31是第3个默尼森数）

代码格式如下：

def prime(num):

...

def monisen(no):

… …

return xxx

# -*- coding:utf-8 -*-
"""
Finding No.? Monisen Number

author：WangSheng
"""

# sqrt：平方根运算
# pow ：幂运算
from math import sqrt,pow

def prime(num):
    """检测num是不是质数"""
    k = sqrt(num)
    for i in range(2,int(k+1)):
        if num % i == 0:
            return False
    
    return True

def monisen(no):
    """找出第no个莫尼森数"""
    n = 0
    num = 2
    while n < no:
        m = pow(2,num) - 1
        if prime(num) == True and prime(m) == True:
            # 只有num和m都为质数时,n才会加一,即n是莫尼森数的序号
            n += 1
        num += 1
    return int(m),num-1

# 输出前五个莫尼森数M 以及对应的质数P
for i in range(1,6):
    print(monisen(i))

 

公众号：Wang的胜宴

关注和我一起学python！！！