OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..94
EXAMPLE
70 is a term since it is a weird number, its aliquot divisors are {1, 2, 5, 7, 10, 14, 35} and 71 = 5 + 7 + 10 + 14 + 35.
MATHEMATICA
q[n_] := Module[{d = Most @ Divisors[n], x, s, c}, If[Plus @@ d <= n, False, s = Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, n + 1}]; c = SeriesCoefficient[s, #] & /@ (n + {0, 1}); c[[1]] == 0 && c[[2]] > 0]]; Select[Range[10000], q]
PROG
(PARI) is(n, d=divisors(n)[^-1], s=vecsum(d))={s>n && !is_A005835(n, d, s) && is_A005835(n+1, d, s)}; \\ using is_A005835() by M. F. Hasler at A005835
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, May 22 2022
STATUS
approved