|
%I #20 Feb 19 2022 11:17:22
%S 6,28,120,496,672,8128,35640,199584,523776,2142720,12999168,33550336,
%T 459818240,1476304896,2836487808,6039429120,6399679104,8589869056,
%U 36639203328,51001180160,137438691328,266653296000,658470384000,2732372020224,6164773235712
%N Numbers m such that the sum of the first k divisors of m, for some k, is equal to the sum of its other divisors.
%C Includes all the perfect numbers (A000396), for them k = d(m) - 1 and the even 3-perfect numbers (A005820), for them k = d(m) - 2 (where d(m) = A000005(m) is the number of divisors of m).
%C 36639203328 is also a term.
%e 6 is a term since its set of divisors, {1, 2, 3, 6}, can be partitioned into two disjoint sets with equal sum, {1, 2, 3} and {6}, such that the first 3 divisors are in the first set.
%t Select[Range[200000], MemberQ[Accumulate[(d = Divisors[#])], (Plus @@ d)/2] &]
%o (Python)
%o from itertools import count, islice
%o from sympy import divisors
%o def A334410_gen(startvalue=1): # generator of terms >= startvalue
%o for n in count(max(startvalue,1)):
%o ds = divisors(n)
%o s = sum(ds)
%o if s % 2 == 0 and any(2*a==s for a in accumulate(ds)):
%o yield n
%o A334410_list = list(islice(A334410_gen(),7)) # _Chai Wah Wu_, Feb 19 2022
%Y Subsequence of A083207.
%Y Cf. A000005, A000396, A005820, A334409.
%K nonn,more
%O 1,1
%A _Amiram Eldar_, Apr 27 2020
%E a(19)-a(25) from _Giovanni Resta_, May 08 2020
|
