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

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