さて，下請け関数をコンパイラ言語で書いたらどうなるかということで，私が一番とっつきやすい FORTRAN（雀百まで踊り忘れず）を使って，以下のようなプログラムを書いた。 function is_prime(x) implicit none integer(8) mx, x, i logical is_prime if (x …

さて，前の記事で Python による素数判定プログラムを評価したが，今回は Python プログラムを Julia に移植した結果を書く。 まず，2 を特別にして，あとは奇数で試し割りをするプログラムである。 Julia は色々な変数タイプが使えるが，符号付き 64 ビット…

原作者が最初に書いた関数。本人も認めるとおり，速度的には超マズイ。 def is_prime000(x: int) -> bool: if x <= 1: return False for i in range(2, x): if x % i == 0: return False return True 試し割りはその数の平方根までで十分。そこまでの数で割…

ツェラーの公式を利用して列挙する def Zeller(y, m, d): w = ['Sat', 'Sun', 'Mon', 'Tue', 'Wed', 'Thu', 'Fri'] if m == 1 or m ==2: m += 12 y -= 1 C = y // 100 Y = y % 100 return w[(d + 26*(m+1)//10 + Y + Y//4 -2*C + C // 4) % 7] for y in rang…