使用 sagemath 的哥德巴赫图答案

【问题标题】：Goldbach graph using sagemath使用 sagemath 的哥德巴赫图
【发布时间】：2020-12-31 11:02:43
【问题描述】：

我正在学习 SageMath（使用 Python 3）并玩弄哥德巴赫猜想。

我写了这个函数（它有效！）：

def Goldbach(n):
    if n % 2 != 0 or n <= 2:
        show("No és parell")
    else:
        for i in srange(n):
            if is_prime(i):
                for j in srange(n):
                    if is_prime(j) and i + j == n:                        
                        a = [i, j]
                        show(a)
    return

现在我正在尝试（不知道）做以下情节：

用r(2k)表示2k的Goldbach分区数，如果k > 1，猜想肯定r(2k) > 0。

我必须做一个点图(k, r(2k))，k > 2。我该怎么办？

【问题讨论】：

标签： graph sage goldbach-conjecture

【解决方案1】：

首先，让我们在Sage 中得到一些更好的实现例程计算 p + q = 2k 的解的数量 r(K)（对于 K > 2 一些偶数整数）， p, q 素数。

当解决方案 (p, q) 和 (q, p) 不同时，我们计算它们。

def r(K):
    if K not in ZZ or K <= 2 or K % 2:
        return None
    if K == 4:
        return 1
    count = 0
    for p in primes(3, K):
        for q in primes(3, K + 1 - p):
            if p + q == K:
                count += 1
    return count

goldbach_points = [(K, r(K)) for K in range(4, 100,2)]
show(points(goldbach_points))

这给出了：

【讨论】：