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.

Previous Showing 31-40 of 47 results. Next

A362082 Numbers k achieving record deficiency via a residue-based measure, M(k) = (k+1)*(1 - zeta(2)/2) - 1 - ( Sum_{j=1..k} k mod j )/k.

Original entry on oeis.org

1, 5, 11, 23, 47, 59, 167, 179, 359, 503, 719, 1439, 5039, 6719, 7559, 15119, 20159, 52919, 75599, 83159, 166319, 415799, 720719, 831599, 1081079, 2162159, 4324319, 5266799, 7900199, 10533599, 18345599, 28274399, 41081039, 136936799, 205405199, 410810399
Offset: 1

Views

Author

Richard Joseph Boland, Apr 17 2023

Keywords

Comments

M(k) = (k+1)*(1 - zeta(2)/2) - 1 - ( Sum_{j=1..k} k mod j )/k is a measure of either abundance (sigma(k) > 2*k), or deficiency (sigma(k) < 2*k), of a positive integer k. The measure follows from the known facts that Sum_{j=1..k} (sigma(j) + k mod j) = k^2 and that the average order of sigma(k)/k is Pi^2/6 = zeta(2) (see derivation below).
M(k) ~ 0 when sigma(k) ~ 2*k and for sufficiently large k, M(k) is positive when k is an abundant number (A005101) and negative when k is a deficient number (A005100).
The terms of this sequence are the deficient k for which M(k) < M(m) for all m < k and may be thought of as "superdeficient", contra-analogous to the superabundant numbers A004394 utilizing sigma(k)/k as the measure of abundance, which is otherwise not particularly meaningful as a deficiency measure.
15119=13*1163 is the first term that is composite and subsequently, up to 1000000000, roughly half of the terms are composite.

Examples

			First few terms with their M(k) measure and factorizations as generated by the Mathematica program:
    1   -0.64493406684822643647   {{1,1}}
    5   -0.73480220054467930942   {{5,1}}
   11   -0.86960440108935861883  {{11,1}}
   23   -1.0000783673961085420   {{23,1}}
   47   -1.0528856894638174541   {{47,1}}
   59   -1.1107338698535727552   {{59,1}}
  167   -1.1984137110594038972  {{167,1}}
  179   -1.2619431113124463216  {{179,1}}
  359   -1.3499704727921791778  {{359,1}}
  503   -1.3722914063892448936  {{503,1}}
  719   -1.4363475145965658088  {{719,1}}
		

Crossrefs

Cf. A362081 (analogous to superabundant A004394).
Cf. A362083 (analogous to A335067, A326393).

Programs

  • Mathematica
    Clear[min, Rp, R, seqtable, M]; min = 1; Rp = 0; seqtable = {};
    Do[R = Rp + 2 k - 1 - DivisorSigma[1, k];
      M = N[(k + 1)*(1 - Zeta[2]/2) - 1 - R/k, 20];
      If[M < min, min = M; Print[k, "   ", min, "   ", FactorInteger[k]];
       AppendTo[seqtable, k]];
      Rp = R, {k, 1, 1000000000}];
    Print[seqtable]
  • PARI
    M(n) = (n+1)*(1 - zeta(2)/2) - 1 - sum(k=2, n, n%k)/n;
    lista(nn) = my(m=+oo, list=List()); for (n=1, nn, my(mm = M(n)); if (mm < m, listput(list, n); m = mm);); Vec(list); \\ Michel Marcus, Apr 21 2023

Formula

Derived starting with lemmas 1-3:
1) Sum_{j=1..k} (sigma(j) + k mod j) = k^2.
2) The average order of sigma(k)/k is Pi^2/6 = zeta(2).
3) R(k) = Sum_{j=1..k} k mod j, so R(k)/k is the average order of (k mod j).
Then:
Sum_{j=1..k} sigma(j) ~ zeta(2)*Sum_{j=1..k} j = zeta(2)*(k^2+k)/2.
R(k)/k ~ k - k*zeta(2)/2 - zeta(2)/2.
0 ~ (k+1)*(1 - zeta(2)/2) - 1 - R(k)/k.
Thus M(k) = (k+1)*(1 - zeta(2)/2) - 1 - R(k)/k is a measure of variance about sigma(k) ~ 2*k corresponding to M(k) ~ 0.

A362083 Numbers k such that, via a residue based measure M(k) (see Comments), k is deficient, k+1 is abundant, and abs(M(k)) + abs(M(k+1)) reaches a new maximum.

Original entry on oeis.org

11, 17, 19, 47, 53, 103, 347, 349, 557, 1663, 1679, 2519, 5039, 10079, 15119, 25199, 27719, 55439, 110879, 166319, 277199, 332639, 554399, 665279, 720719, 1441439, 2162159, 3603599, 4324319, 7207199, 8648639, 10810799, 21621599, 36756719, 61261199, 73513439, 122522399, 147026879
Offset: 1

Views

Author

Richard Joseph Boland, Apr 17 2023

Keywords

Comments

The residue-based quantifier function, M(k), measures either abundance (sigma(k) > 2*k), or deficiency (sigma(k) < 2*k), of a positive integer k. The measure is defined by M(k) = (k+1)*(1 - zeta(2)/2) - 1 - (Sum_{j=1..k} k mod j)/k. It follows from the known facts that Sum_{j=1..k} (sigma(j) + k mod j) = k^2 and that the average order of sigma(k)/k is Pi^2/6 = zeta(2) (see derivation below).
M(k) ~ 0 when sigma(k) ~ 2*k and for sufficiently large k, M(k) is positive when k is an abundant number (A005101) and negative when k is a deficient number (A005100). The terms of this sequence are the deficient k such that k+1 is abundant and abs(M(k)) + abs(M(k+1)) achieves a new maximum, somewhat analogous to A335067 and A326393.

Examples

			The first few terms with measure sums and factorizations generated by the Mathematica program:
0.90610439514731535319   35  {{5,1},{7,1}}   36   {{2,2},{3,2}}
1.1735781643159997761    59  {{59,1}}        60   {{2,2},{3,1},{5,1}}
1.3642976724582397229   119  {{7,1},{17,1}} 120   {{2,3},{3,1},{5,1}}
1.3954100615479538209   179  {{179,1}}      180   {{2,2},{3,2},{5,1}}
1.4600817810807682323   239  {{239,1}}      240   {{2,4},{3,1},{5,1}}
1.6088158511317518390   359  {{359,1}}      360   {{2,3},{3,2},{5,1}}
1.7153941935887132383   719  {{719,1}}      720   {{2,4},{3,2},{5,1}}
1.7851979872921589879   839  {{839,1}}      840   {{2,3},{3,1},{5,1},{7,1}}
		

Crossrefs

Cf. A362081 (analogous to superabundant A004394), A362082 (superdeficient).

Programs

  • Mathematica
    Clear[max, Rp, R, seqtable, Mp, M];max = -1; Rp = 0; Mp = -0.644934066; seqtable = {};
    Do[R = Rp + 2 k - 1 - DivisorSigma[1, k];
     M = N[(k)*(1 - Zeta[2]/2) - 1  - R/k, 20];
     If[DivisorSigma[1, k - 1] < 2 (k - 1) && DivisorSigma[1, k] > 2 k &&
       Abs[Mp] + Abs[M] > max, max = Abs[Mp] + Abs[M];
      Print[max, "   ", k - 1, "   ", FactorInteger[k - 1], "   ", k,
       "   ", FactorInteger[k]]; AppendTo[seqtable, {k - 1, k}]]; Rp = R;
     Mp = M, {k, 2, 1000000000}]; seq = Flatten[seqtable]; Table[seq[[2 j - 1]], {j, 1, Length[seq]/2}]

Formula

Derived starting with lemmas 1-3:
1) Sum_{j=1..k} (sigma(j) + k mod j) = k^2.
2) The average order of sigma(k)/k is Pi^2/6 = zeta(2).
3) R(k) = Sum_{j=1..k} k mod j, so R(k)/k is the average order of (k mod j).
Then:
Sum_{j=1..k} sigma(j) ~ zeta(2)*Sum_{j=1..k} j = zeta(2)*(k^2+k)/2.
R(k)/k ~ k - k*zeta(2)/2 - zeta(2)/2.
0 ~ (k+1)*(1 - zeta(2)/2) - 1 - R(k)/k.
Thus M(k) = (k+1)*(1 - zeta(2)/2) - 1 - R(k)/k is a measure of variance about sigma(k) ~ 2*k corresponding to M(k) ~ 0.

A370634 A135507(n) is the product of the first n terms of this sequence.

Original entry on oeis.org

1, 4, 5, 3, 3, 3, 9, 4, 3, 3, 13, 3, 3, 9, 3, 3, 19, 3, 3, 3, 9, 13, 25, 3, 3, 3, 3, 9, 31, 3, 3, 4, 13, 19, 9, 3, 39, 3, 3, 3, 43, 9, 3, 13, 3, 25, 49, 3, 3, 3, 19, 3, 55, 3, 3, 3, 3, 31, 61, 3, 3, 3, 3, 3, 3, 3, 69, 19, 3, 3, 73, 3, 3, 39, 3, 3, 3, 3, 81, 3
Offset: 1

Views

Author

Michael De Vlieger, May 19 2024

Keywords

Comments

Compactification of A135507 akin to A000705 with respect to A002201.

Examples

			Table of first 20 terms of this sequence and S = A135507.
   n            S(n)  a(n)
  ------------------------
   1              1     1
   2              4     4
   3             20     5
   4             60     3
   5            180     3
   6            540     3
   7           4860     9
   8          19440     4
   9          58320     3
  10         174960     3
  11        2274480    13
  12        6823440     3
  13       20470320     3
  14      184232880     9
  15      552698640     3
  16     1658095920     3
  17    31503822480    19
  18    94511467440     3
  19   283534402320     3
  20   850603206960     3
		

Crossrefs

Programs

  • Mathematica
    nn = 120; j = 1; {1}~Join~Reap[Do[k = 2 j + LCM[j, i]; Sow[k/j]; j = k, {i, 2, nn}] ][[-1, 1]]

Formula

For n > 1, 3 <= a(n) <= n+2.
For p = A001359(i) such that gcd(a(p-1), p) = 1, a(p) = p+2 = A006512(i).

A376687 Numbers that set records in in A376281.

Original entry on oeis.org

24, 96, 120, 240, 360, 480, 840, 1080, 1680, 2160, 2520, 3360, 4320, 5040, 7560, 10080, 15120, 27720, 30240, 55440, 60480, 83160, 110880, 151200, 166320, 221760, 277200, 332640, 498960, 554400, 665280, 831600, 997920, 1330560, 1441440, 1663200, 2162160, 2882880
Offset: 1

Views

Author

Michael De Vlieger, Jan 08 2025

Keywords

Comments

Proper subset of the intersection A025487 and A379336.
There are three kinds of pairs (d, k/d), d | k, such that gcd(d, k/d) does not equal 1, d, or k/d:
Type A: rad(d) does not divide k/d, and rad(k/d) does not divide d (see A379752), where rad = A007947.
Type B: the squarefree kernel of one divisor divides the other but the reverse is not true (see A379772).
Type C: rad(d) = rad(k/d), i.e., d, k/d, and k are coreful (see A379552).
Conjecture: Numbers k that set records in A376281 do not have type C divisor pairs, i.e., those that are coreful but neither divides the other. This, since type C requires k to be powerful and divisible by cubes of 2 distinct primes (i.e., in A376936). Therefore the record is achieved only through large numbers of type A and B.
Since type A divisor pairs are common for composite k in A375055, this sequence is resembles A379752.
Since d and k/d are both composite, this sequence resembles A059992.
This sequence, to a lesser extent A379752, and a greater extent A059992, contains many highly composite numbers. (See plot of S(n) = union of this sequence and A002182 below, and corresponding graphs in respective other sequences.)

Examples

			Let b(n) = A376281(n).
Table showing exponents of prime power factors of a(n) for n = 1..20.
Example: a(5) = 360 = 2^3 * 3^2 * 5, hence we write "3.2.1".
   n    a(n)  Exp.   b(a(n))
  ----------------------------------
   1     24 *   3.1        1   4*6
   2     96     5.1        2   6*16 = 8*12
   3    120 **  3.1.1      3   4*30 = 6*20 = 10*12
   4    240 *   4.1.1      4   6*40 = 8*30 = 10*24 = 12*20
   5    360 **  3.2.1      5   4*90 = 10*36 = 12*30 = 15*24 = 18*20
   6    480     5.1.1      6   6*80 = 8*60 = 10*48 = 12*40 = 16*30 = 20*24
   7    840 *   3.1.1.1    7
   8   1080     3.3.1      9
   9   1680 *   4.1.1.1   10
  10   2160     4.3.1     11
  11   2520 **  3.2.1.1   13
  12   3360     5.1.1.1   14
*  = a(n) is highly composite (in A002182),
** = a(n) is superior highly composite (in both A002182 and A002201).
		

Crossrefs

Programs

  • Mathematica
    (* Load function f at A025487 *)
    r = 0; s = Select[Union@ Flatten@ f[8][[3 ;; -1]], Not@*SquareFreeQ];
    nn = Length[s]; Print[nn]
    Reap[Monitor[
      Do[k = s[[i]];
        If[# > r, r = #; Sow[k]] &@
          Count[Transpose@{#, k/#} &@ #[[2 ;; Ceiling[Length[#]/2]]] &@ Divisors[k],
            _?(And[1 < GCD @@ {##}, Mod[#1, #2] != 0,
               Mod[#2, #1] != 0] & @@ # &)], {i, nn}], i] ][[-1, 1]]

A379753 Numbers that set records in A379752.

Original entry on oeis.org

60, 120, 240, 480, 840, 1260, 1680, 2520, 3360, 5040, 7560, 10080, 15120, 20160, 27720, 36960, 55440, 83160, 110880, 166320, 221760, 277200, 332640, 443520, 498960, 554400, 665280, 720720, 997920, 1081080, 1330560, 1441440, 2162160, 2882880, 3603600, 4324320, 5765760
Offset: 1

Views

Author

Michael De Vlieger, Jan 01 2025

Keywords

Comments

Proper subset of the intersection of A025487 and A375055.
Conjecture: subset of A332785 = A126706 \ A286708.
This sequence seems to be rich in highly composite numbers, the prime shape of a(n) resembles that of highly composite numbers, with long tails of large prime factors with multiplicity 1.
Terms not in A002182 are not all of the form 2^5 * prime(i..j), 1 < i < j, for example, a(24) = 443520 = 2^7 * 3^2 * 5 * 7 * 11.

Examples

			Let b(n) = A379752(n).
Table showing exponents of prime power factors of a(n) for n = 1..12.
Example: a(6) = 1260 = 2^2 * 3^2 * 5 * 7, hence we write "2.2.1.1".
   n      a(n)       Exp.    b(a(n))
  ----------------------------------
   1       60 **   2.1.1        1   6*10
   2      120 **   3.1.1        2   6*20 = 10*12
   3      240 *    4.1.1        3   6*40 = 10*24 = 12*20
   4      480      5.1.1        4   6*80 = 10*48 = 12*40 = 20*24
   5      840 *    3.1.1.1      6   6*140 = 10*84 = 12*70 = 14*60 = 20*42 = 28*30
   6     1260 *    2.2.1.1      7
   7     1680 *    4.1.1.1      9
   8     2520 **   3.2.1.1     11
   9     3360      5.1.1.1     12
  10     5040 **   4.2.1.1     15
  11     7560 *    3.3.1.1     16
  12    10080 *    5.2.1.1     19
*  = a(n) is highly composite (in A002182),
** = a(n) is superior highly composite (in both A002182 and A002201).
		

Crossrefs

Programs

  • Mathematica
    (* Load function f at A025487 *)
    r = 0;
    s = Select[Union@ Flatten@ f[14][[4 ;; -1]], Not@*SquareFreeQ];
    nn = Length[s]; Print[nn];
    Reap[Do[k = s[[i]]; If[# > r, r = #; Sow[k]] &@
      Count[Transpose@ {#, k/#} &@ #[[2 ;; Ceiling[Length[#]/2] ]] &@ Divisors[k],
        _?(And[1 < GCD @@ {##},
           Nor[Divisible[#2, rad[#1]],
               Divisible[#1, rad[#2]] ] ] & @@ # &)], {i, nn}] ][[-1, 1]]

A002497 Numbers N in A002809 such that there is rho > 0 such that for all A > 0, A008475(A)-A008475(N) >= rho*log(A/N).

Original entry on oeis.org

3, 12, 60, 420, 4620, 60060, 180180, 360360, 6126120, 116396280, 2677114440, 77636318760, 2406725881560, 89048857617720, 3651003162326520, 156993135980040360, 313986271960080720, 14757354782123793840, 14757354782123793840, 782139803452561073520, 46146248403701103337680
Offset: 1

Views

Author

Keywords

Comments

The numbers contain the starred entries on pp. 187-190 of Nicolas. It is a subsequence of A002809 by selecting only elements of a set/property "G" (page 150). G contains all N such that a real, strictly positive rho exists such that for all strictly positive integers A we have l(A)-l(N) >= rho*log(A/N). The function l()=A008475() is defined on page 139. - R. J. Mathar, Mar 23 2012
These numbers were named superior l-composite numbers (nombres l-composes superieurs, the function l(n) is A002809) by Massias, in analogy to Ramanujan's superior highly composite numbers (A002201). Deléglise and Nicolas named these numbers l-superchampion numbers. They are used by Deléglise et al. in calculating values of Landau's function g(n) (A000793). - Amiram Eldar, Aug 23 2019

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Extensions

Edited by M. F. Hasler, Mar 29 2015
a(16)-a(21) from the paper by Massias added by Amiram Eldar, Aug 23 2019

A098895 Number of divisors of the n-th superior highly composite number.

Original entry on oeis.org

2, 4, 6, 12, 16, 24, 48, 60, 120, 240, 288, 384, 576, 1152, 2304, 2688, 5376, 8064, 16128, 20160, 40320, 46080, 92160, 184320, 368640, 737280, 983040, 1966080, 3932160, 4423680, 6635520, 13271040, 15925248, 31850496, 63700992, 127401984
Offset: 1

Views

Author

David Terr, Oct 14 2004

Keywords

Comments

Sequence A002201 gives the values of the n-th superior highly composite number N(n) and A000705 gives the values of the (prime) ratio N(n)/N(n-1).

Examples

			a(8)=60 because the eighth superior highly composite number, 5040, has 60 divisors.
		

References

  • S. Ramanujan, Highly composite numbers, Proc. London Math. Soc., 14 (1915), 347-407. Reprinted in Collected Papers, Ed. G. H. Hardy et al., Cambridge 1927; Chelsea, NY, 1962, pp. 78-129. See esp. pp. 87, 115.

Crossrefs

Formula

a(n) = a(n-1) * (1 + 1/k(n)), where k(n) is the p(n)-adic valuation of the n-th superior highly composite number N(n), with p(n) = N(n)/N(n-1) and N(0)=1.

A098896 p(n)-adic valuation of the n-th superior highly composite number N(n), where p(n) = N(n)/N(n-1).

Original entry on oeis.org

1, 1, 2, 1, 3, 2, 1, 4, 1, 1, 5, 3, 2, 1, 1, 6, 1, 2, 1, 4, 1, 7, 1, 1, 1, 1, 3, 1, 1, 8, 2, 1, 5, 1, 1, 1, 1, 2, 1, 1, 9, 1, 1, 1, 1, 3, 1, 1, 2, 1, 1, 1, 1, 6, 4, 1, 1, 2, 10, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 11, 1, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

David Terr, Oct 14 2004

Keywords

Comments

(1+1/a(n)) appears in the denominators of the log arguments of the denominators of the numbers in the table of the reference, pp. 115-117.

Examples

			a(8) = 4 since N(8)=5040 has 2-adic valuation of 4 and N(8)/N(7)=2.
		

Crossrefs

Extensions

More terms from Amiram Eldar, Aug 12 2019

A113069 Number of highly composite numbers between two consecutive superior highly composite numbers.

Original entry on oeis.org

1, 0, 3, 0, 2, 4, 0, 8, 9, 1, 3, 6, 11, 12, 1, 13, 7, 14, 3, 14, 2, 15, 16, 15, 17, 7, 18, 19, 4, 12, 22, 3, 22, 24, 23, 25, 14, 25, 25, 3, 26, 26, 30, 30, 11, 32, 33, 19, 34, 33, 34, 35, 6, 10, 35, 35, 23, 4, 35, 35, 36, 37, 38, 36, 35
Offset: 1

Views

Author

T. D. Noe, Oct 13 2005

Keywords

Comments

The SHC numbers are a subset of the HC numbers. Is there a formula for a(n) that depends on the two consecutive SHC numbers A002201(n) and A002201(n+1)?

Examples

			Example: a(3)=3 because between the SHC numbers 12 and 60 there are three HC numbers: 24, 36 and 48.
		

Crossrefs

Cf. A002182 (highly composite numbers), A002201 (superior highly composite numbers).

A189466 Number of superior highly composite numbers < 10^n.

Original entry on oeis.org

2, 4, 6, 8, 9, 10, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 23, 24, 25, 25, 26, 27, 28, 28, 29, 30, 31, 32, 33, 33, 34, 34, 35, 36, 36, 37, 38, 38, 39, 39, 40, 41, 41, 42, 42, 43, 43, 44, 44, 46, 46, 46, 47, 47, 48, 49, 49, 50, 50, 51, 51, 52, 52, 53
Offset: 1

Views

Author

Krzysztof Ostrowski, Apr 22 2011

Keywords

Comments

Number of superior highly composite numbers (A002201) with at most n digits.

Examples

			a(2) = 4 since there are 4 superior highly composite numbers < 10^2 {2,6,12,60}
		

Crossrefs

Previous Showing 31-40 of 47 results. Next