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-2 of 2 results.

A333872 Numbers at which the sum of the iterated absolute Möbius divisor function (A173557) attains a record.

Original entry on oeis.org

1, 2, 3, 5, 7, 11, 17, 19, 23, 31, 41, 43, 47, 59, 71, 79, 83, 103, 107, 131, 139, 167, 223, 227, 263, 347, 359, 383, 467, 479, 563, 587, 659, 719, 839, 863, 887, 1019, 1163, 1187, 1223, 1259, 1283, 1307, 1319, 1367, 1439, 1823, 1979, 2027, 2039, 2207, 2447, 2879
Offset: 1

Views

Author

Amiram Eldar, Apr 08 2020

Keywords

Comments

Analogous to A181659 with the absolute Möbius divisor function (A173557) instead of the Euler totient function phi (A000010).
The corresponding record values are 0, 1, 3, 5, 9, 15, 17, 21, 37, 39, 45, ... (see the link for more values).

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := p - 1; u[1] = 1; u[n_] := Times @@ (f @@@ FactorInteger[n]); s[n_] := Plus @@ FixedPointList[u, n] - n - 1; seq = {}; smax = -1; Do[s1 = s[n];  If[s1 > smax, smax = s1; AppendTo[seq, n]], {n, 1, 3000}]; seq

A333873 Numbers that equal to the sum of their iterated absolute Möbius divisor function (A173557).

Original entry on oeis.org

3, 5, 17, 257, 413, 611, 1391, 1589, 1903, 2327, 5599, 27959, 29623, 36647, 36983, 38863, 42851, 43919, 46463, 49513, 65537, 76759, 82969, 86567, 88759, 96839, 111179, 116479, 129307, 171191, 184979, 213041, 277619, 301157, 310519, 346151, 362263, 372227, 375167
Offset: 1

Views

Author

Amiram Eldar, Apr 08 2020

Keywords

Comments

A variant of A082897 (perfect totient numbers) in which the absolute Möbius divisor function (A173557) replaces the Euler totient function (A000010).

Examples

			5 is a term since A173557(5) = 4, A173557(4) = 1, and 4 + 1 = 5.
		

Crossrefs

A019434 is a subsequence.

Programs

  • Mathematica
    f[p_, e_] := p - 1; u[1] = 1; u[n_] := Times @@ (f @@@ FactorInteger[n]); Select[Range[10^4], Plus @@ FixedPointList[u, #] == 2*# + 1 &]
Showing 1-2 of 2 results.