1,1

Perfect numbers (A000396) are a subsequence. Perfect numbers are defined as numbers equal to the sum of their proper divisors. The numbers in this sequence are not perfect, but are "permutably perfect" since they are equal to a permutation of the sum of their proper divisors.

Note a(79)=195921 and a(80)=199251 and these are both permutations of each other. The same is true of a(90)=226944 and a(93)=244692.

V. Raman, Table of n, a(n) for n = 1..10000

604 is a term because the divisors of 604 are {1, 2, 4, 151, 302, 604} and 1 + 2 + 4 + 151 + 302 = 460.

16012 is included because it has the proper divisors 1, 2, 4, 4003, 8006. These divisors sum to 12016, which is a rearrangement of the digits of 16012.

(PARI) is_A085844(n)={vecsort(digits(n))==vecsort(digits(sigma(n)-n))} \\ M. F. Hasler, Feb 24 2014

(Python)

from sympy import divisors

def ok(n):

divs = divisors(n)[:-1]

return sorted(str(n)) == sorted(str(sum(divs)))

print(list(filter(ok, range(1, 57000)))) # Michael S. Branicky, Jun 28 2021

Cf. A000396.

Sequence in context: A332465 A335290 A173360 * A325638 A331752 A083387

Adjacent sequences: A085841 A085842 A085843 * A085845 A085846 A085847

base,nonn

Chuck Seggelin, Jul 05 2003; corrected Jul 07 2003

Typo in example corrected by Alan Frank, Jan 20 2010

approved