OFFSET
1,1
COMMENTS
A088826 is a subset of this sequence: it lists the fixed points of the transform n -> sigma(n)-phi(n)-n.
a(29) > 10^9. - Michel Marcus, Mar 06 2014
EXAMPLE
Fixed points: 12, 42, 1242, 6137440, 1385119360, ...
sigma(2010) = 4896, phi(2010) = 528 and 4896 - 528 - 2010 = 2358.
sigma(2358) = 5148, phi(2358) = 780 and 5148 - 780 - 2358 = 2010.
MAPLE
with(numtheory); P:=proc(q)local a, n;
for n from 1 to q do a:=sigma(n)-phi(n)-n;
if a>0 and sigma(a)-phi(a)-a=n then print(n);
fi; od; end: P(10^6);
PROG
(PARI) isok(n) = !isprime(n) && ((x = (sigma(n) - eulerphi(n) - n)) > 0) && (n == (sigma(x) - eulerphi(x) - x)); \\ Michel Marcus, Mar 06 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Feb 20 2014
EXTENSIONS
a(8)-a(28) from Michel Marcus, Mar 06 2014
STATUS
approved