OFFSET
1,2
COMMENTS
Additional terms include 537444352, 2148958208, 137454419968, 35184644718592, 9007202811510784. Are there any terms > 1 not of the form 2^k*p where p is prime and k>0? - David A. Corneth, Jan 08 2020
Terms not of the form 2^k*p do exist, for example 2^15*65713*24194197 and 2^19*1739719*2639431. - Giovanni Resta, Jan 08 2020
a(20) > 10^13. - Giovanni Resta, Jan 14 2020
EXAMPLE
The first term that is not 1 or a single-digit prime is obtained by adding the proper divisors of 14 other than 1 (2,7) to its digits (1,4): (2+7) + (1+4) = 14.
The second such term is 52: the proper divisors of 52 other than 1 (2,4,13,26) and its digits (5,2) sum to (2+4+13+26) + (5+2) = 52.
MATHEMATICA
Select[Range[10^7], DivisorSigma[1, #] - # - If[# == 1, 0, 1] + Plus @@ IntegerDigits[#] == # &] (* Amiram Eldar, Jan 12 2020 *)
PROG
(PARI) is(n) = n == sigma(n)-1-if(n>1, n, 0)+sumdigits(n) \\ Rémy Sigrist, Jan 08 2020
CROSSREFS
KEYWORD
nonn,base,more,less
AUTHOR
Joseph E. Marrow, Jan 08 2020
EXTENSIONS
a(14)-a(16) from Rémy Sigrist, Jan 08 2020
a(17)-a(19) from Giovanni Resta, Jan 14 2020
STATUS
approved