OFFSET
1,1
LINKS
David A. Corneth, Table of n, a(n) for n = 1..10000
EXAMPLE
The distinct primes dividing 24 are 2 and 3, since 24 is factored as 2^3 *3^1. 2 + 3 = 5 is a divisor of 24 + 1 = 25. So 24 is a term of this sequence.
MAPLE
with(numtheory): a:= proc(n) local f: f:=factorset(n); if `mod`(n+1, add(i, i=f))=0 then n end if end proc: seq(a(n), n=2..1200); # Emeric Deutsch, Aug 14 2008
MATHEMATICA
Select[Range[2, 1500], Divisible[#+1, Total[FactorInteger[#][[All, 1]]]]&] (* Harvey P. Dale, Aug 27 2022 *)
PROG
(PARI) is(n) = n > 1 && (n + 1) % vecsum(factor(n)[, 1]) == 0 \\ David A. Corneth, Mar 10 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Aug 07 2008
EXTENSIONS
More terms from Emeric Deutsch, Aug 14 2008
More terms from Max Alekseyev, Mar 10 2009
STATUS
approved