A286325 -id:A286325 - cp's OEIS frontend

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

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

Amiram Eldar, Mar 24 2023

nonn

			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.

Amiram Eldar, Table of n, a(n) for n = 1..313
Jozsef Sandor, On bi-unitary harmonic numbers, arXiv:1105.0294 [math.NT], 2011.

Cf. A188999, A286324, A286325, A361782, A361783.

Similar sequences: A001600, A006087, A361318.

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]

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, ", "))); }

a(n) = A361782(A286325(n)).

A319745 Nonunitary harmonic numbers: numbers such that the harmonic mean of their nonunitary divisors is an integer.

Original entry on oeis.org

4, 9, 12, 18, 24, 25, 45, 49, 54, 60, 112, 121, 126, 150, 168, 169, 270, 289, 294, 336, 361, 529, 560, 594, 637, 726, 841, 961, 1014, 1232, 1369, 1638, 1680, 1681, 1734, 1849, 1984, 2166, 2184, 2209, 2430, 2520, 2688, 2700, 2809, 2850, 3174, 3481, 3721, 3780

Offset: 1

Amiram Eldar, Sep 27 2018

nonn

Includes all the numbers with a single nonunitary divisor. Those with more than one: 12, 18, 24, 45, 54, 60, 112, ...
Supersequence of A064591 (nonunitary perfect numbers).
Ligh & Wall showed that if p, 2p-1 and 2^p-1 are distinct primes (A172461, except for 2), then the following numbers are in the sequence: 6*p^2, p^2*(2p-1), 6*p^2*(2p-1), 2^(p+1)*3*(2^p-1), 2^(p+1)*15*(2^p-1) and 2^(p+1)*(2p-1)*(2^p-1).

Amiram Eldar, Table of n, a(n) for n = 1..10000
Steve Ligh and Charles R. Wall, Functions of Nonunitary Divisors, Fibonacci Quarterly, Vol. 25 (1987), pp. 333-338.

Cf. A001599, A006086, A063947, A064591, A172461, A286325.

nudiv[n_] := Block[{d = Divisors[n]}, Select[d, GCD[#, n/#] > 1 &]]; nhQ[n_]:= Module[ {divs=nudiv[n]}, Length[divs] > 0 && IntegerQ[HarmonicMean[divs]]]; Select[Range[30000], nhQ]

hm(v) = #v/sum(k=1, #v, 1/v[k]);
vnud(n) = select(x->(gcd(x, n/x)!=1), divisors(n));
isok(n) = iferr(denominator(hm(vnud(n))) == 1, E, 0); \\ Michel Marcus, Oct 28 2018

A348964 Exponential harmonic (or e-harmonic) numbers of type 2: numbers k such that the harmonic mean of the exponential divisors of k is an integer.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 12, 13, 14, 15, 17, 18, 19, 21, 22, 23, 26, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 46, 47, 51, 53, 55, 57, 58, 59, 60, 61, 62, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 82, 83, 84, 85, 86, 87, 89, 90, 91, 93, 94

Offset: 1

Amiram Eldar, Nov 05 2021

nonn

Sándor (2006) proved that all the squarefree numbers are e-harmonic of type 2.
Equivalently, numbers k such that A348963(k) | k * A049419(k).
Apparently, most exponential harmonic numbers of type 1 (A348961) are also terms of this sequence. Those that are not exponential harmonic numbers of type 2 are 1936, 5808, 9680, 13552, 17424, 29040, ...

			The squarefree numbers are trivial terms. If k is squarefree, then it has a single exponential divisor, k itself, and thus the harmonic mean of its exponential divisors is also k, which is an integer.
12 is a term since its exponential divisors are 6 and 12, and their harmonic mean, 8, is an integer.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Nicuşor Minculete, Contribuţii la studiul proprietăţilor analitice ale funcţiilor aritmetice - Utilizarea e-divizorilor, Ph.D. thesis, Academia Română, 2012. See section 4.3, pp. 90-94.
József Sándor, On exponentially harmonic numbers, Scientia Magna, Vol. 2, No. 3 (2006), pp. 44-47.
József Sándor, Selected Chapters of Geomety, Analysis and Number Theory, 2005, pp. 141-145.

A005117 and A348965 are subsequences.
Cf. A049419, A322791, A348961, A348963.
Similar sequences: A001599, A006086, A063947, A286325, A319745.

f[p_, e_] := p^e * DivisorSigma[0, e] / DivisorSum[e, p^(e-#) &]; ehQ[1] = True; ehQ[n_] := IntegerQ[Times @@ f @@@ FactorInteger[n]]; Select[Range[100], ehQ]

A361783 Denominators of the harmonic means of the bi-unitary divisors of the positive integers.

Original entry on oeis.org

1, 3, 2, 5, 3, 1, 4, 15, 5, 9, 6, 5, 7, 3, 2, 27, 9, 5, 10, 3, 8, 9, 12, 5, 13, 21, 10, 5, 15, 3, 16, 21, 4, 27, 12, 25, 19, 15, 14, 9, 21, 2, 22, 15, 1, 9, 24, 9, 25, 39, 6, 35, 27, 5, 18, 15, 20, 45, 30, 1, 31, 12, 20, 119, 21, 3, 34, 45, 8, 9, 36, 25, 37, 57

Offset: 1

Amiram Eldar, Mar 24 2023

Amiram Eldar, Table of n, a(n) for n = 1..10000
Jozsef Sandor, On bi-unitary harmonic numbers, arXiv:1105.0294 [math.NT], 2011.

Cf. A188999, A222266, A286324, A286325 (positions of 1's), A361782 (numerators).

Similar sequences: A099378, A103340, A361317.

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_] := Denominator[Times @@ f @@@ FactorInteger[n]]; Array[a, 100]

a(n) = {my(f = factor(n), p, e); denominator(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))))); }

a(n) = denominator(n*A286324(n)/A188999(n)).

A349026 Exponential unitary harmonic numbers: numbers k such that the harmonic mean of the exponential unitary divisors of k is an integer.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 12, 13, 14, 15, 17, 18, 19, 21, 22, 23, 26, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 46, 47, 51, 53, 55, 57, 58, 59, 60, 61, 62, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 82, 83, 84, 85, 86, 87, 89, 90, 91, 93, 94

Offset: 1

Amiram Eldar, Nov 06 2021

nonn

First differs from A348964 at n = 102. a(102) = 144 is not an exponential harmonic number of type 2.
The exponential unitary divisors of n = Product p(i)^e(i) are all the numbers of the form Product p(i)^b(i) where b(i) is a unitary divisor of e(i) (see A278908).
Equivalently, numbers k such that A349025(k) | k * A278908(k).

			The squarefree numbers are trivial terms. If k is squarefree, then it has a single exponential unitary divisor, k itself, and thus the harmonic mean of its exponential unitary divisors is also k, which is an integer.
144 is a term since its exponential unitary divisors are 6, 18, 48 and 144, and their harmonic mean, 16, is an integer.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Nicuşor Minculete, Contribuţii la studiul proprietăţilor analitice ale funcţiilor aritmetice - Utilizarea e-divizorilor, Ph.D. thesis, Academia Română, 2012. See section 4.3, pp. 90-94.

Cf. A278908 (number of exponential unitary divisors), A322857, A322858, A323310, A349025, A349027.

Similar sequences: A001599, A006086, A063947, A286325, A319745, A348964.

f[p_, e_] := p^e * 2^PrimeNu[e] / DivisorSum[e, p^(e - #) &, CoprimeQ[#, e/#] &]; euhQ[1] = True; euhQ[n_] := IntegerQ[Times @@ f @@@ FactorInteger[n]]; Select[Range[100], euhQ]

A362804 Numbers k such that the set of divisors {d | k, BitOr(k, d) = k} has an integer harmonic mean.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 16, 24, 28, 30, 32, 45, 48, 56, 60, 64, 90, 96, 112, 120, 128, 180, 192, 224, 240, 256, 360, 384, 448, 480, 496, 512, 720, 768, 896, 960, 992, 1024, 1440, 1536, 1792, 1920, 1984, 2048, 2880, 3072, 3584, 3840, 3968, 4096, 5760, 6144, 7168, 7680

Offset: 1

Amiram Eldar, May 04 2023

Equivalently, the set of divisors can be defined by {d | k, BitAnd(k, d) = d}.
Analogous to harmonic (or Ore) numbers (A001599) where the divisors d of k are restricted by BitOr(k, d) = k or BitAnd(k, d) = d.
If k is a term then so is 2*k. The primitive terms are in A362805. Thus, this sequence includes all the powers of 2 (A000079), all the numbers of the form 3*2^m and 15*2^m for m >= 1, and all the numbers of the form 7*2^m for m >= 2.
All the even perfect numbers (A000396) are terms: if k = 2^(p-1)*(2^p-1) is a perfect number (where p is a Mersenne exponent, A000043), then the only divisors of k such that BitOr(k, d) = k are 2^(p-1) and k itself, and the harmonic mean of 2^(p-1) and 2^(p-1)*(2^p-1) is 2^p - 1.
Are 1 and 45 the only odd terms in this sequence?

Amiram Eldar, Table of n, a(n) for n = 1..406

Cf. A000043, A000396, A246600, A246601.
Subsequences: A000079, A007283 \ {3}, A005009 \ {7, 14}, A110286 \ {15}, A362805.
Similar sequences: A001599, A006086, A063947, A286325, A319745.

q[n_] := IntegerQ[HarmonicMean[Select[Divisors[n], BitAnd[n, #] == # &]]]; Select[Range[10^4], q]

div(n) = select(x->(bitor(x, n) == n), divisors(n));
is(n) = {my(d = div(n)); denominator(#d/sum(i = 1, #d, 1/d[i])) == 1;}

A335387 Tri-unitary harmonic numbers: numbers k such that the harmonic mean of the tri-unitary divisors of k is an integer.

Original entry on oeis.org

1, 6, 45, 60, 90, 270, 420, 630, 2970, 5460, 8190, 9100, 15925, 27300, 36720, 40950, 46494, 47520, 54600, 81900, 95550, 136500, 163800, 172900, 204750, 232470, 245700, 257040, 332640, 409500, 464940, 491400, 646425, 716625, 790398, 791700, 819000, 900900, 929880

Offset: 1

Amiram Eldar, Jun 04 2020

nonn

Equivalently, numbers k such that A324706(k) | (k * A335385(k)).

Differs from A063947 from n >= 18.

			45 is a term since its tri-unitary divisors are {1, 5, 9, 45} and their harmonic mean, 3, in an integer.

Amiram Eldar, Table of n, a(n) for n = 1..200

A324707 is a subsequence.
Analogous sequences: A001599 (harmonic numbers), A006086 (unitary), A063947 (infinitary), A286325 (bi-unitary), A319745 (nonunitary).
Cf. A324706, A335385.

f1[p_, e_] := If[e == 3 || e == 6, 4, 2]; f2[p_, e_] := If[e == 3, (p^4 - 1)/(p - 1), If[e == 6, (p^8 - 1)/(p^2 - 1), p^e + 1]]; f[p_, e_] := p^e * f1[p, e]/f2[p, e]; tuhQ[1] = True; tuhQ[n_] := IntegerQ[Times @@ (f @@@ FactorInteger[n])]; Select[Range[10^4], tuhQ]

A348918 Noninfinitary harmonic numbers: numbers such that the harmonic mean of their noninfinitary divisors is an integer.

Original entry on oeis.org

4, 9, 12, 18, 25, 45, 49, 60, 96, 112, 121, 126, 150, 169, 289, 294, 336, 361, 448, 486, 529, 540, 560, 600, 637, 672, 726, 841, 961, 1014, 1232, 1344, 1350, 1369, 1638, 1680, 1681, 1734, 1849, 2166, 2209, 2430, 2809, 2850, 3174, 3481, 3721, 3822, 4200, 4320, 4489

Offset: 1

Amiram Eldar, Nov 04 2021

nonn

Includes all the squares of primes (A001248), since they are the numbers with a single noninfinitary divisor.

			12 is a term since its noninfinitary divisors are {2, 6}, and their harmonic mean, 3, is an integer.

Amiram Eldar, Table of n, a(n) for n = 1..5000

Cf. A001248, A001599, A006086, A063947, A286325, A319745.

nidiv[1] = {}; nidiv[n_] := Complement[Divisors[n], Sort@ Flatten@ Outer[Times, Sequence @@ (FactorInteger[n] /. {p_, m_Integer} :> p^Select[Range[0, m], BitOr[m, #] == m &])]]; Select[Range[5000], (d = nidiv[#]) != {} && IntegerQ@ HarmonicMean[d] &]

A349180 Coreful harmonic numbers: nonsquarefree numbers k such that the harmonic mean of the coreful divisors of k is an integer.

Original entry on oeis.org

12, 18, 36, 56, 60, 75, 84, 90, 126, 132, 150, 156, 168, 180, 198, 204, 228, 234, 240, 252, 276, 280, 306, 342, 348, 351, 372, 392, 396, 414, 420, 444, 450, 468, 492, 504, 516, 522, 525, 558, 564, 588, 612, 616, 630, 636, 660, 666, 684, 702, 708, 720, 726, 728

Offset: 1

Amiram Eldar, Nov 09 2021

nonn

A divisor of a number k is coreful if it is divisible by every prime that divides k.

The sequence is restricted to nonsquarefree numbers since the squarefree numbers have a single coreful divisor and thus they trivially have an integer harmonic mean.

			12 is a term since its coreful divisors are 6 and 12 and their harmonic mean, 8, is an integer.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A005117, A013929, A057723, A307958, A307959.

Similar sequences: A001599, A006086, A063947, A286325, A319745, A348964.

rad[n_] := Times @@ FactorInteger[n][[;; , 1]]; corHarmQ[n_] := Module[{r = rad[n], d}, d = Select[Divisors[n], rad[#] == r &]; IntegerQ[HarmonicMean[d]]]; Select[Range[10^3], !SquareFreeQ[#] && corHarmQ[#] &]

A349181 Powerful harmonic numbers: numbers k such that the set of powerful divisors of k that are larger than 1 has more than one element and that the harmonic mean of this set is an integer.

Original entry on oeis.org

100, 300, 700, 1100, 1225, 1300, 1700, 1900, 2100, 2300, 2450, 2900, 3100, 3300, 3675, 3700, 3900, 4100, 4225, 4300, 4700, 5100, 5300, 5700, 5900, 6100, 6700, 6900, 7100, 7300, 7350, 7700, 7900, 8300, 8450, 8700, 8900, 9100, 9300, 9700, 10100, 10300, 10700, 10900

Offset: 1

Amiram Eldar, Nov 09 2021

nonn

Numbers with a single powerful divisor > 1 are A060687 and trivially have an integer harmonic mean.

The least term that is not divisible by 5 (or 25) is a(5446) = 1413721.

			100 is a term since its powerful divisors > 1 are 4, 25 and 100 and their harmonic mean, 10, is an integer.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A001694, A060687.

Similar sequences: A001599, A006086, A063947, A286325, A319745, A348964.

powQ[n_] := Min[FactorInteger[n][[;; , 2]]] > 1; powHarmQ[n_] := Module[{d = Select[Divisors[n], powQ]}, Length[d] > 1 && IntegerQ[HarmonicMean[d]]]; Select[Range[10^4], powHarmQ]

cp's OEIS Frontend

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

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A319745 Nonunitary harmonic numbers: numbers such that the harmonic mean of their nonunitary divisors is an integer.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

A348964 Exponential harmonic (or e-harmonic) numbers of type 2: numbers k such that the harmonic mean of the exponential divisors of k is an integer.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A361783 Denominators of the harmonic means of the bi-unitary divisors of the positive integers.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A349026 Exponential unitary harmonic numbers: numbers k such that the harmonic mean of the exponential unitary divisors of k is an integer.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A362804 Numbers k such that the set of divisors {d | k, BitOr(k, d) = k} has an integer harmonic mean.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

A335387 Tri-unitary harmonic numbers: numbers k such that the harmonic mean of the tri-unitary divisors of k is an integer.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A348918 Noninfinitary harmonic numbers: numbers such that the harmonic mean of their noninfinitary divisors is an integer.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs