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

A361785 Indices of records in the sequence of bi-unitary harmonic means A361782(k)/A361783(k).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 10, 12, 15, 20, 24, 30, 40, 54, 56, 60, 84, 96, 120, 168, 210, 240, 270, 280, 360, 420, 480, 672, 840, 1080, 1320, 1512, 1680, 1890, 2160, 2310, 2520, 3080, 3360, 4320, 5280, 6048, 7392, 7560, 9240, 10920, 11880, 14040, 15120, 18480, 20790
Offset: 1

Views

Author

Amiram Eldar, Mar 24 2023

Keywords

Examples

			The harmonic means of the bi-unitary divisors of the first 6 positive integers are 1 < 4/3 < 3/2 < 8/5 < 5/3 < 2. A361782(7)/A361783(7) = 9/5 < 2, and the next record, A361782(8)/A361783(8) = 32/15, occurs at 8. Therefore, the first 7 terms of this sequence are 1, 2, 3, 4, 5, 6 and 8.

Links

Crossrefs

Cf. A361782, A361783.
Similar sequences: A179971, A348654, A361319.
Other sequences related to records of bi-unitary divisors: A293185, A292983, A292984.

Programs

  • Mathematica
    f[p_, e_] := p^e * If[OddQ[e], (e + 1)*(p - 1)/(p^(e + 1) - 1), e/((p^(e + 1) - 1)/(p - 1) - p^(e/2))]; buhmean[1] = 1; buhmean[n_] := Times @@ f @@@ FactorInteger[n]; seq[kmax_] := Module[{buh, buhmax = 0, s = {}}, Do[buh = buhmean[k]; If[buh > buhmax, buhmax = buh; AppendTo[s, k]], {k, 1, kmax}]; s]; seq[20000]
  • PARI
    buhmean(n) = {my(f = factor(n), p, e); n * prod(i = 1, #f~, p = f[i, 1]; e = f[i, 2];  if(e%2, (e + 1)*(p - 1)/(p^(e + 1) - 1), e/((p^(e + 1) - 1)/(p - 1) - p^(e/2)))); }
lista(kmax) = {my(buh, buhmax=0); for(k = 1, kmax, buh = buhmean(k); if(buh > buhmax, buhmax = buh; print1(k, ", "))); }

A286325 Bi-unitary harmonic numbers.

Original entry on oeis.org

1, 6, 45, 60, 90, 270, 420, 630, 672, 2970, 5460, 8190, 9072, 9100, 10080, 15925, 22680, 22848, 27300, 30240, 40950, 45360, 54600, 81900, 95550, 99792, 136500, 163800, 172900, 204750, 208656, 245700, 249480, 312480, 332640, 342720, 385560, 409500, 472500, 491400
Offset: 1

Views

Author

Michel Marcus, May 07 2017

Keywords

Comments

A number m is a term if the sum of its bi-unitary divisors, A188999(m) divides the product of m by the number of its bi-unitary divisors A286324(m).
Numbers k whose harmonic mean of their bi-unitary divisors, A361782(k)/A361783(k), is an integer. - Amiram Eldar, Mar 24 2023

Links

Crossrefs

Cf. A222266, A188999, A286324, A361782, A361783, A361784.
Cf. A001599 (Ore harmonic), A006086 (unitary harmonic).

Programs

  • Mathematica
    f[p_, e_] := p^e * If[OddQ[e], (e + 1)*(p - 1)/(p^(e + 1) - 1), e/((p^(e + 1) - 1)/(p - 1) - p^(e/2))]; bhQ[n_] := IntegerQ[Times @@ f @@@ FactorInteger[n]]; bhQ[1] = True; Select[Range[10^5], bhQ] (* Amiram Eldar, Mar 24 2023 *)
  • PARI
    udivs(n) = {my(d = divisors(n)); select(x->(gcd(x, n/x)==1), d); }
gcud(n, m) = vecmax(setintersect(udivs(n), udivs(m)));
biudivs(n) = select(x->(gcud(x, n/x)==1), divisors(n));
isok(n) = my(v=biudivs(n)); denominator(n*#v/vecsum(v))==1;

A361782 Numerators of the harmonic means of the bi-unitary divisors of the positive integers.

Original entry on oeis.org

1, 4, 3, 8, 5, 2, 7, 32, 9, 20, 11, 12, 13, 7, 5, 64, 17, 12, 19, 8, 21, 22, 23, 16, 25, 52, 27, 14, 29, 10, 31, 64, 11, 68, 35, 72, 37, 38, 39, 32, 41, 7, 43, 44, 3, 23, 47, 32, 49, 100, 17, 104, 53, 18, 55, 56, 57, 116, 59, 4, 61, 31, 63, 384, 65, 11, 67, 136
Offset: 1

Views

Author

Amiram Eldar, Mar 24 2023

Keywords

Examples

			Fractions begin with 1, 4/3, 3/2, 8/5, 5/3, 2, 7/4, 32/15, 9/5, 20/9, 11/6, 12/5, ...

Links

Crossrefs

Cf. A188999, A222266, A286324, A361783 (denominators).
Similar sequences: A099377, A103339, A361316.

Programs

  • Mathematica
    f[p_, e_] := p^e * If[OddQ[e], (e + 1)*(p - 1)/(p^(e + 1) - 1), e/((p^(e + 1) - 1)/(p - 1) - p^(e/2))]; a[1] = 1; a[n_] := Numerator[Times @@ f @@@ FactorInteger[n]]; Array[a, 100]
  • PARI
    a(n) = {my(f = factor(n), p, e); numerator(n * prod(i = 1, #f~, p = f[i, 1]; e = f[i, 2];  if(e%2, (e + 1)*(p - 1)/(p^(e + 1) - 1), e/((p^(e + 1) - 1)/(p - 1) - p^(e/2))))); }

Formula

a(n) = numerator(n*A286324(n)/A188999(n)).
f(n) = a(n)/A361783(n) is multiplicative with f(p^e) = (e+1)*(p-1)/(p^(e+1)-1) if e is odd, and e/((p^(e+1)-1)/(p-1) - p^(e/2)) if e is even.

A361784 Harmonic means the bi-unitary divisors of the bi-unitary harmonic numbers (A286325).

Original entry on oeis.org

1, 2, 3, 4, 4, 6, 7, 7, 8, 11, 13, 13, 12, 10, 16, 7, 18, 16, 15, 24, 15, 20, 20, 18, 14, 22, 25, 24, 19, 25, 23, 27, 33, 31, 44, 32, 34, 30, 25, 36, 13, 46, 31, 21, 29, 40, 38, 33, 28, 40, 48, 38, 29, 45, 34, 47, 28, 32, 32, 44, 60, 27, 32, 28, 46, 26, 51
Offset: 1

Views

Author

Amiram Eldar, Mar 24 2023

Keywords

Examples

			a(3) = 3 since A286325(3) = 45, the bi-unitary divisors of 45 are 1, 5, 9, and 45, and their harmonic mean is 3.

Links

Crossrefs

Cf. A188999, A286324, A286325, A361782, A361783.
Similar sequences: A001600, A006087, A361318.

Programs

  • Mathematica
    f[p_, e_] := p^e * If[OddQ[e], (e + 1)*(p - 1)/(p^(e + 1) - 1), e/((p^(e + 1) - 1)/(p - 1) - p^(e/2))]; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; Select[Array[s, 10^5], IntegerQ]
  • PARI
    bhmean(n) = {my(f = factor(n), p, e); n * prod(i = 1, #f~, p = f[i, 1]; e = f[i, 2];  if(e%2, (e + 1)*(p - 1)/(p^(e + 1) - 1), e/((p^(e + 1) - 1)/(p - 1) - p^(e/2)))); }
lista(kmax) = {my(bh); for(k = 1, kmax, bh = bhmean(k); if(denominator(bh) == 1, print1(bh, ", "))); }

Formula

a(n) = A361782(A286325(n)).
Showing 1-4 of 4 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.