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 115 results. Next

A281372 Coefficients in q-expansion of (E_2*E_4 - E_6)/720, where E_2, E_4, E_6 are the Eisenstein series shown in A006352, A004009, A013973, respectively.

Original entry on oeis.org

0, 1, 18, 84, 292, 630, 1512, 2408, 4680, 6813, 11340, 14652, 24528, 28574, 43344, 52920, 74896, 83538, 122634, 130340, 183960, 202272, 263736, 279864, 393120, 393775, 514332, 551880, 703136, 707310, 952560, 923552, 1198368, 1230768, 1503684, 1517040, 1989396, 1874198, 2346120, 2400216, 2948400
Offset: 0

Views

Author

N. J. A. Sloane, Feb 04 2017

Keywords

Comments

The q-expansion of the square of this expression is given in A281371.
Multiplicative because A001158 is. - Andrew Howroyd, Jul 23 2018

Crossrefs

Programs

  • Magma
    [0] cat [n*DivisorSigma(3, n): n in [1..50]]; // Vincenzo Librandi, Mar 01 2018
  • Maple
    with(gfun):
    with(numtheory); M:=100;
    E := proc(k) local n, t1; global M;
    t1 := 1-(2*k/bernoulli(k))*add(sigma[k-1](n)*q^n, n=1..M+1);
    series(t1, q, M+1); end;
    e2:=E(2); e4:=E(4); e6:=E(6);
    t1:=series((e2*e4-e6)/720,q,M+1);
    seriestolist(t1);
    # alternative program
    seq(add(sigma[4](d)*phi(n/d), d in divisors(n)), n = 1..100); # Peter Bala, Jan 20 2024
  • Mathematica
    Table[If[n==0, 0, n * DivisorSigma[3, n]], {n, 0, 40}] (* Indranil Ghosh, Mar 11 2017 *)
    terms = 41; Ei[n_] = 1-(2n/BernoulliB[n]) Sum[k^(n-1) x^k/(1-x^k), {k, terms}]; CoefficientList[(Ei[2] Ei[4] - Ei[6])/720 + O[x]^terms, x] (* Jean-François Alcover, Mar 01 2018 *)
  • PARI
    for(n=0, 40, print1(if(n==0, 0, n * sigma(n, 3)), ", ")) \\ Indranil Ghosh, Mar 11 2017
    

Formula

a(n) = A145094(n)/240 for n > 0. - Seiichi Manyama, Feb 04 2017
G.f.: phi_{4, 1}(x) where phi_{r, s}(x) = Sum_{n, m>0} m^r * n^s * x^{m*n}. - Seiichi Manyama, Feb 04 2017
a(n) = n*A001158(n) for n > 0. - Seiichi Manyama, Feb 18 2017
G.f.: x*f'(x), where f(x) = Sum_{k>=1} k^3*x^k/(1 - x^k). - Ilya Gutkovskiy, Aug 31 2017
Sum_{k=1..n} a(k) ~ Pi^4 * n^5 / 450. - Vaclav Kotesovec, May 09 2022
From Amiram Eldar, Oct 30 2023: (Start)
Multiplicative with a(p^e) = p^e * (p^(3*e+3)-1)/(p^3-1).
Dirichlet g.f.: zeta(s-1)*zeta(s-4). (End)
a(n) = Sum_{k = 1..n} sigma_4( gcd(k, n) ) = Sum_{d divides n} sigma_4(d) * phi(n/d). - Peter Bala, Jan 19 2024
a(n) = Sum_{1 <= i, j, k, l <= n} sigma_1( gcd(i, j, k, l, n) ) = Sum_{d divides n} sigma_1(d) * J_4(n/d), where the Jordan totient function J_4(n) = A059377(n). - Peter Bala, Jan 22 2024

A282019 Coefficients in q-expansion of E_2*E_4, where E_2 and E_4 are the Eisenstein series shown in A006352 and A004009, respectively.

Original entry on oeis.org

1, 216, -3672, -62496, -322488, -1121904, -2969568, -6737472, -13678200, -24978312, -43826832, -70620768, -112325472, -166558896, -248342976, -346320576, -491604984, -655461072, -897864696, -1154109600, -1532856528, -1921344768, -2488726944, -3042415296, -3876616800, -4639932504
Offset: 0

Views

Author

N. J. A. Sloane, Feb 05 2017

Keywords

Crossrefs

Programs

  • Maple
    with(numtheory); M:=100;
    E := proc(k) local n, t1; global M;
    t1 := 1-(2*k/bernoulli(k))*add(sigma[k-1](n)*q^n, n=1..M+1);
    series(t1, q, M+1); end;
    e2:=E(2); e4:=E(4); e6:=E(6);
    series(e2*e4,q,M+1);
    seriestolist(%);
  • Mathematica
    terms = 26;
    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}];
    E2[x]*E4[x] + O[x]^terms // CoefficientList[#, x]& (* Jean-François Alcover, Feb 23 2018 *)

A282060 Coefficients in q-expansion of E_4*(E_2*E_4 - E_6)/720, where E_2, E_4, E_6 are the Eisenstein series shown in A006352, A004009, A013973, respectively.

Original entry on oeis.org

0, 1, 258, 6564, 66052, 390630, 1693512, 5764808, 16909320, 43066413, 100782540, 214358892, 433565328, 815730734, 1487320464, 2564095320, 4328785936, 6975757458, 11111134554, 16983563060, 25801892760, 37840199712, 55304594136, 78310985304, 110992776480
Offset: 0

Views

Author

Seiichi Manyama, Feb 05 2017

Keywords

Comments

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

Examples

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

Crossrefs

Cf. A064987 (phi_{2, 1}), A281372 (phi_{4, 1}), A282050 (phi_{6, 1}), this sequence (phi_{8, 1}).
Cf. A006352 (E_2), A004009 (E_4), A013973 (E_6), A282101 (E_2*E_4^2), A013974 (E_4*E_6 = E_10).

Programs

  • Mathematica
    terms = 25;
    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}];
    E4[x]*(E2[x]*E4[x] - E6[x])/720 + O[x]^terms // CoefficientList[#, x]& (* Jean-François Alcover, Feb 26 2018 *)
    Table[n * DivisorSigma[7, n], {n, 0, 24}] (* Amiram Eldar, Sep 06 2023 *)
    nmax = 40; CoefficientList[Series[x*Sum[k^8*x^(k-1)/(1 - x^k)^2, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 01 2025 *)
  • PARI
    a(n) = if(n < 1, 0, n*sigma(n, 7)) \\ Andrew Howroyd, Jul 25 2018

Formula

G.f.: phi_{8, 1}(x) where phi_{r, s}(x) = Sum_{n, m>0} m^r * n^s * x^{m*n}.
a(n) = (A282101(n) - A013974(n))/720. - Seiichi Manyama, Feb 10 2017
If p is a prime, a(p) = p^8 + p = A196288(p). - Seiichi Manyama, Feb 10 2017
a(n) = n*A013955(n) for n > 0. - Seiichi Manyama, Feb 18 2017
Sum_{k=1..n} a(k) ~ zeta(8) * n^9 / 9. - Amiram Eldar, Sep 06 2023
From Amiram Eldar, Oct 30 2023: (Start)
Multiplicative with a(p^e) = p^e * (p^(7*e+7)-1)/(p^7-1).
Dirichlet g.f.: zeta(s-1)*zeta(s-8). (End)
G.f. Sum_{k>=1} k^8*x^(k-1)/(1 - x^k)^2. - Vaclav Kotesovec, Aug 01 2025

A281374 Coefficients in q-expansion of E_2^2, where E_2 is the Eisenstein series shown in A006352.

Original entry on oeis.org

1, -48, 432, 3264, 9456, 21600, 39744, 66432, 105840, 147984, 220320, 281664, 393792, 475104, 646272, 743040, 980592, 1091232, 1432944, 1536960, 1965600, 2118144, 2649024, 2761344, 3516480, 3557040, 4433184, 4594560, 5575296, 5603040, 6998400, 6864384, 8407152, 8494848, 10085472, 9918720, 12319152
Offset: 0

Views

Author

N. J. A. Sloane, Feb 05 2017

Keywords

Crossrefs

Programs

  • Maple
    with(numtheory); M:=100;
    E := proc(k) local n, t1; global M;
    t1 := 1-(2*k/bernoulli(k))*add(sigma[k-1](n)*q^n, n=1..M+1);
    series(t1, q, M+1); end;
    e2:=E(2); e4:=E(4); e6:=E(6);
    series(e2^2,q,M+1);
    seriestolist(%);
  • Mathematica
    terms = 37;
    E2[x_] = 1 - 24*Sum[k*x^k/(1 - x^k), {k, 1, terms}];
    E2[x]^2 + O[x]^terms // CoefficientList[#, x]& (* Jean-François Alcover, Feb 23 2018 *)
    (* or *)
    Join[{1}, Table[240*DivisorSigma[3, n] - 288*n*DivisorSigma[1, n], {n, 1, 50}]] (* Vaclav Kotesovec, Aug 02 2025 *)

Formula

For n>0, a(n) = 240*A001158(n) - 288*n*A000203(n). - Vaclav Kotesovec, Aug 02 2025

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

Original entry on oeis.org

0, 1, 66, 732, 4228, 15630, 48312, 117656, 270600, 533637, 1031580, 1771572, 3094896, 4826822, 7765296, 11441160, 17318416, 24137586, 35220042, 47045900, 66083640, 86124192, 116923752, 148035912, 198079200, 244218775, 318570252, 389021400, 497449568, 594823350
Offset: 0

Views

Author

Seiichi Manyama, Feb 05 2017

Keywords

Comments

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

Examples

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

Crossrefs

Cf. A064987 (phi_{2, 1}), A281372 (phi_{4, 1}), this sequence (phi_{6, 1}), A282060 (phi_{8, 1}).
Cf. A006352 (E_2), A004009 (E_4), A013973 (E_6), A145095 (-q*E'_6), A008410 (E_4^2 = E_8), A282096 (E_2*E_6).

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}];
    (E4[x]^2 - E2[x]*E6[x])/1008 + O[x]^terms // CoefficientList[#, x]& (* Jean-François Alcover, Feb 26 2018 *)
    Table[n * DivisorSigma[5, n], {n, 0, 30}] (* Amiram Eldar, Sep 06 2023 *)
    nmax = 40; CoefficientList[Series[x*Sum[k^6*x^(k-1)/(1 - x^k)^2, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 01 2025 *)
  • PARI
    a(n) = if(n < 1, 0, n * sigma(n, 5)); \\ Andrew Howroyd, Jul 23 2018

Formula

a(n) = A145095(n)/504 for n > 0.
G.f.: phi_{6, 1}(x) where phi_{r, s}(x) = Sum_{n, m>0} m^r * n^s * x^{m*n}.
a(n) = (A008410(n) - A282096(n))/1008. - Seiichi Manyama, Feb 10 2017
If p is a prime, a(p) = p^6 + p = A131472(p). - Seiichi Manyama, Feb 10 2017
a(n) = n*A001160(n) for n > 0. - Seiichi Manyama, Feb 18 2017
Sum_{k=1..n} a(k) ~ zeta(6) * n^7 / 7. - Amiram Eldar, Sep 06 2023
From Amiram Eldar, Oct 30 2023: (Start)
Multiplicative with a(p^e) = p^e * (p^(5*e+5)-1)/(p^5-1).
Dirichlet g.f.: zeta(s-1)*zeta(s-6). (End)
G.f. Sum_{k>=1} k^6*x^(k-1)/(1 - x^k)^2. - Vaclav Kotesovec, Aug 01 2025

A282102 Coefficients in q-expansion of E_2*E_4*E_6, where E_2, E_4, E_6 are the Eisenstein series shown in A006352, A004009, A013973, respectively.

Original entry on oeis.org

1, -288, -129168, -1927296, 65152656, 1535768640, 15223408704, 98001292032, 474055120080, 1870878793824, 6312358836000, 18835985199744, 50831420617152, 126257508465984, 292348744636032, 637474437331200, 1319883180896592, 2610964045674432, 4963491913583664
Offset: 0

Views

Author

Seiichi Manyama, Feb 06 2017

Keywords

Comments

The series expansion of the 12th root of the generating function gives A341801. - Peter Bala, Feb 23 2021

Crossrefs

Cf. A006352 (E_2), A004009 (E_4), A013973 (E_6), A013974 (E_10).
Cf. A281374 (E_2^2), A282019 (E_2*E_4), A282096 (E_2*E_6), A282101 (E_2*E_8), this sequence (E_2*E_10), A341801.

Programs

  • Mathematica
    terms = 19;
    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]*E4[x]*E6[x] + O[x]^terms // CoefficientList[#, x]& (* Jean-François Alcover, Feb 23 2018 *)

A288968 Exponents a(1), a(2), ... such that E_2, 1 - 24*q - 72*q^2 - ... (A006352) is equal to (1-q)^a(1) (1-q^2)^a(2) (1-q^3)^a(3) ... .

Original entry on oeis.org

24, 348, 6424, 129300, 2778648, 62114524, 1428337176, 33527349924, 799482197272, 19302454317660, 470740035601176, 11575875047000596, 286650683468840472, 7140515309818664028, 178783562850377621272, 4496350112540599930692
Offset: 1

Views

Author

Seiichi Manyama, Jun 20 2017

Keywords

Crossrefs

Cf. this sequence (k=2), A110163 (k=4), A288851 (k=6), A288471 (k=8), A289024 (k=10), A288990/A288989 (k=12), A289029 (k=14).
Cf. A006352 (E_2), A008683, A288877 (E_4/E_2), A289635.

Formula

a(n) = 2 + (1/(12*n)) * Sum_{d|n} A008683(n/d) * A288877(d).
a(n) = (1/n) * Sum_{d|n} A008683(n/d) * A289635(d).
a(n) ~ 1 / (n * r^(2*n)), where r = A057823. - Vaclav Kotesovec, Mar 08 2018

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.
		

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)

A282101 Coefficients in q-expansion of E_2*E_4^2, where E_2, E_4 are the Eisenstein series shown in A006352, A004009, respectively.

Original entry on oeis.org

1, 456, 50328, -470496, -21784008, -234371664, -1446514848, -6502690752, -23328111240, -71276388312, -191952331632, -468159788448, -1052750026272, -2212261706256, -4394299104576, -8303419066176, -15060718806024, -26284654025712, -44471780630856
Offset: 0

Views

Author

Seiichi Manyama, Feb 06 2017

Keywords

Crossrefs

Cf. A006352 (E_2), A004009 (E_4), A008410 (E_8).
Cf. A281374 (E_2^2), A282019 (E_2*E_4), A282096 (E_2*E_6), this sequence (E_2*E_8), A282102 (E_2*E_10).

Programs

  • Mathematica
    terms = 19;
    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}];
    E2[x]*E4[x]^2 + O[x]^terms // CoefficientList[#, x]& (* Jean-François Alcover, Feb 23 2018 *)

A282018 Coefficients in q-expansion of E_2^3, where E_2 is the Eisenstein series shown in A006352.

Original entry on oeis.org

1, -72, 1512, -3744, -95544, -473904, -1538784, -3947328, -8597880, -16987176, -30607632, -52030944, -83972448, -129500784, -194056128, -279446976, -397468152, -544155408, -743106744, -978896160, -1296984528, -1654458624, -2139055776, -2661349824, -3370243680, -4106376504, -5113466064
Offset: 0

Views

Author

N. J. A. Sloane, Feb 05 2017

Keywords

Crossrefs

Cf. A006352.

Programs

  • Maple
    with(numtheory); M:=100;
    E := proc(k) local n, t1; global M;
    t1 := 1-(2*k/bernoulli(k))*add(sigma[k-1](n)*q^n, n=1..M+1);
    series(t1, q, M+1); end;
    e2:=E(2); e4:=E(4); e6:=E(6);
    series(e2^3,q,M+1);
    seriestolist(%);
  • Mathematica
    terms = 27;
    E2[x_] = 1 - 24*Sum[k*x^k/(1 - x^k), {k, 1, terms}];
    E2[x]^3 + O[x]^terms // CoefficientList[#, x]& (* Jean-François Alcover, Feb 23 2018 *)
Showing 1-10 of 115 results. Next