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.

Showing 1-3 of 3 results.

A349306 Numbers k that divide A092694(k).

Original entry on oeis.org

1, 8, 16, 27, 32, 54, 64, 81, 108, 125, 128, 162, 216, 243, 250, 256, 324, 343, 432, 486, 500, 512, 625, 648, 686, 729, 864, 972, 1000, 1024, 1029, 1250, 1296, 1331, 1372, 1458, 1728, 1944, 2000, 2048, 2058, 2187, 2197, 2401, 2500, 2592, 2662, 2744, 2916, 3087
Offset: 1

Views

Author

Amiram Eldar, Nov 14 2021

Keywords

Comments

Includes all the numbers of the form 2^k * m, where k >= 0 and m is an odd cubefull number, except for 2 and 4. In particular, includes all the cubefull numbers (A036966).
Terms whose odd part is not cubefull are 1029, 2058, 3087, 4116, 6174, 6591, ...

Examples

			8 is a term since A092694(8) = 8 is divisible by 8.
16 is a term since A092694(16) = 64 is divisible by 16.
		

Crossrefs

A036966 is a subsequence.

Programs

  • Mathematica
    Select[Range[3000], Divisible[Times @@ FixedPointList[EulerPhi, #]/#, #] &]

A092693 Sum of iterated phi(n).

Original entry on oeis.org

0, 1, 3, 3, 7, 3, 9, 7, 9, 7, 17, 7, 19, 9, 15, 15, 31, 9, 27, 15, 19, 17, 39, 15, 35, 19, 27, 19, 47, 15, 45, 31, 35, 31, 39, 19, 55, 27, 39, 31, 71, 19, 61, 35, 39, 39, 85, 31, 61, 35, 63, 39, 91, 27, 71, 39, 55, 47, 105, 31, 91, 45, 55, 63, 79, 35, 101, 63, 79, 39, 109, 39, 111
Offset: 1

Views

Author

T. D. Noe, Mar 04 2004

Keywords

Comments

Iannucci, Moujie and Cohen examine perfect totient numbers: n such that a(n) = n.

Examples

			a(100) = 71 because the iterations of phi (40, 16, 8, 4, 2, 1) sum to 71.
		

Crossrefs

Cf. A003434 (iterations of phi(n) needed to reach 1), A092694 (iterated phi product).
Cf. A082897 and A091847 (perfect totient numbers).

Programs

  • Haskell
    a092693 1 = 0
    a092693 n = (+ 1) $ sum $ takeWhile (/= 1) $ iterate a000010 $ a000010 n
    -- Reinhard Zumkeller, Oct 27 2011
    
  • Mathematica
    nMax=100; a=Table[0, {nMax}]; Do[e=EulerPhi[n]; a[[n]]=e+a[[e]], {n, 2, nMax}]; a (* T. D. Noe *)
    Table[Plus @@ FixedPointList[EulerPhi, n] - (n + 1), {n, 72}] (* Alonso del Arte, Jan 29 2007 *)
  • PARI
    a(n)=my(k);while(n>1,k+=n=eulerphi(n));k \\ Charles R Greathouse IV, Mar 22 2012
    
  • Python
    from sympy import totient
    from math import prod
    def f(n):
        m = n
        while m > 1:
            m = totient(m)
            yield m
    def A092693(n): return sum(f(n)) # Chai Wah Wu, Nov 14 2021

Formula

a(1) = 0, a(n) = phi(n) + a(phi(n))
a(n) = A053478(n) - n. - Vladeta Jovovic, Jul 02 2004
Erdős & Subbarao prove that a(n) ~ phi(n) for almost all n. In particular, a(n) < n for almost all n. The proportion of numbers up to N for which a(n) > n is at most 1/log log log log N. - Charles R Greathouse IV, Mar 22 2012

A291782 Let f_k(n) be the result of applying phi (the Euler totient function A000010) k times to n; a(n) = n*Product_{k=1..oo} f_k(n).

Original entry on oeis.org

1, 2, 6, 8, 40, 12, 84, 64, 108, 80, 880, 96, 1248, 168, 960, 1024, 17408, 216, 4104, 1280, 2016, 1760, 40480, 1536, 32000, 2496, 5832, 2688, 77952, 1920, 59520, 32768, 42240, 34816, 53760, 3456, 127872, 8208, 59904, 40960, 1679360, 4032, 173376, 56320
Offset: 1

Views

Author

N. J. A. Sloane, Sep 02 2017

Keywords

Comments

The logarithmic scatterplot of this sequence shows a banded structure similar to that of A092694. - Rémy Sigrist, Sep 03 2017

Examples

			Under phi, 7 -> 6 -> 2 -> 1, so a(7) = 7*6*2 = 84.
		

Crossrefs

Programs

  • Mathematica
    Table[Times @@ FixedPointList[EulerPhi, n], {n, 44}] (* Michael De Vlieger, Sep 03 2017 *)

Formula

a(n) = n * A092694(n). - Rémy Sigrist, Sep 03 2017
Showing 1-3 of 3 results.