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-5 of 5 results.

A003436 Number of inequivalent labeled Hamiltonian circuits on n-octahedron. Interlacing chords joining 2n points on circle.

Original entry on oeis.org

1, 0, 1, 4, 31, 293, 3326, 44189, 673471, 11588884, 222304897, 4704612119, 108897613826, 2737023412199, 74236203425281, 2161288643251828, 67228358271588991, 2225173863019549229, 78087247031912850686, 2896042595237791161749, 113184512236563589997407
Offset: 0

Views

Author

Keywords

Comments

Also called the relaxed ménage problem (cf. A000179).
a(n) can be seen as a subset of the unordered pairings of the first 2n integers (A001147) with forbidden pairs (1,2n) and (i,i+1) for all i in [1,2n-1] (all adjacent integers modulo 2n). The linear version of this constraint is A000806. - Olivier Gérard, Feb 08 2011
Number of perfect matchings in the complement of C_{2n} where C_{2n} is the cycle graph on 2n vertices. - Andrew Howroyd, Mar 15 2016
Also the number of 2-uniform set partitions of {1...2n} containing no two cyclically successive vertices in the same block. - Gus Wiseman, Feb 27 2019

References

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

Crossrefs

Cf. A003435, A129348. A003437 gives unlabeled case.
First differences of A000806.
Column k=2 of A324428.

Programs

  • Maple
    A003436 := proc(n) local k;
          if n = 0 then 1
        elif n = 1 then 0
        else add( (-1)^k*binomial(n,k)*2*n/(2*n-k)*2^k*(2*n-k)!/2^n/n!,k=0..n) ;
        end if;
    end proc: # R. J. Mathar, Dec 11 2013
    A003436 := n-> `if`(n<2, 1-n, (-1)^n*2*hypergeom([n, -n], [], 1/2)):
    seq(simplify(A003436(n)), n=0..18); # Peter Luschny, Nov 10 2016
  • Mathematica
    a[n_] := (2*n-1)!! * Hypergeometric1F1[-n, 1-2*n, -2]; a[1] = 0;
    Table[a[n], {n, 0, 19}] (* Jean-François Alcover, Apr 05 2013 *)
    twounifll[{}]:={{}};twounifll[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@twounifll[Complement[set,s]]]/@Table[{i,j},{j,If[i==1,Select[set,2<#i+1&]]}];
    Table[Length[twounifll[Range[n]]],{n,0,14,2}] (* Gus Wiseman, Feb 27 2019 *)

Formula

a(n) = A003435(n)/(n!*2^n).
a(n) = 2*n*a(n-1)-2*(n-3)*a(n-2)-a(n-3) for n>4. [Corrected by Vasu Tewari, Apr 11 2010, and by R. J. Mathar, Oct 02 2013]
G.f.: x + ((1-x)/(1+x)) * Sum_{n>=0} A001147(n)*(x/(1+x)^2)^n. - Vladeta Jovovic, Jun 27 2007
a(n) ~ 2^(n+1/2)*n^n/exp(n+1). - Vaclav Kotesovec, Aug 13 2013
a(n) = (-1)^n*2*hypergeom([n, -n], [], 1/2) for n >= 2. - Peter Luschny, Nov 10 2016

Extensions

a(0)=1 prepended by Gus Wiseman, Feb 27 2019

A324428 Number T(n,k) of labeled cyclic chord diagrams with n chords such that every chord has length at least k; triangle T(n,k), n>=1, 1<=k<=n, read by rows.

Original entry on oeis.org

1, 3, 1, 15, 4, 1, 105, 31, 7, 1, 945, 293, 68, 11, 1, 10395, 3326, 837, 159, 18, 1, 135135, 44189, 11863, 2488, 381, 29, 1, 2027025, 673471, 189503, 43169, 7601, 879, 47, 1, 34459425, 11588884, 3377341, 822113, 160784, 23559, 2049, 76, 1, 654729075, 222304897, 66564396, 17066007, 3621067, 607897, 72989, 4788, 123, 1
Offset: 1

Views

Author

Alois P. Heinz, Feb 27 2019

Keywords

Comments

T(n,k) is defined for all n,k >= 0. The triangle contains only the terms with 1 <= k <= n. T(n,0) = A001147(n), T(0,k) = 1, T(n,k) = 0 for k > n > 0.

Examples

			Triangle T(n,k) begins:
        1;
        3,      1;
       15,      4,      1;
      105,     31,      7,     1;
      945,    293,     68,    11,    1;
    10395,   3326,    837,   159,   18,   1;
   135135,  44189,  11863,  2488,  381,  29,  1;
  2027025, 673471, 189503, 43169, 7601, 879, 47, 1;
  ...
		

Crossrefs

T(n,n-1) gives A000204.

Programs

  • Maple
    b:= proc(n, f, m, l, j) option remember; (k-> `if`(n `if`(n=0 or k<2, doublefactorial(2*n-1),
                  b(2*n-k+1, [1$k-1], 0, [0$k-1], k-1)):
    seq(seq(T(n, k), k=1..n), n=1..10);
  • Mathematica
    b[n_, f_List, m_, l_List, j_] := b[n, f, m, l, j] = Function[k, If[n < Total[f] + m + Total[l], 0, If[n == 0, 1, Sum[If[f[[i]] == 0, 0, b[n - 1, ReplacePart[f, i -> 0], m + l[[1]], Append[ReplacePart[l, 1 -> Nothing], 0], Max[0, j - 1]]],{i, Max[1, j + 1], Min[k, n - 1]}] + If[m == 0, 0, m*b[n - 1, f, m - 1 + l[[1]], Append[ReplacePart[l, 1 -> Nothing], 0], Max[0, j - 1]]] + b[n - 1, f, m + l[[1]], Append[ReplacePart[l, 1 -> Nothing], 1], Max[0, j - 1]]]]][Length[l]];
    T[n_, k_] := If[n == 0 || k < 2, 2^(n-1) Pochhammer[3/2, n-1], b[2n-k+1, Table[1, {k-1}], 0, Table[0, {k-1}], k-1]];
    Table[T[n, k], {n, 1, 10}, {k, 1, n}] // Flatten (* Jean-François Alcover, Apr 27 2020, after Alois P. Heinz *)

Formula

T(n,k) = Sum_{j=k..n} A324429(n,j).

A190823 Number of permutations of 2 copies of 1..n introduced in order 1..n with no element equal to another within a distance of 2.

Original entry on oeis.org

1, 0, 0, 1, 10, 99, 1146, 15422, 237135, 4106680, 79154927, 1681383864, 39034539488, 983466451011, 26728184505750, 779476074425297, 24281301468714902, 804688068731837874, 28269541494090294129, 1049450257149017422000, 41050171013933837206545
Offset: 0

Views

Author

R. H. Hardin, May 21 2011

Keywords

Comments

From Gus Wiseman, Feb 27 2019: (Start)
Also the number of 2-uniform set partitions of {1..2n} such that no block has its two vertices differing by less than 3. For example, the a(4) = 10 set partitions are:
{{1,4}, {2,6}, {3,7}, {5,8}}
{{1,4}, {2,7}, {3,6}, {5,8}}
{{1,5}, {2,6}, {3,7}, {4,8}}
{{1,5}, {2,6}, {3,8}, {4,7}}
{{1,5}, {2,7}, {3,6}, {4,8}}
{{1,5}, {2,8}, {3,6}, {4,7}}
{{1,6}, {2,5}, {3,7}, {4,8}}
{{1,6}, {2,5}, {3,8}, {4,7}}
{{1,7}, {2,5}, {3,6}, {4,8}}
{{1,8}, {2,5}, {3,6}, {4,7}}
(End)

Examples

			All solutions for n=4 (read downwards):
  1    1    1    1    1    1    1    1    1    1
  2    2    2    2    2    2    2    2    2    2
  3    3    3    3    3    3    3    3    3    3
  4    4    4    4    1    4    4    1    4    4
  1    1    2    1    4    2    1    4    2    2
  3    3    1    2    2    3    2    3    1    3
  2    4    4    4    3    4    3    2    3    1
  4    2    3    3    4    1    4    4    4    4
		

Crossrefs

Distance of 1 instead of 2 gives |A000806|.
Column k=3 of A293157.
Cf. A000699, A001147 (2-uniform set partitions), A003436, A005493, A011968, A170941, A278990 (distance 2+ version), A306386 (cyclical version).

Programs

  • Magma
    I:=[1,0,0,1,10,99]; [n le 5 select I[n] else 2*n*Self(n-1) -2*(3*n-8)*Self(n-2) +2*(3*n-11)*Self(n-3) -2*(n-5)*Self(n-4) -Self(n-5): n in [1..40]]; // G. C. Greubel, Dec 03 2023
    
  • Mathematica
    a[0]=1; a[1]=0; a[2]=0; a[3]=1; a[4]=10; a[5]=99; a[n_] := a[n] = (2*n+2) a[n-1] - (6*n-10) a[n-2] + (6*n-16) a[n-3] - (2*n-8) a[n-4] - a[n-5]; Array[a, 20, 0] (* based on Sullivan's formula, Giovanni Resta, Mar 20 2017 *)
    dtui[{}]:={{}};dtui[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@dtui[Complement[set,s]]]/@Table[{i,j},{j,Select[set,#>i+2&]}];
    Table[Length[dtui[Range[n]]],{n,0,12,2}] (* Gus Wiseman, Feb 27 2019 *)
  • SageMath
    @CachedFunction
    def a(n): # a = A190823
        if (n<6): return (1,0,0,1,10,99)[n]
        else: return 2*(n+1)*a(n-1) - 2*(3*n-5)*a(n-2) + 2*(3*n-8)*a(n-3) - 2*(n-4)*a(n-4) - a(n-5)
    [a(n) for n in range(41)] # G. C. Greubel, Dec 03 2023

Formula

a(n) = 2*(n+1)*a(n-1) - 2*(3*n-5)*a(n-2) + 2*(3*n-8)*a(n-3) - 2*(n-4)*a(n-4) - a(n-5) (proved). - Everett Sullivan, Mar 16 2017
a(n) ~ 2^(n+1/2) * n^n / exp(n+2), based on Sullivan's formula. - Vaclav Kotesovec, Mar 21 2017

Extensions

a(16)-a(20) (using Everett Sullivan's formula) from Giovanni Resta, Mar 20 2017
a(0)=1 prepended by Alois P. Heinz, Oct 17 2017

A324446 Number of labeled cyclic chord diagrams with n chords such that the minimal chord length equals two.

Original entry on oeis.org

1, 3, 24, 225, 2489, 32326, 483968, 8211543, 155740501, 3265307342, 74995101843, 1872508994356, 50500982610620, 1463062187672336, 45314261742435296, 1494164679669072424, 52257665502536426741, 1932255827699763531474, 75312621088768346098203
Offset: 2

Views

Author

Alois P. Heinz, Feb 28 2019

Keywords

Crossrefs

Column k=2 of A324429.

Formula

a(n) = A003436(n) - A306386(n).

A324447 Number of labeled cyclic chord diagrams with n chords such that the minimal chord length equals three.

Original entry on oeis.org

1, 6, 57, 678, 9375, 146334, 2555228, 49498389, 1054882116, 24545326598, 619415271672, 16854164945997, 491956017936245, 15335174679459134, 508475723159116576, 17870522292434937714, 663618269365250877446, 25964529622925026990819, 1067602511902741954478695
Offset: 3

Views

Author

Alois P. Heinz, Feb 28 2019

Keywords

Crossrefs

Column k=3 of A324429.

Formula

a(n) = A306386(n) - A324430(n).
Showing 1-5 of 5 results.