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.

A007539 a(n) = first n-fold perfect (or n-multiperfect) number.

Original entry on oeis.org

1, 6, 120, 30240, 14182439040, 154345556085770649600, 141310897947438348259849402738485523264343544818565120000
Offset: 1

Views

Author

Keywords

Comments

On the Riemann Hypothesis, a(n) > exp(exp(n / e^gamma)) for n > 3. Unconditionally, a(n) > exp(exp(0.9976n / e^gamma)) for n > 3, where the constant is related to A004394(1000000). - Charles R Greathouse IV, Sep 06 2012
Each of the terms 1, 6, 120, 30240 divides all larger terms <= a(8). See A227765, A227766, ..., A227769. - Jonathan Sondow, Jul 30 2013
Is a(n) < a(n+1)? - Jeppe Stig Nielsen, Jun 16 2015
Equivalently, a(n) is the smallest number k such that sigma(k)/k = n. - Derek Orr, Jun 19 2015
The number of divisors of these terms are: 1, 4, 16, 96, 1920, 110592, 1751777280, 63121588161085440. - Michel Marcus, Jun 20 2015
Given n, let S_n be the sequence of integers k that satisfy numerator(sigma(k)/k) = n. Then a(n) is a member of S_n. In fact a(n) = S_n(i), where the successive values of i are 1, 1, 2, 2, 4, 2, (23, 6, 31, 12, ...), where the terms in parentheses need to be confirmed. - Michel Marcus, Nov 22 2015
The first four terms are the only multiperfect numbers in A025487 among the 1600 initial terms of A007691. Can it be proved that these are the only ones among the whole A007691? See also A349747. - Antti Karttunen, Dec 04 2021

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 22.
  • A. Brousseau, Number Theory Tables. Fibonacci Association, San Jose, CA, 1973, p. 138.
  • R. K. Guy, Unsolved Problems in Number Theory, B2.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    Table[k = 1; While[DivisorSigma[1, k]/k != n, k++]; k, {n, 4}] (* Michael De Vlieger, Jun 20 2015 *)
  • PARI
    a(n)=k=1;while((sigma(k)/k)!=n,k++);k
    vector(4,n,a(n)) \\ Derek Orr, Jun 19 2015

Extensions

More terms sent by Robert G. Wilson v, Nov 30 2000

A349746 Numbers k for which k * gcd(sigma(k), u) is equal to sigma(k) * gcd(k, u), where u is obtained by shifting the prime factorization of k two steps toward larger primes [with u = A003961(A003961(k))].

Original entry on oeis.org

1, 11466, 114660, 411264, 804384, 871416, 4999680, 46332000, 176417280, 378069120, 396168192, 485188704, 709430400, 2004912000, 3921372000, 5600534400, 6128179200, 6956471808, 7556976000, 7746979968, 9904204800, 14092001280, 14182439040, 23423662080, 31998395520
Offset: 1

Views

Author

Antti Karttunen, Nov 30 2021

Keywords

Comments

Sigma preserves both the 2-adic and 3-adic valuation of the terms of this sequence.
All 65 known 5-multiperfect numbers (A046060) are included in this sequence, as well as the smallest 7-multiperfect number, 141310897947438348259849402738485523264343544818565120000 = A007539(7), and probably the majority of other p-multiperfect numbers as well, where p is a prime > 3. However, any term that is in A349747 is not included in this sequence.

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := NextPrime[p, 2]^e; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; q[n_] := n * GCD[(sigma = DivisorSigma[1, n]), (u = s[n])] == sigma * GCD[n, u]; Select[Range[10^6], q] (* Amiram Eldar, Dec 01 2021 *)
  • PARI
    A003961twice(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(1+nextprime(1+f[i, 1]))); factorback(f); };
    isA349746(n) = { my(s=sigma(n),u=A003961twice(n)); (n*gcd(s,u) == (s*gcd(n,u))); };

Formula

For all n >= 1, A007814(A000203(a(n))) = A007814(a(n)) and A007949(A000203(a(n))) = A007949(a(n)). [See comment]

Extensions

a(15)-a(25) from Martin Ehrenstein, Dec 17 2021
Showing 1-2 of 2 results.