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.

User: Martin Rubey

Martin Rubey's wiki page.

Martin Rubey has authored 9 sequences.

A362907 Number of graphs on n unlabeled nodes with treewidth 3.

Original entry on oeis.org

0, 0, 0, 0, 1, 6, 53, 501, 5889, 81786
Offset: 0

Author

Martin Rubey, May 09 2023

Keywords

Examples

			There is a single graph of treewidth 3 on 4 vertices, which is the complete graph.
		

Crossrefs

Column k=3 of A263294.

Programs

  • SageMath
    [sum(1 for g in graphs(n) if g.treewidth() == 3) for n in range(9)]

A362908 Number of graphs on n unlabeled nodes with treewidth 2.

Original entry on oeis.org

0, 0, 0, 1, 4, 17, 72, 323, 1639, 9203
Offset: 0

Author

Martin Rubey, May 09 2023

Keywords

Examples

			There is a single graph of treewidth 2 on 3 vertices, which is the complete graph.
		

Crossrefs

Column 2 of A263294.

Programs

  • SageMath
    [sum(1 for g in graphs(n) if g.treewidth() == 2) for n in range(9)]

A361920 Number of unlabeled ranked posets with n elements.

Original entry on oeis.org

1, 1, 2, 5, 16, 61, 280, 1501, 9394, 68647, 591570, 6108298, 77162708, 1219779207, 24648006828, 647865966973, 22437052221282, 1032905858402302, 63591727342096158, 5258562027225785955, 586001891321599337103, 88241281449605821921186, 17996565026907866304071630
Offset: 0

Author

Martin Rubey, Mar 29 2023

Keywords

Comments

A partially ordered set is ranked if there is a function from the poset elements to the integers such that the function value of a covering element is precisely one larger than the function value of the covered element. This is called graded by some authors.

Examples

			For n=5, A000112(n) - a(n) = 63 - 61 = 2 because we have 2 posets with 5 elements that are not ranked: a<b<c<d  a<e<d  and  a<c<e  a<d  b<d  b<e where < means "is covered by". - _Geoffrey Critzer_, Oct 29 2023
		

Crossrefs

Row sums of A361953.

Programs

  • PARI
    \\ See PARI link in A361953 for program code.
    A361920seq(20) \\ Andrew Howroyd, Apr 01 2023
  • Sage
    sum(1 for P in posets(n) if P.is_ranked())
    

Extensions

Terms a(8) and beyond from Andrew Howroyd, Mar 31 2023

A361921 The number of unlabeled bounded Eulerian posets with n elements.

Original entry on oeis.org

0, 1, 1, 0, 1, 0, 1, 0, 2, 0, 5, 0, 11
Offset: 0

Author

Martin Rubey, Mar 30 2023

Keywords

Comments

A graded partially ordered set is Eulerian if every nontrivial interval has the same number of elements of even rank as of odd rank. It is bounded if it has a unique maximal and a unique minimal element.

Programs

  • Sage
    sum(1 for P in posets(n-2) if (Q := P.with_bounds()).is_graded() and Q.is_eulerian())

A361912 The number of unlabeled graded posets with n elements.

Original entry on oeis.org

1, 1, 2, 4, 10, 28, 93, 354, 1621, 9110, 64801, 595976, 7204091, 115561423, 2473540433, 70853213144, 2720354016419, 140170631441858, 9702605436760235, 903309202327818566, 113234129823368903523, 19137461395401601912043, 4366007821745938984134203
Offset: 0

Author

Martin Rubey, Mar 29 2023

Keywords

Comments

A partially ordered set is graded if all maximal chains have the same length. This is called tiered by some authors.

Crossrefs

Row sums of A361957.
Cf. A000112, A223911 (labeled), A001833, A361920, A361959 (connected).

Programs

  • PARI
    \\ See PARI link in A361957 for program code.
    A361912seq(20) \\ Andrew Howroyd, Apr 03 2023
  • Sage
    sum(1 for P in posets(n) if P.is_graded())
    

Extensions

Terms a(8) and beyond from Andrew Howroyd, Mar 30 2023

A289971 Number of permutations of [n] determined by their antidiagonal sums.

Original entry on oeis.org

1, 1, 2, 4, 9, 20, 49, 114, 277, 665, 1608, 3875
Offset: 0

Author

Martin Rubey, Jul 16 2017

Keywords

Crossrefs

Programs

  • Mathematica
    xray[perm_List] := Module[{P, n = Length[perm]}, P[, ] = 0; Thread[perm -> Range[n]] /. Rule[i_, j_] :> Set[P[i, j], 1]; Table[Sum[P[i - j + 1, j], {j, Max[1, i - n + 1], Min[i, n]}], {i, 1, 2n - 1}]];
    a[n_] := xray /@ Permutations[Range[n]] // Tally // Count[#, {_List, 1}]&;
    Do[Print[n, " ", a[n]], {n, 0, 10}] (* Jean-François Alcover, Feb 28 2020 *)
  • Sage
    def X_ray(pi):
        P = Permutation(pi).to_matrix()
        n = P.nrows()
        return tuple(sum(P[k-1-j][j] for j in range(max(0, k-n), min(k,n)))
                     for k in range(1,2*n))
    @cached_function
    def X_rays(n):
        return sorted(X_ray(pi) for pi in Permutations(n))
    def statistic(pi): return X_rays(pi.size()).count(X_ray(pi))
    [[statistic(pi) for pi in Permutations(n)].count(1) for n in range(7)]

Extensions

a(8)-a(11) from Alois P. Heinz, Jul 24 2017

A269890 Number of conjugacy classes of subgroups of the hyperoctahedral group.

Original entry on oeis.org

2, 8, 33, 193, 953, 7440, 55200, 627187, 7510549
Offset: 1

Author

Martin Rubey, Mar 07 2016

Keywords

Crossrefs

Cf. A000638.

Programs

  • GAP
    Length(ConjugacyClassesSubgroups(CoxeterGroup( "B", n )));
    
  • Magma
    [#Subgroups(CoxeterGroup("B" cat IntegerToString(n))) : n in [1..9]];  // Robin Visser, Aug 09 2023

Extensions

a(6)-a(9) from Robin Visser, Aug 09 2023

A181519 Number of torsion pairs in the cluster category of type A_n up to Auslander-Reiten translation.

Original entry on oeis.org

1, 3, 5, 19, 62, 301, 1413, 7304, 38294, 208052, 1149018, 6466761, 36899604, 213245389, 1245624985, 7345962126, 43688266206, 261791220038, 1579363550250, 9586582997562, 58513327318992, 358957495385684, 2212294939905234
Offset: 3

Author

Martin Rubey, Oct 26 2010

Keywords

Comments

a(n) is also the number of Ptolemy diagrams on n vertices up to rotation.
a(n) is the sum over all polygon dissections up to rotation, where each region of size at least four has weight two.

Examples

			For n=5 there are 5 Ptolemy diagrams up to rotation: the pentagon with no diagonal, the pentagon with one diagonal, the pentagon with two noncrossing diagonals, the pentagon with three diagonals and the pentagon with all five diagonals.
		

Crossrefs

Cf. A181517.

Programs

  • PARI
    seq(n)={my(p=serreverse(x - x^2*(1 + x)/(1 - x) + O(x*x^n)), P(d)=subst(p + O(x^(n\d+1)), x, x^d)); Vec(2*sum(d=1, n, eulerphi(d)/d*log(1/(1-P(d)))) - (P(1)^3 + 2*P(3))/3 - (3*P(1)^2+P(2))/2 - (2 - x)*P(1) - x^2)} \\ Andrew Howroyd, May 09 2023

Formula

G.f.: (2*Sum_{d>=1} phi(d)*log(1/(1-P(y^d)))/d ) - (1/3)*(P(y)^3+2*P(y^3)) - (1/2)*(3*P(y)^2+P(y^2)) - 2*P(y) + y*P(y) - y^2 where y*P(y) - y^2 is the g.f. of A181517. [corrected by Andrew Howroyd, May 09 2023]

A181517 Number of torsion pairs in the cluster category of type A_n.

Original entry on oeis.org

1, 4, 17, 82, 422, 2274, 12665, 72326, 421214, 2492112, 14937210, 90508256, 553492552, 3411758334, 21175624713, 132226234854, 830077057878, 5235817447752, 33166634502334, 210904780742860, 1345806528336772, 8614979593487972, 55307373497626442, 356012579697723084
Offset: 3

Author

Martin Rubey, Oct 26 2010

Keywords

Comments

a(n) is also the number of Ptolemy diagrams on n vertices with distinguished base edge.
a(n) is the sum over all polygon dissections in a polygon with distinguished base edge, where each region of size at least four has weight two.

Examples

			For n=4 there are 4 Ptolemy diagrams: the square with no diagonal, two diagrams with one diagonal, and the square with both diagonals .
		

Crossrefs

Cf. A181519.

Programs

  • PARI
    a(n) = n-=3; sum(i=0, floor((n+1)/2), 2^i*binomial(n+1+i,i)*binomial(2*n+2,n+1-2*i))/(n+2); \\ Michel Marcus, Jan 14 2012; corrected Jun 13 2022
    
  • PARI
    seq(n) = Vec(x*(serreverse(x - x^2*(1 + x)/(1 - x) + O(x^(n+2))) - x)) \\ Andrew Howroyd, May 09 2023

Formula

G.f.: y*P(y) - y^2 where P(y) satisfies P(y) = y + P(y)^2*(1+P(y))/(1-P(y)).