A319526 -id:A319526 - cp's OEIS frontend

A065764 Sum of divisors of square numbers.

1, 7, 13, 31, 31, 91, 57, 127, 121, 217, 133, 403, 183, 399, 403, 511, 307, 847, 381, 961, 741, 931, 553, 1651, 781, 1281, 1093, 1767, 871, 2821, 993, 2047, 1729, 2149, 1767, 3751, 1407, 2667, 2379, 3937, 1723, 5187, 1893, 4123, 3751, 3871, 2257, 6643

Offset: 1

Labos Elemer, Nov 19 2001

Unlike A065765, the sums of divisors of squares give remainders r=1,3,5 modulo 6: sigma(4)==1, sigma(49)==3, sigma(2401)==5 (mod 6). See also A097022.
a(n) is also the number of ordered pairs of positive integers whose LCM is n, (see LeVeque). - Enrique Pérez Herrero, Aug 26 2013
Main diagonal of A319526. - Omar E. Pol, Sep 25 2018
Subsequence of primes is A023195 \ {3}; also, 31 is the only known prime to be twice in the data because 31 = sigma(16) = sigma(25) (see A119598 and Goormaghtigh conjecture link). - Bernard Schott, Jan 17 2021

W. J. LeVeque, Fundamentals of Number Theory, pp. 125 Problem 4, Dover NY 1996.

T. D. Noe, Table of n, a(n) for n=1..10000
Wikipedia, Goormaghtigh conjecture.

Cf. A000203, A000290, A028982, A319526.

a:=List([1..50],n->Sigma(n^2));; Print(a); # Muniru A Asiru, Jan 01 2019

[SumOfDivisors(n^2): n in [1..48]]; // Bruno Berselli, Apr 12 2011

with(numtheory): [sigma(n^2)$n=1..50]; # Muniru A Asiru, Jan 01 2019

Table[Plus@@Divisors[n^2], {n, 48}] (* Alonso del Arte, Feb 24 2012 *)
f[p_, e_] := (p^(2*e + 1) - 1)/(p - 1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 50] (* Amiram Eldar, Sep 10 2020 *)

numlib::sigma(n^2)$ n=1..81 // Zerinvary Lajos, May 13 2008

a(n) = sigma(n^2); \\ Harry J. Smith, Oct 30 2009

from math import prod
from sympy import factorint
def A065764(n): return prod((p**((e<<1)+1)-1)//(p-1) for p,e in factorint(n).items()) # Chai Wah Wu, Oct 25 2023

[sigma(n^2,1)for n in range(1,49)] # Zerinvary Lajos, Jun 13 2009

a(n) = sigma(n^2) = A000203(A000290(n)).
Multiplicative with a(p^e) = (p^(2*e+1)-1)/(p-1). - Vladeta Jovovic, Dec 01 2001
Dirichlet g.f.: zeta(s)*zeta(s-1)*zeta(s-2)/zeta(2*s-2), inverse Mobius transform of A000082. - R. J. Mathar, Mar 06 2011
Dirichlet convolution of A001157 by the absolute terms of A055615. Also the Dirichlet convolution of A048250 by A000290. - R. J. Mathar, Apr 12 2011
a(n) = Sum_{d|n} d*Psi(d), where Psi is A001615. - Enrique Pérez Herrero, Feb 25 2012
a(n) >= (n+1) * sigma(n) - n, where sigma is A000203, equality holds if n is in A000961. - Enrique Pérez Herrero, Apr 21 2012
Sum_{k=1..n} a(k) ~ 5*Zeta(3) * n^3 / Pi^2. - Vaclav Kotesovec, Jan 30 2019
Sum_{k>=1} 1/a(k) = 1.3947708738535614499846243600124612760835313454790187655653356563282177118... - Vaclav Kotesovec, Sep 20 2020

A062731 Sum of divisors of 2*n.

Original entry on oeis.org

3, 7, 12, 15, 18, 28, 24, 31, 39, 42, 36, 60, 42, 56, 72, 63, 54, 91, 60, 90, 96, 84, 72, 124, 93, 98, 120, 120, 90, 168, 96, 127, 144, 126, 144, 195, 114, 140, 168, 186, 126, 224, 132, 180, 234, 168, 144, 252, 171, 217, 216, 210, 162, 280, 216, 248, 240, 210

Offset: 1

Jason Earls, Jul 11 2001

a(n) is also the total number of parts in all partitions of 2*n into equal parts. - Omar E. Pol, Feb 14 2021

N. J. A. Sloane, Table of n, a(n) for n = 1..20000 [First 1000 terms from Harry J. Smith]

Sigma(k*n): A000203 (k=1), A144613 (k=3), A193553 (k=4, even bisection), A283118 (k=5), A224613 (k=6), A283078 (k=7), A283122 (k=8), A283123 (k=9).
Cf. A008438, A074400, A182818, A239052 (odd bisection), A326124 (partial sums), A054784, A215947, A336923, A346870, A346878, A346880, A355750.
Row 2 of A319526. Column & Row 2 of A216626. Row 1 of A355927.
Shallow diagonal (2n,n) of A265652. See also A244658.

[SumOfDivisors(2*n): n in [1..70]]; // Vincenzo Librandi, Oct 31 2014

lst={};Do[AppendTo[lst, DivisorSigma[1, n]], {n, 2, 6!, 2}];lst (* Vladimir Joseph Stephan Orlovsky, Sep 20 2008 *)
DivisorSigma[1,2*Range[60]] (* Harvey P. Dale, Jun 08 2022 *)

numlib::sigma(2*n)$ n=0..81 // Zerinvary Lajos, May 13 2008

PARI
```
vector(66,n,sigma(2*n,1))
```

for (n=1, 1000, write("b062731.txt", n, " ", sigma(2*n)) ) \\ Harry J. Smith, Aug 09 2009

a(n) = A000203(2*n). - R. J. Mathar, Apr 06 2011
a(n) = A000203(n) + A054785(n). - R. J. Mathar, May 19 2020
From Vaclav Kotesovec, Aug 07 2022: (Start)
Dirichlet g.f.: zeta(s) * zeta(s-1) * (3 - 2^(1-s)).
Sum_{k=1..n} a(k) ~ 5 * Pi^2 * n^2 / 24. (End)
From Miles Wilson, Sep 30 2024: (Start)
G.f.: Sum_{k>=1} k*x^(k/gcd(k, 2))/(1 - x^(k/gcd(k, 2))).
G.f.: Sum_{k>=1} k*x^(2*k/(3 + (-1)^k))/(1 - x^(2*k/(3 + (-1)^k))). (End)

Zero removed and offset corrected by Omar E. Pol, Jul 17 2009

A319073 Square array read by antidiagonals upwards: T(n,k) = k*sigma(n), n >= 1, k >= 1.

Original entry on oeis.org

1, 3, 2, 4, 6, 3, 7, 8, 9, 4, 6, 14, 12, 12, 5, 12, 12, 21, 16, 15, 6, 8, 24, 18, 28, 20, 18, 7, 15, 16, 36, 24, 35, 24, 21, 8, 13, 30, 24, 48, 30, 42, 28, 24, 9, 18, 26, 45, 32, 60, 36, 49, 32, 27, 10, 12, 36, 39, 60, 40, 72, 42, 56, 36, 30, 11, 28, 24, 54, 52, 75, 48, 84, 48, 63, 40, 33, 12

Offset: 1

Omar E. Pol, Sep 22 2018

			The corner of the square array begins:
         A000203 A074400 A272027 A239050 A274535 A274536 A319527 A319528
A000027:       1,      2,      3,      4,      5,      6,      7,      8, ...
A008585:       3,      6,      9,     12,     15,     18,     21,     24, ...
A008586:       4,      8,     12,     16,     20,     24,     28,     32, ...
A008589:       7,     14,     21,     28,     35,     42,     49,     56, ...
A008588:       6,     12,     18,     24,     30,     36,     42,     48, ...
A008594:      12,     24,     36,     48,     60,     72,     84,     96, ...
A008590:       8,     16,     24,     32,     40,     48,     56,     64, ...
A008597:      15,     30,     45,     60,     75,     90,    105,    120, ...
A008595:      13,     26,     39,     52,     65,     78,     91,    104, ...
A008600:      18,     36,     54,     72,     90,    108,    126,    144, ...
...

Index entries for sequences related to sigma(n)

Another version of A274824.
Antidiagonal sums give A175254.
Main diagonal gives A064987.
Row n lists the multiples of A000203(n).
Row 1 is A000027.
Initial zeros should be omitted in the following sequences related to the rows of the array:
Row 2-5: A008585, A008586, A008589, A008588.
Rows 6 and 11: A008594.
Rows 7-9: A008590, A008597, A008595.
Rows 10 and 17: A008600.
Rows 12-13: A135628, A008596.
Rows 14, 15 and 23: A008606.
Rows 16 and 25: A135631.
(Note that in the OEIS there are many other sequences that are also rows of this square array.)
Cf. A000203, A237593, A319526.

T:=Flat(List([1..12],n->List([1..n],k->k*Sigma(n-k+1))));; Print(T); # Muniru A Asiru, Jan 01 2019

with(numtheory): T:=(n,k)->k*sigma(n-k+1): seq(seq(T(n,k),k=1..n),n=1..12); # Muniru A Asiru, Jan 01 2019

Table[k DivisorSigma[1, #] &[m - k + 1], {m, 12}, {k, m}] // Flatten (* Michael De Vlieger, Dec 31 2018 *)

cp's OEIS Frontend

A065764 Sum of divisors of square numbers.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

MuPAD

PARI

Python

Sage

Formula

A062731 Sum of divisors of 2*n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

MuPAD

PARI

PARI

Formula

Extensions

A319073 Square array read by antidiagonals upwards: T(n,k) = k*sigma(n), n >= 1, k >= 1.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

GAP

Maple

Mathematica