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 31-40 of 54 results. Next

A326279 Number of labeled n-vertex simple graphs containing either a crossing or a nesting pair of edges.

Original entry on oeis.org

0, 0, 0, 0, 28, 864, 32064, 2094064
Offset: 0

Views

Author

Gus Wiseman, Jun 23 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.

Examples

			The a(4) = 28 edge-sets:
  {13,24}  {12,13,24}  {12,13,14,23}  {12,13,14,23,24}  {12,13,14,23,24,34}
  {14,23}  {12,14,23}  {12,13,14,24}  {12,13,14,23,34}
           {13,14,23}  {12,13,23,24}  {12,13,14,24,34}
           {13,14,24}  {12,13,24,34}  {12,13,23,24,34}
           {13,23,24}  {12,14,23,24}  {12,14,23,24,34}
           {13,24,34}  {12,14,23,34}  {13,14,23,24,34}
           {14,23,24}  {13,14,23,24}
           {14,23,34}  {13,14,23,34}
                       {13,14,24,34}
                       {13,23,24,34}
                       {14,23,24,34}
		

Crossrefs

Crossing and nesting simple graphs are (both) A326210, while non-crossing, non-nesting simple graphs are A326244.

Programs

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

Formula

A006125(n) = a(n) + A326244(n).

A268814 Number of purely crossing partitions of [n].

Original entry on oeis.org

1, 0, 0, 0, 1, 0, 5, 14, 62, 298, 1494, 8140, 47146, 289250, 1873304, 12756416, 91062073, 679616480, 5290206513, 42858740990, 360686972473, 3147670023632, 28439719809159, 265647698228954, 2561823514680235, 25475177517626196, 260922963832247729, 2749617210928715246
Offset: 0

Views

Author

Michel Marcus, Feb 14 2016

Keywords

Comments

For the definition of a purely crossing partition refer to Dykema link (see PC(n) Definition 1.2 and Table 2).
From Gus Wiseman, Feb 23 2019: (Start)
For n >= 1, a set partition of {1,...,n} is purely crossing if it is topologically connected (A099947), has no successive elements in the same block (A000110(n - 1)), and the first and last vertices belong to different blocks (A005493(n - 2)). For example, the a(4) = 1, a(6) = 5, and a(7) = 14 purely crossing set partitions are:
{{13}{24}} {{135}{246}} {{13}{246}{57}}
{{13}{25}{46}} {{13}{257}{46}}
{{14}{25}{36}} {{135}{26}{47}}
{{14}{26}{35}} {{135}{27}{46}}
{{15}{24}{36}} {{136}{24}{57}}
{{136}{25}{47}}
{{14}{257}{36}}
{{14}{26}{357}}
{{146}{25}{37}}
{{146}{27}{35}}
{{15}{246}{37}}
{{15}{247}{36}}
{{16}{24}{357}}
{{16}{247}{35}}
(End)

Examples

			G.f.: A(x) = 1 + x^4 + 5*x^6 + 14*x^7 + 62*x^8 + 298*x^9 + 1494*x^10 + 8140*x^11 + 47146*x^12 +...
		

Crossrefs

Programs

  • Mathematica
    n = 30; F = x*Sum[BellB[k] x^k, {k, 0, n}] + O[x]^n; B = ComposeSeries[1/( InverseSeries[F, w]/w)-1, x/(1+x) + O[x]^n]; A = (B-x)/(1+x); Join[{1}, CoefficientList[A, x] // Rest] (* Jean-François Alcover, Feb 23 2016, adapted from K. J. Dykema's code *)
    intvQ[set_]:=Or[set=={},Sort[set]==Range[Min@@set,Max@@set]];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    Table[Length[Select[sps[Range[n]],And[!MatchQ[#,{_,{_,x_,y_,_},_}/;x+1==y],#=={}||And@@Not/@intvQ/@Union@@@Subsets[#,{1,Length[#]-1}],#=={}||Position[#,1][[1,1]]!=Position[#,n][[1,1]]]&]],{n,0,10}] (* Gus Wiseman, Feb 23 2019 *)
  • PARI
    lista(nn) = {c = x/serreverse(x*serlaplace(exp(exp(x+x*O(x^nn)) -1))); b = subst(c, x, x/(1+x)+ O(x^nn)); vb = Vec(b-1); va = vector(#vb); va[1] = 0; va[2] = 0; for (k=3, #va, va[k] = vb[k] - va[k-1]; ); concat(1, va); }
    
  • PARI
    {a(n) = my(A=1+x^3); for(i=1, n, A = sum(m=0, n, x^m/prod(k=1, m, (1+x)^2*A - k*x +x*O(x^n)) )/(1+x) ); polcoeff( A, n)}
    for(n=0,35,print1(a(n),", ")) \\ Paul D. Hanna, Mar 07 2016
    
  • PARI
    {Stirling2(n, k) = n!*polcoeff(((exp(x+x*O(x^n)) - 1)^k)/k!, n)}
    {Bell(n) = sum(k=0,n, Stirling2(n, k) )}
    {a(n) = my(A=1+x); for(i=1, n, A = sum(m=0, n, Bell(m)*x^m/((1+x +x*O(x^n))^(2*m+1)*A^m)) ); polcoeff(A, n)}
    for(n=0,25,print1(a(n),", ")) \\ Paul D. Hanna, Mar 07 2016

Formula

G.f.: G(x) satisfies B(x) = x + (1 + x)*G(x) where B(x) is the g.f. of A268815 (see A(x) in Dykema link p. 7).
From Paul D. Hanna, Mar 07 2016: (Start)
O.g.f. A(x) satisfies:
(1) A(x) = Sum_{n>=0} A000110(n)*x^n / ((1+x)^(2*n+1) * A(x)^n), where A000110 are the Bell numbers.
(2) A(x) = 1/(1+x) * Sum_{n>=0} x^n / Product_{k=1..n} ((1+x)^2*A(x) - k*x).
(3) A(x) = 1/(1+x - x/((1+x)*A(x) - 1*x/(1+x - x/((1+x)*A(x) - 2*x/(1+x - x/((1+x)*A(x) - 3*x/(1+x - x/((1+x)*A(x) - 4*x/(1+x - x/((1+x)*A(x) -...)))))))))), a continued fraction. (End)

A322402 Triangle read by rows: The number of chord diagrams with n chords and k topologically connected components, 0 <= k <= n.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 4, 6, 5, 0, 27, 36, 28, 14, 0, 248, 310, 225, 120, 42, 0, 2830, 3396, 2332, 1210, 495, 132, 0, 38232, 44604, 29302, 14560, 6006, 2002, 429, 0, 593859, 678696, 430200, 204540, 81900, 28392, 8008, 1430, 0, 10401712, 11701926, 7204821, 3289296, 1263780, 431256, 129948, 31824, 4862
Offset: 0

Views

Author

R. J. Mathar, Dec 06 2018

Keywords

Comments

If all subsets are allowed instead of just pairs (chords), we get A324173. The rightmost column is A000108 (see Riordan). - Gus Wiseman, Feb 27 2019

Examples

			From _Gus Wiseman_, Feb 27 2019: (Start)
Triangle begins:
  1
  0      1
  0      1      2
  0      4      6      5
  0     27     36     28     14
  0    248    310    225    120     42
  0   2830   3396   2332   1210    495    132
  0  38232  44604  29302  14560   6006   2002    429
  0 593859 678696 430200 204540  81900  28392   8008   1430
Row n = 3 counts the following chord diagrams (see link for pictures):
  {{1,3},{2,5},{4,6}}  {{1,2},{3,5},{4,6}}  {{1,2},{3,4},{5,6}}
  {{1,4},{2,5},{3,6}}  {{1,3},{2,4},{5,6}}  {{1,2},{3,6},{4,5}}
  {{1,4},{2,6},{3,5}}  {{1,3},{2,6},{4,5}}  {{1,4},{2,3},{5,6}}
  {{1,5},{2,4},{3,6}}  {{1,5},{2,3},{4,6}}  {{1,6},{2,3},{4,5}}
                       {{1,5},{2,6},{3,4}}  {{1,6},{2,5},{3,4}}
                       {{1,6},{2,4},{3,5}}
(End)
		

Crossrefs

Cf. A000699 (k = 1 column), A001147 (row sums), A000108 (diagonal), A002694 (subdiagonal k = n - 1).

Formula

The g.f. satisfies g(z,w) = 1+w*A000699(w*g^2), where A000699(z) is the g.f. of A000699.

Extensions

Offset changed to 0 by Gus Wiseman, Feb 27 2019

A324323 Regular triangle read by rows where T(n,k) is the number of topologically connected set partitions of {1,...,n} with k blocks, 0 <= k <= n.

Original entry on oeis.org

1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 5, 0, 0, 0, 0, 1, 16, 4, 0, 0, 0, 0, 1, 42, 42, 0, 0, 0, 0, 0, 1, 99, 258, 27, 0, 0, 0, 0, 0, 1, 219, 1222, 465, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Gus Wiseman, Feb 22 2019

Keywords

Comments

A set partition of {1,...,n} 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.

Examples

			Triangle begins:
    1
    0    1
    0    1    0
    0    1    0    0
    0    1    1    0    0
    0    1    5    0    0    0
    0    1   16    4    0    0    0
    0    1   42   42    0    0    0    0
    0    1   99  258   27    0    0    0    0
    0    1  219 1222  465    0    0    0    0    0
Row n = 6 counts the following set partitions:
  {{123456}}  {{1235}{46}}  {{13}{25}{46}}
              {{124}{356}}  {{14}{25}{36}}
              {{1245}{36}}  {{14}{26}{35}}
              {{1246}{35}}  {{15}{24}{36}}
              {{125}{346}}
              {{13}{2456}}
              {{134}{256}}
              {{1345}{26}}
              {{1346}{25}}
              {{135}{246}}
              {{1356}{24}}
              {{136}{245}}
              {{14}{2356}}
              {{145}{236}}
              {{146}{235}}
              {{15}{2346}}
		

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[#]]<=1&&Length[#]==k&]],{n,0,6},{k,0,n}]

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
    				

A326335 Number of set partitions of {1..n} whose nesting blocks are connected.

Original entry on oeis.org

1, 1, 1, 1, 2, 6, 21, 86, 394, 1974, 10696
Offset: 0

Views

Author

Gus Wiseman, Jun 27 2019

Keywords

Comments

Two blocks are nesting if they are of the form {...x,y...}, {...z,t...} where x < z < t < y or z < x < y < t. A set partition has its nesting blocks connected if the graph whose vertices are the blocks and whose edges are nesting pairs of blocks is connected.

Examples

			The a(0) = 1 through a(6) = 21 set partitions:
  {}  {1}  {12}  {123}  {1234}    {12345}    {123456}
                        {14}{23}  {125}{34}  {1236}{45}
                                  {134}{25}  {1245}{36}
                                  {14}{235}  {125}{346}
                                  {145}{23}  {1256}{34}
                                  {15}{234}  {126}{345}
                                             {134}{256}
                                             {1345}{26}
                                             {1346}{25}
                                             {136}{245}
                                             {14}{2356}
                                             {145}{236}
                                             {1456}{23}
                                             {146}{235}
                                             {15}{2346}
                                             {156}{234}
                                             {16}{2345}
                                             {15}{26}{34}
                                             {16}{23}{45}
                                             {16}{24}{35}
                                             {16}{25}{34}
		

Crossrefs

Simple graphs whose nesting blocks are connected are A326330.
Set partitions whose crossing blocks are connected are A099947.
Set partitions whose capturing blocks are connected are A326336.

Programs

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

A326336 Number of set partitions of {1..n} whose capturing blocks are connected.

Original entry on oeis.org

1, 1, 1, 1, 2, 7, 24, 100, 458, 2279, 12270
Offset: 0

Views

Author

Gus Wiseman, Jun 28 2019

Keywords

Comments

Two blocks are capturing if they are of the form {...x...y...}, {...z...t...} where x < z < t < y or z < x < y < t. A set partition has its capturing blocks connected if the graph whose vertices are the blocks and whose edges are capturing pairs of blocks is connected.

Examples

			The a(0) = 1 through a(6) = 24 set partitions:
  {}  {1}  {12}  {123}  {1234}    {12345}    {123456}
                        {14}{23}  {125}{34}  {1236}{45}
                                  {134}{25}  {1245}{36}
                                  {135}{24}  {1246}{35}
                                  {14}{235}  {125}{346}
                                  {145}{23}  {1256}{34}
                                  {15}{234}  {126}{345}
                                             {134}{256}
                                             {1345}{26}
                                             {1346}{25}
                                             {135}{246}
                                             {1356}{24}
                                             {136}{245}
                                             {14}{2356}
                                             {145}{236}
                                             {1456}{23}
                                             {146}{235}
                                             {15}{2346}
                                             {156}{234}
                                             {16}{2345}
                                             {15}{26}{34}
                                             {16}{23}{45}
                                             {16}{24}{35}
                                             {16}{25}{34}
		

Crossrefs

Simple graphs whose capturing blocks are connected are A326330.
Set partitions whose crossing blocks are connected are A099947.
Set partitions whose nesting blocks are connected are A326335.

Programs

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

A268815 Number of purely crossing + partitions of [n].

Original entry on oeis.org

1, 1, 0, 0, 1, 1, 5, 19, 76, 360, 1792, 9634, 55286, 336396, 2162554, 14629720, 103818489, 770678553, 5969822993, 48148947503, 403545713463, 3508356996105, 31587389832791, 294087418038113, 2827471212909189, 28037001032306431, 286398141349873925, 3010540174760962975
Offset: 0

Views

Author

Michel Marcus, Feb 14 2016

Keywords

Comments

For the definition of these special purely crossing partitions refer to Dykema link (see PC+(n) Definition 2.1 and Table 2).
From Gus Wiseman, Feb 23 2019: (Start)
a(n) is the number of topologically connected (A099947) set partitions of {1,...,n} with no successive elements in the same block. For example, the a(4) = 1 through a(7) = 19 set partitions are:
{{13}{24}} {{135}{24}} {{135}{246}} {{1357}{246}}
{{13}{25}{46}} {{13}{246}{57}}
{{14}{25}{36}} {{13}{257}{46}}
{{14}{26}{35}} {{135}{26}{47}}
{{15}{24}{36}} {{135}{27}{46}}
{{136}{24}{57}}
{{136}{25}{47}}
{{137}{25}{46}}
{{14}{257}{36}}
{{14}{26}{357}}
{{146}{25}{37}}
{{146}{27}{35}}
{{147}{25}{36}}
{{147}{26}{35}}
{{15}{246}{37}}
{{15}{247}{36}}
{{157}{24}{36}}
{{16}{24}{357}}
{{16}{247}{35}}
(End)

Examples

			G.f.: A(x) = 1 + x + x^4 + x^5 + 5*x^6 + 19*x^7 + 76*x^8 + 360*x^9 + 1792*x^10 +...
		

Crossrefs

Programs

  • Mathematica
    n = 30; F = x*Sum[BellB[k] x^k, {k, 0, n}] + O[x]^n; B = ComposeSeries[1/( InverseSeries[F, w] /w)-1, x/(1+x) + O[x]^n]; CoefficientList[B, x] // Rest (* Jean-François Alcover, Feb 16 2016, adapted from K. J. Dykema's code *)
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    intvQ[set_]:=Or[set=={},Sort[set]==Range[Min@@set,Max@@set]];
    Table[Length[Select[sps[Range[n]],And[!MatchQ[#,{_,{_,x_,y_,_},_}/;x+1==y],#=={}||And@@Not/@intvQ/@Union@@@Subsets[#,{1,Length[#]-1}]]&]],{n,0,10}] (* Gus Wiseman, Feb 23 2019 *)
  • PARI
    lista(nn) = {c = x/serreverse(x*serlaplace(exp(exp(x+x*O(x^nn)) -1))); b = subst(c, x, x/(1+x) + O(x^nn)); Vec(b);}
    
  • PARI
    {a(n) = my(A=1+x); for(i=1, n, A = sum(m=0, n, x^m/prod(k=1, m, (1+x)*A - k*x +x*O(x^n)) )); polcoeff(A, n)}
    for(n=0,25,print1(a(n),", ")) \\ Paul D. Hanna, Mar 07 2016
    
  • PARI
    {Stirling2(n, k) = n!*polcoeff(((exp(x+x*O(x^n)) - 1)^k)/k!, n)}
    {Bell(n) = sum(k=0,n, Stirling2(n, k) )}
    {a(n) = my(A=1+x); for(i=1, n, A = sum(m=0, n, Bell(m)*x^m/((1+x)*A +x*O(x^n))^m) ); polcoeff(A, n)}
    for(n=0,25,print1(a(n),", ")) \\ Paul D. Hanna, Mar 07 2016

Formula

G.f.: G(x) satisfies C(x) = G(x/1-x) where C(x) is the g.f. of A099947 (see B(x) in Dykema link p. 7).
From Paul D. Hanna, Mar 07 2016: (Start)
O.g.f. A(x) satisfies
(1) A(x) = Sum_{n>=0} A000110(n)*x^n/((1+x)^n*A(x)^n), where A000110 are the Bell numbers.
(2) A(x) = Sum_{n>=0} x^n / Product_{k=1..n} ((1+x)*A(x) - k*x).
(3) A(x) = 1/(1 - x/((1+x)*A(x) - 1*x/(1 - x/((1+x)*A(x) - 2*x/(1 - x/((1+x)*A(x) - 3*x/(1 - x/((1+x)*A(x) - 4*x/(1 - x/((1+x)*A(x) - ... )))))))))), a continued fraction. (End)

A324326 Number of crossing multiset partitions of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 10, 0, 0, 0, 10, 0, 0, 0, 7, 0, 0, 0, 0, 0, 0, 0, 31, 0, 0, 0, 0, 0, 36, 0, 14, 0, 0, 0, 25, 0, 0, 0, 71, 0, 0, 0, 0, 0, 0, 0, 103, 0, 0, 0, 0, 0, 0, 0, 75
Offset: 1

Views

Author

Gus Wiseman, Feb 22 2019

Keywords

Comments

This multiset (row n of A305936) is generally not the same as the multiset of prime indices of n. For example, the prime indices of 12 are {1,1,2}, while a multiset whose multiplicities are {1,1,2} is {1,1,2,3}.
A multiset partition is crossing if it contains two blocks of the form {{...x...y...},{...z...t...}} with x < z < y < t or z < x < t < y.

Examples

			The a(36) = 10 crossing multiset partitions of {1,1,2,2,3,4}:
  {{1,3},{1,2,2,4}}
  {{2,4},{1,1,2,3}}
  {{1,1,3},{2,2,4}}
  {{1,2,3},{1,2,4}}
  {{1},{1,3},{2,2,4}}
  {{1},{2,4},{1,2,3}}
  {{2},{1,3},{1,2,4}}
  {{2},{1,1,3},{2,4}}
  {{1,2},{1,3},{2,4}}
  {{1},{2},{1,3},{2,4}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    croXQ[stn_]:=MatchQ[stn,{_,{_,x_,_,y_,_},_,{_,z_,_,t_,_},_}/;x
    				

Formula

a(n) + A324325(n) = A318284(n).

A125107 Subtract compositions (A011782) from Catalan numbers (A000108).

Original entry on oeis.org

0, 0, 0, 1, 6, 26, 100, 365, 1302, 4606, 16284, 57762, 205964, 738804, 2666248, 9678461, 35324902, 129579254, 477507628, 1767001046, 6563596132, 24465218444, 91480466488, 343055419346, 1289895758716, 4861929624236, 18367319517720, 69533483807140, 263747817532632
Offset: 0

Views

Author

Alford Arnold, Dec 15 2006

Keywords

Comments

Apparently the number of Dyck n-paths with more than half of the path lying between the first and last peaks. - David Scambler, Sep 14 2012
From Peter Luschny, Nov 28 2024: (Start)
A Touchard walk T(n) of length n is, as defined by Dershowitz, "a sequence of n steps, each of which is one of N/S/E/W, such that at each point along the way the number of N-steps that have been taken is never less than the number of S-steps, and are in the end equal."
There are Sum_{k=0..n} binomial(n, k) Touchard walks that have no N/S-steps at all and since by Touchard's identity T(n) = Catalan(n+1), it follows that a(n) = T(n-1) - Sum_{k=0..n-1} binomial(n-1, k) = Catalan(n) - 2^(n-1) for n >= 1. Thus a(n+1) is the number of Touchard walks of length n that have at least one N-step. (End)

Examples

			A000108 begins 1 1 2 5 14 42 132 429 ...
A011782 begins 1 1 2 4  8 16  32  64 ...
so we get .... 0 0 0 1  6 26 100 365 ...
.
The 26 Touchard walks of length 4 that have at least one N-step are:
   NNSS, NSNS, NSEE, NSEW, NSWE, NSWW, NESE, NESW, NWSE,
   NWSW, NEES, NEWS, NWES, NWWS, ENSE, ENSW, WNSE, WNSW,
   ENES, ENWS, WNES, WNWS, EENS, EWNS, WENS, WWNS.
		

Crossrefs

Cf. A000079, A000108, A000110, A011782, A016098, A097805, A091894 (Touchard distribution), A377659 (similar with Motzkin).

Programs

  • Maple
    # From Peter Luschny, Nov 28 2024: (Start)
    a := n -> ifelse(n = 0, 0, binomial(2*n, n)/(n+1) - 2^(n-1)): seq(a(n), n = 0..28);
    # Series expansion:
    gf := (1 - sqrt(1 - 4*x)) / (2*x) - (1 - x) / (1 - 2*x): ser := series(gf, x, 30): seq(coeff(ser, x, n), n = 0..28);
    # Evaluating polynomials:
    p := (n, x) -> ifelse(n = 0, 0, 2^(n-1)*(hypergeom([1 - n/2, 1/2 - n/2], [2], x) - 1)): seq(subs(x = 1, expand(simplify(p(n, x)))), n = 0..28);  # (End)
  • Mathematica
    Table[CatalanNumber[n] - If[n==0, 1, 2^(n-1)], {n, 0, 30}] (* David Scambler, Sep 14 2012 *)
  • Python
    # Generates the walks (for illustration only).
    C = str.count
    def aGen(n: int) -> Generator[str, Any, list[str]]:
        a = [""]
        if n <= 0: return a
        for w in a:
            if len(w) == n - 1:
                if C(w, "N") > 0 and C(w, "N") == C(w, "S"):
                    yield w
            else:
                for j in "NSEW":
                    U = w + j
                    if C(U, "N") >= C(U, "S"):
                        a += [U]
        return a
    for n in range(6): print([w for w in aGen(n)])  # Peter Luschny, Dec 03 2024

Formula

a(n) = A000108(n) - A011782(n).
(n+1)*a(n) + 2*(1-4*n)*a(n-1) + 4*(5*n-7)*a(n-2) + 8*(5-2*n)*a(n-3) = 0. - R. J. Mathar, Aug 10 2013
From Peter Luschny, Nov 28 2024: (Start)
a(n) = [x^n] (1 - sqrt(1 - 4*x)) / (2*x) - (1 - x) / (1 - 2*x).
a(n) = n! * [x^n] (exp(2*x)*(BesselI_{0}(2*x) - BesselI_{1}(2*x) - 1/2) - 1/2).
a(n) = p(n, 1) for n >= 1, where p(n, x) = 2^(n-1)*(hypergeom([1-n/2, (1-n)/2], [2], x) - 1).
a(n) = Sum_{k=0..n} (A091894(n, k) - A097805(n, n-k)). (End)

Extensions

More terms from David Scambler, Sep 14 2012
Previous Showing 31-40 of 54 results. Next