OFFSET
1,1
COMMENTS
a(15) > 10^10. - Amiram Eldar, Jan 30 2019
1.25*10^11 < a(20) <= 13421772796875. - Giovanni Resta, Apr 15 2019
EXAMPLE
(8 - phi(8))’ = (8 - 4)’ = 4’ = 4 = phi(8) ;
(189 - phi(189))’ = (189 - 108)’ = 81’ = 108 = phi(189).
MAPLE
with(numtheory): P:= proc(q) local n, p; for n from 1 to q do
if (n-phi(n))*add(op(2, p)/op(1, p), p=ifactors(n-phi(n))[2])=phi(n) then print(n);
fi; od; end: P(10^9);
MATHEMATICA
Select[Range[10^6], Function[m, Function[k, If[Abs@ k < 2, 0, k Total[#2/#1 & @@@ FactorInteger[Abs@ k]]]][# - m] == m]@ EulerPhi@ # &] (* Michael De Vlieger, May 02 2016, after Michael Somos at A003415 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Paolo P. Lava, May 02 2016
EXTENSIONS
a(14) from Amiram Eldar, Jan 30 2019
a(15)-a(19) from Giovanni Resta, Apr 15 2019
STATUS
approved