OFFSET
1,1
EXAMPLE
n=432, n'=1296, n''=4320 and 4320=1296*3+432.
MAPLE
with(numtheory);
A213241:= proc(i)
local a, b, n, p, pfs;
for n from 1 to i do
pfs:=ifactors(n)[2]; a:=n*add(op(2, p)/op(1, p), p=pfs);
pfs:=ifactors(a)[2]; b:=a*add(op(2, p)/op(1, p), p=pfs);
pfs:=ifactors(b)[2]; c:=b*add(op(2, p)/op(1, p), p=pfs);
if a>0 then if (b mod a)=n then print(n); fi; fi;
od;
end:
A213241(10000000);
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Jun 12 2012
STATUS
approved