def parfaits(n):
    L=[]
    for i in range(2,n+1):
        s=0
        for j in range(1,i//2+1):
            if i%j==0:
                s+=j
                if s==i:
                    L.append(i)       
    return L