OFFSET
1,2
EXAMPLE
For n = 6583 we have tau(n) = 4, n’ = 256, n’’= 1024 and 1024 = 4 * 256.
MAPLE
with(numtheory); P:= proc(q) local a1, a2, n, p;
for n from 1 to q do a1:=n*add(op(2, p)/op(1, p), p=ifactors(n)[2]); ;
a2:=a1*add(op(2, p)/op(1, p), p=ifactors(a1)[2]);
if a2=tau(n)*a1 then print(n); fi; od; end: P(10^6);
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Oct 25 2013
STATUS
approved