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.

Previous Showing 11-20 of 36 results. Next

A336419 Number of divisors d of the n-th superprimorial A006939(n) with distinct prime exponents such that the quotient A006939(n)/d also has distinct prime exponents.

Original entry on oeis.org

1, 2, 4, 10, 24, 64, 184, 536, 1608, 5104, 16448, 55136, 187136, 658624, 2339648, 8618208, 31884640, 121733120, 468209408, 1849540416, 7342849216
Offset: 0

Views

Author

Gus Wiseman, Jul 25 2020

Keywords

Comments

A number has distinct prime exponents iff its prime signature is strict.
The n-th superprimorial or Chernoff number is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).

Examples

			The a(0) = 1 through a(3) = 10 divisors:
  1  2  12  360
-----------------
  1  1   1    1
     2   3    5
         4    8
        12    9
             18
             20
             40
             45
             72
            360
		

Crossrefs

A000110 shifted once to the left dominates this sequence.
A006939 lists superprimorials or Chernoff numbers.
A022915 counts permutations of prime indices of superprimorials.
A130091 lists numbers with distinct prime exponents.
A181796 counts divisors with distinct prime exponents.
A181818 gives products of superprimorials.
A317829 counts factorizations of superprimorials.
A336417 counts perfect-power divisors of superprimorials.

Programs

  • Mathematica
    chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
    Table[Length[Select[Divisors[chern[n]],UnsameQ@@Last/@FactorInteger[#]&&UnsameQ@@Last/@FactorInteger[chern[n]/#]&]],{n,0,6}]
  • PARI
    recurse(n,k,b,d)={if(k>n, 1, sum(i=0, k, if((i==0||!bittest(b,i)) && (i==k||!bittest(d,k-i)), self()(n, k+1, bitor(b, 1<Andrew Howroyd, Aug 30 2020

Extensions

a(10)-a(20) from Andrew Howroyd, Aug 31 2020

A336426 Numbers that cannot be written as a product of superprimorials {2, 12, 360, 75600, ...}.

Original entry on oeis.org

3, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 25, 26, 27, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76
Offset: 1

Views

Author

Gus Wiseman, Jul 26 2020

Keywords

Comments

The n-th superprimorial is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).

Examples

			We have 288 = 2*12*12 so 288 is not in the sequence.
		

Crossrefs

A181818 is the complement.
A336497 is the version for superfactorials.
A001055 counts factorizations.
A006939 lists superprimorials or Chernoff numbers.
A022915 counts permutations of prime indices of superprimorials.
A317829 counts factorizations of superprimorials.
A336417 counts perfect-power divisors of superprimorials.

Programs

  • Mathematica
    chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
    facsusing[s_,n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[facsusing[Select[s,Divisible[n/d,#]&],n/d],Min@@#>=d&],{d,Select[s,Divisible[n,#]&]}]];
    Select[Range[100],facsusing[Array[chern,30],#]=={}&]

A048742 a(n) = n! - (n-th Bell number).

Original entry on oeis.org

0, 0, 0, 1, 9, 68, 517, 4163, 36180, 341733, 3512825, 39238230, 474788003, 6199376363, 86987391878, 1306291409455, 20912309745853, 355604563226196, 6401691628921841, 121639267666626943, 2432850284018404628, 51090467301893283249, 1123996221061869232677
Offset: 0

Views

Author

Keywords

Comments

Number of permutations of [n] which have at least one cycle that has at least one inversion when written with its smallest element in the first position. Example: a(4)=9 because we have (1)(243), (1432), (142)(3), (132)(4), (1342), (1423), (1243), (143)(2) and (1324). - Emeric Deutsch, Apr 29 2008
Number of permutations of [n] having consecutive runs of increasing elements with initial elements in increasing order. a(4) = 9: `124`3, `13`24, `134`2, `14`23, `14`3`2, `2`14`3, `24`3`1, `3`14`2, `4`13`2. - Alois P. Heinz, Apr 27 2016
From Gus Wiseman, Aug 11 2020: (Start)
Also the number of divisors of the superfactorial A006939(n - 1) without distinct prime multiplicities. For example, the a(4) = 9 divisors together with their prime signatures are the following. Note that A076954 can be used here instead of A006939.
6: (1,1)
10: (1,1)
15: (1,1)
30: (1,1,1)
36: (2,2)
60: (2,1,1)
90: (1,2,1)
120: (3,1,1)
180: (2,2,1)
(End)

Crossrefs

A000110 lists Bell numbers.
A000142 lists factorial numbers.
A006939 lists superprimorials or Chernoff numbers.
A181796 counts divisors with distinct prime multiplicities.
A336414 counts divisors of n! with distinct prime multiplicities.

Programs

Formula

a(n) = A000142(n) - A000110(n).
E.g.f.: 1/(1-x) - exp(exp(x)-1). - Alois P. Heinz, Apr 27 2016

A336107 Number of permutations of the prime indices of n with at least one non-singleton run, or non-separations.

Original entry on oeis.org

0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 2, 0, 0, 0, 1, 0, 2, 0, 2, 0, 0, 0, 4, 1, 0, 1, 2, 0, 0, 0, 1, 0, 0, 0, 4, 0, 0, 0, 4, 0, 0, 0, 2, 2, 0, 0, 5, 1, 2, 0, 2, 0, 4, 0, 4, 0, 0, 0, 6, 0, 0, 2, 1, 0, 0, 0, 2, 0, 0, 0, 9, 0, 0, 2, 2, 0, 0, 0, 5, 1, 0, 0, 6, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Sep 03 2020

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
A separation (or Carlitz composition) of a multiset is a permutation with no adjacent equal parts.

Examples

			The a(n) non-separations for n = 12, 36, 60, 72, 180, 420:
  (11)  (112)  (1122)  (1123)  (11122)  (11223)  (11234)
        (211)  (1221)  (1132)  (11212)  (11232)  (11243)
               (2112)  (2113)  (11221)  (11322)  (11324)
               (2211)  (2311)  (12112)  (12213)  (11342)
                       (3112)  (12211)  (12231)  (11423)
                       (3211)  (21112)  (13122)  (11432)
                               (21121)  (13221)  (21134)
                               (21211)  (21123)  (21143)
                               (22111)  (21132)  (23114)
                                        (22113)  (23411)
                                        (22131)  (24113)
                                        (22311)  (24311)
                                        (23112)  (31124)
                                        (23211)  (31142)
                                        (31122)  (32114)
                                        (31221)  (32411)
                                        (32112)  (34112)
                                        (32211)  (34211)
                                                 (41123)
                                                 (41132)
                                                 (42113)
                                                 (42311)
                                                 (43112)
                                                 (43211)
		

Crossrefs

A005117 lists positions of zeros, with complement A013929.
A008480 counts permutations of prime indices, ranked by A333221.
A003242 and A335452 count separations, ranked by A333489.
A325535 counts inseparable partitions, ranked by A335448.
A325534 counts separable partitions, ranked by A335433.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Select[Permutations[primeMS[n]],MatchQ[#,{_,x_,x_,_}]&]],{n,100}]

Formula

a(n) = A008480(n) - A335452(n).
a(A000961(n)) = 0 if n is in A027883, otherwise 1.
a(A005117(n)) = 0.
a(n!) = A335459(n).
a(A006939(n)) = A022915(n).

A337070 Number of strict chains of divisors starting with the superprimorial A006939(n).

Original entry on oeis.org

1, 2, 16, 1208, 1383936, 32718467072, 20166949856488576, 391322675415566237681536
Offset: 0

Views

Author

Gus Wiseman, Aug 15 2020

Keywords

Comments

The n-th superprimorial is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).

Examples

			The a(0) = 1 through a(2) = 16 chains:
  1  2    12
     2/1  12/1
          12/2
          12/3
          12/4
          12/6
          12/2/1
          12/3/1
          12/4/1
          12/4/2
          12/6/1
          12/6/2
          12/6/3
          12/4/2/1
          12/6/2/1
          12/6/3/1
		

Crossrefs

A022915 is the maximal case.
A076954 can be used instead of A006939 (cf. A307895, A325337).
A336571 is the case with distinct prime multiplicities.
A336941 is the case ending with 1.
A337071 is the version for factorials.
A000005 counts divisors.
A000142 counts divisors of superprimorials.
A006939 lists superprimorials or Chernoff numbers.
A067824 counts chains of divisors starting with n.
A074206 counts chains of divisors from n to 1.
A253249 counts chains of divisors.
A317829 counts factorizations of superprimorials.

Programs

  • Mathematica
    chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
    chnsc[n_]:=If[n==1,{{1}},Prepend[Join@@Table[Prepend[#,n]&/@chnsc[d],{d,Most[Divisors[n]]}],{n}]];
    Table[Length[chnsc[chern[n]]],{n,0,3}]

Formula

a(n) = 2*A336941(n) for n > 0.
a(n) = A067824(A006939(n)).

A095149 Triangle read by rows: Aitken's array (A011971) but with a leading diagonal before it given by the Bell numbers (A000110), 1, 1, 2, 5, 15, 52, ...

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 5, 2, 3, 5, 15, 5, 7, 10, 15, 52, 15, 20, 27, 37, 52, 203, 52, 67, 87, 114, 151, 203, 877, 203, 255, 322, 409, 523, 674, 877, 4140, 877, 1080, 1335, 1657, 2066, 2589, 3263, 4140, 21147, 4140, 5017, 6097, 7432, 9089, 11155, 13744, 17007, 21147
Offset: 0

Views

Author

Gary W. Adamson, May 30 2004

Keywords

Comments

Or, prefix Aitken's array (A011971) with a leading diagonal of 0's and take the differences of each row to get the new triangle.
With offset 1, triangle read by rows: T(n,k) is the number of partitions of the set {1,2,...,n} in which k is the largest entry in the block containing 1 (1 <= k <= n). - Emeric Deutsch, Oct 29 2006
Row term sums = the Bell numbers starting with A000110(1): 1, 2, 5, 15, ...
The k-th term in the n-th row is the number of permutations of length n starting with k and avoiding the dashed pattern 23-1. Equivalently, the number of permutations of length n ending with k and avoiding 1-32. - Andrew Baxter, Jun 13 2011
From Gus Wiseman, Aug 11 2020: (Start)
Conjecture: Also the number of divisors d with distinct prime multiplicities of the superprimorial A006939(n) that are of the form d = m * 2^k where m is odd. For example, row n = 4 counts the following divisors:
1 2 4 8 16
3 18 12 24 48
5 50 20 40 80
7 54 28 56 112
9 1350 108 72 144
25 540 200 400
27 756 360 432
45 504 720
63 600 1008
75 1400 1200
135 2160
175 2800
189 3024
675 10800
4725 75600
Equivalently, T(n,k) is the number of length-n vectors 0 <= v_i <= i whose nonzero values are distinct and such that v_n = k.
Crossrefs:
A008278 is the version counted by omega A001221.
A336420 is the version counted by Omega A001222.
A006939 lists superprimorials or Chernoff numbers.
A008302 counts divisors of superprimorials by Omega.
A022915 counts permutations of prime indices of superprimorials.
A098859 counts partitions with distinct multiplicities.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
(End)

Examples

			Triangle starts:
   1;
   1,  1;
   2,  1,  2;
   5,  2,  3,  5;
  15,  5,  7, 10, 15;
  52, 15, 20, 27, 37, 52;
From _Gus Wiseman_, Aug 11 2020: (Start)
Row n = 3 counts the following set partitions (described in Emeric Deutsch's comment above):
  {1}{234}      {12}{34}    {123}{4}    {1234}
  {1}{2}{34}    {12}{3}{4}  {13}{24}    {124}{3}
  {1}{23}{4}                {13}{2}{4}  {134}{2}
  {1}{24}{3}                            {14}{23}
  {1}{2}{3}{4}                          {14}{2}{3}
(End)
		

Crossrefs

Programs

  • Maple
    with(combinat): T:=proc(n,k) if k=1 then bell(n-1) elif k=2 and n>=2 then bell(n-2) elif k<=n then add(binomial(k-2,i)*bell(n-2-i),i=0..k-2) else 0 fi end: matrix(8,8,T): for n from 1 to 11 do seq(T(n,k),k=1..n) od; # yields sequence in triangular form
    Q[1]:=t*s: for n from 2 to 11 do Q[n]:=expand(t^n*subs(t=1,Q[n-1])+s*diff(Q[n-1],s)-Q[n-1]+s*Q[n-1]) od: for n from 1 to 11 do P[n]:=sort(subs(s=1,Q[n])) od: for n from 1 to 11 do seq(coeff(P[n],t,k),k=1..n) od; # yields sequence in triangular form - Emeric Deutsch, Oct 29 2006
    A011971 := proc(n,k) option remember ; if k = 0 then if n=0 then 1; else A011971(n-1,n-1) ; fi ; else A011971(n,k-1)+A011971(n-1,k-1) ; fi ; end: A000110 := proc(n) option remember; if n<=1 then 1 ; else add( binomial(n-1,i)*A000110(n-1-i),i=0..n-1) ; fi ; end: A095149 := proc(n,k) option remember ; if k = 0 then A000110(n) ; else A011971(n-1,k-1) ; fi ; end: for n from 0 to 11 do for k from 0 to n do printf("%d, ",A095149(n,k)) ; od ; od ; # R. J. Mathar, Feb 05 2007
    # alternative Maple program:
    b:= proc(n, m, k) option remember; `if`(n=0, 1, add(
          b(n-1, max(j, m), max(k-1, -1)), j=`if`(k=0, m+1, 1..m+1)))
        end:
    T:= (n, k)-> b(n, 0, n-k):
    seq(seq(T(n, k), k=0..n), n=0..10);  # Alois P. Heinz, Dec 20 2018
  • Mathematica
    nmax = 10; t[n_, 1] = t[n_, n_] = BellB[n-1]; t[n_, 2] = BellB[n-2]; t[n_, k_] /; n >= k >= 3 := t[n, k] = t[n, k-1] + t[n-1, k-1]; Flatten[ Table[ t[n, k], {n, 1, nmax}, {k, 1, n}]] (* Jean-François Alcover, Nov 15 2011, from formula with offset 1 *)
  • Python
    # requires Python 3.2 or higher.
    from itertools import accumulate
    A095149_list, blist = [1,1,1], [1]
    for _ in range(2*10**2):
        b = blist[-1]
        blist = list(accumulate([b]+blist))
        A095149_list += [blist[-1]]+ blist
    # Chai Wah Wu, Sep 02 2014, updated Chai Wah Wu, Sep 20 2014

Formula

With offset 1, T(n,1) = T(n,n) = T(n+1,2) = B(n-1) = A000110(n-1) (the Bell numbers). T(n,k) = T(n,k-1) + T(n-1,k-1) for n >= k >= 3. T(n,n-1) = B(n-1) - B(n-2) = A005493(n-3) for n >= 3 (B(q) are the Bell numbers A000110). T(n,k) = A011971(n-2,k-2) for n >= k >= 2. In other words, deleting the first row and first column we obtain Aitken's array A011971 (also called Bell or Pierce triangle; offset in A011971 is 0). - Emeric Deutsch, Oct 29 2006
T(n,1) = B(n-1); T(n,2) = B(n-2) for n >= 2; T(n,k) = Sum_{i=0..k-2} binomial(k-2,i)*B(n-2-i) for 3 <= k <= n, where B(q) are the Bell numbers (A000110). Generating polynomial of row n is P[n](t) = Q[n](t,1), where Q[n](t,s) = t^n*Q[n-1](1,s) + s*dQ[n-1](t,s)/ds + (s-1) Q[n-1](t,s); Q[1](t,s) = ts. - Emeric Deutsch, Oct 29 2006

Extensions

Corrected and extended by R. J. Mathar, Feb 05 2007
Entry revised by N. J. A. Sloane, Jun 01 2005, Jun 16 2007

A336421 Number of ways to choose a divisor of a divisor, both having distinct prime exponents, of the n-th superprimorial number A006939(n).

Original entry on oeis.org

1, 3, 13, 76, 571, 5309, 59341, 780149
Offset: 0

Views

Author

Gus Wiseman, Jul 25 2020

Keywords

Comments

A number has distinct prime exponents iff its prime signature is strict.
The n-th superprimorial or Chernoff number is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).

Examples

			The a(2) = 13 ways:
  12/1/1  12/2/1  12/3/1  12/4/1  12/12/1
          12/2/2  12/3/3  12/4/2  12/12/2
                          12/4/4  12/12/3
                                  12/12/4
                                  12/12/12
		

Crossrefs

A000258 shifted once to the left is dominated by this sequence.
A336422 is the generalization to non-superprimorials.
A000110 counts divisors of superprimorials with distinct prime exponents.
A006939 lists superprimorials or Chernoff numbers.
A008302 counts divisors of superprimorials by bigomega.
A022915 counts permutations of prime indices of superprimorials.
A076954 can be used instead of A006939.
A130091 lists numbers with distinct prime exponents.
A181796 counts divisors with distinct prime exponents.
A181818 gives products of superprimorials.
A317829 counts factorizations of superprimorials.

Programs

  • Mathematica
    chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
    strsig[n_]:=UnsameQ@@Last/@FactorInteger[n];
    Table[Total[Cases[Divisors[chern[n]],d_?strsig:>Count[Divisors[d],e_?strsig]]],{n,0,5}]

A336941 Number of strict chains of divisors starting with the superprimorial A006939(n) and ending with 1.

Original entry on oeis.org

1, 1, 8, 604, 691968, 16359233536, 10083474928244288, 195661337707783118840768, 139988400203593571474134024847360, 4231553868972506381329450624389969130848256, 6090860257621637852755610879241895108657182173073604608, 464479854191019594417264488167571483344961210693790188774166838214656
Offset: 0

Views

Author

Gus Wiseman, Aug 13 2020

Keywords

Examples

			The a(2) = 8 chains:
  12/1
  12/2/1
  12/3/1
  12/4/1
  12/6/1
  12/4/2/1
  12/6/2/1
  12/6/3/1
		

Crossrefs

A022915 is the maximal case.
A076954 can be used instead of A006939.
A336571 is the case with distinct prime multiplicities.
A336942 is the case using members of A130091.
A337070 is the version ending with any divisor of A006939(n).
A000005 counts divisors.
A074206 counts chains of divisors from n to 1.
A006939 lists superprimorials or Chernoff numbers.
A067824 counts divisor chains starting with n.
A181818 gives products of superprimorials, with complement A336426.
A253249 counts chains of divisors.
A317829 counts factorizations of superprimorials.
A336423 counts chains using A130091, with maximal case A336569.

Programs

  • Mathematica
    chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
    chns[n_]:=If[n==1,1,Sum[chns[d],{d,Most[Divisors[n]]}]];
    Table[chns[chern[n]],{n,0,3}]
  • PARI
    a(n)={my(sig=vector(n,i,i), m=vecsum(sig)); sum(k=0, m, prod(i=1, #sig, binomial(sig[i]+k-1, k-1))*sum(r=k, m, binomial(r,k)*(-1)^(r-k)))} \\ Andrew Howroyd, Aug 30 2020

Formula

a(n) = A337070(n)/2 for n > 0.
a(n) = A074206(A006939(n)).

Extensions

Terms a(8) and beyond from Andrew Howroyd, Aug 30 2020

A336496 Products of superfactorials (A000178).

Original entry on oeis.org

1, 2, 4, 8, 12, 16, 24, 32, 48, 64, 96, 128, 144, 192, 256, 288, 384, 512, 576, 768, 1024, 1152, 1536, 1728, 2048, 2304, 3072, 3456, 4096, 4608, 6144, 6912, 8192, 9216, 12288, 13824, 16384, 18432, 20736, 24576, 27648, 32768, 34560, 36864, 41472, 49152, 55296
Offset: 1

Views

Author

Gus Wiseman, Aug 03 2020

Keywords

Comments

First differs from A317804 in having 34560, which is the first term with more than two distinct prime factors.

Examples

			The sequence of terms together with their prime indices begins:
    1: {}
    2: {1}
    4: {1,1}
    8: {1,1,1}
   12: {1,1,2}
   16: {1,1,1,1}
   24: {1,1,1,2}
   32: {1,1,1,1,1}
   48: {1,1,1,1,2}
   64: {1,1,1,1,1,1}
   96: {1,1,1,1,1,2}
  128: {1,1,1,1,1,1,1}
  144: {1,1,1,1,2,2}
  192: {1,1,1,1,1,1,2}
  256: {1,1,1,1,1,1,1,1}
  288: {1,1,1,1,1,2,2}
  384: {1,1,1,1,1,1,1,2}
  512: {1,1,1,1,1,1,1,1,1}
		

Crossrefs

A001013 is the version for factorials, with complement A093373.
A181818 is the version for superprimorials, with complement A336426.
A336497 is the complement.
A000178 lists superfactorials.
A001055 counts factorizations.
A006939 lists superprimorials or Chernoff numbers.
A049711 is the minimum prime multiplicity in A000178.
A174605 is the maximum prime multiplicity in A000178.
A303279 counts prime factors of superfactorials.
A317829 counts factorizations of superprimorials.
A322583 counts factorizations into factorials.
A325509 counts factorizations of factorials into factorials.

Programs

  • Mathematica
    supfac[n_]:=Product[k!,{k,n}];
    facsusing[s_,n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facsusing[Select[s,Divisible[n/d,#]&],n/d],Min@@#>=d&]],{d,Select[s,Divisible[n,#]&]}]];
    Select[Range[1000],facsusing[Rest[Array[supfac,30]],#]!={}&]

A208437 Triangular array read by rows: T(n,k) is the number of set partitions of {1,2,...,n} that have exactly k distinct block sizes.

Original entry on oeis.org

1, 2, 2, 3, 5, 10, 2, 50, 27, 116, 60, 2, 560, 315, 142, 1730, 2268, 282, 6123, 14742, 1073, 30122, 72180, 12600, 2, 116908, 464640, 97020, 32034, 507277, 2676366, 997920, 2, 2492737, 16400098, 8751600, 136853, 15328119, 94209206, 81225144, 1527528, 56182092, 673282610, 614128515, 37837800
Offset: 1

Views

Author

Geoffrey Critzer, Feb 26 2012

Keywords

Comments

Column 1 = A038041.
Column 2 = A088142.
Column 3 = A133118.
Row sums = A000110 (Bell numbers).
Row n has floor([sqrt(1+8n)-1]/2) terms (number of terms increases by one at each triangular number). - Franklin T. Adams-Watters, Feb 26 2012

Examples

			:    1;
:    2;
:    2,      3;
:    5,     10;
:    2,     50;
:   27,    116,     60;
:    2,    560,    315;
:  142,   1730,   2268;
:  282,   6123,  14742;
: 1073,  30122,  72180,   12600;
		

Crossrefs

Programs

  • Maple
    with(combinat):
    b:= proc(n, i) option remember; expand(`if`(n=0, 1,
          `if`(i<1, 0, add(multinomial(n, n-i*j, i$j)/j!*
          b(n-i*j, i-1)*`if`(j=0, 1, x), j=0..n/i))))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=1..degree(p)))(b(n$2)):
    seq(T(n), n=1..16);  # Alois P. Heinz, Aug 21 2014
  • Mathematica
    nn = 15; p = Product[1 + y (Exp[x^i/i!] - 1), {i, 1, nn}];f[list_] := Select[list, # > 0 &];
    Map[f, Drop[ Range[0, nn]! CoefficientList[Series[p, {x, 0, nn}], {x, y}], 1]] // Flatten

Formula

E.g.f.: Product_{i>=1} 1 + y *(exp(x^i/i!)-1).
T(n*(n+1)/2,n) = A022915(n). - Alois P. Heinz, Apr 08 2016
Previous Showing 11-20 of 36 results. Next