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.

Showing 1-10 of 15 results. Next

A035116 a(n) = tau(n)^2, where tau(n) = A000005(n).

Original entry on oeis.org

1, 4, 4, 9, 4, 16, 4, 16, 9, 16, 4, 36, 4, 16, 16, 25, 4, 36, 4, 36, 16, 16, 4, 64, 9, 16, 16, 36, 4, 64, 4, 36, 16, 16, 16, 81, 4, 16, 16, 64, 4, 64, 4, 36, 36, 16, 4, 100, 9, 36, 16, 36, 4, 64, 16, 64, 16, 16, 4, 144, 4, 16
Offset: 1

Views

Author

Keywords

References

  • G. H. Hardy, Ramanujan: twelve lectures on subjects suggested by his life and work, Cambridge, University Press, 1940, p. 59.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, Theorem 304.

Crossrefs

Programs

Formula

Dirichlet g.f.: zeta(s)^4/zeta(2s).
tau(n)^2 = Sum_{d|n} tau(d^2), Dirichlet convolution of A048691 and A000012 (i.e.: inverse Mobius transform of A048691).
Multiplicative with a(p^e) = (e+1)^2. - Vladeta Jovovic, Dec 03 2001
G.f.: Sum_{n>=1} A000005(n^2)*x^n/(1-x^n). - Mircea Merca, Feb 25 2014
a(n) = A066446(n) + A184389(n). - Reinhard Zumkeller, Sep 08 2015
Let b(n), n > 0, be the Dirichlet inverse of a(n). Then b(n) is multiplicative with b(p^e) = (-1)^e*(Sum_{i=0..e} binomial(3,i)) for prime p and e >= 0, where binomial(n,k)=0 if n < k; abs(b(n)) is multiplicative and has the Dirichlet g.f.: (zeta(s))^4/(zeta(2*s))^3. - Werner Schulte, Feb 07 2021

Extensions

Additional comments from Vladeta Jovovic, Apr 29 2001

A057434 a(n) = Sum_{k=1..n} phi(k)^2.

Original entry on oeis.org

1, 2, 6, 10, 26, 30, 66, 82, 118, 134, 234, 250, 394, 430, 494, 558, 814, 850, 1174, 1238, 1382, 1482, 1966, 2030, 2430, 2574, 2898, 3042, 3826, 3890, 4790, 5046, 5446, 5702, 6278, 6422, 7718, 8042, 8618, 8874, 10474, 10618, 12382, 12782
Offset: 1

Views

Author

N. J. A. Sloane, Sep 08 2000

Keywords

Comments

Partial sums of A127473. - R. J. Mathar, Sep 29 2008

Crossrefs

Programs

  • Mathematica
    FoldList[Plus, 1, EulerPhi[Range[2, 50]]^2] (* Ivan Neretin, May 30 2015 *)
  • PARI
    a(n) = sum(k=1, n, eulerphi(k)^2); \\ Michel Marcus, Dec 20 2015

Formula

We can derive an asymptotic formula from a general formula given in the reference, namely: a(n) = C*n^3 + O(log(x)^(4/3)log(log(x))^(8/3)) where C = (1/3)/zeta(2)^2*Product_{p prime}(1+1/(p-1)/(p+1)^2) = 0.142749835225698(...). - Benoit Cloitre, Dec 22 2015
a(n) ~ c * n^3 / 3, where c = A065464 = Product_{primes p} (1 - 2/p^2 + 1/p^3) = 0.4282495056770944402187657075818235461212985133559361440319... - Vaclav Kotesovec, Dec 18 2019

A072379 Sum_{k<=n} (sigma(k)^2), where sigma(k) denotes the sum of the divisors of k A000203.

Original entry on oeis.org

1, 10, 26, 75, 111, 255, 319, 544, 713, 1037, 1181, 1965, 2161, 2737, 3313, 4274, 4598, 6119, 6519, 8283, 9307, 10603, 11179, 14779, 15740, 17504, 19104, 22240, 23140, 28324, 29348, 33317, 35621, 38537, 40841, 49122, 50566, 54166, 57302
Offset: 1

Views

Author

Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Jul 20 2002

Keywords

Crossrefs

Programs

  • Maple
    A072379 := proc(n)
        add( numtheory[sigma](k)^2,k=0..n) ;
    end proc:
    seq(A072379(n),n=1..80) ; # R. J. Mathar, Jul 09 2024
  • Mathematica
    Accumulate[Table[DivisorSigma[1, k]^2, {k, 1, 50}]] (* Vaclav Kotesovec, Sep 10 2018 *)
  • PARI
    a(n) = sum(k=1, n, sigma(k)^2) \\ Michel Marcus, Jun 20 2013

Formula

Ramanujan's asymptotic formula: (5/6)*Zeta(3)*n^3+O(n^2*log(n)^2)

A061503 a(n) = Sum_{k=1..n} tau(k^2), where tau is the number of divisors function A000005.

Original entry on oeis.org

1, 4, 7, 12, 15, 24, 27, 34, 39, 48, 51, 66, 69, 78, 87, 96, 99, 114, 117, 132, 141, 150, 153, 174, 179, 188, 195, 210, 213, 240, 243, 254, 263, 272, 281, 306, 309, 318, 327, 348, 351, 378, 381, 396, 411, 420, 423, 450, 455, 470, 479, 494, 497
Offset: 1

Views

Author

N. J. A. Sloane, Jun 14 2001

Keywords

Comments

a(n) is the number of pairs of positive integers <= n with their LCM <= n. - Andrew Howroyd, Sep 01 2019

References

  • Mentioned by Steven Finch in a posting to the Number Theory List (NMBRTHRY(AT)LISTSERV.NODAK.EDU), Jun 13 2001.

Crossrefs

Cf. A000005, A061502. Partial sums of A048691.

Programs

  • GAP
    List([1..60],n->Sum([1..n],k->Tau(k^2))); # Muniru A Asiru, Mar 09 2019
    
  • Maple
    with(numtheory): a:=n->add(tau(k^2),k=1..n): seq(a(n),n=1..60); # Muniru A Asiru, Mar 09 2019
  • Mathematica
    DivisorSigma[0, Range[60]^2] // Accumulate (* Jean-François Alcover, Nov 25 2013 *)
  • PARI
    for (n=1, 1024, write("b061503.txt", n, " ", sum(k=1, n, numdiv(k^2)))) \\ Harry J. Smith, Jul 23 2009
    
  • PARI
    t=0;v=vector(60,n,t+=numdiv(n^2)) \\ Charles R Greathouse IV, Nov 08 2012
    
  • Python
    from math import prod
    from sympy import factorint
    def A061503(n): return sum(prod(2*e+1 for e in factorint(k).values()) for k in range(1,n+1)) # Chai Wah Wu, May 10 2022
  • Sage
    def A061503(n) :
        tau = sloane.A000005
        return add(tau(k^2) for k in (1..n))
    [ A061503(i) for i in (1..19)] # Peter Luschny, Sep 15 2012
    

Formula

a(n) = Sum_{j=1..n^2} floor(n/A019554(j)). - Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Jul 20 2002
a(n) = Sum_{i=1..n} 2^omega(i) * floor(n/i). - Enrique Pérez Herrero, Sep 15 2012
a(n) ~ 3/Pi^2 * n log^2 n. - Charles R Greathouse IV, Nov 08 2012
a(n) ~ 3*n/Pi^2 * (log(n)^2 + log(n)*(-2 + 6*g - 24*z/Pi^2) + 2 - 6*g + 6*g^2 - 6*sg1 + 288*z^2/Pi^4 - 24*(-z + 3*g*z + z2)/ Pi^2), where g is the Euler-Mascheroni constant A001620, sg1 is the first Stieltjes constant (see A082633), z = Zeta'(2) (see A073002), z2 = Zeta''(2) = A201994. - Vaclav Kotesovec, Jan 30 2019
a(n) = Sum_{k=1..n} A064608(floor(n/k)). - Daniel Suteu, Mar 09 2019

Extensions

Name corrected by Peter Luschny, Sep 15 2012

A318755 a(n) = Sum_{k=1..n} tau(k)^3, where tau is A000005.

Original entry on oeis.org

1, 9, 17, 44, 52, 116, 124, 188, 215, 279, 287, 503, 511, 575, 639, 764, 772, 988, 996, 1212, 1276, 1340, 1348, 1860, 1887, 1951, 2015, 2231, 2239, 2751, 2759, 2975, 3039, 3103, 3167, 3896, 3904, 3968, 4032, 4544, 4552, 5064, 5072, 5288, 5504, 5568, 5576
Offset: 1

Views

Author

Vaclav Kotesovec, Sep 02 2018

Keywords

Crossrefs

Programs

  • Mathematica
    Accumulate[DivisorSigma[0, Range[50]]^3]
  • PARI
    a(n) = sum(k=1, n, numdiv(k)^3); \\ Michel Marcus, Sep 03 2018

Formula

a(n) ~ n * (A1*log(n)^7 + A2*log(n)^6 + A3*log(n)^5 + A4*log(n)^4 + A5*log(n)^3 + A6*log(n)^2 + A7*log(n) + A8) [Ramanujan, 1916, formula (8)].
From Vaclav Kotesovec, Mar 12 2023: (Start)
Let f(s) = Product_{p prime} (1 - 9/p^(2*s) + 16/p^(3*s) - 9/p^(4*s) + 1/p^(6*s)), then
A1 = f(1)/5040 = 0.0000097860463451190658257888710490039661018239924009134296302566263529129...
A2 = ((8*gamma - 1)*f(1) + f'(1)) / 720 = 0.0007019997226174095261771358653540021199703406583347258622085873074052900...
A3 = (2 * f(1) * (1 - 8*gamma + 28*gamma^2 - 8*sg1) + 2*(8*gamma - 1)*f'(1) + f''(1)) / 240 = 0.0171707557268638504150726777646428533953516776541779590118582753709080243...
A4 = (6*f(1)*(-1 - 28*gamma^2 + 56*gamma^3 + gamma*(8 - 56*sg1) + 8*sg1 + 4*sg2) + 6*(1 - 8*gamma + 28*gamma^2 - 8*sg1)*f'(1) + (24*gamma - 3)*f''(1) + f'''(1)) / 144 = 0.1758477246705824231478998937203303065702508974398264386862202155788...,
where f(1) = Product_{p prime} (1 - 9/p^2 + 16/p^3 - 9/p^4 + 1/p^6) = 0.0493216735794000917619759100869799891531929217006036853364933968186814900...,
f'(1) = f(1) * Sum_{p prime} 6*(3*p + 1) * log(p) / ((p-1) * (p^2 + 4*p + 1)) = 0.3270075329904166293296173488834535949530448497141635531152019426434776932...,
f''(1) = f'(1)^2 / f(1) + f(1) * Sum_{p prime} -36 * p^2 * (p+1)^2 * log(p)^2 / ((p-1)^2 * (p^2 + 4*p + 1)^2) = 1.1340946589859924227356699847227569935993284591079455746283572890834872890...,
f'''(1) = 3*f'(1)*f''(1)/f(1) - 2*f'(1)^3/f(1)^2 + f(1) * Sum_{p prime} 72*p^2 * (p^5 + 3*p^4 + 8*p^3 + 8*p^2 + 3*p+ 1) * log(p)^3 / ((p-1)^3 * (p^2+ 4*p + 1)^3) = -1.3447542210274297874241826540796632006263184659735145444999327537246287...,
gamma is the Euler-Mascheroni constant A001620 and sg1, sg2 are the Stieltjes constants, see A082633 and A086279.
Approximate values of other constants:
A5 = 0.7626157870664479996781152281270580148665443022014605423466363134512...
A6 = 1.3720912878905940866975369743071441424192833481004753922122458993040...
A7 = 1.1416118168318711437057727816148048057614284471759625288073915723140...
A8 = 0.2618221765943171424958051160111945242076019991649774700610674747694...
(End)

A074789 Partial sums of usigma(n)^2: square of the sum of unitary divisors of n.

Original entry on oeis.org

1, 10, 26, 51, 87, 231, 295, 376, 476, 800, 944, 1344, 1540, 2116, 2692, 2981, 3305, 4205, 4605, 5505, 6529, 7825, 8401, 9697, 10373, 12137, 12921, 14521, 15421, 20605, 21629, 22718, 25022, 27938, 30242, 32742, 34186, 37786, 40922, 43838
Offset: 1

Views

Author

Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Sep 07 2002

Keywords

Crossrefs

Programs

  • Mathematica
    Accumulate[Table[DivisorSum[n, # &, CoprimeQ[#, n/#] &]^2, {n, 1, 50}]] (* Vaclav Kotesovec, Feb 04 2019 *)
  • PARI
    A034448(n) = {my(f = factor(n)); prod(i=1, #f~, 1 + f[i, 1]^f[i, 2]);}
    lista(nmax) = {my(s = 0); for(n = 1, nmax, s += A034448(n)^2; print1(s, ", "));} \\ Amiram Eldar, Jul 24 2024

Formula

a(n) = Sum_{k=1..n} usigma(k)^2 = Sum_{k=1..n} A034448(k)^2.
Asymptotic expression: a(n) = Sum_{k<=n} usigma(k)^2 = (zeta(2)*zeta(3)*alpha_1/3)*n^3 + O(n^2*log(n)^4), alpha_1 = Product_{p prime} (1+1/p^2-2/p^3-1/p^4-2/p^5+3/p^6), zeta(2) = A013661 and zeta(3) = A002117.
alpha_1 = 1.001619936509160661474009830789... . - Amiram Eldar, Jul 24 2024

A245074 Decimal expansion of B, the coefficient of n*log(n)^2 in the asymptotic formula of Ramanujan for Sum_{k=1..n} (d(k)^2), where d(k) is the number of distinct divisors of k.

Original entry on oeis.org

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

Views

Author

Jean-François Alcover, Jul 11 2014

Keywords

Comments

The coefficient of n*log(n)^3 in the same asymptotic formula is A = 1/Pi^2.

Examples

			0.744341276391456640439006036785694615691377808839427047585292094877364...
		

References

  • Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, Section Sierpinski's Constant, p. 124.

Crossrefs

Programs

  • Mathematica
    B = (12*EulerGamma - 3)/Pi^2 - (36/Pi^4)*Zeta'[2]; RealDigits[B, 10, 103] // First
  • PARI
    (4*Euler-1)/(2*zeta(2)) - zeta'(2)/zeta(2)^2 \\ Amiram Eldar, Apr 27 2025

Formula

B = (12*gamma - 3)/Pi^2 - (36/Pi^4)*zeta'(2).

A319090 Decimal expansion of C, the coefficient of n*log(n) in the asymptotic formula of Ramanujan for Sum_{k=1..n} (d(k)^2), where d(k) is the number of distinct divisors of k.

Original entry on oeis.org

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

Views

Author

Vaclav Kotesovec, Sep 10 2018

Keywords

Examples

			0.823265208269485020156816453947090406301273270321142250892524579208530395971755...
		

Crossrefs

Programs

  • Mathematica
    36*EulerGamma^2/Pi^2 - (288*Zeta'[2]/Pi^4 + 24/Pi^2)*EulerGamma + (864*Zeta'[2]^2/Pi^6 + 72*Zeta'[2]/Pi^4 - 72/Pi^4*Zeta''[2] + 6/Pi^2) - 24*StieltjesGamma[1]/Pi^2

Formula

C = 36*gamma^2/Pi^2 - (288*z1/Pi^4 + 24/Pi^2)*gamma + (864*z1^2/Pi^6 + 72*z1/Pi^4 - 72/Pi^4*z2 + 6/Pi^2) - 24*g1/Pi^2, where gamma is the Euler-Mascheroni constant A001620, z1 = Zeta'(2) = A073002, z2 = Zeta''(2) = A201994 and g1 is the first Stieltjes constant, see A082633.

A319091 Decimal expansion of D, the coefficient of n in the asymptotic formula of Ramanujan for Sum_{k=1..n} (d(k)^2), where d(k) is the number of distinct divisors of k.

Original entry on oeis.org

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

Views

Author

Vaclav Kotesovec, Sep 10 2018

Keywords

Examples

			0.4603233722587214303937620863844189747632149035387392240584250348445902629324...
		

Crossrefs

Programs

  • Mathematica
    24*EulerGamma^3/Pi^2 - (432*Zeta'[2] /Pi^4+ 36/Pi^2)*EulerGamma^2 + (3456*Zeta'[2]^2/Pi^6 + 288*(Zeta'[2]-Zeta''[2])/Pi^4 + 24/Pi^2 - 72*StieltjesGamma[1]/Pi^2)*EulerGamma + StieltjesGamma[1]*(288*Zeta'[2]/Pi^4 + 24/Pi^2)-10368*Zeta'[2]^3/Pi^8 - 864*Zeta'[2]^2/Pi^6 + 1728*Zeta''[2] * Zeta'[2]/Pi^6 + 72*(Zeta''[2]-Zeta'[2])/Pi^4 - 48*Zeta'''[2]/Pi^4 + (12*StieltjesGamma[2] - 6)/Pi^2

Formula

D = 24*gamma^3/Pi^2 - (432*z1 /Pi^4+ 36/Pi^2)*gamma^2 + (3456*z1^2/Pi^6 + 288*(z1-z2)/Pi^4 + 24/Pi^2 - 72*g1/Pi^2)*gamma + g1*(288*z1/Pi^4 + 24/Pi^2)-10368*z1^3/Pi^8 - 864*z1^2/Pi^6 + 1728*z2*z1/Pi^6 + 72*(z2-z1)/Pi^4- 48*z3/Pi^4 + (12*g2-6)/Pi^2, where gamma is the Euler-Mascheroni constant A001620, z1 = Zeta'(2) = A073002, z2 = Zeta''(2) = A201994, z3 = Zeta'''(2) = A201995 and g1, g2 are the Stieltjes constants, see A082633 and A086279.

A320896 a(n) = Sum_{k=1..n} k * tau(k)^2, where tau is A000005.

Original entry on oeis.org

1, 9, 21, 57, 77, 173, 201, 329, 410, 570, 614, 1046, 1098, 1322, 1562, 1962, 2030, 2678, 2754, 3474, 3810, 4162, 4254, 5790, 6015, 6431, 6863, 7871, 7987, 9907, 10031, 11183, 11711, 12255, 12815, 15731, 15879, 16487, 17111, 19671, 19835, 22523, 22695, 24279
Offset: 1

Views

Author

Vaclav Kotesovec, Oct 23 2018

Keywords

Crossrefs

Programs

  • Mathematica
    Accumulate[Table[k*DivisorSigma[0, k]^2, {k, 1, 50}]]
  • PARI
    a(n) = sum(k=1, n, k*numdiv(k)^2); \\ Michel Marcus, Oct 23 2018

Formula

a(n) ~ n^2 * (3*(Pi^6*(-1 - 24*g^2 + 32*g^3 + g*(8 - 96*s1) + 16*s1 + 16*s2) - 13824*z1^3 + 576*Pi^2*z1*((-1 + 8*g)*z1 + 4*z2) - 8*Pi^4*(3*(1 - 8*g + 24*g^2 - 16*s1)*z1 - 6*z2 + 48*g*z2 + 8*z3)) + 6*(Pi^6*(1 - 8*g + 24*g^2 - 16*s1) + 576*Pi^2*z1^2 - 24*Pi^4*(-z1 + 8*g*z1 + 2*z2))*log(n) + 6*((-1 + 8*g)*Pi^6 - 24*Pi^4*z1)*log(n)^2 + 4*Pi^6*log(n)^3) / (8*Pi^8), where g is the Euler-Mascheroni constant A001620, z1 = Zeta'(2) = A073002, z2 = Zeta''(2) = A201994, z3 = Zeta'''(2) = A201995 and s1, s2 are the Stieltjes constants, see A082633 and A086279.
Showing 1-10 of 15 results. Next