This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A053211 #21 Mar 13 2023 06:04:55
%S A053211 2,4,3,8,5,9,16,7,32,27,11,25,64,13,81,128,17,49,19,256,23,125,243,29,
%T A053211 31,512,121,37,41,43,1024,729,169,47,343,53,625,59,61,2048,67,289,71,
%U A053211 73,79,2187,361,83,89,4096,97,101,103,107,109,529,113,1331,3125,127
%N A053211 Cototients of consecutive pure powers of primes.
%C A053211 Cototients of prime powers do not remain always prime powers, but are primes if their exponent is 2.
%H A053211 Michael De Vlieger, <a href="/A053211/b053211.txt">Table of n, a(n) for n = 1..16384</a>
%H A053211 Michael De Vlieger, <a href="/A053211/a053211.png">Log log scatterplot of a(n)</a> n = 1..2^20, showing even a(n) in blue, 3 | a(n) in green, and prime a(n) in red, else black.
%F A053211 a(n) = A051953(A025475(n+1)) = cototient(p^k) = p^(k-1).
%e A053211 The 10th pure power of prime (but not a prime) is 81, so a(10) = 81 - EulerPhi(81) = 81 - 54 = 27. For n=p^2, a(n)=p.
%t A053211 Map[# - EulerPhi@ # &, Select[Range[16200], And[! PrimeQ@ #, PrimePowerQ@ #] &]] (* _Michael De Vlieger_, Jun 11 2018 *)
%t A053211 With[{nn = 2^14}, Map[Times @@ Map[#1^(#2 - 1) & @@ FactorInteger[#][[1]]] &, Select[Union@ Flatten@ Table[a^2*b^3, {b, nn^(1/3)}, {a, Sqrt[nn/b^3]}], PrimePowerQ] ] ] (* _Michael De Vlieger_, Mar 11 2023 *)
%Y A053211 Cf. A000010, A051953, A001248, A002618, A036689, A053650, A053191, A053192, A246547.
%K A053211 nonn
%O A053211 1,1
%A A053211 _Labos Elemer_, Mar 03 2000