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

A060043 Triangle T(n,k), n >= 1, k >= 1, of generalized sum of divisors function, read by rows.

Original entry on oeis.org

1, 3, 1, 4, 3, 7, 9, 6, 1, 15, 12, 3, 30, 8, 9, 45, 15, 22, 67, 13, 1, 42, 99, 18, 3, 81, 135, 12, 9, 140, 175, 28, 22, 231, 231, 14, 51, 351, 306, 24, 1, 97, 551, 354, 24, 3, 188, 783, 465, 31, 9, 330, 1134, 540, 18, 22, 568, 1546, 681, 39, 51, 918, 2142, 765, 20
Offset: 1

Views

Author

N. J. A. Sloane, Mar 19 2001

Keywords

Comments

Lengths of rows are 1 1 2 2 2 3 3 3 3 ... (A003056).

Examples

			Triangle turned on its side begins:
1 3 4 7 6 12  8 15 13 18 ...
    1 3 9 15 30 45 67 99 ...
           1  3  9 22 42 ...
                       1 ...
For example, T(6,2) = 15.
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          b(n, i-1)+add(expand(b(n-i*j, i-1)*j*x), j=1..n/i)))
        end:
    T:= n-> (p-> seq(coeff(p, x, degree(p)-i), i=0..degree(p)-1))(b(n$2)):
    seq(T(n), n=1..20);  # Alois P. Heinz, Jul 21 2025
  • Mathematica
    Clear[diag, m]; nmax = 19; kmax = Floor[(Sqrt[8*nmax+1]-1)/2]; m[0] = 0; diag[k_] := diag[k] = Sum[q^(Sum[m[i], {i, 1, k}])/(Times @@ (1 - q^Array[m, k]))^2, Sequence @@ Table[{m[j], m[j-1]+1, nmax}, {j, 1, k}] // Evaluate] + O[q]^(nmax+1) // CoefficientList[#, q]&; Table[ Select[ Table[diag[k][[j+1]], {k, 1, kmax}], IntegerQ[#] && # > 0&] // Reverse, {j, 1, nmax}] // Flatten (* Jean-François Alcover, Jul 18 2017 *)

Formula

T(n, 1) = sum of divisors of n (A000203), T(n, k) = sum of s_1*s_2*...*s_k where s_1, s_2, ..., s_k are such that s_1*m_1 + s_2*m_2 + ... + s_k*m_k = n and the sum is over all such k-partitions of n.
G.f. for k-th diagonal (the k-th row of the sideways triangle shown in the example): Sum_{ m_1 < m_2 < ... < m_k} q^(m_1+m_2+...+m_k)/((1-q^m_1)*(1-q^m_2)*...*(1-q^m_k))^2 = Sum_n T(n, k)*q^n.
G.f. for k-th diagonal: (-1)^k * (1/(2*k+1)) * ( Sum_{j>=k} (-1)^j * (2*j+1) * binomial(j+k,2*k) * q^(j*(j+1)/2) ) / ( Sum_{j>=0} (-1)^j * (2*j+1) * q^(j*(j+1)/2) ). - Seiichi Manyama, Sep 15 2023

Extensions

More terms from Naohiro Nomoto, Jan 24 2002

A385001 Irregular triangle read by rows: T(n,k) is the number of partitions of n with k designated summands, n >= 0, 0 <= k <= A003056(n).

Original entry on oeis.org

1, 0, 1, 0, 3, 0, 4, 1, 0, 7, 3, 0, 6, 9, 0, 12, 15, 1, 0, 8, 30, 3, 0, 15, 45, 9, 0, 13, 67, 22, 0, 18, 99, 42, 1, 0, 12, 135, 81, 3, 0, 28, 175, 140, 9, 0, 14, 231, 231, 22, 0, 24, 306, 351, 51, 0, 24, 354, 551, 97, 1, 0, 31, 465, 783, 188, 3, 0, 18, 540, 1134, 330, 9
Offset: 0

Views

Author

Omar E. Pol, Jul 17 2025

Keywords

Comments

The divisor function sigma_1(n) = A000203(n) is also the number of partitions of n with only one designated summand, n >= 1.
When part i has multiplicity j > 0 exactly one part i is "designated".
The length of the row n is A002024(n+1) = 1 + A003056(n), hence the first positive element in column k is in the row A000217(k).
Alternating row sums give A329157.
Columns converge to A000716.
This triangle equals A060043 with reversed rows and an additional column 0.

Examples

			Triangle begins:
--------------------------------------------
   n\k:   0    1     2     3     4    5   6
--------------------------------------------
   0 |    1;
   1 |    0,   1;
   2 |    0,   3;
   3 |    0,   4,    1;
   4 |    0,   7,    3;
   5 |    0,   6,    9;
   6 |    0,  12,   15,    1;
   7 |    0,   8,   30,    3;
   8 |    0,  15,   45,    9;
   9 |    0,  13,   67,   22;
  10 |    0,  18,   99,   42,    1;
  11 |    0,  12,  135,   81,    3;
  12 |    0,  28,  175,  140,    9;
  13 |    0,  14,  231,  231,   22;
  14 |    0,  24,  306,  351,   51;
  15 |    0,  24,  354,  551,   97,   1;
  16 |    0,  31,  465,  783,  188,   3;
  17 |    0,  18,  540, 1134,  330,   9;
  18 |    0,  39,  681, 1546,  568,  22;
  19 |    0,  20,  765, 2142,  918,  51;
  20 |    0,  42,  945, 2835, 1452, 108;
  21 |    0,  32, 1040, 3758, 2233, 208,  1;
  ...
For n = 6 and k = 1 there are 12 partitions of 6 with only one designated summand as shown below:
   6'
   3'+ 3
   3 + 3'
   2'+ 2 + 2
   2 + 2'+ 2
   2 + 2 + 2'
   1'+ 1 + 1 + 1 + 1 + 1
   1 + 1'+ 1 + 1 + 1 + 1
   1 + 1 + 1'+ 1 + 1 + 1
   1 + 1 + 1 + 1'+ 1 + 1
   1 + 1 + 1 + 1 + 1'+ 1
   1 + 1 + 1 + 1 + 1 + 1'
So T(6,1) = 12, the same as A000203(6) = 12.
.
For n = 6 and k = 2 there are 15 partitions of 6 with two designated summands as shown below:
   5'+ 1'
   4'+ 2'
   4'+ 1'+ 1
   4'+ 1 + 1'
   3'+ 1'+ 1 + 1
   3'+ 1 + 1'+ 1
   3'+ 1 + 1 + 1'
   2'+ 2 + 1'+ 1
   2'+ 2 + 1 + 1'
   2 + 2'+ 1'+ 1
   2 + 2'+ 1 + 1'
   2'+ 1'+ 1 + 1 + 1
   2'+ 1 + 1'+ 1 + 1
   2'+ 1 + 1 + 1'+ 1
   2'+ 1 + 1 + 1 + 1'
So T(6,2) = 15, the same as A002127(6) = 15.
.
For n = 6 and k = 3 there is only one partition of 6 with three designated summands as shown below:
   3'+ 2'+ 1'
So T(6,3) = 1, the same as A002128(6) = 1.
There are 28 partitions of 6 with designated summands, so A077285(6) = 28.
.
		

Crossrefs

Columns: A000007 (k=0), A000203 (k=1), A002127 (k=2), A002128 (k=3), A365664 (k=4), A365665 (k=5), A384926 (k=6).
Row sums give A077285.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          b(n, i-1)+add(expand(b(n-i*j, i-1)*j*x), j=1..n/i)))
        end:
    T:= n-> (p-> seq(coeff(p,x,i), i=0..degree(p)))(b(n$2)):
    seq(T(n), n=0..20);  # Alois P. Heinz, Jul 18 2025

Formula

From Alois P. Heinz, Jul 18 2025: (Start)
Sum_{k>=1} k * T(n,k) = A293421(n).
T(A000096(n),n) = A000716(n). (End)
G.f.: Product_{i>0} 1 + (y*x^i)/(1 - x^i)^2. - John Tyler Rascoe, Jul 23 2025
Conjecture: For fixed k >= 1, Sum_{j=1..n} T(j,k) ~ Pi^(2*k) * n^(2*k) / ((2*k)! * (2*k+1)!). - Vaclav Kotesovec, Aug 01 2025

A365665 Expansion of Sum_{0

Original entry on oeis.org

1, 3, 9, 22, 51, 108, 208, 390, 693, 1193, 1977, 3195, 4995, 7722, 11583, 17164, 24882, 35685, 50205, 70083, 96300, 131101, 176358, 235377, 310651, 407352, 529074, 682750, 874038, 1112085, 1405521, 1766259, 2206413, 2741431, 3389052, 4168089, 5103450, 6218469
Offset: 15

Views

Author

Seiichi Manyama, Sep 15 2023

Keywords

Comments

Number of partitions of n with five designated summands (when part i has multiplicity j > 0 exactly one part i is "designated"). For example: a(16) = 3 because there are three partitions of 16 with five designated summands: [6'+ 4'+ 3'+ 2'+ 1'], [5'+ 4'+ 3'+ 2'+ 1'+ 1], [5'+ 4'+ 3'+ 2'+ 1 + 1']. - Omar E. Pol, Jul 29 2025

Crossrefs

A diagonal of A060043.
Column k=5 of A385001.
Cf. A384926.

Programs

  • Mathematica
    nmax = 60; Drop[CoefficientList[Series[-1/11 * Sum[(-1)^k*(2*k + 1)*Binomial[k + 5, 10]*x^(k*(k + 1)/2), {k, 5, nmax}]/Sum[(-1)^k*(2*k + 1)*x^(k*(k + 1)/2), {k, 0, nmax}], {x, 0, nmax}], x], 15] (* Vaclav Kotesovec, Jul 29 2025 *)
    (* or *)
    Table[(10679/17203200 - 1571*n/774144 + 133*n^2/92160 - n^3/3072 + n^4/46080) * DivisorSigma[1, n] + (1571/1548288 - 133*n/122880 + 3*n^2/10240 - n^3/46080) * DivisorSigma[3, n] + (133/1228800 - n/20480 + n^2/215040) * DivisorSigma[5, n] + (1/516096 - n/3096576) * DivisorSigma[7, n] + DivisorSigma[9, n]/154828800, {n, 15, 60}] (* Vaclav Kotesovec, Jul 29 2025 *)

Formula

G.f.: -(1/11) * ( Sum_{k>=5} (-1)^k * (2*k+1) * binomial(k+5,10) * q^(k*(k+1)/2) ) / ( Sum_{k>=0} (-1)^k * (2*k+1) * q^(k*(k+1)/2) ).
From Vaclav Kotesovec, Jul 29 2025: (Start)
a(n) = (10679/17203200 - 1571*n/774144 + 133*n^2/92160 - n^3/3072 + n^4/46080)*sigma(n) + (1571/1548288 - 133*n/122880 + 3*n^2/10240 - n^3/46080)*sigma_3(n) + (133/1228800 - n/20480 + n^2/215040)*sigma_5(n) + (1/516096 - n/3096576)*sigma_7(n) + sigma_9(n)/154828800.
Sum_{k=1..n} a(k) ~ Pi^10 * n^10 / 144850083840000.
(End)

A384926 Number of partitions of n with six designated summands.

Original entry on oeis.org

1, 3, 9, 22, 51, 108, 221, 414, 765, 1344, 2310, 3834, 6248, 9894, 15408, 23550, 35394, 52353, 76402, 109959, 156366, 219850, 305796, 421281, 574568, 777234, 1042083, 1387037, 1831362, 2402595, 3128995, 4051797, 5211639, 6668490, 8482089, 10737063, 13516615
Offset: 21

Views

Author

Omar E. Pol, Jul 23 2025

Keywords

Examples

			21 has only one partition with six designated summands: [6'+ 5'+ 4'+ 3'+ 2'+ 1'], so a(21) = 1.
22 has three partitions with six designated summands: [7'+ 5'+ 4'+ 3'+ 2'+ 1'], [6'+ 5'+ 4'+ 3'+ 2'+ 1'+ 1], [6'+ 5'+ 4'+ 3'+ 2'+ 1 + 1'], so a(22) = 3.
		

Crossrefs

Column k=6 of A385001.
Columns of A385001 converge to A000716.
Other columns k of A385001 are A000007 (k=0), A000203 (k=1), A002127 (k=2), A002128 (k=3), A365664 (k=4), A365665 (k=5).

Programs

  • Maple
    b:= proc(n, i) option remember; series(`if`(n=0, 1, `if`(i<1, 0,
          b(n, i-1)+add(b(n-i*j, i-1)*j*x, j=1..n/i))), x, 7)
        end:
    a:= n-> coeff(b(n$2), x, 6):
    seq(a(n), n=21..57);  # Alois P. Heinz, Jul 23 2025
  • Mathematica
    nmax=60; Drop[CoefficientList[Series[1/13 * Sum[(-1)^k*(2*k + 1)*Binomial[k + 6, 12]*x^(k*(k + 1)/2), {k, 6, nmax}]/Sum[(-1)^k*(2*k + 1)*x^(k*(k + 1)/2), {k, 0, nmax}], {x, 0, nmax}], x] , 21] (* Vaclav Kotesovec, Jul 29 2025 *)

Formula

A000716(n) >= a(21+n) with equality only for n <= 6.
Sum_{k=1..n} a(k) ~ Pi^12 * n^12 / (12! * 13!). - Vaclav Kotesovec, Aug 01 2025

Extensions

More terms from Alois P. Heinz, Jul 23 2025

A374931 Expansion of Sum_{1<=i<=j<=k<=l} q^(i+j+k+l)/( (1-q^i)*(1-q^j)*(1-q^k)*(1-q^l) )^2.

Original entry on oeis.org

1, 9, 44, 156, 450, 1121, 2508, 5139, 9867, 17831, 30888, 51117, 82212, 127458, 193809, 285702, 415701, 588797, 826803, 1133352, 1545986, 2062203, 2744811, 3579155, 4665375, 5964543, 7638983, 9600392, 12104136, 14992496, 18642996, 22787829, 28006341, 33830697, 41125050, 49183310, 59202486
Offset: 4

Views

Author

Seiichi Manyama, Jul 24 2024

Keywords

Crossrefs

Showing 1-5 of 5 results.