技術メモ

役に立てる技術的な何か、時々自分用の覚書。幅広く色々なことに興味があります。

2017-06-22から1日間の記事一覧

pythonで完全数と友愛数と婚約数を求める

以前の記事の続きとして、せっかくなので友愛数・完全数・婚約数を求めてみたい。 約数関数(再掲) divisor(n)として約数関数を使うので、再掲しておく。 sympyを使う。 import sympy def divisor(n): factors = sympy.factorint(n) rst = 1 for i,j in fac…

sympyで素因数分解ができるからオイラー関数と約数関数を書いてみた

pythonのライブラリの一つsympyを使えば、簡単に素因数分解ができるということを知った。 import sympy sympy.factorint(1000) #{2: 3, 5: 3} ちなみに、因数分解も簡単にできる(◎.◎)!! 凄い。 import sympy x = sympy.Symbol('x') eqn = x**2 - 3*x + 2 pri…