OFFSET
1,2
COMMENTS
A246914 gives the primes in this sequence.
LINKS
Martin Møller Skarbiniks Pedersen, Table of n, a(n) for n = 1..1000
EXAMPLE
Number 26 (with sigma(26) = 42) is in sequence because sigma(26+sigma(26)) = sigma(68) = 126 = 3*42.
MAPLE
with(numtheory): A246910:=n->`if`(sigma(n+sigma(n)) = 3*sigma(n), n, NULL): seq(A246910(n), n=1..5000); # Wesley Ivan Hurt, Sep 07 2014
PROG
(Magma) [n:n in[1..10000] | SumOfDivisors(n+SumOfDivisors(n)) eq 3*SumOfDivisors(n)]
(PARI)
for(n=1, 10^4, if(sigma(n+sigma(n))==3*sigma(n), print1(n, ", "))) \\ Derek Orr, Sep 07 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Sep 07 2014
STATUS
approved