A247551 -id:A247551 - cp's OEIS frontend

A005651 Sum of multinomial coefficients (n_1+n_2+...)!/(n_1!n_2!...) where (n_1, n_2, ...) runs over all integer partitions of n.

1, 1, 3, 10, 47, 246, 1602, 11481, 95503, 871030, 8879558, 98329551, 1191578522, 15543026747, 218668538441, 3285749117475, 52700813279423, 896697825211142, 16160442591627990, 307183340680888755, 6147451460222703502, 129125045333789172825, 2841626597871149750951

Offset: 0

Simon Plouffe

This is the total number of hierarchies of n labeled elements arranged on 1 to n levels. A distribution of elements onto levels is "hierarchical" if a level l+1 contains <= elements than level l. Thus for n=4 the arrangement {1,2}:{3}{4} is not allowed. See also A140585. Examples: Let the colon ":" separate two consecutive levels l and l+1. Then n=2 --> 3: {1}{2}, {1}:{2}, {2}:{1}, n=3 --> 10: {1}{2}{3}, {1}{2}:{3}, {3}{1}:{2}, {2}{3}:{1}, {1}:{2}:{3}, {3}:{1}:{2}, {2}:{3}:{1}, {1}:{3}:{2}, {2}:{1}:{3}, {3}:{2}:{1}. - Thomas Wieder, May 17 2008
n identical objects are painted by dipping them into a long row of cans of paint of distinct colors. Begining with the first can and not skipping any cans k, 1<=k<=n, objects are dipped (painted) and not more objects are dipped into any subsequent can than were dipped into the previous can. The painted objects are then linearly ordered. - Geoffrey Critzer, Jun 08 2009
a(n) is the number of partitions of n where each part i is marked with a word of length i over an n-ary alphabet whose letters appear in alphabetical order and all n letters occur exactly once in the partition. a(3) = 10: 3abc, 2ab1c, 2ac1b, 2bc1a, 1a1b1c, 1a1c1b, 1b1a1c, 1b1c1a, 1c1a1b, 1c1b1a. - Alois P. Heinz, Aug 30 2015
Also the number of ordered set partitions of {1,...,n} with weakly decreasing block sizes. - Gus Wiseman, Sep 03 2018
The parity of a(n) is that of A000110(A000120(n)), so a(n) is even if and only if A000120(n) == 2 (mod 3). - Álvar Ibeas, Aug 11 2020

			For n=3, say the first three cans in the row contain red, white, and blue paint respectively. The objects can be painted r,r,r or r,r,w or r,w,b and then linearly ordered in 1 + 3 + 6 = 10 ways. - _Geoffrey Critzer_, Jun 08 2009
From _Gus Wiseman_, Sep 03 2018: (Start)
The a(3) = 10 ordered set partitions with weakly decreasing block sizes:
  {{1},{2},{3}}
  {{1},{3},{2}}
  {{2},{1},{3}}
  {{2},{3},{1}}
  {{3},{1},{2}}
  {{3},{2},{1}}
  {{2,3},{1}}
  {{1,2},{3}}
  {{1,3},{2}}
  {{1,2,3}}
(End)

Abramowitz and Stegun, Handbook, p. 831, column labeled "M_1".
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 126.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Alois P. Heinz, Table of n, a(n) for n = 0..450 (first 101 terms from T. D. Noe)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
M. E. Hoffman, Updown categories: Generating functions and universal covers, arXiv preprint arXiv:1207.1705 [math.CO], 2012.
A. Knopfmacher, A. M. Odlyzko, B. Pittel, L. B. Richmond, D. Stark, G. Szekeres, and N. C. Wormald, The Asymptotic Number of Set Partitions with Unequal Block Sizes, The Electronic Journal of Combinatorics, 6 (1999), R2.
S. Schreiber & N. J. A. Sloane, Correspondence, 1980.

Main diagonal of: A226873, A261719, A309973.
Row sums of: A226874, A262071, A327803.
Column k=1 of A309951.
Column k=0 of A327801.
Cf. A000041, A000110, A000258, A000670, A007837, A008277, A008480, A036038, A140585, A178682, A212855, A247551, A300335, A318762.

A005651b := proc(k) add( d/(d!)^(k/d),d=numtheory[divisors](k)) ; end proc:
A005651 := proc(n) option remember; local k ; if n <= 1 then 1; else (n-1)!*add(A005651b(k)*procname(n-k)/(n-k)!, k=1..n) ; end if; end proc:
seq(A005651(k), k=0..10) ; # R. J. Mathar, Jan 03 2011
# second Maple program:
b:= proc(n, i) option remember; `if`(n=0 or i=1, n!,
      b(n, i-1) +binomial(n, i)*b(n-i, min(n-i, i)))
    end:
a:= n-> b(n$2):
seq(a(n), n=0..25);  # Alois P. Heinz, Aug 29 2015, Dec 12 2016

Table[Total[n!/Map[Function[n, Apply[Times, n! ]], IntegerPartitions[n]]], {n, 0, 20}] (* Geoffrey Critzer, Jun 08 2009 *)
Table[Total[Apply[Multinomial, IntegerPartitions[n], {1}]], {n, 0, 20}] (* Jean-François Alcover and Olivier Gérard, Sep 11 2014 *)
b[n_, i_, t_] := b[n, i, t] = If[t==1, 1/n!, Sum[b[n-j, j, t-1]/j!, {j, i, n/t}]]; a[n_] := If[n==0, 1, n!*b[n, 0, n]]; Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Nov 20 2015, after Alois P. Heinz *)

a(m,n):=if n=m then 1 else sum(binomial(n,k)*a(k,n-k),k,m,(n/2))+1;
makelist(a(1,n),n,0,17); /* Vladimir Kruchinin, Sep 06 2014 */

a(n)=my(N=n!,s);forpart(x=n,s+=N/prod(i=1,#x,x[i]!));s \\ Charles R Greathouse IV, May 01 2015

{ my(n=25); Vec(serlaplace(prod(k=1, n, 1/(1-x^k/k!) + O(x*x^n)))) } \\ Andrew Howroyd, Dec 20 2017

E.g.f.: 1 / Product (1 - x^k/k!).
a(n) = Sum_{k=1..n} (n-1)!/(n-k)!*b(k)*a(n-k), where b(k) = Sum_{d divides k} d*d!^(-k/d). - Vladeta Jovovic, Oct 14 2002
a(n) ~ c * n!, where c = Product_{k>=2} 1/(1-1/k!) = A247551 = 2.52947747207915264... . - Vaclav Kotesovec, May 09 2014
a(n) = S(n,1), where S(n,m) = sum(k=m..n/2 , binomial(n,k)*S(n-k,k))+1, S(n,n)=1, S(n,m)=0 for nVladimir Kruchinin, Sep 06 2014
E.g.f.: exp(Sum_{k>=1} Sum_{j>=1} x^(j*k)/(k*(j!)^k)). - Ilya Gutkovskiy, Jun 18 2018

More terms from Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Mar 29 2003

A035341 Sum of ordered factorizations over all prime signatures with n factors.

Original entry on oeis.org

1, 1, 5, 25, 173, 1297, 12225, 124997, 1492765, 19452389, 284145077, 4500039733, 78159312233, 1460072616929, 29459406350773, 634783708448137, 14613962109584749, 356957383060502945, 9241222160142506097, 252390723655315856437, 7260629936987794508973

Offset: 0

Alford Arnold

Let f(n) = number of ordered factorizations of n (A074206(n)); a(n) = sum of f(k) over all terms k in A025487 that have n factors.
When the unordered spectrum A035310 is so ordered the sequences A000041 A000070 ...A035098 A000110 yield A000079 A001792 ... A005649 A000670 respectively.
Row sums of A095705. - David Wasserman, Feb 22 2008
From Ludovic Schwob, Sep 23 2023: (Start)
a(n) is the number of nonnegative integer matrices with sum of entries equal to n and no zero rows or columns, with weakly decreasing row sums. The a(3) = 25 matrices:
  [1 1 1] [1 2] [2 1] [3]
  .
  [1 1] [1 1] [1 1 0] [1 0 1] [0 1 1] [2] [0 2] [2 0]
  [1 0] [0 1] [0 0 1] [0 1 0] [1 0 0] [1] [1 0] [0 1]
  .
  [1] [1 0] [0 1] [1 0] [0 1] [1 0 0] [1 0 0] [0 1] [1 0]
  [1] [1 0] [0 1] [0 1] [1 0] [0 1 0] [0 0 1] [1 0] [0 1]
  [1] [0 1] [1 0] [1 0] [0 1] [0 0 1] [0 1 0] [1 0] [0 1]
  .
  [0 1 0] [0 1 0] [0 0 1] [0 0 1]
  [1 0 0] [0 0 1] [1 0 0] [0 1 0]
  [0 0 1] [1 0 0] [0 1 0] [1 0 0] (End)

			a(3) = 25 because there are 3 terms in A025487 with 3 factors, namely 8, 12, 30; and f(8)=4, f(12)=8, f(30)=13 and 4+8+13 = 25.

Alois P. Heinz, Table of n, a(n) for n = 0..250 (first 36 terms from David Wasserman)
Eric Weisstein's World of Mathematics, Perfect Partition

Cf. A005651, A025487, A035310, A255906, A365961.

Row sums of A261719.

b:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<1, 0,
      add(b(n-i*j, i-1, k)*binomial(i+k-1, k-1)^j, j=0..n/i)))
    end:
a:= n->add(add(b(n$2, k-i)*(-1)^i*binomial(k, i), i=0..k), k=0..n):
seq(a(n), n=0..25);  # Alois P. Heinz, Aug 29 2015

b[n_, i_, k_] := b[n, i, k] = If[n == 0, 1, If[i < 1, 0, Sum[b[n - i*j, i - 1, k]*If[j == 0, 1, Binomial[i + k - 1, k - 1]^j], {j, 0, n/i}]]];
a[n_] := Sum[Sum[b[n, n, k-i]*(-1)^i*Binomial[k, i], {i, 0, k}], {k, 0, n}];
Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Oct 26 2015, after Alois P. Heinz, updated Dec 15 2020 *)

R(n,k)=Vec(-1 + 1/prod(j=1, n, 1 - binomial(k+j-1,j)*x^j + O(x*x^n)))
seq(n) = {concat([1], sum(k=1, n, R(n, k)*sum(r=k, n, binomial(r, k)*(-1)^(r-k)) ))} \\ Andrew Howroyd, Sep 23 2023

a(n) ~ c * n! / log(2)^n, where c = 1/(2*log(2)) * Product_{k>=2} 1/(1-1/k!) = A247551 / (2*log(2)) = 1.8246323... . - Vaclav Kotesovec, Jan 21 2017

More terms from Erich Friedman.

More terms from David Wasserman, Feb 22 2008

A140585 Total number of all hierarchical orderings for all set partitions of n.

Original entry on oeis.org

1, 4, 20, 129, 1012, 9341, 99213, 1191392, 15958404, 235939211, 3817327362, 67103292438, 1273789853650, 25973844914959, 566329335460917, 13150556885604115, 324045146807055210, 8446201774570017379, 232198473069120178475, 6715304449424099384968

Offset: 1

Thomas Wieder, May 17 2008

nonn

			We are considering all set partitions of the n-set {1,2,3,...,n}.
For each such set partition we examine all possible hierarchical arrangements of the subsets. A hierarchy is a distribution of elements (sets in the present case) onto levels.
A distribution onto levels is "hierarchical" if a level L+1 contains at most as many elements as level L. Thus for n=4 the arrangement {1,2}:{3}{4} is not allowed.
Let the colon ":" separate two consecutive levels L and L+1.
n=2 --> 1+3=4
{1,2} {1}{2}
{1}:{2}
{2}:{1}
-----------------------
n=3 --> 1+9+10=20
1*1 3*3=9 1*10
{1,2,3} {1,2}{3} {1}{2}{3}
{1,3}{2}
{2,3}{1} {1}{2}:{3}
{3}{1}:{2}
{1,2}:{3} {2}{3}:{1}
{1,3}:{2}
{2,3}:{1} {1}:{2}:{3}
{3}:{1}:{2}
{3}:{1,2} {2}:{3}:{1}
{2}:{1,3} {1}:{3}:{2}
{1}:{2,3} {2}:{1}:{3}
{3}:{2}:{1}
-----------------------
n=4 --> 1+12+9+60+47=129
1*1 4*3=12 3*3=9 6*10=60 1*47
{1,2,3,4} {1,2,3}{4} {1,2}{3,4} {1,2}{3}{4} {1}{2}{3}{4}
{1,2,4}{3} {1,3}{2,4} {1,2}{3}:{4}
{1,3,4}{2} {1,4}{2,3} {1,2}{4}:{3} {1}{2}:{3}:{4}
{2,3,4}{1} {1}{2}:{3,4} {1}{3}:{2}:{4}
{1,2}:{3,4} {1,2}:{3}:{4} {1}{4}:{2}:{3}
{1,2,3}:{4} {1,3}:{2,4} {1,2}:{4}:{3} {1}{2}:{4}:{3}
{1,2,4}:{3} {1,4}:{2,3} {1}:{2}:{3,4} {1}{3}:{4}:{2}
{1,3,4}:{2} {3,4}:{1,2} {2}:{1}:{3,4} {1}{4}:{3}:{2}
{2,3,4}:{1} {2,4}:{1,3} {1}:{3,4}:{2}
{2,3}:{1,4} {2}:{3,4}:{1} {2}{3}:{1}:{4}
{4}:{1,2,3} {2}{4}:{1}:{3}
{3}:{1,2,4} likewise for: {2}{3}:{4}:{1}
{2}:{1,3,4} {3,4}{1}{2} {2}{4}:{3}:{1}
{1}:{2,3,4} {2,4}{1}{3}
{2,3}{1}{4} {3}{4}:{1}:{2}
{1,4}{2}{3} {3}{4}:{2}:{1}
{1,3}{2}{4}
{1}{2}:{3}{4}
{1}{3}:{2}{4}
{1}{4}:{2}{3}
{2}{3}:{1}{4}
{2}{4}:{1}{3}
{3}{4}:{1}{2}
{2}{3}{4}:{1}
{1}{3}{4}:{2}
{1}{2}{4}:{3}
{1}{2}{3}:{4}
{1}:{2}:{3}:{4}
+23 permutations

Vaclav Kotesovec, Table of n, a(n) for n = 1..420 (first 160 terms from Alois P. Heinz)
Thomas Wieder, The number of certain rankings and hierarchies formed from labeled or unlabeled elements and sets, Applied Mathematical Sciences, vol. 3, 2009, no. 55, 2707 - 2724. [_Thomas Wieder_, Nov 14 2009]

Cf. A005651, A075729, A034691, A083355, A109186, A000670.

A140585 := proc(n::integer) local k, Result; Result:=0: for k from 1 to n do Result:=Result+stirling2(n,k)*A005651(k); end do; lprint(Result); end proc;
E.g.f.: series(1/mul(1-(exp(x)-1)^k/k!,k=1..10),x=0,10). # Thomas Wieder, Sep 04 2008
# second Maple program:
with(numtheory): b:= proc(k) option remember; add(d/d!^(k/d), d=divisors(k)) end: c:= proc(n) option remember; `if`(n=0, 1, add((n-1)!/ (n-k)!* b(k)* c(n-k), k=1..n)) end: a:= n-> add(Stirling2(n, k) *c(k), k=1..n): seq(a(n), n=1..30); # Alois P. Heinz, Oct 10 2008

Table[n!*SeriesCoefficient[1/Product[(1-(E^x-1)^k/k!),{k,1,n}],{x,0,n}],{n,1,20}] (* Vaclav Kotesovec, Sep 03 2014 *)

Stirling transform of A005651 = Sum of multinomial coefficients: a(n) = Sum_{i=1..n} S2(n,k) A005651(k).
E.g.f.: 1/Product_{k>=1} (1 - (exp(x)-1)^k/k!). - Thomas Wieder, Sep 04 2008
a(n) ~ c * n! / (log(2))^n, where c = 1/(2*log(2)) * Product_{k>=2} 1/(1-1/k!) = A247551 / (2*log(2)) = 1.82463230250447246267598544320244231645906135137... . - Vaclav Kotesovec, Sep 04 2014, updated Jan 21 2017

More terms from Alois P. Heinz, Oct 10 2008

A209668 a(n) = count of monomials, of degree k = n, in the complete homogeneous symmetric polynomials h(mu,k) summed over all partitions mu of n.

Original entry on oeis.org

1, 1, 7, 55, 631, 8001, 130453, 2323483, 48916087, 1129559068, 29442232007, 835245785452, 26113646252773, 880685234758941, 32191922753658129, 1259701078978200555, 52802268925363689079, 2352843030410455053891, 111343906794849929711260, 5567596199767400904172045

Offset: 0

Wouter Meeussen, Mar 11 2012

nonn

a(n) is the number of partitions of n where each part i is marked with a word of length i over an n-ary alphabet whose letters appear in alphabetical order. a(2) = 7: 2aa, 2ab, 2bb, 1a1a, 1a1b, 1b1a, 1b1b. - Alois P. Heinz, Aug 30 2015

Alois P. Heinz, Table of n, a(n) for n = 0..386
Wikipedia, Symmetric Polynomials

Cf. A209664-A209673.
Main diagonal of A209666 and A261718.
Cf. A261783.

b:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<1, 0,
      add(b(n-i*j, i-1, k)*binomial(i+k-1, k-1)^j, j=0..n/i)))
    end:
a:= n-> b(n$3):
seq(a(n), n=0..25);  # Alois P. Heinz, Aug 29 2015

h[n_, v_] := Tr@ Apply[Times, Table[Subscript[x, j], {j, v}]^# & /@ Compositions[n, v], {1}]; h[par_?PartitionQ, v_] := Times @@ (h[#, v] & /@ par); Tr /@ Table[(h[#, l] & /@ Partitions[l]) /. Subscript[x, _] -> 1, {l, 10}]
b[n_, i_, k_] := b[n, i, k] = If[n==0, 1, If[i<1, 0, Sum[b[n-i*j, i-1, k] * Binomial[i+k-1, k-1]^j, {j, 0, n/i}]]]; a[n_] := b[n, n, n]; Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Jan 15 2016, after Alois P. Heinz *)

a(n) ~ c * n^n, where c = A247551 = Product_{k>=2} 1/(1-1/k!) = 2.529477472... . - Vaclav Kotesovec, Nov 15 2016

a(n) = [x^n] Product_{k>=1} 1 / (1 - binomial(k+n-1,n-1)*x^k). - Ilya Gutkovskiy, May 09 2021

a(0)=1 prepended and a(11)-a(19) added by Alois P. Heinz, Aug 29 2015

A072605 Number of necklaces with n beads over an n-ary alphabet {a1,a2,...,an} such that #(w,a1) >= #(w,a2) >= ... >= #(w,ak) >= 0, where #(w,x) counts the letters x in word w.

Original entry on oeis.org

1, 1, 2, 4, 13, 50, 270, 1641, 11945, 96784, 887982, 8939051, 99298354, 1195617443, 15619182139, 219049941201, 3293800835940, 52746930894774, 897802366250126, 16167544246362567, 307372573011579188, 6148811682561390279, 129164845357784003661

Offset: 0

Wouter Meeussen, Aug 06 2002

Alois P. Heinz, Table of n, a(n) for n = 0..200
Frank Ruskey, Necklaces, Lyndon words, De Bruijn sequences, etc.
F. Ruskey, Necklaces, Lyndon words, De Bruijn sequences, etc. [Cached copy, with permission, pdf format only]
Eric Weisstein's world of Mathematics, Necklaces

Cf. A005651, A052307, A247551, A261531, A261599, A261600.

neck[li:{__Integer}] := Module[{n, d}, n=Plus@@li; d=n-First[li]; Fold[ #1+(EulerPhi[ #2]*(n/#2)!)/Times@@((li/#2)!)&, 0, Divisors[GCD@@li]]/n]; Table[ Plus@@(neck /@ IntegerPartitions[n]), {n, 24}]

a(n)={if(n==0, 1, my(p=prod(k=1, n, 1/(1-x^k/k!) + O(x*x^n))); sumdiv(n, d, eulerphi(n/d)*d!*polcoeff(p,d))/n)} \\ Andrew Howroyd, Dec 20 2017

a(n) = (1/n) * Sum_{d|n} phi(n/d) * A005651(d) for n > 0. - Andrew Howroyd, Sep 25 2017
See Mathematica line.
a(n) ~ c * (n-1)!, where c = Product_{k>=2} 1/(1-1/k!) = A247551 = 2.52947747207915264818011615... . - Vaclav Kotesovec, Aug 27 2015

a(0)=1 prepended by Alois P. Heinz, Aug 23 2015

Name changed by Andrew Howroyd, Sep 25 2017

A261600 Number of primitive (aperiodic, or Lyndon) necklaces with n beads such that beads of a largest subset have label 0, beads of a largest remaining subset have label 1, and so on.

Original entry on oeis.org

1, 1, 1, 3, 11, 49, 265, 1640, 11932, 96780, 887931, 8939050, 99298073, 1195617442, 15619180497, 219049941148, 3293800823995, 52746930894773, 897802366153076, 16167544246362566, 307372573010691195, 6148811682561388635, 129164845357775064609

Offset: 0

Alois P. Heinz, Aug 27 2015

nonn

			a(3) = 3: 001, 012, 021.
a(4) = 11: 0001, 0011, 0012, 0021, 0102, 0123, 0132, 0213, 0231, 0312, 0321.

Alois P. Heinz, Table of n, a(n) for n = 0..300
F. Ruskey, Necklaces, Lyndon words, De Bruijn sequences, etc.
F. Ruskey, Necklaces, Lyndon words, De Bruijn sequences, etc. [Cached copy, with permission, pdf format only]
Eric Weisstein's World of Mathematics, Necklace
Wikipedia, Lyndon word
Wikipedia, Necklace (combinatorics)
Index entries for sequences related to necklaces

Cf. A072605, A247551, A261531, A261599.

with(numtheory):
b:= proc(n, i, g, d, j) option remember; `if`(g>0 and gn, 0, binomial(n/j, i/j)*b(n-i, i, igcd(i, g), d, j)))))
    end:
a:= n-> `if`(n=0, 1, add(add((f-> `if`(f=0, 0, f*b(n$2, 0, d, j)))(
                     mobius(j)), j=divisors(d)), d=divisors(n))/n):
seq(a(n), n=0..25);

b[n_, i_, g_, d_, j_] := b[n, i, g, d, j] = If[g>0 && gn, 0, Binomial[n/j, i/j]*b[n-i, i, GCD[i, g], d, j]]]]]; a[n_] := If[n==0, 1, Sum[Sum[ Function[f, If[f==0, 0, f*b[n, n, 0, d, j]]][MoebiusMu[j]], {j, Divisors[ d]}], {d, Divisors[n]}]/n]; Table[a[n], {n, 0, 25}] (* Jean-François Alcover, Feb 22 2017, translated from Maple *)

a(n) ~ c * (n-1)!, where c = Product_{k>=2} 1/(1-1/k!) = A247551 = 2.52947747207915264818011615... . - Vaclav Kotesovec, Aug 27 2015

A282529 Decimal expansion of Product_{n>=2} (1-1/n!).

Original entry on oeis.org

3, 9, 5, 3, 3, 8, 5, 6, 7, 3, 6, 7, 4, 4, 5, 5, 6, 6, 0, 3, 2, 3, 5, 6, 2, 0, 0, 4, 3, 1, 1, 8, 0, 6, 1, 3, 0, 2, 6, 6, 5, 6, 4, 8, 6, 7, 5, 8, 7, 2, 0, 1, 3, 2, 5, 0, 1, 1, 9, 7, 3, 9, 0, 1, 5, 1, 7, 4, 1, 6, 2, 0, 3, 8, 4, 7, 7, 5, 6, 1, 2, 6, 7, 0, 9, 2, 0, 0, 7, 7, 9, 3, 8, 9, 1, 0, 2, 2, 7, 0, 1, 1, 5, 8, 6, 8, 5, 4, 9, 6, 4, 0, 8, 4, 6, 7

Offset: 0

Lucian Craciun, Feb 22 2017

Conjectured to be irrational, transcendental and normal, none have been shown.

			0.39533856736744556603235620043118061302665648675872...

Lucian Craciun, Table of n, a(n) for n = 0..9900
StackExchange, Closed form for Product_{n>=2} (1-1/n!), Nov 11 2016.
Wolfram Alpha, Product_{n>=2} (1-1/n!).

Cf. A238695, A247551.

evalf(product(1-1/n!, n = 2..infinity), 100); # Lucian Craciun, Feb 22 2017

RealDigits[Product[1-1/n!, {n, 2, 70}], 10, 100][[1]] (* Robert G. Wilson v, Mar 19 2014 *)

prodinf(n=3, 1-1/gamma(n)) \\ Charles R Greathouse IV, Nov 13 2014

Equals 1/A247551. - Amiram Eldar, Sep 15 2023

A320566 Expansion of e.g.f. exp(x) * Product_{k>=1} 1/(1 - x^k/k!).

Original entry on oeis.org

1, 2, 6, 23, 110, 617, 4035, 29927, 249926, 2316317, 23674841, 264329177, 3207278255, 42011308653, 591460307157, 8905905152798, 142897741683846, 2433947385964373, 43873382718719949, 834402502632550589, 16699964488044322205, 350869837371828862607, 7721899536993122262447

Offset: 0

Ilya Gutkovskiy, Oct 15 2018

nonn

Binomial transform of A005651.

Vaclav Kotesovec, Table of n, a(n) for n = 0..440
N. J. A. Sloane, Transforms

Cf. A005651, A320567, A327870.

Row sums of A327801.

seq(coeff(series(factorial(n)*exp(x)*mul((1-x^k/factorial(k))^(-1),k=1..n),x,n+1), x, n), n = 0 .. 22); # Muniru A Asiru, Oct 15 2018

nmax = 22; CoefficientList[Series[Exp[x] Product[1/(1 - x^k/k!), {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]!
nmax = 22; CoefficientList[Series[Exp[x + Sum[Sum[x^(j k)/(k (j!)^k), {j, 1, nmax}], {k, 1, nmax}]], {x, 0, nmax}], x] Range[0, nmax]!
Table[Sum[Binomial[n, k] Total[Apply[Multinomial, IntegerPartitions[k], {1}]], {k, 0, n}], {n, 0, 22}]

E.g.f.: exp(x + Sum_{k>=1} Sum_{j>=1} x^(j*k)/(k*(j!)^k)).
a(n) = Sum_{k=0..n} binomial(n,k)*A005651(k).
a(n) ~ exp(1) * A247551 * n!. - Vaclav Kotesovec, Jul 21 2019

A183239 G.f.: exp( Sum_{n>=1} A005651(n)*x^n/n ), where A005651 gives the sums of multinomial coefficients.

Original entry on oeis.org

1, 1, 2, 5, 17, 69, 352, 2077, 14505, 114354, 1023839, 10130051, 110878314, 1320375213, 17086334702, 237832320231, 3552995476517, 56590659564489, 958653346775294, 17192978984630744, 325681548343314833, 6494280460641306608

Offset: 0

Paul D. Hanna, Jan 03 2011

nonn

			G.f.: A(x) = 1 + x + 2*x^2 + 5*x^3 + 17*x^4 + 69*x^5 + 352*x^6 +...
log(A(x)) = x + 3*x^2/2 + 10*x^3/3 + 47*x^4/4 + 246*x^5/5 + 1602*x^6/6 + 11481*x^7/7 + 95503*x^8/8 +...+ A005651(n)*x^n/n +...

Vaclav Kotesovec, Table of n, a(n) for n = 0..420

Cf. A005651, A183238, A183241.

{a(n)=polcoeff(exp(intformal(1/x*(-1+serlaplace(1/prod(k=1,n+1,1-x^k/k!+O(x^(n+2))))))),n)}

a(n) ~ c * (n-1)!, where c = Product_{k>=2} 1/(1-1/k!) = A247551 = 2.52947747207915264... . - Vaclav Kotesovec, Feb 19 2015

A292713 Number of multisets of nonempty words with a total of n letters over n-ary alphabet such that within each word every letter of the alphabet is at least as frequent as the subsequent alphabet letter.

Original entry on oeis.org

1, 1, 4, 14, 67, 343, 2151, 14900, 119259, 1055520, 10465854, 113479756, 1350508150, 17373376892, 241576630993, 3596468789967, 57232276979726, 967517444008250, 17339617861447844, 328037083000497867, 6537494747743375847, 136820214583596515519

Offset: 0

Alois P. Heinz, Sep 21 2017

nonn

			a(0) = 1: {}.
a(1) = 1: {a}.
a(2) = 4: {aa}, {ab}, {ba}, {a,a}.
a(3) = 14: {aaa}, {aab}, {aba}, {baa}, {abc}, {acb}, {bac}, {bca}, {cab}, {cba}, {aa,a}, {ab,a}, {ba,a}, {a,a,a}.

Alois P. Heinz, Table of n, a(n) for n = 0..450

Main diagonal of A292712.
Row sums of A319495.
Cf. A226873, A292796.

b:= proc(n, i, t) option remember; `if`(t=1, 1/n!,
      add(b(n-j, j, t-1)/j!, j=i..n/t))
    end:
g:= (n, k)-> `if`(k=0, `if`(n=0, 1, 0), n!*b(n, 0, k)):
A:= proc(n, k) option remember; `if`(n=0, 1, add(add(d*
      g(d, k), d=numtheory[divisors](j))*A(n-j, k), j=1..n)/n)
    end:
a:= n-> A(n$2):
seq(a(n), n=0..25);

b[n_, i_, t_] := b[n, i, t] = If[t == 1, 1/n!, Sum[b[n - j, j, t - 1]/j!, {j, i, n/t}]];
g[n_, k_] := If[k == 0, If[n == 0, 1, 0], n!*b[n, 0, k]];
A[n_, k_] := A[n, k] = If[n == 0, 1, Sum[Sum[d*g[d, k], {d, Divisors[j]}]* A[n - j, k], {j, 1, n}]/n];
a[n_] := A[n, n];
a /@ Range[0, 25] (* Jean-François Alcover, Dec 19 2020, after Alois P. Heinz *)

a(n) = [x^n] Product_{j=1..n} 1/(1-x^j)^A226873(j,n).
a(n) = A292712(n,n).
a(n) ~ c * n!, where c = A247551 = 2.5294774720791526... - Vaclav Kotesovec, Oct 05 2017

cp's OEIS Frontend

A005651 Sum of multinomial coefficients (n_1+n_2+...)!/(n_1!*n_2!*...) where (n_1, n_2, ...) runs over all integer partitions of n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Maxima

PARI

PARI

Formula

Extensions

A035341 Sum of ordered factorizations over all prime signatures with n factors.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A140585 Total number of all hierarchical orderings for all set partitions of n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A209668 a(n) = count of monomials, of degree k = n, in the complete homogeneous symmetric polynomials h(mu,k) summed over all partitions mu of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A072605 Number of necklaces with n beads over an n-ary alphabet {a1,a2,...,an} such that #(w,a1) >= #(w,a2) >= ... >= #(w,ak) >= 0, where #(w,x) counts the letters x in word w.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A261600 Number of primitive (aperiodic, or Lyndon) necklaces with n beads such that beads of a largest subset have label 0, beads of a largest remaining subset have label 1, and so on.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A282529 Decimal expansion of Product_{n>=2} (1-1/n!).

Views

A005651 Sum of multinomial coefficients (n_1+n_2+...)!/(n_1!n_2!...) where (n_1, n_2, ...) runs over all integer partitions of n.