The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A217715 Numbers equal to the Euler totient function of their arithmetic derivative: k = phi(k'). 1
 16, 30, 32, 48, 54, 64, 80, 120, 176, 288, 368, 432, 464, 656, 848, 858, 864, 1328, 1424, 1722, 1808, 1944, 2096, 2768, 2864, 2916, 3056, 3728, 3824, 4016, 4496, 4688, 5744, 5832, 6704, 6896, 7088, 7856, 8144, 9488, 10256, 10448, 10544, 10928, 11504, 11888 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,1 COMMENTS If p is a Sophie Germain prime (A005384) then m = 16*p is a term. Indeed: m' = (16*p)' = 32*p + 16 = 16*(2*p + 1) and phi(m') = phi(32*p + 16) = phi(16*(2*p + 1)) = 8*phi(2*p + 1) = 8*2*p = m for odd p. If p = 2 then m = 16*2 = 32 is a term. - Marius A. Burtea, Apr 10 2022 LINKS Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..200 from Paolo P. Lava) EXAMPLE For k=368, k'=752 and phi(752)=368. MAPLE with(numtheory); A217715:= proc(q) local n, p; for n from 1 to q do if phi(n*add(op(2, p)/op(1, p), p=ifactors(n)[2]))=n then print(n); fi; od; end: A217715(10^6); MATHEMATICA aQ[1]=1; aQ[n_] := EulerPhi[n * Total[#2/#1 & @@@ FactorInteger[n]]] == n; Select[Range[10000], aQ] (* Amiram Eldar, Jul 11 2019 *) PROG (Magma) f:=func; [k:k in [2..12000]| k eq EulerPhi(Floor(f(k)))]; // Marius A. Burtea, Apr 09 2022 CROSSREFS Cf. A000010, A003415, A005384, A166374. Sequence in context: A372980 A228664 A359935 * A064634 A255265 A183372 Adjacent sequences: A217712 A217713 A217714 * A217716 A217717 A217718 KEYWORD nonn AUTHOR Paolo P. Lava, Mar 21 2013 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified September 16 07:48 EDT 2024. Contains 375959 sequences. (Running on oeis4.)