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

A274634 a(n) = n!*A003436(n).

Original entry on oeis.org

0, 2, 24, 744, 35160, 2394720, 222712560, 27154350720, 4205374225920, 806700010233600, 187793061031699200, 52162131258836121600, 17043501717850146739200, 6471785359791584459827200, 2826261760629911644744704000, 1406604814631643298586923008000
Offset: 1

Views

Author

N. J. A. Sloane, Jul 08 2016

Keywords

Crossrefs

Cf. A003436.

Programs

  • Mathematica
    a[n_]:= n! (2 n - 1)!! Hypergeometric1F1[-n, 1 - 2 n, - 2]; a[1] = 0; Table[a[n], {n, 1, 20}] (* Vincenzo Librandi, Aug 18 2016 *)

Formula

a(n) ~ 2*sqrt(Pi*n)*2^n*n^(2*n)/exp(2*n+1). - Ilya Gutkovskiy, Aug 18 2016

A000858 Duplicate of A003436.

Original entry on oeis.org

0, 1, 4, 31, 293, 3326, 44189, 673471, 11588884, 222304897, 4704612119
Offset: 1

Views

Author

Keywords

A000806 Bessel polynomial y_n(-1).

Original entry on oeis.org

1, 0, 1, -5, 36, -329, 3655, -47844, 721315, -12310199, 234615096, -4939227215, 113836841041, -2850860253240, 77087063678521, -2238375706930349, 69466733978519340, -2294640596998068569, 80381887628910919255, -2976424482866702081004, 116160936719430292078411
Offset: 0

Views

Author

Keywords

Comments

a(n) can be seen as a subset of the unordered pairings of the first 2n integers (A001147) with forbidden pairs (i,i+1) for all i in [1,2n-1] (all adjacent integers). The circular version of this constraint is A003436. - Olivier Gérard, Feb 08 2011
|a(n)| is the number of perfect matchings in the complement of P_{2n} where P_{2n} is the path graph on 2n vertices. - Andrew Howroyd, Mar 15 2016
The unsigned version of these numbers now has its own entry: see A278990. - N. J. A. Sloane, Dec 07 2016

Examples

			For n=3, the a(3) = 5 solutions are (14) (25) (36), (14) (26) (35), (15) (24) (36), (16) (24) (35), (13) (25) (46) excluding 10 other possible pairings.
G.f. = 1 + x^2 - 5*x^3 + 36*x^4 - 329*x^5 + 3655*x^6 - 47844*x^7 + ...
		

References

  • G. Kreweras and Y. Poupard, Sur les partitions en paires d'un ensemble fini totalement ordonné, Publications de l'Institut de Statistique de l'Université de Paris, 23 (1978), 57-74.
  • J. Riordan, Combinatorial Identities, Wiley, 1968, p. 77.
  • 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

Polynomial coefficients are in A001498. Cf. A003436.

Programs

  • Magma
    I:=[0,1]; [1] cat [n le 2 select I[n] else (1-2*n)*Self(n-1)+Self(n-2): n in [1..30]]; // Vincenzo Librandi, Apr 19 2015
  • Maple
    A000806 := proc(n) option remember; if n<=1 then 1-n else (1-2*n)*procname(n-1)+procname(n-2); fi; end proc;
    a := n -> hypergeom([n+1,-n],[],1/2): seq(simplify(a(n)),n=0..20); # Peter Luschny, Nov 10 2016
  • Mathematica
    a[n_] := a[n] = (-2n+1)*a[n-1] + a[n-2]; a[0] = 1; a[1] = 0; Table[a[n], {n, 0, 19}] (* Jean-François Alcover, Nov 29 2011, after T. D. Noe *)
    Table[Sum[Binomial[n, i]*(2*n-i)!/2^(n-i)*(-1)^(n-i)/n!, {i, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Aug 07 2013 *)
    a[ n_] := With[ {m = If[ n<0, -n-1, n]}, (-1)^m (2 m - 1)!! Hypergeometric1F1[ -m, -2 m, -2] ]; (* Michael Somos, Jan 27 2014 *)
    a[ n_] := With[ {m = If[ n<0, -n-1, n]}, Sum[ (-1)^(m - i) (2 m - i)! / (2^(m - i) i! (m - i)!), {i, 0, m}] ]; (* Michael Somos, Jan 27 2014 *)
    a[ n_] := With[ {m = If[ n<0, -n-1, n]}, If[ m<1, 1, (-1)^m Numerator @ FromContinuedFraction[ Table[ (-1)^Quotient[k, 2] If[ OddQ[k], k, 1], {k, 2 m}] ] ] ]; (* Michael Somos, Jan 27 2014 *)
    Table[(-1)^n (2 n - 1)!! Hypergeometric1F1[-n, -2 n, -2], {n, 0, 20}] (* Eric W. Weisstein, Nov 14 2018 *)
  • PARI
    {a(n) = if( n<0, n = -n-1); sum(k=0, n, (2*n-k)! / (k! * (n-k)!) * (-1/2)^(n-k) )}; /* Michael Somos, Apr 02 2007 */
    
  • PARI
    {a(n) = local(A); if( n<0, n = -n-1); A = sqrt(1 + 2*x + x * O(x^n)); n! * polcoeff( exp(A-1) / A, n)}; /* Michael Somos, Apr 02 2007 */
    
  • PARI
    {a(n) = local(A); if( n<0, n = -n-1); n+=2; -(-1)^n * n! * polcoeff( serreverse( sum(k=1, n, k^(k-2)* x^k / k!, x * O(x^n))), n)}; /* Michael Somos, Apr 02 2007 */
    
  • PARI
    {a(n) = if( n<0, n=-n-1); contfracpnqn( vector( 2*n, k, (-1)^(k\2) * if( k%2, k, 1))) [1,1] }; /* Michael Somos, Jan 27 2014 */
    

Formula

E.g.f.: exp(sqrt(1 + 2*x) - 1) / sqrt(1 + 2*x). - Michael Somos, Feb 16 2002
D-finite with recurrence a(n) = (-2*n+1)*a(n-1) + a(n-2). - T. D. Noe, Oct 26 2006
If y = x + Sum_{k>1} A000272(k) * x^k/k!, then y = x + Sum{k>1} a(k-2) * (-y)^k/k!. - Michael Somos, Sep 07 2005
a(-1-n) = a(n). - Michael Somos, Apr 02 2007
a(n) = Sum_{m=0..n} A001498(n,m)*(-1)^m, n>=0 (alternating row sums of Bessel triangle).
E.g.f. for unsigned version: -exp(sqrt(1-2*x)-1). - Karol A. Penson, Mar 20 2010 [gives -1, 1, 0, 1, 5, 36, 329, ... ]
E.g.f. for unsigned version: 1/(sqrt(1-2*x))*exp(sqrt(1-2*x)-1). - Sergei N. Gladkovskii, Jul 03 2012
G.f.: 1/G(0) where G(k) = 1 - x + x*(2*k+1)/(1 - x + 2*x*(k+1)/G(k+1)); (continued fraction, 2-step). - Sergei N. Gladkovskii, Jul 10 2012
G.f.: 1+x/U(0) where U(k) = 1 - x + x*(k+1)/U(k+1) ; (continued fraction, Euler's 1st kind, 1-step). - Sergei N. Gladkovskii, Oct 06 2012
a(n) = BesselK[n+1/2,-1]/BesselK[5/2,-1]. - Vaclav Kotesovec, Aug 07 2013
|a(n)| ~ 2^(n+1/2)*n^n/exp(n+1). - Vaclav Kotesovec, Aug 07 2013
0 = a(n) * (a(n+2)) + a(n+1) * (-a(n+1) + 2*a(n+2) + a(n+3)) + a(n+2) * (-a(n+2)) for all n in Z. - Michael Somos, Jan 27 2014
a(n) = -i*(BesselK[3/2,1]*BesselI[n+3/2,-1] - BesselI[3/2,-1]*BesselK[n+3/2,1]), n>=0 for unsigned version - G. C. Greubel , Apr 19 2015
a(n) = hypergeom( [n+1, -n], [], 1/2). - Peter Luschny, Nov 10 2016
From G. C. Greubel, Aug 16 2017: (Start)
a(n) = (1/2)_{n} * (-2)^n * hypergeometric1f1(-n; -2*n; -2).
G.f.: (1/(1-t))*hypergeometric2f0(1, 1/2; -; -2*t/(1-t)^2). (End)

A278990 Number of loopless linear chord diagrams with n chords.

Original entry on oeis.org

1, 0, 1, 5, 36, 329, 3655, 47844, 721315, 12310199, 234615096, 4939227215, 113836841041, 2850860253240, 77087063678521, 2238375706930349, 69466733978519340, 2294640596998068569, 80381887628910919255, 2976424482866702081004, 116160936719430292078411
Offset: 0

Views

Author

N. J. A. Sloane, Dec 07 2016

Keywords

Comments

See the signed version of these numbers, A000806, for much more information about these numbers.
From Gus Wiseman, Feb 27 2019: (Start)
Also the number of 2-uniform set partitions of {1..2n} containing no two successive vertices in the same block. For example, the a(3) = 5 set partitions are:
{{1,3},{2,5},{4,6}}
{{1,4},{2,5},{3,6}}
{{1,4},{2,6},{3,5}}
{{1,5},{2,4},{3,6}}
{{1,6},{2,4},{3,5}}
(End)
From Gus Wiseman, Jul 05 2020: (Start)
Also the number of permutations of the multiset {1,1,2,2,...,n,n} with no two consecutive terms equal and where the first i appears before the first j for i < j. For example, the a(3) = 5 permutations are the following.
(1,2,3,1,2,3)
(1,2,3,1,3,2)
(1,2,3,2,1,3)
(1,2,3,2,3,1)
(1,2,1,3,2,3)
(End)

Crossrefs

Column k=0 of A079267.
Column k=2 of A293157.
Row n=2 of A322013.
Cf. A000110, A000699 (topologically connected 2-uniform), A000806, A001147 (2-uniform), A003436 (cyclical version), A005493, A170941, A190823 (distance 3+ version), A322402, A324011, A324172.
Anti-run compositions are A003242.
Separable partitions are A325534.
Other sequences involving the multiset {1,1,2,2,...,n,n}: A001147, A007717, A020555, A094574, A316972.

Programs

  • Magma
    [n le 2 select 2-n else (2*n-3)*Self(n-1) + Self(n-2): n in [1..30]]; // G. C. Greubel, Sep 26 2023
    
  • Mathematica
    RecurrenceTable[{a[n]== (2n-1)a[n-1] +a[n-2], a[0]==1, a[1]==0}, a, {n,0,20}] (* Vaclav Kotesovec, Sep 15 2017 *)
    FullSimplify[Table[-I*(BesselI[1/2+n,-1] BesselK[3/2,1] - BesselI[3/2,-1] BesselK[1/2+ n,1]), {n,0,20}]] (* Vaclav Kotesovec, Sep 15 2017 *)
    Table[(2 n-1)!! Hypergeometric1F1[-n,-2 n,-2], {n,0,20}] (* Eric W. Weisstein, Nov 14 2018 *)
    Table[Sqrt[2/Pi]/E ((-1)^n Pi BesselI[1/2+n,1] +BesselK[1/2+n,1]), {n,0,20}] // FunctionExpand // FullSimplify (* Eric W. Weisstein, Nov 14 2018 *)
    twouniflin[{}]:={{}};twouniflin[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@twouniflin[Complement[set,s]]]/@Table[{i,j},{j,Select[set,#>i+1&]}];
    Table[Length[twouniflin[Range[n]]],{n,0,14,2}] (* Gus Wiseman, Feb 27 2019 *)
  • PARI
    seq(N) = {
      my(a = vector(N)); a[1] = 0; a[2] = 1;
      for (n = 3, N, a[n] = (2*n-1)*a[n-1] + a[n-2]);
      concat(1, a);
    };
    seq(20) \\ Gheorghe Coserea, Dec 09 2016
    
  • SageMath
    def A278990_list(prec):
        P. = PowerSeriesRing(QQ, prec)
        return P( exp(-1+sqrt(1-2*x))/sqrt(1-2*x) ).egf_to_ogf().list()
    A278990_list(30) # G. C. Greubel, Sep 26 2023

Formula

From Gheorghe Coserea, Dec 09 2016: (Start)
D-finite with recurrence a(n) = (2*n-1)*a(n-1) + a(n-2), with a(0) = 1, a(1) = 0.
E.g.f. y satisfies: 0 = (1-2*x)*y'' - 3*y' - y.
a(n) - a(n-1) = A003436(n) for all n >= 2. (End)
From Vaclav Kotesovec, Sep 15 2017: (Start)
a(n) = sqrt(2)*exp(-1)*(BesselK(1/2 + n, 1)/sqrt(Pi) - i*sqrt(Pi)*BesselI(1/2 + n, -1)), where i is the imaginary unit.
a(n) ~ 2^(n+1/2) * n^n / exp(n+1). (End)
a(n) = A114938(n)/n! - Gus Wiseman, Jul 05 2020 (from Alexander Burstein's formula at A114938).
From G. C. Greubel, Sep 26 2023: (Start)
a(n) = (-1)^n * (i/e)*Sqrt(2/Pi) * BesselK(n + 1/2, -1).
G.f.: sqrt(Pi/(2*x)) * exp(-(1+x)^2/(2*x)) * Erfi((1+x)/sqrt(2*x)).
E.g.f.: exp(-1 + sqrt(1-2*x))/sqrt(1-2*x). (End)

Extensions

a(0)=1 prepended by Gheorghe Coserea, Dec 09 2016

A003437 Number of unlabeled Hamiltonian circuits on n-octahedron (cross polytope); also number of circular chord diagrams with n chords, modulo symmetries.

Original entry on oeis.org

0, 1, 2, 7, 29, 196, 1788, 21994, 326115, 5578431, 107026037, 2269254616, 52638064494, 1325663757897, 36021577975918, 1050443713185782, 32723148860301935, 1084545122297249077, 38105823782987999742, 1414806404051118314077
Offset: 1

Views

Author

Keywords

References

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

Crossrefs

See A003436 for labeled case.
See also A278990, A007474.

Programs

  • Mathematica
    nn = 20; M = Array[0&, {2nn, 2nn}];
    Mget[n_, k_] := Which[n < 0, 0, n==0, 1, n==1, 1-Mod[k, 2], n==2, k - Mod[k, 2], True, M[[n, k]]];
    Mset[n_, k_, v_] := (M[[n, k]] = v);
    Minit = Module[{tmp = 0}, For[n = 3, n <= 2nn, n++, For[k = 1, k <= 2nn, k++, tmp = If[OddQ[k], k(n-1) Mget[n-2, k] + Mget[n-4, k], Mget[n-1, k] + k(n-1) Mget[n-2, k] - Mget[n-3, k] + Mget[n-4, k]]; Mset[n, k, tmp]]]];
    A007474[n_] := Sum[EulerPhi[d] (Mget[2n/d, d] - Mget[2n/d - 2, d]), {d, Divisors[2n]}]/(2n);
    a[n_] := A007474[n]/2 + (Mget[n, 2] - Mget[n-1, 2] + Mget[n-2, 2])/4;
    Array[a, nn] (* Jean-François Alcover, Aug 12 2018, after Gheorghe Coserea *)
  • PARI
    N = 20; M = matrix(2*N, 2*N);
    Mget(n,k) = { if (n<0, 0, n==0, 1, n==1, 1-(k%2), n==2, k-(k%2), M[n,k]) };
    Mset(n,k,v) = { M[n,k] = v;};
    Minit() = {
      my(tmp = 0);
      for (n=3, 2*N, for(k=1, 2*N,
        tmp = if (k%2, k*(n-1) * Mget(n-2, k) + Mget(n-4, k),
        Mget(n-1, k) + k*(n-1) * Mget(n-2, k) - Mget(n-3, k) + Mget(n-4, k));
        Mset(n, k, tmp)));
    };
    Minit();
    A007474(n) = sumdiv(2*n, d, eulerphi(d) * (Mget(2*n/d, d) - Mget(2*n/d-2, d)))/(2*n);
    a(n) = A007474(n)/2 + (Mget(n,2) - Mget(n-1,2) + Mget(n-2,2))/4;
    vector(N, n, a(n))  \\ Gheorghe Coserea, Dec 10 2016

Formula

a(n) ~ 2^(n-3/2) * n^(n-1) / exp(n+1). - Vaclav Kotesovec, Dec 10 2016

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).

A324461 Number of simple graphs with n vertices and distinct rotations.

Original entry on oeis.org

1, 1, 0, 6, 48, 1020, 32232, 2097144, 268369920, 68719472640, 35184338533920, 36028797018963936, 73786976226114539520, 302231454903657293676480, 2475880078570197599844819072, 40564819207303340847860140736640, 1329227995784915854457062986570792960
Offset: 0

Views

Author

Gus Wiseman, Feb 28 2019

Keywords

Comments

A simple graph with n vertices has distinct rotations if all n rotations of its vertex set act on the edge set to give distinct graphs. These are different from aperiodic graphs and acyclic graphs but are similar to aperiodic sequences (A000740) and aperiodic arrays (A323867).

Crossrefs

Programs

  • Mathematica
    rotgra[g_,m_]:=Sort[Sort/@(g/.k_Integer:>If[k==m,1,k+1])];
    Table[Length[Select[Subsets[Subsets[Range[n],{2}]],UnsameQ@@Table[Nest[rotgra[#,n]&,#,j],{j,n}]&]],{n,0,5}]
  • PARI
    a(n)={if(n==0, 1, sumdiv(n, d, moebius(d)*2^(n*(n/d-1)/2 + n*(d\2)/d)))} \\ Andrew Howroyd, Aug 15 2019
    
  • Python
    from sympy import mobius, divisors
    def A324461(n): return sum(mobius(m:=n//d)<<(n*(d-1)>>1)+d*(m>>1) for d in divisors(n,generator=True)) if n else 1 # Chai Wah Wu, Jul 03 2024

Formula

a(n > 0) = A306715(n) * n.
a(n) = Sum_{d|n} mu(d)*2^(n*(n/d-1)/2 + n*floor(d/2)/d) for n > 0. - Andrew Howroyd, Aug 15 2019

Extensions

Terms a(7) and beyond from Andrew Howroyd, Aug 15 2019

A129348 Number of (directed) Hamiltonian circuits in the cocktail party graph of order n.

Original entry on oeis.org

0, 2, 32, 1488, 112512, 12771840, 2036229120, 434469611520, 119619533537280, 41303040523960320, 17481826772405452800, 8902337068174698086400, 5370014079716477003366400, 3786918976243761421064601600, 3087031512410698159166482022400, 2880726660365605475506018320384000
Offset: 1

Views

Author

Eric W. Weisstein, Apr 10 2007

Keywords

Comments

Also, the number of ways (up to rotations) to seat n married couples at a circular table with no spouses next to each other. Cf. A007060, A193639. - Geoffrey Critzer, Feb 09 2014
The cocktail party graph may also be called the n-octahedron, n-orthoplex or n-dimensional cross polytope. - Andrew Howroyd, May 14 2017

Crossrefs

Programs

  • Maple
    a:= proc(n) option remember; `if`(n<3, n*(n-1),
         ((136*n^3-608*n^2+762*n-470) *a(n-1)
           +4*(n-2)*(14*n^2+29*n-193) *a(n-2)
           -80*(n-2)*(n-3)*(n-4) *a(n-3)) /(34*n-101))
        end:
    seq(a(n), n=0..20);  # Alois P. Heinz, Feb 09 2014
  • Mathematica
    Prepend[Table[Sum[(-1)^i Binomial[n, i] (2n - 1 - i)! 2^i, {i, 0, n}], {n, 2, 16}], 0] (* Geoffrey Critzer, Feb 09 2014 *)
    Table[Piecewise[{{(-1 + 2 n)! Hypergeometric1F1[-n, 1 - 2 n, -2],
        n > 1}}], {n, 16}] (* Eric W. Weisstein, Mar 29 2014 *)
  • PARI
    { A129348(n) = sum(m=0,n-1, sum(k=1,n-m, (-1)^k * binomial(n-1,m) * binomial(n-m-1,k-1) * 2^(k-1) * ([0,k-1,2*(n-m-k);1,k-2,2*(n-m-k);1,k-1,2*(n-m-k-1)]^(2*n))[1,1] ) + sum(k=0,n-m, (-1)^k * binomial(n-1,m) * binomial(n-m-1,k) * 2^(k-1) * ([0,k,2*(n-m-k-1);2,k-1,2*(n-m-k-1);2,k,2*(n-m-k-2)]^(2*n))[1,1] ) ) } \\ Max Alekseyev, Dec 22 2013

Formula

For n>=2, a(n) = Sum_{k=0..n} (-1)^k*binomial(n,k)*(2*n-1-k)!*2^k. - Geoffrey Critzer, Feb 09 2014
Recurrence (for n>=4): (2*n-3)*a(n) = 2*(n-1)*(4*n^2 - 8*n + 5)*a(n-1) + 4*(n-2)*(n-1)*(2*n-1)*a(n-2). - Vaclav Kotesovec, Feb 09 2014
a(n) ~ sqrt(Pi) * 2^(2*n) * n^(2*n-1/2) / exp(2*n+1). - Vaclav Kotesovec, Feb 09 2014
For n>=2, a(n) = (-1 + 2 n)! Hypergeometric1F1[-n, 1 - 2 n, -2]. - Eric W. Weisstein, Mar 29 2014
a(n) = A003435(n) / (2*n) = A003436(n) * (n-1)! * 2^(n-1). - Andrew Howroyd, May 14 2017

Extensions

Terms a(6) onward from Max Alekseyev, Nov 10 2007

A348815 a(n) = number of chord labeled loopless diagrams by number of K_4.

Original entry on oeis.org

0, 1, 134, 75843, 83002866, 158861646466, 490294453324924, 2292204611710892971, 15459367618357013402267, 144663877588996810362218074, 1819753109993633276315632934129, 29976383544377113242613349012354566, 632574848906117234957565158900144038734
Offset: 1

Views

Author

Michael De Vlieger, Nov 01 2021

Keywords

Crossrefs

Extensions

a(11) onwards from Andrew Howroyd, Feb 05 2024

A369923 Array read by antidiagonals: A(n,k) is the number of permutations of n copies of 1..k with values introduced in order and without cyclically adjacent elements equal.

Original entry on oeis.org

0, 1, 0, 1, 1, 0, 1, 4, 1, 0, 1, 31, 22, 1, 0, 1, 293, 1415, 134, 1, 0, 1, 3326, 140343, 75843, 866, 1, 0, 1, 44189, 20167651, 83002866, 4446741, 5812, 1, 0, 1, 673471, 3980871156, 158861646466, 55279816356, 276154969, 40048, 1, 0
Offset: 1

Views

Author

Andrew Howroyd, Feb 05 2024

Keywords

Comments

Also, T(n,k) is the number of generalized chord labeled loopless diagrams with k parts of K_n. See the Krasko reference for a full definition.

Examples

			Array begins:
n\k| 1 2    3         4              5                    6 ...
---+-----------------------------------------------------------
 1 | 0 1    1         1              1                    1 ...
 2 | 0 1    4        31            293                 3326 ...
 3 | 0 1   22      1415         140343             20167651 ...
 4 | 0 1  134     75843       83002866         158861646466 ...
 5 | 0 1  866   4446741    55279816356     1450728060971387 ...
 6 | 0 1 5812 276154969 39738077935264 14571371516350429940 ...
 ...
		

Crossrefs

Column 3 is A197657, column 4 appears to be A209183(n)/2.
Cf. A322013 (without linearly adjacent elements equal), A322093.

Programs

  • Mathematica
    T[n_, k_] := If[k == 1, 0, Expand[(-1)^(k (n + 1))/(k - 1)! n Hypergeometric1F1[1 - n, 2, x]^k x^(k - 1)] /. x^p_ :> p!] (* Eric W. Weisstein, Feb 20 2025 *)
  • PARI
    \\ compare with A322013.
    q(n, x) = sum(i=1, n, (-1)^(n-i) * binomial(n-1, n-i) * x^i/i!)
    T(n, k) = if(k > 1, subst(serlaplace(n*q(n, x)^k/x), x, 1)/(k-1)!, 0)
Showing 1-10 of 29 results. Next