OFFSET
1,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000 (first 200 terms from Paolo P. Lava)
EXAMPLE
Arithmetic derivative of 56 is 92 and phi(92) = 44, while phi(56) = 24 and the arithmetic derivative of 24 is 44.
MAPLE
with(numtheory): P:=proc(q) local a, b, c, n, p; for n from 1 to q do
a:=n*add(op(2, p)/op(1, p), p=ifactors(n)[2]); b:=phi(n);
c:=b*add(op(2, p)/op(1, p), p=ifactors(b)[2]);
if phi(a)=c then print(n); fi; od; end: P(10^9);
MATHEMATICA
f[n_] := If[Abs@ n < 2, 0, n Total[#2/#1 & @@@ FactorInteger@ Abs@ n]]; Select[Range@ 12000, EulerPhi@ f@ # == f@ EulerPhi@ # &] (* Michael De Vlieger, Aug 07 2015, after Michael Somos at A003415 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Paolo P. Lava, Aug 06 2015
STATUS
approved