A235796 -id:A235796 - cp's OEIS frontend

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.

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

Richard Joseph Boland, Apr 17 2023

nonn

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.

			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}}

Jeffrey C. Lagarias, An Elementary Problem Equivalent to the Riemmann Hypothesis, arXiv:math/0008177 [math.NT], 2000-2001; Amer. Math. Monthly, 109 (2002), 534-543.

Cf. A362081 (analogous to superabundant A004394).
Cf. A362083 (analogous to A335067, A326393).
Cf. A004490, A002201, A005100, A005101, A004125, A024916, A000290, A120444, A235796, A000396, A000079.

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]

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

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

Richard Joseph Boland, Apr 17 2023

nonn

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.

			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}}

Jeffrey C. Lagarias, An Elementary Problem Equivalent to the Riemmann Hypothesis, arXiv:math/0008177 [math.NT], 2000-2001; Amer. Math. Monthly, 109 (2002), 534-543.

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

Cf. A335067, A326393, A004490, A002201, A326393, A005100, A005101, A004125, A024916, A000290, A120444, A235796, A000396, A000079.

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}]

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.

A244824 Sum of all divisors of all positive integers <= 2^n.

Original entry on oeis.org

1, 4, 15, 56, 220, 857, 3403, 13535, 54077, 215900, 862954, 3450545, 13802279, 55201838, 220792018, 883134861, 3532518195, 14129951284, 56519699688, 226078355122, 904312961284, 3617249936000, 14468996179294, 57875977567596, 231503907383054, 926015589350438

Offset: 0

Omar E. Pol, Jul 06 2014

nonn

Has a symmetric representation, the same as A024916.

			For n = 2 the sum of all divisors of all positive integers <= 4 is [1] + [1+2] + [1+3] + [1+2+4] = 1 + 3 + 4 + 7 = 15, so a(2) = 15.

Chai Wah Wu, Table of n, a(n) for n = 0..75

Cf. A000079, A000203, A024916, A196020, A235796, A236104, A237270, A237593.

a(n) = sum(k=1, 2^n, k*floor(2^n/k) ) \\ Jens Kruse Andersen, Jul 26 2014

from math import isqrt
def A244824(n): return -(s:=isqrt(m:=1<>1 # Chai Wah Wu, Oct 23 2023

a(n) = A024916(A000079(n)).

a(n) ~ 2^(2*n-2) * Pi^2/3. - Vaclav Kotesovec, Oct 23 2023

More terms from Jens Kruse Andersen, Jul 26 2014

A380231 Alternating row sums of triangle A237591.

Original entry on oeis.org

1, 2, 1, 2, 1, 4, 3, 4, 5, 4, 3, 6, 5, 4, 7, 8, 7, 8, 7, 10, 9, 8, 7, 10, 11, 10, 9, 12, 11, 14, 13, 14, 13, 12, 15, 16, 15, 14, 13, 16, 15, 18, 17, 16, 19, 18, 17, 20, 21, 22, 21, 20, 19, 22, 21, 24, 23, 22, 21, 24, 23, 22, 25, 26, 25, 28, 27, 26, 25, 28, 27, 32, 31, 30, 29, 28, 31, 30, 29

Offset: 1

Omar E. Pol, Jan 17 2025

nonn

Consider the symmetric Dyck path in the first quadrant of the square grid described in the n-th row of A237593. Let C = (A240542(n), A240542(n)) be the middle point of the Dyck path.
a(n) is also the coordinate on the x axis of the point (a(n),n) and also the coordinate on the y axis of the point (n,a(n)) such that the middle point of the line segment [(a(n),n),(n,a(n))] coincides with the middle point C of the symmetric Dyck path.
The three line segments [(a(n),n),C], [(n,a(n)),C] and [(n,n),C] have the same length.
For n > 2 the points (n,n), C and (a(n),n) are the vertices of a virtual isosceles right triangle.
For n > 2 the points (n,n), C and (n,a(n)) are the vertices of a virtual isosceles right triangle.
For n > 2 the points (a(n),n), (n,n) and (n,a(n)) are the vertices of a virtual isosceles right triangle.

			For n = 14 the 14th row of A237591 is [8, 3, 1, 2] hence the alternating row sum is 8 - 3 + 1 - 2 = 4, so a(14) = 4.
On the other hand the 14th row of A237593 is the 14th row of A237591 together with the 14 th row of A237591 in reverse order as follows: [8, 3, 1, 2, 2, 1, 3, 8].
Then with the terms of the 14th row of A237593 we can draw a Dyck path in the first quadrant of the square grid as shown below:
.
         (y axis)
          .
          .
          .    (4,14)              (14,14)
          ._ _ _ . _ _ _ _            .
          .               |
          .               |
          .               |_
          .                 |
          .                 |_ _
          .                C    |_ _ _
          .                           |
          .                           |
          .                           |
          .                           |
          .                           . (14,4)
          .                           |
          .                           |
          . . . . . . . . . . . . . . | . . . (x axis)
        (0,0)
.
In the example the point C is the point (9,9).
The three line segments [(4,14),(9,9)], [(14,4),(9,9)] and [(14,14),(9,9)] have the same length.
The points (14,14), (9,9) and (4,14) are the vertices of a virtual isosceles right triangle.
The points (14,14), (9,9) and (14,4) are the vertices of a virtual isosceles right triangle.
The points (4,14), (14,14) and (14,4) are the vertices of a virtual isosceles right triangle.

Other alternating row sums (ARS) related to the Dyck paths of A237593 and the stepped pyramid described in A245092 are as follows:
ARS of A237593 give A000004.
ARS of A196020 give A000203.
ARS of A252117 give A000203.
ARS of A271343 give A000593.
ARS of A231347 give A001065.
ARS of A236112 give A004125.
ARS of A236104 give A024916.
ARS of A249120 give A024916.
ARS of A271344 give A033879.
ARS of A231345 give A033880.
ARS of A239313 give A048050.
ARS of A237048 give A067742.
ARS of A236106 give A074400.
ARS of A235794 give A120444.
ARS of A266537 give A146076.
ARS of A236540 give A153485.
ARS of A262612 give A175254.
ARS of A353690 give A175254.
ARS of A239446 give A235796.
ARS of A239662 give A239050.
ARS of A235791 give A240542.
ARS of A272026 give A272027.
ARS of A211343 give A336305.
Cf. A221529, A237270, A237271, A237591, A262626, A322141.

row235791(n) = vector((sqrtint(8*n+1)-1)\2, i, 1+(n-(i*(i+1)/2))\i);
a(n) = my(orow = concat(row235791(n), 0)); vecsum(vector(#orow-1, i, (-1)^(i+1)*(orow[i] - orow[i+1]))); \\ Michel Marcus, Apr 13 2025

a(n) = 2*A240542(n) - n.
a(n) = n - 2*A322141(n).
a(n) = A240542(n) - A322141(n).

cp's OEIS Frontend

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.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

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.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A244824 Sum of all divisors of all positive integers <= 2^n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Python

Formula

Extensions

A380231 Alternating row sums of triangle A237591.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

PARI

Formula