【寻找第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!!!

 

 

 

阅读更多

没有更多推荐了,返回首页