コラッツ予想がとけたらいいな2

自分の考察を書いていきます。

Pythonとともに

・弱いゴールドバッハ予想 - Wikipedia を見る。

5 より大きい奇数は 3 個の素数の和で表せる。 という予想

        • -

2002年、廖明哲と王天沢は e^3100 ≒ 2*10^1346 より大きい奇数については弱いゴールドバッハ予想が成り立つことを証明した。

コンピュータで解けるんじゃね?(あさはか

Pythonを初めてさわってみる。
Pythonでやったこと
素数判定
・並列計算
・メモ化

ソースコード

#coding:utf8
import random #乱数を使いたい
from multiprocessing import Pool #並列処理したい
import time #時間計測したい
import itertools #直積を使いたい
import sys #システムパラメータをとる
import math #logを使いたい

def wgold(emin, emax):

    tm1 = time.clock()

    print "%d %d" % (emin, emax)

    #素数のリストを作る
    primes = filter( is_prime, range(3, 100000, 2) )

    #偶数のリストを作る
    po = Pool()
    func_args = []
    for even in range(emin, emax+1, 2):
        func_args.append( (chkp, even, primes) )
    print "--chkp start--"
    rst = po.map(argwrapper, func_args)
    print "--chkp end--"

    tm2 = time.clock()

    print tm2-tm1


def chkp(even, primes):
    if even%100000==0: print even

    for p in primes:
        if even-p < 0:
            print "counter exsample %d" % even
            return True
        else:
            pass
        if is_prime(even-p) == True:
            return False
        else:
            pass
    print "counter exsample %d" % even
    return True

・実行結果