技術メモ

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