Pyhonでの完全数とほぼ完全数の判別

Question

前提・実現したいこと

ここに質問の内容を詳しく書いてください。
Python3で複数の整数の完全数とほぼ完全数の判別を行うプログラムを作成しています
入力値は
3
28
16
777
です。3は入力する整数の個数を示しています。

実現したい内容は完全数ならperfect、ほぼ完全数ならnearly、どちらでもなければnoと出力したいです。

発生している問題・エラーメッセージ

[28.0, 14.0, 7.0, 4.0, 2.0, 1.0]
[28.0, 14.0, 7.0, 4.0, 2.0, 1.0, 16.0, 8.0, 4.0, 2.0, 1.0]
[28.0, 14.0, 7.0, 4.0, 2.0, 1.0, 16.0, 8.0, 4.0, 2.0, 1.0, 777.0, 259.0, 111.0, 37.0, 21.0, 7.0, 3.0, 1.0]
以上のような結果が出力され、毎回28の約数等が加算されるためうまくいきません。

該当のソースコード

num = int(input())
yakusu = 0
amari = 0
list = []
total = 0
s = [int(input()) for i in range(num)]
for i in range(0,num):
    for j in range(1,s[i]+1):
        yakusu = s[i] / j
        amari = s[i] % j
        if amari == 0:
            list.append(yakusu)
    print(list)     

試したこと

変数listにも番地をつけようとしたができなかった

補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

Answer 1

以下で、完全数、擬似完全数を判定できます。約数の列挙はsympyを使いました。それぞれ関数の名前をわかりやすく書きましたので、説明は不要と思います。1〜99の判定を例として実行させています。

なお、擬似完全数の判定は、約数の部分集合の列挙にO(n^2)もしくはO(2^n)の計算量になると思いますので、numが大きいと時間がかかります。このあたりは、工夫のしどころだと思います。

import sympy

def divisors_without(num):
	return sympy.divisors(num)[:-1]

def is_perfect(num):
	return num == sum(divisors_without(num))

def is_semiperfect(num):
	for l in enum_sublists(divisors_without(num)):
		if num == sum(l):
			return True
	return False

def enum_sublists(l):
	lists = [[]]
	for i in range(len(l)): 
		orig = lists[:] 
		new = l[i] 
		for j in range(len(lists)): 
			lists[j] = lists[j] + [new] 
		lists = orig + lists
	return lists

for num in range(1, 100):
	if is_perfect(num):
		classified = 'perfect'
	elif is_semiperfect(num):
		classified = 'nearly'
	else:
		classified = 'no'
	print(num, ': ', classified)

Answer 2

以下でどうでしょうか。（ほぼ完全数については考えてません）

• 変数名listはpythonの型名なのでyakusu_lstに変更 & for文の中で毎回リセット
• 辞書dctに対象の数と、約数リスト＆完全数判定を入れ込んでいく
• 2段目のfor文がs[i]+1になっていて、そのままだと自分自身を約数リストに入れてしまうので+1を消した
• yakusuが名前に反して「対象の数を約数で割った値」であり、それを約数リストに追加していたので、appendの引数をint(j)に変更

num = int(input())
s = [int(input()) for i in range(num)]
dct = {}
for i in range(0,num):
    yakusu_lst = []
    typ = "no"
    for j in range(1,s[i]):
        yakusu = s[i] / j
        amari = s[i] % j
        if amari == 0:
            yakusu_lst.append(int(j))
    if sum(yakusu_lst) == s[i]:
        typ = "perfect"
    dct[s[i]] = (yakusu_lst,typ)
print(dct)
# {3: ([1], 'no'), 6: ([1, 2, 3], 'perfect'), 28: ([1, 2, 4, 7, 14], 'perfect'), 777: ([1, 3, 7, 21, 37, 111, 259], 'no')}