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

A330973 Least positive integer with exactly n factorizations into factors > 1, and 0 if no such number exists.

Original entry on oeis.org

1, 4, 8, 12, 16, 0, 24, 0, 36, 0, 60, 48, 0, 0, 128, 72, 0, 0, 96, 0, 120, 256, 0, 0, 0, 180, 0, 0, 144, 192, 216, 0, 0, 0, 0, 420, 0, 240, 0, 0, 0, 1024, 0, 0, 384, 0, 288, 0, 0, 0, 0, 360, 0, 0, 0, 2048, 432, 0, 0, 0, 0, 0, 0, 480, 0, 900, 768, 0, 0, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Jan 06 2020

Keywords

Crossrefs

All nonzero terms belong to A025487.
Includes all highly factorable numbers A033833.
Factorizations are A001055, with image A045782.
The version without zeros is A045783.
The sorted version is A330972.
The strict version is A330974.
Positions of zeros are A330976.

Programs

  • Mathematica
    nn=10;
    fam[n_]:=fam[n]=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[fam[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    nds=Length/@Array[fam[#]&,2^nn];
    Table[If[#=={},0,#[[1,1]]]&[Position[nds,i]],{i,nn}]

Extensions

More terms from Jinyuan Wang, Jul 07 2021

A346520 Number A(n,k) of partitions of the (n+k)-multiset {0,...,0,1,2,...,k} with n 0's into distinct multisets; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 5, 5, 3, 2, 15, 15, 9, 5, 2, 52, 52, 31, 16, 7, 3, 203, 203, 120, 59, 25, 10, 4, 877, 877, 514, 244, 100, 38, 14, 5, 4140, 4140, 2407, 1112, 442, 161, 56, 19, 6, 21147, 21147, 12205, 5516, 2134, 750, 249, 80, 25, 8, 115975, 115975, 66491, 29505, 11147, 3799, 1213, 372, 111, 33, 10
Offset: 0

Views

Author

Alois P. Heinz, Jul 21 2021

Keywords

Comments

Also number A(n,k) of factorizations of 2^n * Product_{i=1..k} prime(i+1) into distinct factors; A(3,1) = 5: 2*3*4, 4*6, 3*8, 2*12, 24; A(1,2) = 5: 2*3*5, 5*6, 3*10, 2*15, 30.

Examples

			A(2,2) = 9: 00|1|2, 001|2, 1|002, 0|01|2, 0|1|02, 01|02, 00|12, 0|012, 0012.
Square array A(n,k) begins:
  1,  1,   2,   5,   15,    52,   203,    877,    4140, ...
  1,  2,   5,  15,   52,   203,   877,   4140,   21147, ...
  1,  3,   9,  31,  120,   514,  2407,  12205,   66491, ...
  2,  5,  16,  59,  244,  1112,  5516,  29505,  168938, ...
  2,  7,  25, 100,  442,  2134, 11147,  62505,  373832, ...
  3, 10,  38, 161,  750,  3799, 20739, 121141,  752681, ...
  4, 14,  56, 249, 1213,  6404, 36332, 220000, 1413937, ...
  5, 19,  80, 372, 1887, 10340, 60727, 379831, 2516880, ...
  6, 25, 111, 539, 2840, 16108, 97666, 629346, 4288933, ...
  ...
		

Crossrefs

Main diagonal gives A346519.
Antidiagonal sums give A346521.

Programs

  • Maple
    g:= proc(n) option remember; `if`(n=0, 1, add(g(n-j)*add(
         `if`(d::odd, d, 0), d=numtheory[divisors](j)), j=1..n)/n)
        end:
    s:= proc(n) option remember; expand(`if`(n=0, 1,
          x*add(s(n-j)*binomial(n-1, j-1), j=1..n)))
        end:
    S:= proc(n, k) option remember; coeff(s(n), x, k) end:
    b:= proc(n, i) option remember; `if`(n=0, 1,
         `if`(i=0, g(n), add(b(n-j, i-1), j=0..n)))
        end:
    A:= (n, k)-> add(S(k, j)*b(n, j), j=0..k):
    seq(seq(A(n, d-n), n=0..d), d=0..12);
  • Mathematica
    g[n_] := g[n] = If[n == 0, 1, Sum[g[n - j]*Sum[If[OddQ[d], d, 0], {d, Divisors[j]}], {j, 1, n}]/n];
    s[n_] := s[n] = Expand[If[n == 0, 1, x*Sum[s[n - j]*Binomial[n - 1, j - 1], {j, 1, n}]]];
    S[n_, k_] := S[n, k] = Coefficient[s[n], x, k];
    b[n_, i_] := b[n, i] = If[n == 0, 1, If[i == 0, g[n], Sum[b[n - j, i - 1], {j, 0, n}]]];
    A[n_, k_] := Sum[S[k, j]*b[n, j], {j, 0, k}];
    Table[Table[A[n, d - n], {n, 0, d}], {d, 0, 12}] // Flatten (* Jean-François Alcover, Jul 31 2021, after Alois P. Heinz *)

Formula

A(n,k) = A045778(A000079(n)*A070826(k+1)).
A(n,k) = Sum_{j=0..k} Stirling2(k,j)*Sum_{i=0..n} binomial(j+i-1,i)*A000009(n-i).

A035098 Near-Bell numbers: partitions of an n-multiset with multiplicities 1, 1, 1, ..., 1, 2.

Original entry on oeis.org

1, 2, 4, 11, 36, 135, 566, 2610, 13082, 70631, 407846, 2504071, 16268302, 111378678, 800751152, 6027000007, 47363985248, 387710909055, 3298841940510, 29119488623294, 266213358298590, 2516654856419723, 24566795704844210
Offset: 1

Views

Author

Keywords

Comments

A035098 and A000070 are near the two ends of a spectrum. Another way to look at A000070 is as the number of partitions of an n-multiset with multiplicities n-1, 1.
The very ends are the number of partitions and the Stirling numbers of the second kind, which count the n-multiset partitions with multiplicities n and 1,1,1,...,1, respectively.
Intermediate sequences are the number of ways of partitioning an n-multiset with multiplicities some partition of n.

Examples

			a(3)=4 because there are 4 ways to partition the multiset {1,2,2} (with multiplicities {1,2}): {{1,2,2}} {{1,2},{2}} {{1},{2,2}} {{1},{2},{2}}.
		

References

  • D. E. Knuth, The Art of Computer Programming, Vol. 4A, Table A-1, page 778. - N. J. A. Sloane, Dec 30 2018

Crossrefs

Row sums of A241500.
Column 1 of array in A322765.
Row n=2 of A346426.

Programs

  • Maple
    with(combinat): a:= n-> floor(1/2*(bell(n-2)+bell(n-1)+bell(n))): seq(a(n), n=1..25); # Zerinvary Lajos, Oct 07 2007
  • Mathematica
    f[n_] := Sum[ StirlingS2[n, k] ((k + 1) (k + 2)/2 + 1), {k, 0, n}]; Array[f, 22, 0]
    f[n_] := (BellB[n] + BellB[n + 1] + BellB[n + 2])/2; Array[f, 22, 0]
    Range[0, 21]! CoefficientList[ Series[ (1 + Exp@ x)^2/2 Exp[ Exp@ x - 1], {x, 0, 21}], x] (* 3 variants by Robert G. Wilson v, Jan 13 2011 *)
    Join[{1},Total[#]/2&/@Partition[BellB[Range[0,30]],3,1]] (* Harvey P. Dale, Jan 02 2019 *)

Formula

Sum_{k=0..n} Stirling2(n, k)*((k+1)*(k+2)/2+1). E.g.f.: 1/2*(1+exp(x))^2*exp(exp(x)-1). (1/2)*(Bell(n)+Bell(n+1)+Bell(n+2)). - Vladeta Jovovic, Sep 23 2003 [for offset -1]
a(n) ~ Bell(n)/2 * (1 + LambertW(n)/n). - Vaclav Kotesovec, Jul 28 2021

Extensions

More terms from Vladeta Jovovic, Sep 23 2003

A346500 Number A(n,k) of partitions of the (n+k)-multiset {1,2,...,n,1,2,...,k}; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 2, 2, 2, 5, 4, 4, 5, 15, 11, 9, 11, 15, 52, 36, 26, 26, 36, 52, 203, 135, 92, 66, 92, 135, 203, 877, 566, 371, 249, 249, 371, 566, 877, 4140, 2610, 1663, 1075, 712, 1075, 1663, 2610, 4140, 21147, 13082, 8155, 5133, 3274, 3274, 5133, 8155, 13082, 21147
Offset: 0

Views

Author

Alois P. Heinz, Jul 20 2021

Keywords

Comments

Also number A(n,k) of factorizations of Product_{i=1..n} prime(i) * Product_{i=1..k} prime(i); A(2,2) = 9: 2*2*3*3, 3*3*4, 6*6, 2*3*6, 4*9, 2*2*9, 3*12, 2*18, 36.

Examples

			A(2,2) = 9: 1122, 11|22, 12|12, 1|122, 112|2, 11|2|2, 1|1|22, 1|12|2, 1|1|2|2.
Square array A(n,k) begins:
    1,    1,    2,     5,    15,     52,     203,     877, ...
    1,    2,    4,    11,    36,    135,     566,    2610, ...
    2,    4,    9,    26,    92,    371,    1663,    8155, ...
    5,   11,   26,    66,   249,   1075,    5133,   26683, ...
   15,   36,   92,   249,   712,   3274,   16601,   91226, ...
   52,  135,  371,  1075,  3274,  10457,   56135,  325269, ...
  203,  566, 1663,  5133, 16601,  56135,  198091, 1207433, ...
  877, 2610, 8155, 26683, 91226, 325269, 1207433, 4659138, ...
  ...
		

Crossrefs

Columns (or rows) k=0-10 give: A000110, A035098, A322764, A322768, A346881, A346882, A346883, A346884, A346885, A346886, A346887.
Main diagonal gives A020555.
First upper (or lower) diagonal gives A322766.
Second upper (or lower) diagonal gives A322767.
Antidiagonal sums give A346490.
A(2n,n) gives A322769.

Programs

  • Maple
    g:= proc(n, k) option remember; uses numtheory; `if`(n>k, 0, 1)+
         `if`(isprime(n), 0, add(`if`(d>k or max(factorset(n/d))>d, 0,
            g(n/d, d)), d=divisors(n) minus {1, n}))
        end:
    p:= proc(n) option remember; `if`(n=0, 1, p(n-1)*ithprime(n)) end:
    A:= (n, k)-> g(p(n)*p(k)$2):
    seq(seq(A(n, d-n), n=0..d), d=0..10);
    # second Maple program:
    b:= proc(n) option remember; `if`(n=0, 1,
          add(b(n-j)*binomial(n-1, j-1), j=1..n))
        end:
    A:= proc(n, k) option remember; `if`(n
    				
  • Mathematica
    b[n_] := b[n] = If[n == 0, 1, Sum[b[n-j] Binomial[n-1, j-1], {j, 1, n}]];
    A[n_, k_] := A[n, k] = If[n < k, A[k, n],
         If[k == 0, b[n], (A[n + 1, k - 1] + Sum[A[n - k + j, j]*
         Binomial[k - 1, j], {j, 0, k - 1}] + A[n, k - 1])/2]];
    Table[Table[A[n, d - n], {n, 0, d}], {d, 0, 10}] // Flatten (* Jean-François Alcover, Aug 18 2021, after Alois P. Heinz's second program *)

Formula

A(n,k) = A001055(A002110(n)*A002110(k)).
A(n,k) = A(k,n).
A(n,k) = A322765(abs(n-k),min(n,k)).

A082775 Convolution of natural numbers >= 2 and the partition numbers (A000041).

Original entry on oeis.org

2, 5, 11, 21, 38, 64, 105, 165, 254, 381, 562, 813, 1162, 1636, 2279, 3139, 4285, 5794, 7776, 10353, 13694, 17992, 23502, 30520, 39433, 50687, 64855, 82607, 104785, 132375, 166608, 208921, 261090, 325196, 403779, 499818, 616928, 759335, 932135
Offset: 2

Views

Author

Alford Arnold, May 22 2003

Keywords

Comments

Contribution from George Beck, Jan 08 2011: (Start)
The number of multiset partitions of the n-multiset M={0,0,...,0,1,2} (with n-2 zeros) is sum_{k=0..(n-2)}( (n-k) * p(k) ) where p(k) is the number of partitions of k.
Proof:
For each k = 0, 1, ..., n-2, partition k zeros and add the remaining n-k-2 zeros to the block {1, 2}, to give p(k) partitions.
For each k, partition k zeros and add the remaining n-k-2 zeros to the two blocks {1} and {2} in all possible 1 + n-k-2 ways, which gives (1 + n-k-2) * p(k) partitions.
Together, the number of partitions of M is sum_{k=0..n-2}( (n-k) * p(k) ). (End)
A082775 is the special case of A126442 with n-k = 2.

Examples

			a(7) = 64 because (7,5,3,2,1,1) dot (2,3,4,5,6,7) = 14+15+12+10+6+7= 64.
		

Crossrefs

Column k=2 of A346426.

Programs

  • Mathematica
    f[n_] := Sum[(n - k) PartitionsP[k], {k, 0, n - 2}]; Array[f, 39, 2]

Formula

a(n) = a(n-1) + A000041(n) + A000070(n) for n>1. - Alford Arnold, Dec 10 2007
a(n) = n*A000070(n-2) - A182738(n-2) for n>2. - Vaclav Kotesovec, Jun 23 2015
a(n) ~ sqrt(3) * exp(Pi*sqrt(2*n/3)) / (2*Pi^2). - Vaclav Kotesovec, Jun 23 2015

Extensions

More terms from Ray Chandler, Oct 11 2003

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

A169587 The total number of ways of partitioning the multiset {1,1,1,2,3,...,n-2}.

Original entry on oeis.org

3, 7, 21, 74, 296, 1315, 6393, 33645, 190085, 1145246, 7318338, 49376293, 350384315, 2606467211, 20266981269, 164306340566, 1385709542808, 12133083103491, 110095025916745, 1033601910417425, 10024991744613469, 100316367530768074, 1034373400144455266
Offset: 3

Views

Author

Martin Griffiths, Dec 02 2009

Keywords

Examples

			The partitions of {1,1,1,2} are {{1},{1},{1},{2}}, {{1,1},{1},{2}}, {{1,2},{1},{1}}, {{1,1},{1,2}}, {{1,1,1},{2}}, {{1,1,2},{1}} and {{1,1,1,2}}, so a(4)=7.
		

Crossrefs

This is related to A000110, A035098 and A169588.
Row n=3 of A346426.
Cf. A346813.

Programs

  • Mathematica
    Table[(BellB[n] + 3 BellB[n - 1] + 5 BellB[n - 2] + 2 BellB[n - 3])/ 6, {n, 3, 23}]

Formula

For n>=3, a(n)=(Bell(n)+3Bell(n-1)+5Bell(n-2)+2Bell(n-3))/6, where Bell(n) is the n-th Bell number (the Bell numbers are given in A000110).
E.g.f.: (e^(3x)+6e^(2x)+9e^x+2)(e^(e^x-1))/6.

A169588 The total number of ways of partitioning the multiset {1,1,1,1,2,3,...,n-3}.

Original entry on oeis.org

5, 12, 38, 141, 592, 2752, 13960, 76464, 448603, 2801054, 18516832, 129034659, 944356507, 7235605732, 57879020756, 482189616711, 4174720731316, 37489711726834, 348592657600818, 3350919079643612, 33252861484374737, 340209759518479300, 3584240435109146792
Offset: 4

Views

Author

Martin Griffiths, Dec 02 2009

Keywords

Crossrefs

This is related to A000110, A035098 and A169587.
Row n=4 of A346426.
Cf. A346814.

Programs

  • Mathematica
    Table[(BellB[n] + 6 BellB[n - 1] + 17 BellB[n - 2] + 20 BellB[n - 3] + 21 BellB[n - 4])/24, {n, 4, 23}]

Formula

For n>=4, a(n)=(Bell(n)+6Bell(n-1)+17Bell(n-2)+20Bell(n-3)+21Bell(n-4))/24, where Bell(n) is the n-th Bell number (the Bell numbers are given in A000110). e.g.f. (e^(4x)+12e^(3x)+42e^(2x)+44e^x+21)(e^(e^x-1))/24.

A346428 Total number of partitions of all n-multisets {0,...,0,1,2,...,j} for 0 <= j <= n.

Original entry on oeis.org

1, 2, 6, 17, 53, 180, 683, 2866, 13219, 66307, 358532, 2074229, 12761831, 83086064, 570017222, 4106269668, 30965072776, 243778358992, 1998878586251, 17034471643814, 150591119435358, 1378657063570498, 13050460812585580, 127553991370245410, 1285578058726241427
Offset: 0

Views

Author

Alois P. Heinz, Jul 16 2021

Keywords

Comments

Also total number of factorizations of 2^(n-j) * Product_{i=1..j} prime(i+1) for 0 <= j <= n; a(2) = 6: 2*2, 4, 2*3, 6, 3*5, 15; a(3) = 17: 2*2*2, 2*4, 8, 2*2*3, 3*4, 2*6, 12, 2*3*5, 5*6, 3*10, 2*15, 30, 3*5*7, 7*15, 5*21, 3*35, 105.

Examples

			a(2) = 6: 00, 0|0, 01, 0|1, 12, 1|2.
a(3) = 17: 000, 0|00, 0|0|0, 001, 00|1, 0|01, 0|0|1, 012, 0|12, 02|1, 01|2, 0|1|2, 123, 1|23, 13|2, 12|3, 1|2|3.
		

Crossrefs

Antidiagonal sums of A346426.

Programs

  • Maple
    s:= proc(n) option remember; expand(`if`(n=0, 1,
          x*add(s(n-j)*binomial(n-1, j-1), j=1..n)))
        end:
    S:= proc(n, k) option remember; coeff(s(n), x, k) end:
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i=0,
          combinat[numbpart](n), add(b(n-j, i-1), j=0..n)))
        end:
    a:= n-> add(add(S(n-i, j)*b(i, j), j=0..n-i), i=0..n):
    seq(a(n), n=0..25);
  • Mathematica
    s[n_] := s[n] = Expand[If[n == 0, 1,
         x*Sum[s[n - j]*Binomial[n - 1, j - 1], {j, 1, n}]]];
    S[n_, k_] := S[n, k] = Coefficient[s[n], x, k];
    b[n_, i_] := b[n, i] = If[n == 0, 1, If[i == 0,
         PartitionsP[n], Sum[b[n - j, i - 1], {j, 0, n}]]];
    a[n_] := Sum[Sum[S[n - i, j]*b[i, j], {j, 0, n - i}], {i, 0, n}];
    Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Mar 12 2022, after Alois P. Heinz *)

Formula

a(n) = Sum_{j=0..n} A346426(n-j,j).

A346424 Number of partitions of the 2n-multiset {0,...,0,1,2,...,n}.

Original entry on oeis.org

1, 2, 11, 74, 592, 5317, 52902, 572402, 6670707, 83025806, 1096662664, 15292076689, 224145880470, 3440981816071, 55153081768896, 920494136057715, 15959177281931953, 286834809549486462, 5334308665713522860, 102476857445135062727, 2030589375575413246579
Offset: 0

Views

Author

Alois P. Heinz, Jul 16 2021

Keywords

Comments

Also number of factorizations of 2^n * Product_{i=1..n} prime(i+1); a(2) = 11: 2*2*3*5, 3*4*5, 2*5*6, 6*10, 2*3*10, 5*12, 4*15, 2*2*15, 3*20, 2*30, 60.

Examples

			a(0) = 1: {}.
a(1) = 2: 01, 0|1.
a(2) = 11: 00|1|2, 001|2, 1|002, 0|0|1|2, 0|01|2, 0|1|02, 01|02, 00|12, 0|0|12, 0|012, 0012.
		

Crossrefs

Programs

  • Maple
    s:= proc(n) option remember; expand(`if`(n=0, 1,
          x*add(s(n-j)*binomial(n-1, j-1), j=1..n)))
        end:
    S:= proc(n, k) option remember; coeff(s(n), x, k) end:
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i=0,
          combinat[numbpart](n), add(b(n-j, i-1), j=0..n)))
        end:
    a:= n-> add(S(n, j)*b(n, j), j=0..n):
    seq(a(n), n=0..21);
  • Mathematica
    s[n_] := s[n] = Expand[If[n == 0, 1,
         x*Sum[s[n - j]*Binomial[n - 1, j - 1], {j, 1, n}]]];
    S[n_, k_] := S[n, k] = Coefficient[s[n], x, k];
    b[n_, i_] := b[n, i] = If[n == 0, 1, If[i == 0,
         PartitionsP[n], Sum[b[n - j, i - 1], {j, 0, n}]]];
    a[n_] := Sum[S[n, j]*b[n, j], {j, 0, n}];
    Table[a[n], {n, 0, 21}] (* Jean-François Alcover, Apr 06 2022, after Alois P. Heinz *)

Formula

a(n) = A001055(A000079(n)*A070826(n+1)).
a(n) = Sum_{j=0..n} A048993(n,j)*A292508(n,j+1).
a(n) = A346426(n,n).
Showing 1-10 of 23 results. Next