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

A000699 Number of irreducible chord diagrams with 2n nodes.

Original entry on oeis.org

1, 1, 1, 4, 27, 248, 2830, 38232, 593859, 10401712, 202601898, 4342263000, 101551822350, 2573779506192, 70282204726396, 2057490936366320, 64291032462761955, 2136017303903513184, 75197869250518812754, 2796475872605709079512, 109549714522464120960474, 4509302910783496963256400, 194584224274515194731540740
Offset: 0

Views

Author

Keywords

Comments

Perturbation expansion in quantum field theory: spinor case in 4 spacetime dimensions.
a(n)*2^(-n) is the coefficient of the x^(2*n-1) term in the series reversal of the asymptotic expansion of 2*DawsonF(x) = sqrt(Pi)*exp(-x^2)*erfi(x) for x -> inf. - Vladimir Reshetnikov, Apr 23 2016
The September 2018 talk by Noam Zeilberger (see link to video) connects three topics (planar maps, Tamari lattices, lambda calculus) and eight sequences: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827. - N. J. A. Sloane, Sep 17 2018
A set partition is topologically connected if the graph whose vertices are the blocks and whose edges are crossing pairs of blocks is connected, where two blocks cross each other if they are of the form {{...x...y...},{...z...t...}} for some x < z < y < t or z < x < t < y. Then a(n) is the number of topologically connected 2-uniform set partitions of {1...2n}. See my links for examples. - Gus Wiseman, Feb 23 2019
From Julien Courtiel, Oct 09 2024: (Start)
a(n) is the number of rooted bridgeless combinatorial maps with n edges (genus is not fixed). A map is bridgeless if it has no edge whose removal disconnects the graph. For example, for n=2, there are 4 bridgeless maps with 2 edges: 2 planar maps with 1 vertex (either two consecutive loops, or two nested loops), 1 toric map with 1 vertex, and 1 planar map with 2 vertices connected by a double edge.
Also, a(n) is the number of trees with n edges equipped with a binary tubing. A tube is a connected subgraph. A binary tubing of a tree is a nested set collection S of tubes such that 1. S contains the tube of all vertices 2. Every tube of S is either reduced to one vertex, or it can be can partitioned by 2 tubes of S.
(End)

Examples

			a(31)=627625976637472254550352492162870816129760 was computed using Kreimer's Hopf algebra of rooted trees. It subsumes 2.6*10^21 terms in quantum field theory.
G.f.: A(x) = 1 + x + x^2 + 4*x^3 + 27*x^4 + 248*x^5 + 2830*x^6 +...
where d/dx (A(x) - 1)^2/x = 1 + 4*x + 27*x^2 + 248*x^3 + 2830*x^4 +...
		

References

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

Crossrefs

Sequences mentioned in the Noam Zeilberger 2018 video: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827.
Cf. A004300, A051862, A212273. Column sums of A232223. First column of A322402.

Programs

  • Maple
    A000699 := proc(n)
        option remember;
        if n <= 1 then
            1;
        else
            add((2*i-1)*procname(i)*procname(n-i),i=1..n-1) ;
        end if;
    end proc:
    seq(A000699(n),n=0..30) ; # R. J. Mathar, Jun 12 2018
  • Mathematica
    terms = 22; A[] = 0; Do[A[x] = x + x^2 * D[A[x]^2/x, x] + O[x]^(terms+1) // Normal, terms]; CoefficientList[A[x], x] // Rest (* Jean-François Alcover, Apr 06 2012, after Paul D. Hanna, updated Jan 11 2018 *)
    a = ConstantArray[0,20]; a[[1]]=1; Do[a[[n]] = (n-1)*Sum[a[[i]]*a[[n-i]],{i,1,n-1}],{n,2,20}]; a (* Vaclav Kotesovec, Feb 22 2014 *)
    Module[{max = 20, s}, s = InverseSeries[ComplexExpand[Re[Series[2 DawsonF[x], {x, Infinity, 2 max + 1}]]]]; Table[SeriesCoefficient[s, 2 n - 1] 2^n, {n, 1, max}]] (* Vladimir Reshetnikov, Apr 23 2016 *)
  • PARI
    {a(n)=local(A=1+x*O(x^n)); for(i=1, n, A=1+x+x^2*deriv((A-1)^2/x)+x*O(x^n)); polcoeff(A, n)} \\ Paul D. Hanna, Dec 31 2010 [Modified to include a(0) = 1. - Paul D. Hanna, Nov 06 2020]
    
  • PARI
    {a(n) = my(A); A = 1+O(x) ; for( i=0, n, A = 1+x + (A-1)*(2*x*A' - A + 1)); polcoeff(A, n)}; /* Michael Somos, May 12 2012 [Modified to include a(0) = 1. - Paul D. Hanna, Nov 06 2020] */
    
  • PARI
    seq(N) = {
      my(a = vector(N)); a[1] = 1;
      for (n=2, N, a[n] = sum(k=1, n-1, (2*k-1)*a[k]*a[n-k])); a;
    };
    seq(22)  \\ Gheorghe Coserea, Jan 22 2017
    
  • PARI
    seq(n)={my(g=serlaplace(1 / sqrt(1 - 2*x + O(x*x^n)))); Vec(sqrt((x/serreverse( x*g^2 ))))} \\ Andrew Howroyd, Nov 21 2024
    
  • Python
    def A000699_list(n):
        list = [1, 1] + [0] * (n - 1)
        for i in range(2, n + 1):
            list[i] = (i - 1) * sum(list[j] * list[i - j] for j in range(1, i))
        return list
    print(A000699_list(22)) # M. Eren Kesim, Jun 23 2021

Formula

a(n) = (n-1)*Sum_{i=1..n-1} a(i)*a(n-i) for n > 1, with a(1) = a(0) = 1. [Modified to include a(0) = 1. - Paul D. Hanna, Nov 06 2020]
A212273(n) = n * a(n). - Michael Somos, May 12 2012
G.f. satisfies: A(x) = 1 + x + x^2*[d/dx (A(x) - 1)^2/x]. - Paul D. Hanna, Dec 31 2010 [Modified to include a(0) = 1. - Paul D. Hanna, Nov 06 2020]
a(n) ~ n^n * 2^(n+1/2) / exp(n+1) * (1 - 31/(24*n) - 2207/(1152*n^2) - 3085547/(414720*n^3) - 1842851707/(39813120*n^4) - ...). - Vaclav Kotesovec, Feb 22 2014, extended Oct 23 2017
G.f. A(x) satisfies: 1 = A(x) - x/(A(x) - 2*x/(A(x) - 3*x/(A(x) - 4*x/(A(x) - 5*x/(A(x) - ...))))), a continued fraction relation. - Paul D. Hanna, Nov 04 2020
G.f. A(x) satisfies: A(x*B(x)^2) = B(x) where B(x) is the g.f. of A001147. - Andrew Howroyd, Nov 21 2024

Extensions

More terms from David Broadhurst, Dec 14 1999
Inserted "chord" in definition. - N. J. A. Sloane, Jan 19 2017
Added a(0)=1. - N. J. A. Sloane, Nov 05 2020
Modified formulas slightly to include a(0) = 1. - Paul D. Hanna, Nov 06 2020

A007297 Number of connected graphs on n labeled nodes on a circle with straight-line edges that don't cross.

Original entry on oeis.org

1, 1, 4, 23, 156, 1162, 9192, 75819, 644908, 5616182, 49826712, 448771622, 4092553752, 37714212564, 350658882768, 3285490743987, 30989950019532, 294031964658430, 2804331954047160, 26870823304476690, 258548658860327880
Offset: 1

Views

Author

Keywords

Comments

Apart from the initial 1, reversion of g.f. for A162395 (squares with signs): see A263843.

Examples

			G.f. = x*(1 + x + 4*x^2 + 23*x^3 + 156*x^4 + 1162*x^5 + 9192*x^6 + 75819*x^7 + ...).
		

References

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

Crossrefs

Cf. A162395, A000290. 4th row of A107111. Row sums of A089434.
See A263843 for a variant.
Cf. A000108 (non-crossing set partitions), A001006, A001187, A054726 (non-crossing graphs), A054921, A099947, A194560, A293510, A323818, A324167, A324169, A324173.

Programs

  • Maple
    A007297:=proc(n) if n = 1 then 1 else add(binomial(3*n - 3, n + j)*binomial(j - 1, j - n + 1), j = n - 1 .. 2*n - 3)/(n - 1); fi; end;
  • Mathematica
    CoefficientList[ InverseSeries[ Series[(x-x^2)/(1+x)^3, {x, 0, 20}], x], x] // Rest (* From Jean-François Alcover, May 19 2011, after PARI prog. *)
    Table[Binomial[3n, 2n+1] Hypergeometric2F1[1-n, n, 2n+2, -1]/n, {n, 1, 20}] (* Vladimir Reshetnikov, Oct 25 2015 *)
  • PARI
    a(n)=if(n<0,0,polcoeff(serreverse((x-x^2)/(1+x)^3+O(x^(n+2))),n+1)) /* Ralf Stephan */

Formula

Apart from initial term, g.f. is the series reversion of (x-x^2)/(1+x)^3 (A162395). See A263843. - Vladimir Kruchinin, Feb 08 2013
G.f.: (g-z)/z, where g=-1/3+(2/3)*sqrt(1+9z)*sin((1/3)*arcsin((2+27z+54z^2)/2/(1+9*z)^(3/2))). - Emeric Deutsch, Dec 02 2002
a(n) = (1/n)*Sum_{k=0..n} binomial(3n, n-k-1)*binomial(n+k-1, k). - Paul Barry, May 11 2005
a(n) = 4^(n-1)*(Gamma(3*n/2-1)/Gamma(n/2+1)/Gamma(n) -Gamma((3*n-1)/2)/ Gamma( (n+1)/2)/Gamma(n+1)). - Mark van Hoeij, Aug 27 2005, adapted to offset Feb 21 2020 by R. J. Mathar
a(n) = 4^n * binomial(3*n/2, n/2) / (9*n-6) - 4^(n-1) * binomial(3*(n-1)/2, (n-1)/2 ) / n. - Mark van Hoeij, Aug 27 2005, adapted to offset Feb 21 2020 by R. J. Mathar
D-finite with recurrence: n*(n-1)*(3*n-4)*a(n) +36*(n-1)*a(n-1) -12*(3*n-8)*(3*n-1)*(3*n-7)*a(n-2)=0. - Mark van Hoeij, Aug 27 2005, adapted to offset Feb 21 2020 by R. J. Mathar
a(n) = (1/n)*Sum_{k=0..n} C(3n, k)*C(2n-k-2, n-1). - Paul Barry, Sep 27 2005
a(n) ~ (2-sqrt(3)) * 6^n * 3^(n/2) / (sqrt(2*Pi) * n^(3/2)). - Vaclav Kotesovec, Mar 17 2014
a(n) = binomial(3*n,2*n+1)*hypergeom([1-n,n], [2*n+2], -1)/n. - Vladimir Reshetnikov, Oct 25 2015
a(n) = 2*A078531(n) - A085614(n+1). - Vladimir Reshetnikov, Apr 24 2016

Extensions

Better description from Philippe Flajolet, Apr 20 2000
More terms from James Sellers, Aug 21 2000
Definition revised and initial a(1)=1 added by N. J. A. Sloane, Nov 05 2015 at the suggestion of Axel Boldt. Some of the formulas may now need to be adjusted slightly.

A324170 Numbers whose multiset multisystem (A302242) is crossing.

Original entry on oeis.org

2117, 3973, 4234, 4843, 5183, 5249, 5891, 6351, 6757, 7181, 7801, 7946, 8249, 8468, 8903, 9193, 9686, 9727, 10019, 10063, 10366, 10498, 10585, 11051, 11513, 11567, 11782, 11857, 11919, 12557, 12629, 12702, 12851, 13021, 13193, 13459, 13514, 13631, 14123, 14362
Offset: 1

Views

Author

Gus Wiseman, Feb 17 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 of n is obtained by taking the multiset of prime indices of each prime index of n.
A multiset of multisets is crossing if it contains a 2-element submultiset of the form {{...x...y...}, {...z...t...}} where x < z < y < t or z < x < t < y.

Examples

			The sequence of terms together with their multiset multisystems begins:
  2117: {{1,3},{2,4}}
  3973: {{1,3},{2,5}}
  4234: {{},{1,3},{2,4}}
  4843: {{1,3},{2,6}}
  5183: {{1,1,3},{2,4}}
  5249: {{1,3},{1,2,4}}
  5891: {{1,4},{2,5}}
  6351: {{1},{1,3},{2,4}}
  6757: {{1,3},{2,7}}
  7181: {{1,4},{2,6}}
  7801: {{1,3},{2,8}}
  7946: {{},{1,3},{2,5}}
  8249: {{2,4},{1,2,3}}
  8468: {{},{},{1,3},{2,4}}
  8903: {{1,3},{2,2,4}}
  9193: {{1,3},{1,2,5}}
  9686: {{},{1,3},{2,6}}
  9727: {{1,1,3},{2,5}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    croXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x
    				

A324169 Number of labeled graphs covering the vertex set {1,...,n} with no crossing edges.

Original entry on oeis.org

1, 0, 1, 4, 25, 176, 1353, 11012, 93329, 815104, 7285489, 66324644, 612863337, 5733381616, 54195878137, 516852285668, 4966883732129, 48049936644736, 467566946973537, 4573486005681092, 44942852084894777, 443484037981300144, 4392621673072766505
Offset: 0

Views

Author

Gus Wiseman, Feb 17 2019

Keywords

Comments

Two edges {x,y}, {z,t} are crossing if either x < z < y < t or z < x < t < y. If the vertices are arranged in a circle, this is equivalent to crossing of chords.
Covering means there are no isolated vertices.

Crossrefs

Cf. A000108, A000124, A001006, A001764, A003465, A007297 (connected case), A016098, A054726 (non-crossing graphs), A099947, A306438.

Programs

  • Mathematica
    nn=8;
    croXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x
    				
  • PARI
    seq(n)=Vec((2 + 7*x + 3*x^2 - x*sqrt(1 - 10*x - 7*x^2 + O(x^n)))/(2*(1 + x)^3)) \\ Andrew Howroyd, Jan 20 2023

Formula

Inverse binomial transform of A054726.
G.f.: (2 + 7*x + 3*x^2 - x*sqrt(1 - 10*x - 7*x^2))/(2*(1 + x)^3). - Andrew Howroyd, Jan 20 2023

A324173 Regular triangle read by rows where T(n,k) is the number of set partitions of {1,...,n} with k topologically connected components.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 3, 1, 0, 2, 6, 6, 1, 0, 6, 15, 20, 10, 1, 0, 21, 51, 65, 50, 15, 1, 0, 85, 203, 252, 210, 105, 21, 1, 0, 385, 912, 1120, 938, 560, 196, 28, 1, 0, 1907, 4527, 5520, 4620, 2898, 1302, 336, 36, 1, 0, 10205, 24370, 29700, 24780, 15792, 7812, 2730, 540, 45, 1
Offset: 0

Views

Author

Gus Wiseman, Feb 17 2019

Keywords

Comments

A set partition is crossing if it contains a pair of blocks of the form {{...x...y...}, {...z...t...}} where x < z < y < t or z < x < t < y.
The topologically connected components of a set partition correspond to the blocks of its minimal non-crossing coarsening.

Examples

			Triangle begins:
     1
     0     1
     0     1     1
     0     1     3     1
     0     2     6     6     1
     0     6    15    20    10     1
     0    21    51    65    50    15     1
     0    85   203   252   210   105    21     1
     0   385   912  1120   938   560   196    28     1
     0  1907  4527  5520  4620  2898  1302   336    36     1
     0 10205 24370 29700 24780 15792  7812  2730   540    45     1
Row n = 4 counts the following set partitions:
  {{1234}}    {{1}{234}}  {{1}{2}{34}}  {{1}{2}{3}{4}}
  {{13}{24}}  {{12}{34}}  {{1}{23}{4}}
              {{123}{4}}  {{12}{3}{4}}
              {{124}{3}}  {{1}{24}{3}}
              {{134}{2}}  {{13}{2}{4}}
              {{14}{23}}  {{14}{2}{3}}
		

Crossrefs

Programs

  • Mathematica
    croXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x0]&]},If[c=={},s,csm[Sort[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    crosscmpts[stn_]:=csm[Union[Subsets[stn,{1}],Select[Subsets[stn,{2}],croXQ]]];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    Table[Length[Select[sps[Range[n]],Length[crosscmpts[#]]==k&]],{n,0,8},{k,0,n}]

A136653 G.f.: A(x) satisfies: coefficient of x^n in A(x)^(n+1)/(n+1) = 2^(n*(n-1)/2).

Original entry on oeis.org

1, 1, 1, 4, 39, 748, 27162, 1880872, 252273611, 66358216668, 34506398937158, 35644762692112792, 73356520492898454022, 301274559225693420690360, 2471654510727312089903896948, 40527708183358718551543295827536, 1328579216048284168977214446788083699
Offset: 0

Views

Author

Paul D. Hanna, Jan 15 2008

Keywords

Comments

a(n) is the number of graphs on vertices 1,...,n such that, when these vertices are arranged counterclockwise around a circle and edges are drawn as straight line segments, the resulting diagram is connected. - Jonathan Novak (j2novak(AT)math.uwaterloo.ca), Apr 30 2010
In this interpretation, both intersecting (set theoretically) and crossing (topologically) edges are considered connected. - Gus Wiseman, Feb 23 2019

Examples

			G.f.: A(x) = 1 + x + x^2 + 4*x^3 + 39*x^4 + 748*x^5 + 27162*x^6 +...
Let F(x) = 1 + x + 2*x^2 + 8*x^3 + 64*x^4 + 1024*x^5 +...+ 2^(n*(n-1)/2)*x^n +..
then A(x) = F(x/A(x)), A(x*F(x)) = F(x).
Coefficient of x^n in A(x)^(n+1)/(n+1) = 2^(n*(n-1)/2),
as can be seen by the main diagonal in the array of
coefficients in the initial powers of A(x):
A^1: [(1), 1, 1, 4, 39, 748, 27162, 1880872, 252273611,...;
A^2: [1, (2), 3, 10, 87, 1582, 55914, 3817876, 508370795,...;
A^3: [1, 3, (6), 19, 147, 2517, 86398, 5813550, 768378627,...;
A^4: [1, 4, 10, (32), 223, 3572, 118778, 7870640, 1032387787,...;
A^5: [1, 5, 15, 50, (320), 4771, 153245, 9992130, 1300492845,...;
A^6: [1, 6, 21, 74, 444, (6144), 190023, 12181278, 1572792585,...;
A^7: [1, 7, 28, 105, 602, 7728, (229376), 14441659, 1849390375,...;
A^8: [1, 8, 36, 144, 802, 9568, 271616, (16777216), 2130394591,...;
A^9: [1, 9, 45, 192, 1053, 11718, 317112, 19192320, (2415919104),...;
dividing each diagonal term in row n by (n+1) gives 2^(n*(n-1)/2).
The diagonal above the main diagonal gives coefficients of l.g.f.:
log(F(x)) = x + 3*x^2/2 + 19*x^3/3 + 223*x^4/4 + 4771*x^5/5 +...
		

Crossrefs

Programs

  • Mathematica
    max = 15; s = x*Sum[2^(k*(k-1)/2)*x^k, {k, 0, max}] + O[x]^(max+2); x/InverseSeries[s] + O[x]^(max+1) // CoefficientList[#, x]& (* Jean-François Alcover, Sep 03 2017 *)
    croXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x0]&]},If[c=={},s,csm[Sort[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    bicmpts[stn_]:=csm[Union[Subsets[stn,{1}],Select[Subsets[stn,{2}],Intersection@@#!={}&],Select[Subsets[stn,{2}],croXQ]]];
    Table[Length[Select[Subsets[Subsets[Range[n],{2}]],And[Union@@#==Range[n],Length[bicmpts[#]]<=1]&]],{n,0,5}] (* Gus Wiseman, Feb 23 2019 *)
  • PARI
    a(n)=polcoeff(x/serreverse(x*sum(k=0,n,2^(k*(k-1)/2)*x^k +x*O(x^n))),n)

Formula

G.f.: A(x) = x/Series_Reversion( x*Sum_{k=0..n} 2^(k(k-1)/2)*x^k ).
Equals the free cumulant sequence corresponding to A006125. - Jonathan Novak (j2novak(AT)math.uwaterloo.ca), Apr 30 2010

Extensions

Name changed and part of prior name moved to formula section by Paul D. Hanna, Sep 19 2013

A324167 Number of non-crossing antichain covers of {1,...,n}.

Original entry on oeis.org

1, 1, 2, 9, 67, 633, 6763, 77766, 938957, 11739033, 150649945, 1973059212, 26265513030, 354344889798, 4833929879517, 66568517557803, 924166526830701, 12920482325488761, 181750521972603049, 2570566932237176232, 36532394627404815308, 521439507533582646156
Offset: 0

Views

Author

Gus Wiseman, Feb 17 2019

Keywords

Comments

An antichain is non-crossing if no pair of distinct parts is of the form {{...x...y...}, {...z...t...}} where x < z < y < t or z < x < t < y.

Examples

			The a(3) = 9 antichains:
  {{1,2,3}}
  {{1},{2,3}}
  {{2},{1,3}}
  {{3},{1,2}}
  {{1,2},{1,3}}
  {{1,2},{2,3}}
  {{1,3},{2,3}}
  {{1},{2},{3}}
  {{1,2},{1,3},{2,3}}
		

Crossrefs

Cf. A000108, A000124, A000372 (antichains), A001006, A006126 (antichain covers), A014466, A048143, A054726 (non-crossing graphs), A099947, A261005, A283877, A306438.
Cf. A324166, A324168, A324169, A324170, A324171, A324173, A359984 (no singletons).

Programs

  • Mathematica
    nn=6;
    croXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x
    				
  • PARI
    seq(n)={my(f=O(1)); for(n=2, n, f = 1 + (4*x + x^2)*f^2 - 3*x^2*(1 + x)*f^3); Vec(subst(x*(1 + x^2*f^2 - 3*x^3*f^3), x, x/(1-x))/x) } \\ Andrew Howroyd, Jan 20 2023

Formula

Inverse binomial transform of A324168.
Binomial transform of A359984. - Andrew Howroyd, Jan 20 2023

Extensions

Terms a(9) and beyond from Andrew Howroyd, Jan 20 2023

A326293 Number of non-nesting, topologically connected simple graphs with vertices {1..n}.

Original entry on oeis.org

1, 1, 2, 4, 8, 27, 192, 1750
Offset: 0

Views

Author

Gus Wiseman, Jun 29 2019

Keywords

Comments

Two edges {a,b}, {c,d} are crossing if a < c < b < d or c < a < d < b, and nesting if a < c < d < b or c < a < b < d. A graph with positive integer vertices is topologically connected if the graph whose vertices are the edges and whose edges are crossing pairs of edges is connected.

Crossrefs

The inverse binomial transform is the covering case A326349.
Topologically connected simple graphs are A324328.
Non-crossing simple graphs are A054726.
Topologically connected set partitions are A099947.

Programs

  • Mathematica
    croXQ[eds_]:=MatchQ[eds,{_,{x_,y_},_,{z_,t_},_}/;x_,{x_,y_},_,{z_,t_},_}/;x0]&]},If[c=={},s,csm[Sort[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
    Table[Length[Select[Subsets[Subsets[Range[n],{2}]],!nesXQ[#]&&Length[csm[Union[Subsets[#,{1}],Select[Subsets[#,{2}],croXQ]]]]<=1&]],{n,0,5}]

A324166 Number of totally crossing set partitions of {1,...,n}.

Original entry on oeis.org

1, 1, 1, 1, 2, 6, 18, 57, 207, 842, 3673, 17062, 84897
Offset: 0

Views

Author

Gus Wiseman, Feb 17 2019

Keywords

Comments

A set partition is totally crossing if every pair of distinct blocks is of the form {{...x...y...}, {...z...t...}} for some x < z < y < t or z < x < t < y.

Examples

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

Crossrefs

Cf. A000108 (non-crossing partitions), A000110, A000296, A002662, A016098 (crossing partitions), A054726, A099947 (topologically connected partitions), A305854, A306006, A306418, A306438, A319752.

Programs

  • Mathematica
    nn=6;
    nonXQ[stn_]:=!MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x
    				

A324171 Number of non-crossing multiset partitions of normal multisets of size n.

Original entry on oeis.org

1, 1, 4, 16, 75, 378, 2042, 11489, 66697
Offset: 0

Views

Author

Gus Wiseman, Feb 17 2019

Keywords

Comments

A multiset is normal if its union is an initial interval of positive integers.
A multiset partition is crossing if it has a 2-element submultiset of the form {{...x...y...}, {...z...t...}} where x < z < y < t or z < x < t < y.

Examples

			The A255906(5) - a(5) = 22 crossing multiset partitions:
  {{13}{124}}  {{1}{13}{24}}
  {{13}{224}}  {{1}{24}{35}}
  {{13}{234}}  {{2}{13}{24}}
  {{13}{244}}  {{2}{14}{35}}
  {{13}{245}}  {{3}{13}{24}}
  {{14}{235}}  {{3}{14}{25}}
  {{24}{113}}  {{4}{13}{24}}
  {{24}{123}}  {{4}{13}{25}}
  {{24}{133}}  {{5}{13}{24}}
  {{24}{134}}
  {{24}{135}}
  {{25}{134}}
  {{35}{124}}
		

Crossrefs

Cf. A000108 (non-crossing set partitions), A000124, A001006, A001055, A001263, A007297, A054726 (non-crossing graphs), A099947, A194560, A255906 (multiset partitions of normal multisets), A306438.

Programs

  • Mathematica
    nonXQ[stn_]:=!MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    allnorm[n_]:=If[n<=0,{{}},Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1]];
    Table[Sum[Length[Select[mps[m],nonXQ]],{m,allnorm[n]}],{n,0,8}]
Showing 1-10 of 35 results. Next