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

A218276 Convolution of level 2 of the divisor function.

Original entry on oeis.org

0, 0, 1, 3, 7, 16, 22, 45, 49, 100, 95, 178, 161, 304, 250, 465, 372, 676, 525, 952, 720, 1280, 946, 1702, 1217, 2156, 1570, 2764, 1925, 3376, 2360, 4185, 2912, 4944, 3404, 6121, 4047, 6960, 4858, 8344, 5530, 9600, 6391, 11246, 7513, 12496, 8372, 14926, 9486
Offset: 1

Views

Author

Michel Marcus, Oct 25 2012

Keywords

Comments

Belongs to the family of convolution sums: Sum_{m < n*N} sigma(n)*sigma(n - N*m).
Named W2(n) by S. Alaca and K. S. Williams.
The convolution sum: Sum_{m < n} sigma(n)*sigma(n - m) = W1(n) is A000385(n+1).

Crossrefs

Programs

  • Maple
    with(numtheory): seq((1/48)*(22*sigma[3](n) - 2*sigma[3](2*n) + 5*sigma(n) - sigma(2*n) - 24*n*sigma(n) + 6*n*sigma(2*n)),n=1..60); # Ridouane Oudra, Feb 23 2021
  • Mathematica
    Table[Sum[DivisorSigma[1, k]*DivisorSigma[1, n - 2*k], {k, 1, Floor[(n - 1)/2]}], {n, 1, 50}] (* G. C. Greubel, Dec 24 2016 *)
  • PARI
    lista(nn) = {for (i=1, nn, s = sum(m=1, floor((i-1)/2), sigma(m)*sigma(i-2*m)); print1(s , ", "););}
    
  • PARI
    lista(nn) = {for (i=1, nn, v = sigma(i,3)/12 - i*sigma(i)/8 + sigma(i)/24;if (i%2 == 0, v += sigma(i/2,3)/3 - i*sigma(i/2)/4 + sigma(i/2)/24); print1(v , ", "););}

Formula

a(n) = Sum_{m < 2*n} sigma(n)*sigma(n - 2*m).
a(n) = sigma_3(n)/12 + sigma_3(n/2)/3 - n*sigma(n)/8 - n*sigma(n/2)/4 + sigma(n)/24 + sigma(n/2)/24.
a(n) = (1/48)*(22*sigma_3(n) - 2*sigma_3(2*n) + 5*sigma(n) - sigma(2*n) - 24*n*sigma(n) + 6*n*sigma(2*n)). - Ridouane Oudra, Feb 23 2021

A218277 Convolution of level 3 of the divisor function.

Original entry on oeis.org

0, 0, 0, 1, 3, 4, 10, 15, 24, 33, 45, 65, 77, 102, 143, 155, 180, 268, 255, 315, 434, 435, 462, 695, 593, 735, 960, 918, 945, 1437, 1160, 1395, 1825, 1692, 1668, 2549, 1995, 2385, 3073, 2775, 2730, 4190, 3157, 3747, 4739, 4290, 4140, 6355, 4686, 5523, 7044
Offset: 1

Views

Author

Michel Marcus, Oct 25 2012

Keywords

Comments

Named W3(n) by S. Alaca and K. S. Williams.

Crossrefs

Programs

  • Maple
    f:= n -> add(numtheory:-sigma(m)*numtheory:-sigma(n-3*m),m=1..floor((n-1)/3)):
    map(f, [$1..50]); # Robert Israel, Jun 28 2018
    with(numtheory): seq((1/72)*(31*sigma[3](n) - sigma[3](3*n) + 7*sigma(n) - sigma(3*n) - 30*n*sigma(n) + 6*n*sigma(3*n)), n=1..50); # Ridouane Oudra, Mar 21 2021
  • Mathematica
    a[n_] := Sum[DivisorSigma[1, m] DivisorSigma[1, n-3m], {m, 1, (n-1)/3}];
    Array[a, 50] (* Jean-François Alcover, Sep 19 2018 *)
  • PARI
    lista(n) = {for (i=1, n, s = sum(m=1, floor((i-1)/3), sigma(m)*sigma(i-3*m)); print1(s , ", "););}
    
  • PARI
    lista(n) = {for (i=1, n, v = sigma(i,3)/24 - i*sigma(i)/12 + sigma(i)/24;if (i%3 == 0, v += 3*sigma(i/3,3)/8 - i*sigma(i/3)/4 + sigma(i/3)/24); print1(v , ", "););}

Formula

a(n) = Sum_{m<3n} sigma(m)*sigma(n-3*m).
a(n) = sigma3(n)/24 - n*sigma(n)/12 + sigma(n)/24 + 3*sigma3(n/3)/8 - n*sigma(n/3)/4 + sigma(n/3)/24.
a(n) = (1/72)*(31*sigma_3(n) - sigma_3(3*n) + 7*sigma(n) - sigma(3*n) - 30*n*sigma(n) + 6*n*sigma(3*n)). - Ridouane Oudra, Mar 21 2021

A346193 Convolution of level 5 of the divisor function.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 3, 4, 7, 6, 15, 17, 27, 34, 36, 52, 64, 75, 91, 102, 122, 155, 169, 193, 228, 263, 276, 326, 349, 415, 430, 500, 520, 620, 681, 727, 741, 881, 880, 1090, 1020, 1192, 1178, 1375, 1513, 1590, 1557, 1809, 1756, 2274, 2024, 2323, 2245, 2626, 2865
Offset: 1

Views

Author

Amiram Eldar, Jul 09 2021

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_] := Sum[DivisorSigma[1, k] * DivisorSigma[1, n - 5*k], {k, 1, (n - 1)/5}]; Array[a, 100]
    (* or *)
    c[n_] := SeriesCoefficient[q * (QPochhammer[q] * QPochhammer[q^5])^4, {q, 0, n}]; a[n_] := 5 * DivisorSigma[3, n]/312 + If[Divisible[n, 5], 125 * DivisorSigma[3, n/5]/312, 0] - n * DivisorSigma[1, n]/20 - If[Divisible[n, 5], n * DivisorSigma[1, n/5]/4, 0] + DivisorSigma[1, n]/24 + If[Divisible[n, 5], DivisorSigma[1, n/5]/24, 0] - c[n]/130; Array[a, 100]

Formula

a(n) = Sum_{k < n/5} sigma(k) * sigma(n-5*k).
a(n) = 5*sigma_3(n)/312 + 125*sigma_3(n/5)/312 + (1/24 - n/20)*sigma(n) + (1/24 - n/4)*sigma(n/5) - c_5(n)/130, where sigma_3(n/5) = sigma(n/5) = 0 if n is not divisible by 5, and c_5(n) is the coefficient of q^n in the expansion of (eta(q) * eta(q^5))^4 (A030210).
Showing 1-3 of 3 results.