-p}+\cdots+2^{3p}+2^{2p}+2^p+1)$となり、メルセンヌ素数に$2^{k-1}$を掛けた値が完全数になるのは次の計算ができるためです。整数nが完全数であるためnの約数の合計:Sとすると...2023.09.21Pythonある整数の約数に対するオイラー関数の和がある整数になるPythonでオイラー関数を計算するオイラー関数は、ある整数に対し、その整数より小さくかつ互いに素(最大公約数=1)となる整数の個数を計算します。オイラー関数の前提として、Pythonでは次の通り最大公約数を計算することができます。#最大公約数を計算す...2021.08.31PythonPythonで小数型の数値を2進数、さらには16進数に変換するPythonで小数型の数値を2進数に変換するために関数を作ってきましたが、この結果が本当に正しいかは全くわかりません。そこで、小数を2進数に変換する関数やメソッドを探してみましたが現時点で見つかっていません。そのかわり、16進数に変換すること、次...2021.06.27Python2進数での丸め処理をする関数Pythonで浮動小数点の計算をしていて、桁あふれの処理をするとき、あふれたbitを丸める必要があります。このときに丸める方法ための関数を作成します。このことにより、本当に細かい話ですが小数点以下のおかしな振る舞いがなぜ起ってしまうのかがわかりま...2021.06.22Pythonバーゼル問題と円周率arctanのほかにも、バーゼル問題やさらには素数から円周率を計算する方法があります。\dfrac{1}{n^{2}}=1+\frac{1}{2^2}+\frac{1}{3^2}+\frac{1}{4^2}+\frac{1}{5^2}+\frac{1}{6^2}\cdots=\frac{\pi^2}{...2021.06.14Pythonマチンの公式でどれくらいの精度で円周率を計算できるか?次のようなグレゴリー級数で、$\pi/4$を計算することができます。\sum_{n=0}^{\infty}\frac{(-1)^n}{2n+1}x^{2n+1}=x-\frac{x^3}{3}+\frac{x^5}{5}-\frac{x^7}{7}+\frac{x^9}{9}\cdots$そこで、この計算で100項までループした...2021.06.13Python円周率10万桁これから円周率の精度計算をするにあたり、正確な円周率10万桁をネットで検索してみました。調べて2つのサイトが使えそうです。\pi/4=1$から円周率を計算する微分、積分などの知識を使って円周率を計算します。はじめに、単純なところから。\frac{\pi}{4}=1$$\pi/4$は45度なので、$tan=1$なのはすぐにわかります。そこでtanの逆数arctanを考えます。1$ここで$acrt...2021.06.11PythonPythonでネイピア(Napier)数を計算するPythonでネイピア数を計算します。いってもSymPyモジュールを使えば、好きな桁数だけ計算できます。sympyprint(sympy.E.evalf(100))2.718281828459045235360287471352662497757247093699959574966967627724...2021.06.10PythonPythonで平方根を計算するPythonを使って平方根を計算してみます。いっても、x**(0.5)とかsqrt関数で計算しようという話ではありません。頭の体操で、手計算に近い方法で計算しようという趣旨です。