实验1--哥德巴赫猜想

最新推荐文章于 2021-11-25 21:53:40 发布

(☆_☆)々@

最新推荐文章于 2021-11-25 21:53:40 发布

阅读量859

点赞数

分类专栏：算法作业记录

本文链接：https://blog.csdn.net/qq_44242419/article/details/109395472

版权

算法作业记录专栏收录该内容

2 篇文章 0 订阅

订阅专栏

2020.10.30星期五

实验一：算法时间性能的经验分析

验证哥德巴赫猜想。统计其关键语句的执行次数。并绘制规模-执行次数散点图。
（Goldbach Conjecture）猜想：即任一大于2的偶数都可写成两个质数之和。请验证这个对于
较大的偶数都是成立的。
算法：goldbach(n)
描述：算法验证对于小于等于n的偶数，歌德巴赫猜想都是成立的。注意，goldbach(200)并非
验证200可以拆分成两个质数的和，而是验证对于所有≤n的偶数，都可以拆分成两个质数的
和。
输入：整数n
输出：1表示成立；0表示猜想有误
注1：用excel绘制散点图，excel可以打开csv格式的文件。
注2：csv(comma seperated values)是如下格式的文件。如下所示，n是第一列，times是第二
列。你可以用python直接输出该文本文件。例如，400,1221，表示为了验证歌德巴赫猜想对于小
于等于400的偶数都成立，其关键语句执行了1221次。

一、结果如下

二、源代码

import numpy as np
import csv


# 验证一个正整数是否是素数
def is_prime(num):
    count = 0
    if num < 1 or num % 2 == 0:
        return count, False
    for i in range(2, num, 2):
        if i <= np.sqrt(num):
            if num % i == 0:
                return count, False
            count = count + 1
    return count, True


def goldbach_conjecture(n):
    count = 0
    for i in range(2, n):
        a_prime = is_prime(i)
        b_prime = is_prime(n-i)
        count = count + a_prime[0] + b_prime[0]
        if a_prime[1] and b_prime[1]:
            return [i, n-i], count, True
    return None, count, False


file = open('goldbachConjecture.csv', 'w', newline='', encoding='utf-8')
csv_writer = csv.writer(file)
csv_writer.writerow(['n', 'times'])
for i in range(200, 5001, 50):
    count = 0
    for j in range(2, i+1, 2):
        count = count + goldbach_conjecture(j)[1]
    csv_writer.writerow([i, count])

三、数据分析：

数据分析图片

(☆_☆)々@

关注

0
点赞
踩
1

收藏

觉得还不错? 一键收藏
0
评论
实验1--哥德巴赫猜想

2020.10.30星期五实验一：算法时间性能的经验分析验证哥德巴赫猜想。统计其关键语句的执行次数。并绘制规模-执行次数散点图。（Goldbach Conjecture）猜想：即任一大于2的偶数都可写成两个质数之和。请验证这个对于较大的偶数都是成立的。算法：goldbach(n)描述：算法验证对于小于等于n的偶数，歌德巴赫猜想都是成立的。注意，goldbach(200)并非验证200可以拆分成两个质数的和，而是验证对于所有≤n的偶数，都可以拆分成两个质数的和。输入：整数n输出：1表示成立
复制链接

扫一扫