A354801 n^2 minus the sum of all aliquot divisors of all positive integers <= n.

1, 3, 7, 11, 19, 24, 36, 44, 57, 68, 88, 95, 119, 136, 156, 172, 204, 218, 254, 271, 301, 330, 374, 385, 428, 463, 503, 530, 586, 603, 663, 695, 745, 792, 848, 864, 936, 989, 1049, 1078, 1158, 1187, 1271, 1318, 1374, 1439, 1531, 1550, 1639, 1695, 1775, 1832, 1936, 1977, 2069, 2116

Offset

1,2

Comments

After the Dyck paths described in A237593 we can see that a(n) has a symmetric representation as follows: a(n) is the sum of the areas of two polygons. In the fourth quadrant of the infinite square grid the first polygon has a vertex at (0,0) and its area is equal to A000217(n). The second polygon appears if n >= 3 and it has a vertex at (n,-n) and its area is equal to A004125(n). So the area of the arrowhead-shaped polygon is equal to A153485(n). See the illustration of initial terms in the Links section.

Links

Table of n, a(n) for n=1..56.

Omar E. Pol, Illustration of initial terms

Formula

a(n) = A000217(n) + A004125(n).

a(n) = A000290(n) - A153485(n).

a(n) = A024916(n) + A004125(n) - A153485(n).

Mathematica

Accumulate[Table[3*n - 1 - DivisorSigma[1, n], {n, 1, 60}]] (* Amiram Eldar, Jun 12 2022 *)

Prog

(PARI) a(n) = n^2 - sum(k=1, n, sigma(k)-k); \\ Michel Marcus, Jun 13 2022
(Python)
from math import isqrt
def A354801(n): return n*(3*n+1)+(s:=isqrt(n))**2*(s+1)-sum((q:=n//k)*((k<<1)+q+1) for k in range(1, s+1))>>1 # Chai Wah Wu, Oct 22 2023

Crossrefs

Partial sums of A353190.

Cf. A000217, A000290, A001065, A004125, A024916, A153485, A236104, A237591, A237593.

Sequence in context: A135932 A231847 A105876 * A141101 A236632 A098379

Adjacent sequences: A354798 A354799 A354800 * A354802 A354803 A354804

Keyword

nonn,easy

Author

Omar E. Pol, Jun 06 2022

Status

approved