欧拉项目（python练习）problem 46

最新推荐文章于 2024-02-13 19:39:50 发布

lu123star456

最新推荐文章于 2024-02-13 19:39:50 发布

阅读量874

点赞数

分类专栏：欧拉项目文章标签： python 欧拉项目编程

本文链接：https://blog.csdn.net/u014399656/article/details/27235259

版权

欧拉项目专栏收录该内容

4 篇文章 0 订阅

订阅专栏

本文探讨了Python在欧拉项目中的应用，具体聚焦于问题46，通过`decide_whether_prime`函数找寻特定素数，最终得出运行结果5777。

摘要由CSDN通过智能技术生成

# -*- coding:utf8 -*-

import sys
from math import sqrt

import decide_whether_prime

#题目设置
'''Christian Goldbach 提出每个奇合数都可以写作一个质数与一个平方数的二倍之和：

9 = 7 + 2×12
15 = 7 + 2×22
21 = 3 + 2×32
25 = 7 + 2×32
27 = 19 + 2×22
33 = 31 + 2×12

但是这个推测是错误的。

最小的不能写作一个质数与一个平方数的二倍之和的奇合数是多少？'''

#解题思路：首先创建一个方法，判断一个奇合数能否被以题目中的方式分解，之后从9开始，对所有的奇合数进行判断，就利用刚才创建的方法，直到找到这个不能被分解的奇合数为止

class MinimumUndividedNumber():                                               #定义了一个找到最小的不能分解的奇合数的类

    def __init__(self):
        pass
      
    def judge_whether_divided(self , odd_composite_number):                   #创建了一个方法，判断奇合数是否能被分解，能则返回1，否则返回0
        can_be_divided = 0
        for number in range(3 , odd_composite_number , 2):
            if decide_whether_prime.decide_whether_prime(number) == 1:
                if sqrt((odd_composite_number - number) / 2) % 1 == 0:
                    can_be_divided = 1
                    break
        return can_be_divided
                    
    def find_minimum_number(self):                                            #主方法，找到最小的不能分解的奇合数，从9开始，逐个判断
        odd_number = 9
        while(True):
            if decide_whether_prime.decide_whether_prime(odd_number) == 0:
                if self.judge_whether_divided(odd_number) == 0:
                    return odd_number
            odd_number = odd_number + 2
            
if __name__ == "__main__":
    minimum_undivided_number = MinimumUndividedNumber()
    print minimum_undivided_number.find_minimum_number()

decide_whether_prime:

# -*- coding:utf8 -*-

import math

def decide_whether_prime(number):
    is_prime = 1
    if number == 1:
        is_prime = 0
    if number == 2 :
        pass
    else :
        for divisor in range(2 , int(math.sqrt(number)) + 1):
            if number%divisor == 0:
                is_prime = 0
                break
    return is_prime

运行结果：5777