OFFSET
1,1
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..300
EXAMPLE
phi(57)=36. 57'=22 and sigma(22)=36
phi(1012505)=725760. 1012505'=310156 and sigma(310156)=725760
MAPLE
with(numtheory);
P:=proc(i)
local f, n, p, pfs;
for n from 1 by 1 to i do
pfs:=ifactors(n)[2];
f:=n*add(op(2, p)/op(1, p), p=pfs);
if phi(n)=sigma(f) then print(n); fi;
od;
end:
P(1000000)
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, May 17 2011
STATUS
approved