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

A126442 Triangular array t read by rows: t(0,k) is p(k), the number of partitions of the k-multiset {0,0,...,0} with k zeros. For 0 <= n < k, t(n, k) is the number of partitions of the k-multiset {0, 0, ..., 0, 1, 2, 3, ..., k-n} with n zeros.

Original entry on oeis.org

1, 2, 2, 3, 4, 5, 5, 7, 11, 15, 7, 12, 21, 36, 52, 11, 19, 38, 74, 135, 203, 15, 30, 64, 141, 296, 566, 877, 22, 45, 105, 250, 592, 1315, 2610, 4140, 30, 67, 165, 426, 1098, 2752, 6393, 13082, 21147, 42, 97, 254, 696, 1940, 5317, 13960, 33645, 70631, 115975
Offset: 1

Views

Author

Alford Arnold, Jan 28 2007

Keywords

Comments

First in a series of triangular arrays which comprise subsequences of A096443(n).
The second array begins 9 16 26 29 52 92 47 98 198 371 and when the arrays are aligned as illustrated in triangle A126441 with p(n) values they sum to A035310 which counts unordered multisets.
Let t(n, k) be the number of ways to partition the k-multiset {0,0,...,0,1,2,3,4,...,k-n} with n zeros, 0 <= n < k. Then t(n, k) = sum_i = 0..k j = 0..n S(n, j) C(i, j) p(k - n - i), where S(n, j) are Stirling numbers of the second kind, C(i, j) are the number of compositions of i distinct objects into j parts, and p is the integer partition function.
To see this, partition [n] into j blocks; there are S(n, j) partitions. For such a partition x and for each i, there are C(i, j) ways to distribute i zeros into x, because the blocks of x are all distinct. There are p(k-n-i) ways to partition the remaining k-n-i zeros. Multiplying and summing gives the result. - George Beck, Jan 10 2011
Values are also part of A096443, A129306 and A249620. Columns are also columns of the last one of these irregular triangles. See "Partitions_of_multisets" link. - Tilman Piesk, Nov 09 2014

Examples

			This first array includes only the hook cases. A096443(9,14,16) correspond to partitions [2,2], [3,2] and [2,2,1] so these values do not appear in A126442.
The array begins:
1
2 2
3 4 5
5 7 11 15
7 12 21 36 52
		

Crossrefs

Programs

  • Mathematica
    (* The triangle is flattened to a sequence. *)
    t[n_, k_] := Sum[StirlingS2[n, j] * Binomial[-1 + i + j, i] * PartitionsP[k - n - i], {j, 0, n}, {i, 0, k - n}]; Table[ t[n, k], {k, 10}, {n, 0, k - 1}] // Flatten (* George Beck, Jan 10 2011 *)

Extensions

Definition clarified by George Beck, Jan 11 2011

A093802 Number of distinct factorizations of 105*2^n.

Original entry on oeis.org

5, 15, 36, 74, 141, 250, 426, 696, 1106, 1711, 2593, 3852, 5635, 8118, 11548, 16231, 22577, 31092, 42447, 57464, 77213, 103009, 136529, 179830, 235514, 306751, 397506, 512607, 658030, 841020, 1070490, 1357195, 1714274, 2157539, 2706174, 3383187, 4216358
Offset: 0

Views

Author

Alford Arnold, May 19 2004

Keywords

Examples

			105*A000079 is 105, 210, 420, 840, 1680, 3360, ... and there are 15 distinct factorizations of 210 so a(1) = 15.
a(0) = 5: 105*2^0 = 105 = 3*5*7 = 3*35 = 5*21 = 7*15. - _Alois P. Heinz_, May 26 2013
		

Crossrefs

Similar sequences: 45*A000079 => A002763, [1, 3, 9, 27, 81, 243...]*A000079 => A054225, 1*A002110 => A000110, 2*A002110 => A035098, A000142 => A076716.
Column k=3 of A346426.

Programs

  • Maple
    with(numtheory):
    b:= proc(n, k) option remember;
          `if`(n>k, 0, 1) +`if`(isprime(n), 0,
          add(`if`(d>k, 0, b(n/d, d)), d=divisors(n) minus {1, n}))
        end:
    a:= n-> b((105*2^n)$2):
    seq(a(n), n=0..50);  # Alois P. Heinz, May 26 2013
  • Mathematica
    b[n_, k_] := b[n, k] = If[n > k, 0, 1] + If[PrimeQ[n], 0,
         Sum[If[d > k, 0, b[n/d, d]], {d, Divisors[n][[2;;-2]]}]];
    a[n_] := b[105*2^n, 105*2^n];
    Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Jul 15 2021, after Alois P. Heinz *)

Extensions

2 more terms from Alford Arnold, Aug 29 2007
Corrected offset and extended beyond a(7) by Alois P. Heinz, May 26 2013
Showing 1-2 of 2 results.