OFFSET
1,1
EXAMPLE
phi(555) = 288 and 288' = 912 = sigma(555);
phi(2691) = 1584 and 1584' = 4368 = sigma(2691).
MAPLE
with(numtheory): P:=proc(q) local a, n, p; for n from 1 to q do
a:=phi(n)*add(op(2, p)/op(1, p), p=ifactors(phi(n))[2]);
if a=sigma(n) then print(n); fi; od; end: P(10^9);
MATHEMATICA
d[1]=0; d[n_] := Total[n / Divide@@@ FactorInteger@ n]; Select[Range[10^5], DivisorSigma[1, #] == d@ EulerPhi@ # &] (* Giovanni Resta, Aug 21 2015 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Paolo P. Lava, Aug 18 2015
STATUS
approved