login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A001157 sigma_2(n): sum of squares of divisors of n.
(Formerly M3799 N1551)
304

%I M3799 N1551

%S 1,5,10,21,26,50,50,85,91,130,122,210,170,250,260,341,290,455,362,546,

%T 500,610,530,850,651,850,820,1050,842,1300,962,1365,1220,1450,1300,

%U 1911,1370,1810,1700,2210,1682,2500,1850,2562,2366,2650,2210,3410,2451,3255

%N sigma_2(n): sum of squares of divisors of n.

%C If the canonical factorization of n into prime powers is the product of p^e(p) then sigma_k(n) = Product_p ((p^((e(p)+1)*k))-1)/(p^k-1).

%C sigma_2(n) is the sum of the squares of the divisors of n.

%C Sum_{d|n} 1/d^k is equal to sigma_k(n)/n^k. So sequences A017665-A017712 also give the numerators and denominators of sigma_k(n)/n^k for k = 1..24. The power sums sigma_k(n) are in sequences A000203 (k=1), A001157-A001160 (k=2,3,4,5), A013954-A013972 for k = 6,7,...,24. - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 05 2001.

%C Row sums of triangles A134575 and A134559. - _Gary W. Adamson_, Nov 02 2007

%C Also sum of square divisors of n^2. - _Michel Marcus_, Jan 14 2014

%C Conjecture: For each k = 2,3,..., all the rational numbers sigma_k(n)/n^k = Sum_{d|n} 1/d^k (n = 1,2,3,...) have pairwise distinct fractional parts. - _Zhi-Wei Sun_, Oct 15 2015

%C 5 is the only prime entry in the sequence. - _Drake Thomas_, Dec 18 2016

%C 4*a(n) = sum of squares of even divisors of 2*n. - _Wolfdieter Lang_, Jan 07 2017

%D M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 827.

%D T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 38.

%D D. M. Bressoud, Proofs and Confirmations, Camb. Univ. Press, 1999; p. 11.

%D P. A. MacMahon, The connexion between the sum of the squares of the divisors and the number of partitions of a given number, Messenger Math., 54 (1924), 113-116. Collected Papers, MIT Press, 1978, Vol. I, pp. 1364-1367. See Table I. The entry 53 should be 50. - _N. J. A. Sloane_, May 21 2014

%D N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).

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

%H T. D. Noe, <a href="/A001157/b001157.txt">Table of n, a(n) for n = 1..10000</a>

%H M. Abramowitz and I. A. Stegun, eds., <a href="http://www.convertit.com/Go/ConvertIt/Reference/AMS55.ASP">Handbook of Mathematical Functions</a>, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].

%H George E. Andrews, <a href="http://dx.doi.org/10.1007/BF01608779">Stacked lattice boxes</a>, Ann. Comb. 3 (1999), 115-130. See (2.3).

%H Y. Puri and T. Ward, <a href="http://www.cs.uwaterloo.ca/journals/JIS/VOL4/WARD/short.html">Arithmetic and growth of periodic orbits</a>, J. Integer Seqs., Vol. 4 (2001), #01.2.1.

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/DivisorFunction.html">Divisor Function</a>

%H <a href="/index/Si#SIGMAN">Index entries for sequences related to sigma(n)</a>

%H <a href="/index/Cor#core">Index entries for "core" sequences</a>

%F G.f.: Sum_{k>0} k^2 x^k/(1-x^k). Dirichlet g.f.: zeta(s)*zeta(s-2). - _Michael Somos_, Apr 05 2003

%F Multiplicative with a(p^e) = (p^(2e+2)-1)/(p^2-1). - _David W. Wilson_, Aug 01 2001

%F G.f. for sigma_k(n): Sum_{m>0} m^k*x^m/(1-x^m). - _Vladeta Jovovic_, Oct 18 2002

%F L.g.f.: -log(Product_{j>=1} (1-x^j)^j) = Sum_{n>=1} a(n)/n*x^n. - _Joerg Arndt_, Feb 04 2011

%F Equals A127093 * [1, 2, 3, ...]. - _Gary W. Adamson_, May 10 2007

%F Equals A051731 * [1, 4, 9, 16, 25, ...]. A051731 * [1/1, 1/2, 1/3, 1/4,...] = [1/1, 5/4, 10/9, 21/16, 26/25, ...]. - _Gary W. Adamson_, Nov 02 2007

%F Row sums of triangle A134841. - _Gary W. Adamson_, Nov 12 2007

%F a(n) = A035316(n^2). - _Michel Marcus_, Jan 14 2014

%F Conjecture: a(n) = sigma(n^2*rad(n))/sigma(rad(n)), where sigma = A000203 and rad = A007947. - _Velin Yanev_, Aug 20 2017

%F G.f.: Sum_{k>=1} x^k*(1 + x^k)/(1 - x^k)^3. - _Ilya Gutkovskiy_, Oct 24 2018

%p with(numtheory); A001157 := n->sigma[2](n); [seq(sigma[2](n), n=1..100)];

%t Table[DivisorSigma[2, n], {n, 1, 50}] (* _Stefan Steinerberger_, Mar 24 2006 *)

%t DivisorSigma[2,Range[50]] (* _Harvey P. Dale_, Aug 22 2016 *)

%o (PARI) a(n)=if(n<1,0,sigma(n,2))

%o (PARI) a(n)=if(n<1,0,direuler(p=2,n,1/(1-X)/(1-p^2*X))[n])

%o (PARI) a(n)=if(n<1,0,n*polcoeff(sum(k=1,n,x^k/(x^k-1)^2/k,x*O(x^n)),n)) /* _Michael Somos_, Jan 29 2005 */

%o (PARI) N=99; q='q+O('q^N); Vec(sum(n=1,N,n^2*q^n/(1-q^n))) /* _Joerg Arndt_, Feb 04 2011 */

%o (PARI) a(n) = sumdiv(n^2, d, issquare(d)*d); \\ _Michel Marcus_, Jan 14 2014

%o (Sage) [sigma(n,2)for n in xrange(1,51)] # _Zerinvary Lajos_, Jun 04 2009

%o (Maxima) makelist(divsum(n,2),n,1,20); \\ _Emanuele Munarini_, Mar 26 2011

%o (Haskell)

%o a001157 n = s n 1 1 a000040_list where

%o s 1 1 y _ = y

%o s m x y ps'@(p:ps)

%o | m `mod` p == 0 = s (m `div` p) (x * p^2) y ps'

%o | x > 1 = s m 1 (y * (x * p^2 - 1) `div` (p^2 - 1)) ps

%o | otherwise = s m 1 y ps

%o -- _Reinhard Zumkeller_, Jul 10 2011

%o (MAGMA) [DivisorSigma(2,n): n in [1..50]]; // _Bruno Berselli_, Apr 10 2013

%Y Cf. A000005, A000203, A001158, A001159, A053807, A064602, A127093, A134841, A048250.

%Y Cf. also A192794, A082063 (gcd(a(n),n) and its largest prime factor); A179931, A192795 (gcd(a(n),A000203(n)) and largest prime factor).

%Y Main diagonal of the array in A242639.

%K nonn,core,nice,easy,mult

%O 1,2

%A _N. J. A. Sloane_, _R. K. Guy_

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 18 07:19 EDT 2019. Contains 328146 sequences. (Running on oeis4.)