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 A332737 #5 Feb 22 2020 08:33:39 %S A332737 289,2329,4369,4913,18769,21331,35209,66049,128881,197143,258121, %T A332737 281929,516961,739903,971203,1762249,1942663,2070721,2898703,2952673, %U A332737 3820819,4142881,8288641,16773619,16843009,16974593,20229241,21762361,32472241,132575071,187903693 %N A332737 Composite terms of A181659, where the sum of the iterated totient function attains a record. %C A332737 Most of the terms of A181659 are primes. Out of the first 10^4 terms of A181659 only 28 are composites. %C A332737 The indices of the terms of this sequence in A181659 are 30, 73, 93, 99, 154, 161, 191, 236, 286, 316, ... %C A332737 The corresponding record values (terms of A126106) are 527, 4223, 8191, 8847, 35527, 39423, 67583, 131327, 246869, 376559, 493739, 550911, 1009981, 1466879, 1884671, 3442687, 3819519, 4089245, 5707263, 5791743, 7444991, 8178491, 16464253, 33260031, 33554431, 33718527, 39989247, 42809067, 63932219, 263382015, 372697723. %t A332737 s[n_] := Plus @@ FixedPointList[EulerPhi, n] - n - 1; seq={}; smax = 1; Do[s1 =s[n]; If[s1 >smax, smax = s1; If[CompositeQ[n], AppendTo[seq, n]]], {n, 1, 5000}]; seq %Y A332737 Cf. A092693, A126106, A181659. %K A332737 nonn %O A332737 1,1 %A A332737 _Amiram Eldar_, Feb 21 2020