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

A326237 Number of non-nesting digraphs with vertices {1..n}, where two edges (a,b), (c,d) are nesting if a < c and b > d or a > c and b < d.

Original entry on oeis.org

1, 2, 12, 104, 1008, 10272, 107712, 1150592
Offset: 0

Views

Author

Gus Wiseman, Jun 19 2019

Keywords

Comments

These are digraphs with the property that, if the edges are listed in lexicographic order, the sequence of targets is weakly increasing. For example, the digraph with lexicographically ordered edge set {(1,2),(2,1),(3,1),(3,2)} is nesting because the targets are (2,1,1,2), a sequence that is not weakly increasing.
Also the number of non-semicrossing digraphs with vertices {1..n}, where two edges (a,b), (c,d) are semicrossing if a < c and b < d or a > c and b > d. For example, the a(2) = 4 non-semicrossing digraph edge-sets are:
{}
{11}
{12}
{21}
{22}
{11,12}
{11,21}
{12,21}
{12,22}
{21,22}
{11,12,21}
{12,21,22}
Apparently a duplicate of A152254. - R. J. Mathar, Jul 12 2019

Examples

			The a(2) = 12 non-nesting digraph edge-sets:
  {}
  {11}
  {12}
  {21}
  {22}
  {11,12}
  {11,21}
  {11,22}
  {12,22}
  {21,22}
  {11,12,22}
  {11,21,22}
		

Crossrefs

Nesting digraphs are A326209.
Non-nesting set partitions are A000108.
Non-capturing set partitions are A054391.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Tuples[Range[n],2]],OrderedQ[Last/@#]&]],{n,4}]

Formula

A002416(n) = a(n) + A326209(n).

A033452 "STIRLING" transform of squares A000290.

Original entry on oeis.org

0, 1, 5, 22, 99, 471, 2386, 12867, 73681, 446620, 2856457, 19217243, 135610448, 1001159901, 7714225057, 61904585510, 516347066551, 4468588592739, 40058673825258, 371421499686007, 3556976106133821, 35138574378189700, 357654857584636597, 3746672593640388775
Offset: 0

Views

Author

Keywords

Comments

If an integer N is squarefree and has n+2 distinct prime factors then a(n) is the number of product signs needed to write the factorizations of N, so a(n)=A076277(N). - Floor van Lamoen, Oct 17 2002
Convolved with powers of 2 = A058681: (1, 7, 36, 171, 813, ...). Cf. triangle A180338. - Gary W. Adamson, Aug 28 2010

Examples

			G.f. = x + 5*x^2 + 22*x^3 + 99*x^4 + 471*x^5 + 2386*x^6 + 12867*x^7 + 73681*x^8 + ...
		

Crossrefs

Partial sums of A005494.
Cf. A180338.

Programs

  • Maple
    a := n -> add(Stirling2(n, j)*j^2, j=0..n): seq(a(n), n=0..20); # Zerinvary Lajos, Apr 18 2007
    # second Maple program:
    b:= proc(n, m) option remember;
         `if`(n=0, m^2, m*b(n-1, m)+b(n-1, m+1))
        end:
    a:= n-> b(n, 0):
    seq(a(n), n=0..23);  # Alois P. Heinz, Aug 04 2021
  • Mathematica
    max = 20; Clear[g]; g[max + 2] = 1; g[k_] := g[k] = 2 - 1/(1 - k*x)/(1 - x/(x - 1/g[k + 1])); gf = 1/x + 1/x^2 - g[0]/x^2; CoefficientList[ Series[gf, {x, 0, max}], x] (* Jean-François Alcover, Jan 24 2013, after Sergei N. Gladkovskii *)
  • PARI
    {a(n) = if( n<0, 0, n! * polcoeff( (exp(x + x * O(x^n)) - 1) * exp( exp(x + x * O(x^n)) - 1 + x), n))}; /* Michael Somos, Mar 28 2012 */

Formula

Representation as an infinite series: a(n) = (Sum_{k>=1} k^n*k*(k-2)/k!)/exp(1), n >= 1. This is a Dobinski-type summation formula. - Karol A. Penson, Mar 21 2002
a(n) = A005493(n) - A000110(n+1). - Floor van Lamoen and Christian Bower, Oct 16 2002. (n^2 has e.g.f.: e^x * (x^2+x), a(n) thus has e.g.f: e^(e^x-1) * ( (e^x-1)^2 + (e^x-1) ) which simplifies to e^(e^x-1) * (e^2x - e^x). A005493 has e.g.f.: e^(e^x+2x-1), A000110 has e.g.f.: e^(e^x-1), A000110(n+1) has as e.g.f.: derivative of A000110 which is e^(e^x+x-1).) [corrected by Georg Fischer, Jun 17 2020]
a(n) = Bell(n+2) - 2*Bell(n+1). - Vladeta Jovovic, Jul 28 2003
G.f.: sum{k>=0, k^2*x^k/prod[l=1..k, 1-lx]}. - Ralf Stephan, Apr 18 2004
E.g.f.: exp( exp(x) - 1 + x) * (exp(x) - 1). - Michael Somos, Mar 28 2012
a(n) = A123158(n,3). - Philippe Deléham, Oct 06 2006
G.f.: G(0)/x -1/x, where G(k) = 1 - x^2*(k+1)/( x^2*(k+1) - (2*x+x*k-1)*(3*x+x*k-1)/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Feb 25 2014

A326246 Number of crossing, capturing set partitions of {1..n}.

Original entry on oeis.org

0, 0, 0, 0, 0, 3, 37, 307, 2173, 14344, 92402, 596688
Offset: 0

Views

Author

Gus Wiseman, Jun 20 2019

Keywords

Comments

A set partition is crossing if it has two blocks of the form {...x...y...}, {...z...t...} where x < z < y < t or z < x < t < y, and capturing if it has two blocks of the form {...x...y...}, {...z...t...} where x < z < t < y or z < x < y < t. Capturing is a weaker condition than nesting, so for example {{1,3,5},{2,4}} is capturing but not nesting.

Examples

			The a(5) = 3 set partitions:
  {{1,3,4},{2,5}}
  {{1,3,5},{2,4}}
  {{1,4},{2,3,5}}
		

Crossrefs

MM-numbers of crossing, capturing multiset partitions are A326259.
Crossing set partitions are A016098.
Capturing set partitions are A326243.
Crossing, nesting set partitions are A326248.
Crossing, non-capturing set partitions are A326245.
Non-crossing, capturing set partitions are A122880 (conjecture).

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    croXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x
    				

A326249 Number of capturing set partitions of {1..n} that are not nesting.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 9, 55, 283, 1324, 5838, 24744
Offset: 0

Views

Author

Gus Wiseman, Jun 20 2019

Keywords

Comments

Capturing is a weaker condition than nesting. A set partition is capturing if it has two blocks of the form {...x...y...}, {...z...t...} where x < z < t < y or z < x < y < t, and nesting if it has two blocks of the form {...x,y...}, {...z,t...} where x < z < t < y or z < x < y < t. For example, {{1,3,5},{2,4}} is capturing but not nesting, so is counted under a(5).

Examples

			The a(6) = 9 set partitions:
  {{1},{2,4,6},{3,5}}
  {{1,3,5},{2,4},{6}}
  {{1,3,6},{2,4},{5}}
  {{1,3,6},{2,5},{4}}
  {{1,4,6},{2},{3,5}}
  {{1,4,6},{2,5},{3}}
  {{1,3,5},{2,4,6}}
  {{1,2,4,6},{3,5}}
  {{1,3,5,6},{2,4}}
		

Crossrefs

MM-numbers of capturing, non-nesting multiset partitions are A326260.
Nesting set partitions are A016098.
Capturing set partitions are A326243.
Non-crossing, nesting set partitions are A122880 (conjectured).

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    capXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x_,{_,x_,y_,_},_,{_,z_,t_,_},_}/;x
    				

A005465 Number of n-dimensional hypotheses allowing for conditional independence.

Original entry on oeis.org

0, 0, 1, 10, 70, 431, 2534, 14820, 88267, 542912, 3475978, 23253693, 162723444, 1190464900, 9092400633, 72370378750, 599168889634, 5150536258735, 45891028609826, 423144495659912, 4031842435506171, 39645279656283820, 401806832058661334, 4192631368792015237, 44992655908959220440
Offset: 0

Views

Author

Keywords

Comments

From Petros Hadjicostas, Oct 10 2019: (Start)
Mallows (1979) comments that I. J. Good did not consider all kinds of independence between n random variables in deriving his formula for the e.g.f. of a(n). Unfortunately, the paper by Good (1975), where the proof of the e.g.f. was published, is not easily accessible.
We give a sketch of a proof of the formula below (using only the ideas of independence that I. J. Good considered). Given r.v.'s X_1, X_2, ..., X_n, we choose s of them on which to condition (where s = 0, 1, 2, ..., n-2). By "condition", we mean that we condition on every possible s-tuple of values of those chosen variables. This can be done in C(n, s) ways.
Note that we can only condition on up to n-2 variables, because we need at least two variables to define any kind of independence: conditional (s >= 1) or unconditional (s = 0). Thus, 2 <= n-s <= n.
From the remaining n-s variables, we choose t of them (where 2 <= t <= n-s) on which we will define (or test) independence. [According to Mallows (1979), this is the only kind of independence Good (1975) considers.] There are C(n-s, t) ways to choose the t variables on which to define (or test) independence.
Now, there are Bell(t) - 1 = A058692(t) ways to partition the set of t chosen variables into one or more subsets, say {S_1, ..., S_r} (the order of the subsets is not important). Here |S_i| >= 1 and (S_1 union S_2 union ... union S_r) equals the t chosen variables. Thus, there are Bell(t) - 1 ways to factor the joint p.d.f. of the chosen t variables into the product of the joint marginal p.d.f.'s of the variables in S_i (i = 1, ..., r). [By "joint marginal p.d.f." for group S_i we mean the joint p.d.f of all the variables in group S_i.] Each such factorization defines a different kind of independence of the chosen variables (conditional on the original chosen s variables).
Thus, in the sense of Good (1975, 1979), there are a(n) = Sum_{s = 0..n-2} Sum_{t = 2..n-s} C(n, s) * C(n-s, t) * (Bell(t) - 1) kinds of independence among the variables.
Letting k = n-s, we get a(n) = Sum_{k = 2..n} Sum_{t = 2..k} C(n, n-k) * C(k, t) * (Bell(t) - 1) = Sum_{k = 2..n} Sum_{t = 2..k} C(n, k) * C(k, t) * (Bell(t) - 1).
The second formula for a(n) follows from the fact that Sum_{m = 2..k} binomial(k,m) * (Bell(m) - 1) = A058681(k) = Bell(k+1) - 2^k.
Counting all kinds of independence, as suggested by Mallows (1979), seems to be a very difficult task. For example, for n = 3, he finds a(3) = 17 kinds of independence rather than 10.
(End)

Examples

			From _Petros Hadjicostas_, Oct 10 2019: (Start)
For two r.v.'s X and Y, there is one kind of independence: f(x,y) = f(x)*f(y) (where f denotes a p.d.f., joint or marginal). Thus, a(2) = 1.
For three r.v.'s X, Y, and Z, we have
(i) 3 pairwise independence relations (f(x,y) = f(x)*f(y), or f(x,z) = f(x)*f(z), or f(y,z) = f(y)*f(z));
(ii) a factorization of the form f(x,y,z) = f(x)*f(y)*f(z);
(iii) 3 factorizations of the form f(x,y,z) = f(x,y)*f(z), or f(x,y,z) = f(y,z)*f(x), or f(x,y,z) = f(x,z)*f(y); [e.g. f(x,y,z) = f(x,y)*f(z) means random vector (X,Y) is jointly independent of variable Z]
(iv) 3 conditional factorizations f(x,y|z) = f(x|z)*f(y|z), or f(y,z|x) = f(y|x)*f(z|x), or f(x,z|y) = f(x|y)*f(z|y).
Hence, a(3) = 3 + 1 + 3 + 3 = 10. (See Mallows (1979) for some kinds of independence not considered by Good (1975, 1976).)
For four variables X, Y, Z, W, we have several cases:
(i) If we condition on a single variable, then we have 3 + 1 + 3 = 7  conditional independence cases (cases (i), (ii), and (iii) for n = 3) --> a total of C(4,1)*7 = 28.
(ii) If we condition on 2 variables, we have 1 kind of independence; i.e., a total of C(4,2)*1 = 6.
(iii) If we do not condition on any variables, we may consider:
(A) the independence of every two r.v.'s --> C(4,2) = 6 cases.
(B) the independence of every three variables: 1 + 3 = 4 factorizations (see the unconditional cases (ii) and (iii) above for the case n=3) --> a total of C(4,3)*4 = 16 factorizations.
(C) the factorizations f(x,y,z,w) = f(x)*f(y)*f(z)*f(w), or = f(x,y,z)*f(z), or = f(x,z,w)*f(y), or = f(x,y,w)*f(z), or f(y,z,w)*f(x), or = f(x,y)*f(z,w), or = f(x,z)*f(y,w), or = f(x,w)*f(y,z), or = f(x,w)*f(y)*f(z), etc. --> a total of 15-1 = 14 factorizations.
Hence, a(4) = 28 + 6 + 6 + 16 + 14 = 70.
(End)
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Magma
    [Bell(n+2) -Bell(n+1) -3^n: n in [0..30]]; // G. C. Greubel, Feb 23 2022
    
  • Mathematica
    With[{nn=30},CoefficientList[Series[Exp[Exp[x]+2x-1]-Exp[3x],{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Nov 04 2015 *)
  • Sage
    [bell_number(n+2) -bell_number(n+1) -3^n for n in (0..30)] # G. C. Greubel, Feb 23 2022

Formula

E.g.f.: exp(exp(x)+2*x-1) - exp(3*x).
From Petros Hadjicostas, Oct 10 2019: (Start)
a(n) = Sum_{k = 2..n} Sum_{m = 2..k} binomial(n,k) * binomial(k,m) * (Bell(m) - 1), where Bell(m) = A000110(m) and Bell(m) - 1 = A058692(m).
a(n) = Sum_{k = 2..n} (Bell(k+1) - 2^k) * binomial(n,k) = Sum_{k = 2..n} A058681(k)*binomial(n,k). (End)
From G. C. Greubel, Feb 23 2022: (Start)
a(n) = Sum_{j=0..n} ( binomial(n,j)*2^j*Bell(n-j) ) - 3^n [Good, Iranian J. Sci. Tech., pg. 80, eq (10)].
a(n) = Bell(n+2) - Bell(n+1) - 3^n. (End)

Extensions

More terms from N. J. A. Sloane, Jun 26 2015

A326259 MM-numbers of crossing, capturing multiset partitions (with empty parts allowed).

Original entry on oeis.org

8903, 15167, 16717, 17806, 18647, 20329, 20453, 21797, 22489, 25607, 26709, 27649, 29551, 30334, 31373, 32741, 33434, 34691, 35177, 35612, 35821, 37091, 37133, 37294, 37969, 38243, 39493, 40658, 40906, 41449, 42011, 42949, 43594, 43817, 43873, 44515, 44861
Offset: 1

Views

Author

Gus Wiseman, Jun 22 2019

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. The multiset multisystem with MM-number n is obtained by taking the multiset of prime indices of each prime index of n.
A multiset partition is crossing if it has two blocks of the form {...x...y...}, {...z...t...} where x < z < y < t or z < x < t < y. It is capturing if it has two blocks of the form {...x...y...} and {...z...t...} where x < z and y > t or x > z and y < t. Capturing is a weaker condition than nesting, so for example {{1,3,5},{2,4}} is capturing but not nesting.

Examples

			The sequence of terms together with their multiset multisystems begins:
   8903: {{1,3},{2,2,4}}
  15167: {{1,3},{2,2,5}}
  16717: {{2,4},{1,3,3}}
  17806: {{},{1,3},{2,2,4}}
  18647: {{1,3},{2,2,6}}
  20329: {{1,3},{1,2,2,4}}
  20453: {{1,2,3},{1,2,4}}
  21797: {{1,1,3},{2,2,4}}
  22489: {{1,4},{2,2,5}}
  25607: {{1,3},{2,2,7}}
  26709: {{1},{1,3},{2,2,4}}
  27649: {{1,4},{2,2,6}}
  29551: {{1,3},{2,2,8}}
  30334: {{},{1,3},{2,2,5}}
  31373: {{2,5},{1,3,3}}
  32741: {{1,3},{2,2,2,4}}
  33434: {{},{2,4},{1,3,3}}
  34691: {{1,2,3},{2,2,4}}
  35177: {{1,3},{1,2,2,5}}
  35612: {{},{},{1,3},{2,2,4}}
		

Crossrefs

Crossing set partitions are A000108.
Capturing set partitions are A326243.
Crossing, capturing set partitions are A326246.
MM-numbers of crossing multiset partitions are A324170.
MM-numbers of nesting multiset partitions are A326256.
MM-numbers of capturing multiset partitions are A326255.
MM-numbers of unsortable multiset partitions are A326258.

Programs

  • Mathematica
    croXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;xTable[PrimePi[p],{k}]]]];
    Select[Range[100000],capXQ[primeMS/@primeMS[#]]&&croXQ[primeMS/@primeMS[#]]&]

A058669 Triangle T(n,k) read by rows, giving number of matroids of rank k on n labeled points (n >= 0, 0 <= k <= n).

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 7, 7, 1, 1, 15, 36, 15, 1, 1, 31, 171, 171, 31, 1, 1, 63, 813, 2053, 813, 63, 1, 1, 127, 4012, 33442, 33442, 4012, 127, 1, 1, 255, 20891, 1022217, 8520812, 1022217, 20891, 255, 1, 1, 511
Offset: 0

Views

Author

N. J. A. Sloane, Dec 30 2000

Keywords

Examples

			Triangle T(n,k) (with rows n >= 0 and columns k >= 0) begins as follows:
  1;
  1,   1;
  1,   3,     1;
  1,   7,     7,       1;
  1,  15,    36,      15,       1;
  1,  31,   171,     171,      31,       1;
  1,  63,   813,    2053,     813,      63,     1;
  1, 127,  4012,   33442,   33442,    4012,   127,   1;
  1, 255, 20891, 1022217, 8520812, 1022217, 20891, 255, 1;
  ...
		

Crossrefs

Row sums give A058673.
Columns include (truncated versions of) A000012 (k=0), A000225 (k=1), A058681 (k=2), A058687 (k=3).

Formula

From Petros Hadjicostas, Oct 10 2019: (Start)
T(n,0) = 1 for n >= 0.
T(n,1) = 2^n - 1 for n >= 1. [Dukes (2004), Theorem 2.1 (ii).]
T(n,2) = Bell(n+1) - 2^n = A000110(n+1) - A000079(n) for n >= 2. [Dukes (2004), Theorem 2.1 (ii).]
T(n,k) = Sum_{m = k..n} binomial(n,m) * A058711(m,k) for n >= k. [Dukes (2004), see the equations before Theorem 2.1.]
(End)

A326245 Number of crossing, non-capturing set partitions of {1..n}.

Original entry on oeis.org

0, 0, 0, 0, 1, 7, 34, 141, 537, 1941, 6777, 23096, 77340
Offset: 0

Views

Author

Gus Wiseman, Jun 20 2019

Keywords

Comments

A set partition is crossing if it has two blocks of the form {...x...y...}, {...z...t...} where x < z < y < t or z < x < t < y, and capturing if it has two blocks of the form {...x...y...}, {...z...t...} where x < z < t < y or z < x < y < t. Capturing is a weaker condition than nesting, so for example {{1,3,5},{2,4}} is capturing but not nesting.

Examples

			The a(4) = 1 and a(5) = 7 set partitions:
  {{1,3},{2,4}}  {{1,2,4},{3,5}}
                 {{1,3},{2,4,5}}
                 {{1},{2,4},{3,5}}
                 {{1,3},{2,4},{5}}
                 {{1,3},{2,5},{4}}
                 {{1,4},{2},{3,5}}
                 {{1,4},{2,5},{3}}
		

Crossrefs

Crossing set partitions are A016098.
Non-capturing set partitions are A054391.
Crossing, capturing set partitions are A326246.

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    croXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x
    				

A340264 T(n, k) = Sum_{j=0..k} binomial(n, k - j)*Stirling2(n - k + j, j). Triangle read by rows, 0 <= k <= n.

Original entry on oeis.org

1, 0, 2, 0, 1, 4, 0, 1, 6, 8, 0, 1, 11, 24, 16, 0, 1, 20, 70, 80, 32, 0, 1, 37, 195, 340, 240, 64, 0, 1, 70, 539, 1330, 1400, 672, 128, 0, 1, 135, 1498, 5033, 7280, 5152, 1792, 256, 0, 1, 264, 4204, 18816, 35826, 34272, 17472, 4608, 512
Offset: 0

Views

Author

Peter Luschny, Jan 08 2021

Keywords

Comments

A006905(n) = Sum_{k=0..n} A001035(k) * T(n, k). - Michael Somos, Jul 18 2021
T(n, k) is the number of idempotent relations R on [n] containing exactly k strongly connected components such that the following conditional statement holds for all x, y in [n]: If x, y are in distinct strongly connected components of R then (x, y) is not in R. - Geoffrey Critzer, Jan 10 2024

Examples

			[0] 1;
[1] 0, 2;
[2] 0, 1,   4;
[3] 0, 1,   6,    8;
[4] 0, 1,  11,   24,    16;
[5] 0, 1,  20,   70,    80,    32;
[6] 0, 1,  37,  195,   340,   240,    64;
[7] 0, 1,  70,  539,  1330,  1400,   672,   128;
[8] 0, 1, 135, 1498,  5033,  7280,  5152,  1792,  256;
[9] 0, 1, 264, 4204, 18816, 35826, 34272, 17472, 4608, 512;
		

Crossrefs

Sum of row(n) is A000110(n+1).
Sum of row(n) - 2^n is A058681(n).
Alternating sum of row(n) is A109747(n).

Programs

  • Maple
    egf := exp(t*(exp(-x) - x - 1));
    ser := series(egf, x, 22):
    p := n -> coeff(ser, x, n);
    seq(seq((-1)^n*n!*coeff(p(n), t, k), k=0..n), n = 0..10);
    # Alternative:
    T := (n, k) -> add(binomial(n, k - j)*Stirling2(n - k + j, j), j=0..k):
    seq(seq(T(n, k), k = 0..n), n=0..9); # Peter Luschny, Feb 09 2021
  • Mathematica
    T[ n_, k_] := Sum[ Binomial[n, k-j] StirlingS2[n-k+j, j], {j, 0 ,k}]; (* Michael Somos, Jul 18 2021 *)
  • PARI
    T(n, k) = sum(j=0, k, binomial(n, j)*stirling(n-j, k-j, 2)); /* Michael Somos, Jul 18 2021 */

Formula

T(n, k) = (-1)^n * n! * [t^k] [x^n] exp(t*(exp(-x) - x - 1)).
n-th row polynomial R(n,x) = exp(-x)*Sum_{k >= 0} (x + k)^n * x^k/k! = Sum_{k = 0..n} binomial(n,k)*Bell(k,x)*x^(n-k), where Bell(n,x) denotes the n-th Bell polynomial. - Peter Bala, Jan 13 2022

Extensions

New name from Peter Luschny, Feb 09 2021

A326260 MM-numbers of capturing, non-nesting multiset partitions (with empty parts allowed).

Original entry on oeis.org

2599, 4163, 5198, 6463, 6893, 7291, 7797, 8326, 8507, 9131, 9959, 10396, 10649, 11041, 11639, 12489, 12811, 12926, 12995, 13786, 14237, 14582, 14899, 15157, 15594, 16123, 16403, 16652, 17014, 17063, 17089, 17141, 18101, 18193, 18262, 18643, 18659, 19337, 19389
Offset: 1

Views

Author

Gus Wiseman, Jun 22 2019

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. The multiset multisystem with MM-number n is obtained by taking the multiset of prime indices of each prime index of n.
A set partition is capturing if it has two blocks of the form {...x...y...} and {...z...t...} where x < z and y > t or x > z and y < t. It is nesting if it has two blocks of the form {...x,y...} and {...z,t...} where x < z and y > t or x > z and y < t. Capturing is a weaker condition than nesting, so for example {{1,3,5},{2,4}} is capturing but not nesting.

Examples

			The sequence of terms together with their multiset multisystems begins:
   2599: {{2,2},{1,2,3}}
   4163: {{2,2},{1,2,4}}
   5198: {{},{2,2},{1,2,3}}
   6463: {{2,2},{1,1,2,3}}
   6893: {{1,2,2},{1,2,3}}
   7291: {{2,2},{1,2,5}}
   7797: {{1},{2,2},{1,2,3}}
   8326: {{},{2,2},{1,2,4}}
   8507: {{2,3},{1,2,4}}
   9131: {{2,2},{1,2,6}}
   9959: {{2,2},{1,1,2,4}}
  10396: {{},{},{2,2},{1,2,3}}
  10649: {{2,2},{1,2,2,3}}
  11041: {{1,2,2},{1,2,4}}
  11639: {{2,2,2},{1,2,3}}
  12489: {{1},{2,2},{1,2,4}}
  12811: {{2,2},{1,2,7}}
  12926: {{},{2,2},{1,1,2,3}}
  12995: {{2},{2,2},{1,2,3}}
  13786: {{},{1,2,2},{1,2,3}}
		

Crossrefs

Non-nesting set partitions are A000108.
Capturing set partitions are A326243.
Capturing, non-nesting set partitions are A326249.
MM-numbers of nesting multiset partitions are A326256.
MM-numbers of capturing multiset partitions are A326255.

Programs

  • Mathematica
    capXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x_,{_,x_,y_,_},_,{_,z_,t_,_},_}/;xTable[PrimePi[p],{k}]]]];
    Select[Range[10000],!nesXQ[primeMS/@primeMS[#]]&&capXQ[primeMS/@primeMS[#]]&]
Previous Showing 11-20 of 26 results. Next