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-8 of 8 results.

A282097 Coefficients in q-expansion of (3*E_2*E_4 - 2*E_6 - E_2^3)/1728, where E_2, E_4, E_6 are the Eisenstein series shown in A006352, A004009, A013973, respectively.

Original entry on oeis.org

0, 1, 12, 36, 112, 150, 432, 392, 960, 1053, 1800, 1452, 4032, 2366, 4704, 5400, 7936, 5202, 12636, 7220, 16800, 14112, 17424, 12696, 34560, 19375, 28392, 29160, 43904, 25230, 64800, 30752, 64512, 52272, 62424, 58800, 117936, 52022, 86640, 85176, 144000, 70602
Offset: 0

Views

Author

Seiichi Manyama, Feb 06 2017

Keywords

Comments

Multiplicative because A000203 is. - Andrew Howroyd, Jul 25 2018

Examples

			a(6) = 1^3*6^2 + 2^3*3^2 + 3^3*2^2 + 6^3*1^2 = 432.

Links

Crossrefs

Cf. this sequence (phi_{3, 2}), A282099 (phi_{5, 2}).
Cf. A006352 (E_2), A004009 (E_4), A013973 (E_6), A282018 (E_2^3), A282019 (E_2*E_4).
Cf. A000203 (sigma(n)), A064987 (n*sigma(n)), this sequence (n^2*sigma(n)), A282211 (n^3*sigma(n)).
Cf. A222171.

Programs

  • Magma
    [0] cat [n^2*DivisorSigma(1, n): n in [1..50]]; // Vincenzo Librandi, Mar 01 2018
  • Mathematica
    a[0]=0;a[n_]:=(n^2)*DivisorSigma[1,n];Table[a[n],{n,0,41}] (* Indranil Ghosh, Feb 21 2017 *)
terms = 42; Ei[n_] = 1-(2n/BernoulliB[n]) Sum[k^(n-1) x^k/(1-x^k), {k, terms}]; CoefficientList[(3*Ei[2]*Ei[4] - 2*Ei[6] - Ei[2]^3)/1728 + O[x]^terms, x] (* Jean-François Alcover, Mar 01 2018 *)
  • PARI
    a(n) = if (n==0, 0, n^2*sigma(n)); \\ Michel Marcus, Feb 21 2017

Formula

a(n) = (3*A282019(n) - 2*A013973(n) - A282018(n))/1728.
G.f.: phi_{3, 2}(x) where phi_{r, s}(x) = Sum_{n, m>0} m^r * n^s * x^{m*n}.
a(n) = n^2*A000203(n) for n > 0. - Seiichi Manyama, Feb 19 2017
G.f.: Sum_{k>=1} k^3*x^k*(1 + x^k)/(1 - x^k)^3. - Ilya Gutkovskiy, May 02 2018
From Amiram Eldar, Oct 30 2023: (Start)
Multiplicative with a(p^e) = p^(2*e) * (p^(e+1)-1)/(p-1).
Dirichlet g.f.: zeta(s-2)*zeta(s-3).
Sum_{k=1..n} a(k) ~ (Pi^2/24) * n^4. (End)
From Peter Bala, Jan 22 2024: (Start)
a(n) = Sum_{1 <= i, j, k <= n} sigma_2( gcd(i, j, k, n) ).
a(n) = Sum_{1 <= i, j <= n} sigma_3( gcd(i, j, n) ).
a(n) = Sum_{d divides n} sigma_2(d) * J_3(n/d) = Sum_{d divides n} sigma_3(d) * J_2(n/d), where the Jordan totient functions J_2(n) = A007434(n) and J_3(n) = A059376(n). (End)

A280022 Expansion of phi_{5, 4}(x) where phi_{r, s}(x) = Sum_{n, m>0} m^r * n^s * x^{m*n}.

Original entry on oeis.org

0, 1, 48, 324, 1792, 3750, 15552, 19208, 61440, 85293, 180000, 175692, 580608, 399854, 921984, 1215000, 2031616, 1503378, 4094064, 2606420, 6720000, 6223392, 8433216, 6716184, 19906560, 12109375, 19192992, 21257640, 34420736, 21218430, 58320000, 29552672
Offset: 0

Views

Author

Seiichi Manyama, Feb 22 2017

Keywords

Comments

Multiplicative because A000203 is. - Andrew Howroyd, Jul 23 2018

Links

Crossrefs

Cf. this sequence (phi_{5, 4}), A280025 (phi_{7, 4}).
Cf. A282101 (E_2*E_4^2), A282595 (E_2^2*E_6), A282586 (E_2^3*E_4), A013974 (E_4*E_6 = E_10), A282431 (E_2^5).
Cf. A000203 (sigma(n)), A064987 (n*sigma(n)), A282097 (n^2*sigma(n)), A282211 (n^3*sigma(n)), this sequence (n^4*sigma(n)).
Cf. A353908.

Programs

  • Mathematica
    Table[n^4 * DivisorSigma[1, n], {n, 0, 32}] (* Amiram Eldar, Oct 31 2023 *)
nmax = 30; CoefficientList[Series[Sum[k^4*x^k*(1 + 26*x^k + 66*x^(2*k) + 26*x^(3*k) + x^(4*k))/(1 - x^k)^6, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 02 2025 *)
  • PARI
    a(n) = if(n < 1, 0, n^4 * sigma(n)); \\ Andrew Howroyd, Jul 23 2018

Formula

a(n) = n^4*A000203(n) for n > 0.
a(n) = (15*A282101(n) - 20*A282595(n) + 10*A282586(n) - 4*A013974(n) - A282431(n))/20736.
Sum_{k=1..n} a(k) ~ c * n^6, where c = Pi^2/36 = 0.274155... (A353908). - Amiram Eldar, Dec 08 2022
From Amiram Eldar, Oct 31 2023: (Start)
Multiplicative with a(p^e) = p^(4*e) * (p^(e+1)-1)/(p-1).
Dirichlet g.f.: zeta(s-4)*zeta(s-5). (End)
G.f.: Sum_{k>=1} k^4*x^k*(1 + 26*x^k + 66*x^(2*k) + 26*x^(3*k) + x^(4*k))/(1 - x^k)^6. - Vaclav Kotesovec, Aug 02 2025

A282213 Coefficients in q-expansion of (E_2^3*E_4 - 3*E_2^2*E_6 + 3*E_2*E_4^2 - E_4*E_6)/3456, where E_2, E_4, E_6 are the Eisenstein series shown in A006352, A004009, A013973, respectively.

Original entry on oeis.org

0, 1, 72, 756, 4672, 15750, 54432, 117992, 299520, 551853, 1134000, 1772892, 3532032, 4829006, 8495424, 11907000, 19173376, 24142482, 39733416, 47052740, 73584000, 89201952, 127648224, 148048056, 226437120, 246109375, 347688432, 402320520, 551258624, 594847710
Offset: 0

Views

Author

Seiichi Manyama, Feb 09 2017

Keywords

Comments

Multiplicative because A001158 is. - Andrew Howroyd, Jul 25 2018

Examples

			a(6) = 1^6*6^3 + 2^6*3^3 + 3^6*2^3 + 6^6*1^3 = 54432.

Links

Crossrefs

Cf. A282211 (phi_{4, 3}), this sequence (phi_{6, 3}).
Cf. A006352 (E_2), A004009 (E_4), A013973 (E_6), A282586 (E_2^3*E_4), A282595 (E_2^2*E_6), A282101 (E_2*E_4^2), A013974 (E_4*E_6 = E_10).
Cf. A001158 (sigma_3(n)), A281372 (n*sigma_3(n)), A282099 (n^2*sigma_3(n)), this sequence (n^3*sigma_3(n))
Cf. A013662.

Programs

  • Mathematica
    terms = 30;
E2[x_] = 1 - 24*Sum[k*x^k/(1 - x^k), {k, 1, terms}];
E4[x_] = 1 + 240*Sum[k^3*x^k/(1 - x^k), {k, 1, terms}];
E6[x_] = 1 - 504*Sum[k^5*x^k/(1 - x^k), {k, 1, terms}];
(E2[x]^3*E4[x] - 3 E2[x]^2*E6[x] + 3 E2[x] E4[x]^2 - E4[x] E6[x])/3456 + O[x]^terms // CoefficientList[#, x]&
(* or: *)
Table[n^3*DivisorSigma[3, n], {n, 0, terms-1}] (* Jean-François Alcover, Feb 27 2018 *)
nmax = 30; CoefficientList[Series[Sum[k^6*x^k*(x^(2*k) + 4*x^k + 1)/(1 - x^k)^4, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 02 2025 *)
  • PARI
    a(n) = if (n, n^3*sigma(n, 3), 0); \\ Michel Marcus, Feb 27 2018

Formula

G.f.: phi_{6, 3}(x) where phi_{r, s}(x) = Sum_{n, m>0} m^r * n^s * x^{m*n}.
a(n) = (A282586(n) - 3*A282595(n) + 3*A282101(n) - A013974(n))/3456. - Seiichi Manyama, Feb 19 2017
a(n) = n^3*A001158(n) for n > 0. - Seiichi Manyama, Feb 19 2017
Sum_{k=1..n} a(k) ~ zeta(4) * n^7 / 7. - Amiram Eldar, Sep 06 2023
From Amiram Eldar, Oct 31 2023: (Start)
Multiplicative with a(p^e) = p^(3*e) * (p^(3*e+3)-1)/(p^3-1).
Dirichlet g.f.: zeta(s-3)*zeta(s-6). (End)
G.f.: Sum_{k>=1} k^6*x^k*(x^(2*k) + 4*x^k + 1)/(1 - x^k)^4. - Vaclav Kotesovec, Aug 02 2025

A282781 Expansion of phi_{8, 3}(x) where phi_{r, s}(x) = Sum_{n, m>0} m^r * n^s * x^{m*n}.

Original entry on oeis.org

0, 1, 264, 6588, 67648, 390750, 1739232, 5765144, 17318400, 43224597, 103158000, 214360212, 445665024, 815732918, 1521998016, 2574261000, 4433514496, 6975762354, 11411293608, 16983569900, 26433456000, 37980768672, 56591095968, 78310997448
Offset: 0

Views

Author

Seiichi Manyama, Feb 21 2017

Keywords

Comments

Multiplicative because A001160 is. - Andrew Howroyd, Jul 25 2018

Links

Crossrefs

Cf. A282211 (phi_{4, 3}), A282213 (phi_{6, 3}), this sequence (phi_{8, 3}).
Cf. A282752 (E_2^2*E_4^2), A282780 (E_2^3*E_6), A282102 (E_2*E_4*E_6), A008411 (E_4^3), A280869 (E_6^2).
Cf. A001160 (sigma_5(n)), A282050 (n*sigma_5(n)), A282751 (n^2*sigma_5(n)), this sequence (n^3*sigma_5(n)).
Cf. A013664.

Programs

  • Mathematica
    a[0]=0;a[n_]:=(n^3)*DivisorSigma[5,n];Table[a[n],{n,0,23}] (* Indranil Ghosh, Feb 21 2017 *)
nmax = 30; CoefficientList[Series[Sum[k^8*x^k*(x^(2*k) + 4*x^k + 1)/(1 - x^k)^4, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 02 2025 *)
  • PARI
    a(n) = if (n==0, 0, n^3*sigma(n, 5)); \\ Michel Marcus, Feb 21 2017

Formula

a(n) = n^3*A001160(n) for n > 0.
a(n) = (6*A282752(n) - 2*A282780(n) - 6*A282102(n) + A008411(n) + A280869(n))/5184.
Sum_{k=1..n} a(k) ~ zeta(6) * n^9 / 9. - Amiram Eldar, Sep 06 2023
From Amiram Eldar, Oct 31 2023: (Start)
Multiplicative with a(p^e) = p^(3*e) * (p^(5*e+5)-1)/(p^5-1).
Dirichlet g.f.: zeta(s-3)*zeta(s-8). (End)
G.f.: Sum_{k>=1} k^8*x^k*(x^(2*k) + 4*x^k + 1)/(1 - x^k)^4. - Vaclav Kotesovec, Aug 02 2025

A386781 a(n) = n^3*sigma_7(n).

Original entry on oeis.org

0, 1, 1032, 59076, 1056832, 9765750, 60966432, 282475592, 1082196480, 3488379453, 10078254000, 25937425932, 62433407232, 137858494046, 291514810944, 576921447000, 1108169199616, 2015993905362, 3600007595496, 6131066264660, 10320757104000, 16687528072992, 26767423561824
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 02 2025

Keywords

Links

Crossrefs

Cf. A282211, A386746, A282213, A386748, A282781, A386780, A386782.
Cf. A013955, A282060, A282753.
Cf. A386813, A282549, A282792, A058550, A282576.

Programs

  • Magma
    [0] cat [n^3*DivisorSigma(7, n): n in [1..35]]; // Vincenzo Librandi, Aug 04 2025
  • Mathematica
    Table[n^3*DivisorSigma[7, n], {n, 0, 30}]
(* or *)
nmax = 30; CoefficientList[Series[Sum[k^10*x^k*(x^(2*k) + 4*x^k + 1)/(1 - x^k)^4, {k, 1, nmax}], {x, 0, nmax}], x]
(* or *)
terms = 30; E2[x_] = 1 - 24*Sum[k*x^k/(1 - x^k), {k, 1, terms}]; E4[x_] = 1 + 240*Sum[k^3*x^k/(1 - x^k), {k, 1, terms}]; E6[x_] = 1 - 504*Sum[k^5*x^k/(1 - x^k), {k, 1, terms}]; CoefficientList[Series[(3*E2[x]^3*E4[x]^2 + 5*E2[x]*E4[x]^3 - 9*E2[x]^2*E4[x]*E6[x] - 3*E4[x]^2*E6[x] + 4*E2[x]*E6[x]^2)/3456, {x, 0, terms}], x]

Formula

G.f.: Sum_{k>=1} k^10*x^k*(x^(2*k) + 4*x^k + 1)/(1 - x^k)^4.
a(n) = (3*A386813(n) + 5*A282549(n) - 9*A282792(n) - 3*A058550(n) + 4*A282576(n))/3456.
a(n) = n^3*A013955(n).
Dirichlet g.f.: zeta(s-3)*zeta(s-10). - R. J. Mathar, Aug 03 2025

A386780 a(n) = n^3*sigma_6(n).

Original entry on oeis.org

0, 1, 520, 19710, 266304, 1953250, 10249200, 40353950, 136348160, 387952659, 1015690000, 2357949022, 5248851840, 10604501570, 20984054000, 38498557500, 69810262016, 118587881410, 201735382680, 322687704638, 520158288000, 795376354500, 1226133491440, 1801152673630
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 02 2025

Keywords

Links

Crossrefs

Cf. A282211, A386746, A282213, A386748, A282781, A386781, A386782.
Cf. A013954, A386750, A386777.

Programs

  • Magma
    [0] cat [n^3*DivisorSigma(6, n): n in [1..35]]; // Vincenzo Librandi, Aug 04 2025
  • Mathematica
    Table[n^3*DivisorSigma[6, n], {n, 0, 30}]
nmax = 30; CoefficientList[Series[Sum[k^9*x^k*(x^(2*k) + 4*x^k + 1)/(1 - x^k)^4, {k, 1, nmax}], {x, 0, nmax}], x]

Formula

G.f.: Sum_{k>=1} k^9*x^k*(x^(2*k) + 4*x^k + 1)/(1 - x^k)^4.
a(n) = n^3*A013954(n).
Dirichlet g.f.: zeta(s-3)*zeta(s-9). - R. J. Mathar, Aug 03 2025

A386782 a(n) = n^3*sigma_8(n).

Original entry on oeis.org

0, 1, 2056, 177174, 4210752, 48828250, 364269744, 1977327086, 8623620608, 31385843307, 100390882000, 285311671942, 746035774848, 1792160396234, 4065384488816, 8651096365500, 17661175009280, 34271896312546, 64529293839192, 116490258905078, 205603651344000, 350330949134964
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 02 2025

Keywords

Links

Crossrefs

Cf. A282211, A386746, A282213, A386748, A282781, A386780, A386781.
Cf. A013956, A386751, A386778.

Programs

  • Magma
    [0] cat [n^3*DivisorSigma(8, n): n in [1..35]]; // Vincenzo Librandi, Aug 04 2025
  • Mathematica
    Table[n^3*DivisorSigma[8, n], {n, 0, 30}]
nmax = 30; CoefficientList[Series[Sum[k^11*x^k*(x^(2*k) + 4*x^k + 1)/(1 - x^k)^4, {k, 1, nmax}], {x, 0, nmax}], x]

Formula

G.f.: Sum_{k>=1} k^11*x^k*(x^(2*k) + 4*x^k + 1)/(1 - x^k)^4.
a(n) = n^3*A013956(n).
Dirichlet g.f.: zeta(s-3)*zeta(s-11). - R. J. Mathar, Aug 03 2025

A364194 a(n) = Sum_{k=1..n} k^3*sigma(k), where sigma is A000203.

Original entry on oeis.org

1, 25, 133, 581, 1331, 3923, 6667, 14347, 23824, 41824, 57796, 106180, 136938, 202794, 283794, 410770, 499204, 726652, 863832, 1199832, 1496184, 1879512, 2171520, 3000960, 3485335, 4223527, 5010847, 6240159, 6971829, 8915829, 9869141, 11933525, 13658501
Offset: 1

Views

Author

Seiichi Manyama, Oct 20 2023

Keywords

Links

Crossrefs

Partial sums of A282211.
Cf. A319086, A320895, A356125, A364269.
Cf. A000203, A000537, A143128.

Programs

  • Mathematica
    Accumulate[Table[n^3*DivisorSigma[1, n], {n, 1, 33}]] (* Amiram Eldar, Oct 20 2023 *)
  • PARI
    f(n, m) = (subst(bernpol(m+1, x), x, n+1)-subst(bernpol(m+1, x), x, 0))/(m+1);
a(n, s=3, t=1) = sum(k=1, n, k^(s+t)*f(n\k, s));
  • Python
    def A364194(n): return sum((k**2*(m:=n//k)*(m+1)>>1)**2 for k in range(1,n+1)) # Chai Wah Wu, Oct 20 2023
  • Python
    from math import isqrt
def A364194(n): return ((((s:=isqrt(n))*(s + 1))**3*(2*s+1)*(1-3*s*(s+1))>>1) + sum((q:=n//k)*(q+1)*k**3*(q*(15*k+q*(15*k+12*q+18)+2)-2) for k in range(1,s+1)))//60 # Chai Wah Wu, Oct 21 2023

Formula

a(n) = Sum_{k=1..n} k^4 * A000537(floor(n/k)).
a(n) ~ (zeta(2)/5) * n^5. - Amiram Eldar, Oct 20 2023
Showing 1-8 of 8 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.