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.

A247077 Composite numbers for which the harmonic mean of proper divisors is an integer.

Original entry on oeis.org

1645, 88473, 63626653506
Offset: 1

Views

Author

Daniel Lignon, Nov 17 2014

Keywords

Comments

Of course, for all prime numbers the harmonic mean of proper divisors is an integer.
a(4) >= 2*10^11. - Hiroaki Yamanouchi, Nov 20 2014
Conjecture: all terms are of the form m*(sigma(m)-1) where sigma(m)-1 is prime. - Chai Wah Wu, Dec 15 2020
a(4) <= 22351741783447265625. - Daniel Suteu, Dec 16 2020
From Chai Wah Wu, Feb 04 2021: (Start)
Other terms of the sequence of the form m*(sigma(m)-1) correspond to the following values of m:
3 * 5^143
3 * 5^623
3 * 5^1423
5 * 7^127
5 * 7^6595
101 * 103^25
(End)
Equivalently, composite numbers k such that sigma(k)-1 divides k*(tau(k)-1), where sigma = A000203 and tau = A000005. - Daniel Suteu, Feb 05 2021

Examples

			The proper divisors of 1645 are [1,5,7,35,47,235,329] and their harmonic mean is 7/(1/1 + 1/5 + 1/7 + 1/35 + 1/47 + 1/235 + 1/329) = 5.
		

Crossrefs

Cf. A001599 for harmonic mean of all divisors and A247078 for harmonic mean of nontrivial divisors.

Programs

  • Mathematica
    Select[Range[2,100000],(IntegerQ[HarmonicMean[Most[Divisors[#]]]] && Not[PrimeQ[#]])&] (* Daniel Lignon, Nov 17 2014 *)
  • PARI
    lista(nn) = forcomposite (n=2, nn, my(d=divisors(n)); if (denominator((#d-1)/sum(i=1, #d-1, 1/d[i])) == 1, print1(n, ", "))); \\ Michel Marcus, Nov 17 2014
    
  • PARI
    isok(n) = n > 1 && !isprime(n) && (n*(numdiv(n)-1)) % (sigma(n)-1) == 0; \\ Daniel Suteu, Feb 05 2021

Extensions

a(3) from Hiroaki Yamanouchi, Nov 20 2014

A335267 Composite numbers whose harmonic mean of their divisors that are larger than 1 is an integer.

Original entry on oeis.org

6, 15, 28, 30, 91, 117, 135, 252, 270, 496, 703, 864, 936, 1891, 1989, 2295, 2701, 4284, 4590, 5733, 8128, 8432, 12403, 18721, 19872, 21528, 38503, 41580, 49141, 51319, 56896, 79003, 88831, 104653, 121920, 146611, 188191, 218791, 226801, 235053, 269011, 286903
Offset: 1

Views

Author

Amiram Eldar, May 29 2020

Keywords

Comments

The primes are excluded from this sequence since they are trivial terms.
The corresponding harmonic means are 3, 5, 5, 5, 13, 9, 9, 9, 9, 9, 37, ...
Equivalently, composite numbers m such that (sigma(m)-m) | m*(tau(m)-1), or A001065(m) | A168014(m).
The semiprimes terms of this sequence are of the form p*q where p and q = 2*p - 1 are primes (A129521).
If m is a k-perfect numbers, k = 2, 3, ... (i.e., sigma(m) = k*m), then sigma(m)-m = (k-1)*m. If (k-1)*m | m*(tau(m)-1) then (k-1) | (tau(m)-1). If k is odd then tau(m) is also odd, so m is a square, and sigma(m) is odd. Since m | sigma(m) this means that m is also odd. Since there is no known odd multiply-perfect number except for 1 (A007691), there are no known k-perfect numbers with odd k in this sequence.
The perfect numbers (k=2, A000396) are terms: if m is a perfect number then sigma(m)-m = m.
The 4-perfect number (k=4, A027687) m are terms if 3 | (tau(m)-1). Of the first 36 terms of A027687 there are 8 such terms, the first is A027687(26).
The 6-perfect number (k=6, A046061) m are terms if 5 | (tau(m)-1). Of the first 245 terms of A046061 there are 20 such terms, the first is A046061(19).
Hemiperfect numbers that are terms of this sequence include A055153(i) for i = 10, 18 and 20, A141645(21), and A159271(i) for i = 97 and 103.

Examples

			6 is a term since its divisors other than 1 are 2, 3 and 6, and their harmonic mean, 3/(1/2 + 1/3 + 1/6) = 3, is an integer.
		

Crossrefs

A000396 and A129521 are subsequences.
Similar sequences: A001599, A247077, A247078.
Cf. A000005 (tau), A000203 (sigma).

Programs

  • Mathematica
    Select[Range[10^6], CompositeQ[#] && Divisible[# * (DivisorSigma[0, #] - 1), DivisorSigma[1, #] - #] &]
    Select[Range[287000],CompositeQ[#]&&IntegerQ[HarmonicMean[ Rest[ Divisors[ #]]]]&] (* Harvey P. Dale, Jan 21 2021 *)

A335269 Numbers for which the harmonic mean of the nontrivial unitary divisors is an integer.

Original entry on oeis.org

228, 345, 1645, 2120, 4025, 4386, 4977, 7725, 8041, 13026, 23881, 24157, 24336, 51925, 88473, 115957, 150161, 169893, 229177, 255041, 278721, 322592, 342637, 377201, 490725, 538625, 656937, 1497517, 1566981, 2132021, 3256261, 3847001, 4646101, 5054221, 5524897
Offset: 1

Views

Author

Amiram Eldar, May 29 2020

Keywords

Comments

A number m is a term if the set {d|m ; d > 1, d < m, gcd(d, m/d) = 1} is nonempty and the harmonic mean its members is an integer.
The corresponding harmonic means are 8, 9, 15, 16, 25, 12, 21, 15, 33, 12, ...
Equivalently, numbers m such that omega(m) > 1 and (usigma(m)-m-1) | m*(2^omega(m)-2), where usigma is the sum of unitary divisors (A034448), and 2^omega(m)-2 = A034444(m)-2 = A087893 (m) is the number of the nontrivial unitary divisors of m.
The squarefree terms of A247078 are also terms of this sequence.

Examples

			228 is a term since the harmonic mean of its nontrivial unitary divisors, {3, 4, 12, 19, 57, 76} is 8 which is an integer.
		

Crossrefs

The unitary version of A247078.

Programs

  • Mathematica
    usigma[1] = 1; usigma[n_] := Times @@ (1 + Power @@@ FactorInteger[n]); Select[Range[10^6], (omega = PrimeNu[#]) > 1 && Divisible[#*(2^omega - 2), usigma[#] - # - 1] &]

A247079 Numbers for which the harmonic mean of nontrivial divisors is an integer and which are not a square of prime numbers.

Original entry on oeis.org

345, 1050, 1645, 4386, 6489, 8041, 13026, 23881, 88473, 115957, 255041, 342637, 377201, 1497517, 2132021, 2428489, 3256261, 3847001, 4114285, 4646101, 5054221, 6816865, 7218641, 7587901, 13384681, 14872837, 17897737, 20901553, 23807821, 25863653, 28207957
Offset: 1

Views

Author

Daniel Lignon, Nov 17 2014

Keywords

Comments

That's the numbers which are in A247078 and not in A001248.
a(149) >= 2*10^11. - Hiroaki Yamanouchi, Nov 20 2014

Examples

			345 is not the square of a prime number and the nontrivial divisors of 345 are [3,5,15,23,69,115]. Their harmonic mean is 6/(1/3+1/5+1/15+1/23+1/69+1/115)=9.
		

Crossrefs

Cf. A001599 (harmonic numbers), A247078.

Programs

  • PARI
    isok(n) = !(issquare(n) && isprime(sqrtint(n))) && (d=divisors(n)) && (#d > 2) && (denominator((#d-2)/sum(i=2, #d-1, 1/d[i])) == 1); \\ Michel Marcus, Nov 17 2014

Extensions

a(15)-a(24) from Michel Marcus, Nov 17 2014
Showing 1-4 of 4 results.