cp's OEIS Frontend

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.

A330400 Numbers at which the sum of the iterated unitary totient function (A329153) attains a record.

This page as a plain text file.
%I A330400 #16 Feb 26 2020 00:30:21
%S A330400 1,2,3,4,5,8,9,16,17,29,32,41,45,53,64,73,83,85,101,113,125,128,137,
%T A330400 153,187,197,233,257,389,401,512,577,641,677,685,703,773,901,929,977,
%U A330400 1153,1193,1493,1537,1553,1657,1697,2047,2048,2313,2897,3089,3137,3593,4001
%N A330400 Numbers at which the sum of the iterated unitary totient function (A329153) attains a record.
%C A330400 Analogous to A181659 with the unitary totient function (A047994) instead of the Euler totient function phi (A000010).
%C A330400 The corresponding record values are 0, 1, 3, 6, 10, 16, 24, 39, 55, 70, 85, ... (see the link for more values).
%H A330400 Amiram Eldar, <a href="/A330400/b330400.txt">Table of n, a(n) for n = 1..637</a> (terms below 2*10^9)
%H A330400 Amiram Eldar, <a href="/A330400/a330400.txt">Table of n, a(n), A329153(a(n)) for n = 1..637</a>
%t A330400 uphi[1] = 1; uphi[n_] := Times @@ (-1 + Power @@@ FactorInteger[n]); s[n_] := Plus @@ FixedPointList[uphi, n] - n - 1; seq = {}; smax = -1; Do[s1 = s[n]; If[s1 > smax, smax = s1; AppendTo[seq, n]], {n, 1, 5000}]; seq
%Y A330400 Cf. A047994, A181659, A329153.
%K A330400 nonn
%O A330400 1,2
%A A330400 _Amiram Eldar_, Feb 25 2020