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

A006134 a(n) = Sum_{k=0..n} binomial(2*k,k).

Original entry on oeis.org

1, 3, 9, 29, 99, 351, 1275, 4707, 17577, 66197, 250953, 956385, 3660541, 14061141, 54177741, 209295261, 810375651, 3143981871, 12219117171, 47564380971, 185410909791, 723668784231, 2827767747951, 11061198475551, 43308802158651, 169719408596403, 665637941544507
Offset: 0

Views

Author

Keywords

Comments

The expression a(n) = B^n*Sum_{ k=0..n } binomial(2*k,k)/B^k gives A006134 for B=1, A082590 (B=2), A132310 (B=3), A002457 (B=4), A144635 (B=5). - N. J. A. Sloane, Jan 21 2009
T(n+1,1) from table A045912 of characteristic polynomial of negative Pascal matrix. - Michael Somos, Jul 24 2002
p divides a((p-3)/2) for p=11, 13, 23, 37, 47, 59, 61, 71, 73, 83, 97, 107, 109, 131, 157, 167, ...: A097933. Also primes congruent to {1, 2, 3, 11} mod 12 or primes p such that 3 is a square mod p (excluding 2 and 3) A038874. - Alexander Adamchuk, Jul 05 2006
Partial sums of the even central binomial coefficients. For p prime >=5, a(p-1) = 1 or -1 (mod p) according as p = 1 or -1 (mod 3) (see Pan and Sun link). - David Callan, Nov 29 2007
First column of triangle A187887. - Michel Marcus, Jun 23 2013
From Gus Wiseman, Apr 20 2023: (Start)
Also the number of nonempty subsets of {1,...,2n+1} with median n+1, where the median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length). The odd/even-length cases are A000984 and A006134(n-1). For example, the a(0) = 1 through a(2) = 9 subsets are:
{1} {2} {3}
{1,3} {1,5}
{1,2,3} {2,4}
{1,3,4}
{1,3,5}
{2,3,4}
{2,3,5}
{1,2,4,5}
{1,2,3,4,5}
Alternatively, a(n-1) is the number of nonempty subsets of {1,...,2n-1} with median n.
(End)

Examples

			1 + 3*x + 9*x^2 + 29*x^3 + 99*x^4 + 351*x^5 + 1275*x^6 + 4707*x^7 + 17577*x^8 + ...
		

References

  • Marko Petkovsek, Herbert Wilf and Doron Zeilberger, A=B, A K Peters, 1996, p. 22.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000984 (first differences), A097933, A038874, A132310.
Equals A066796 + 1.
Odd bisection of A100066.
Row sums of A361654 (also column k = 2).
A007318 counts subsets by length, A231147 by median, A013580 by integer median.
A359893 and A359901 count partitions by median.

Programs

  • MATLAB
    n=10; x=pascal(n); trace(x)
    
  • Magma
    &cat[ [&+[ Binomial(2*k, k): k in [0..n]]]: n in [0..30]]; // Vincenzo Librandi, Aug 13 2015
  • Maple
    A006134 := proc(n) sum(binomial(2*k,k),k=0..n); end;
    a := n -> -binomial(2*(n+1),n+1)*hypergeom([1,n+3/2],[n+2], 4) - I/sqrt(3):
    seq(simplify(a(n)), n=0..24); # Peter Luschny, Oct 29 2015
    # third program:
    A006134 := series(exp(2*x)*BesselI(0, 2*x) + exp(x)*int(BesselI(0, 2*x)*exp(x), x), x = 0, 25):
    seq(n!*coeff(A006134, x, n), n=0..24); # Mélika Tebni, Feb 27 2024
  • Mathematica
    Table[Sum[((2k)!/(k!)^2),{k,0,n}], {n,0,50}] (* Alexander Adamchuk, Jul 05 2006 *)
    a[ n_] := (4/3) Binomial[ 2 n, n] Hypergeometric2F1[ 1/2, 1, -n + 1/2, -1/3] (* Michael Somos, Jun 20 2012 *)
    Accumulate[Table[Binomial[2n,n],{n,0,30}]] (* Harvey P. Dale, Jan 11 2015 *)
    CoefficientList[Series[1/((1 - x) Sqrt[1 - 4 x]), {x, 0, 33}], x] (* Vincenzo Librandi, Aug 13 2015 *)
  • Maxima
    makelist(sum(binomial(2*k,k),k,0,n),n,0,12); /* Emanuele Munarini, Mar 15 2011 */
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( charpoly( matrix( n+1, n+1, i, j, -binomial( i+j-2, i-1))), 1))} \\ Michael Somos, Jul 10 2002
    
  • PARI
    {a(n)=binomial(2*n,n)*sum(k=0,2*n,(-1)^k*polcoeff((1+x+x^2)^n,k)/binomial(2*n,k))} \\ Paul D. Hanna, Aug 21 2007
    
  • PARI
    my(x='x+O('x^100)); Vec(1/((1-x)*sqrt(1-4*x))) \\ Altug Alkan, Oct 29 2015
    

Formula

From Alexander Adamchuk, Jul 05 2006: (Start)
a(n) = Sum_{k=0..n} (2k)!/(k!)^2.
a(n) = A066796(n) + 1, n>0. (End)
G.f.: 1/((1-x)*sqrt(1-4*x)).
D-finite with recurrence: (n+2)*a(n+2) - (5*n+8)*a(n+1) + 2*(2*n+3)*a(n) = 0. - Emanuele Munarini, Mar 15 2011
a(n) = C(2n,n) * Sum_{k=0..2n} (-1)^k*trinomial(n,k)/C(2n,k) where trinomial(n,k) = [x^k] (1 + x + x^2)^n. E.g. a(2) = C(4,2)*(1/1 - 2/4 + 3/6 - 2/4 + 1/1) = 6*(3/2) = 9 ; a(3) = C(6,3)*(1/1 - 3/6 + 6/15 - 7/20 + 6/15 - 3/6 + 1/1) = 20*(29/20) = 29. - Paul D. Hanna, Aug 21 2007
From Alzhekeyev Ascar M, Jan 19 2012: (Start)
a(n) = Sum_{ k=0..n } b(k)*binomial(n+k,k), where b(k)=0 for n-k == 2 (mod 3), b(k)=1 for n-k == 0 or 1 (mod 6), and b(k)=-1 for n-k== 3 or 4 (mod 6).
a(n) = Sum_{ k=0..n-1 } c(k)*binomial(2n,k) + binomial(2n,n), where c(k)=0 for n-k == 0 (mod 3), c(k)=1 for n-k== 1 (mod 3), and c(k)=-1 for n-k==2 (mod 3). (End)
a(n) ~ 2^(2*n+2)/(3*sqrt(Pi*n)). - Vaclav Kotesovec, Nov 06 2012
G.f.: G(0)/2/(1-x), where G(k)= 1 + 1/(1 - 2*x*(2*k+1)/(2*x*(2*k+1) + (k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 24 2013
G.f.: G(0)/(1-x), where G(k)= 1 + 4*x*(4*k+1)/( (4*k+2) - x*(4*k+2)*(4*k+3)/(x*(4*k+3) + (k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 26 2013
a(n) = Sum_{k = 0..n} binomial(n+1,k+1)*A002426(k). - Peter Bala, Oct 29 2015
a(n) = -binomial(2*(n+1),n+1)*hypergeom([1,n+3/2],[n+2], 4) - i/sqrt(3). - Peter Luschny, Oct 29 2015
a(n) = binomial(2*n, n)*hypergeom([1,-n], [1/2-n], 1/4). - Peter Luschny, Mar 16 2016
From Gus Wiseman, Apr 20 2023: (Start)
a(n+1) - a(n) = A000984(n).
a(n) = A013580(2n+1,n+1) (conjectured).
a(n) = 2*A024718(n) - 1.
a(n) = A100066(2n+1).
a(n) = A231147(2n+1,n+1) (conjectured). (End)
a(n) = Sum_{k=0..floor(n/3)} 3^(n-3*k) * binomial(n-k,2*k) * binomial(2*k,k) (Sawhney, 2017). - Amiram Eldar, Feb 24 2024
From Mélika Tebni, Feb 27 2024: (Start)
Limit_{n -> oo} a(n) / A281593(n) = 2.
E.g.f.: exp(2*x)*BesselI(0,2*x) + exp(x)*integral( BesselI(0,2*x)*exp(x) ) dx. (End)
a(n) = [(x*y)^n] 1/((1 - (x + y))*(1 - x*y)). - Stefano Spezia, Feb 16 2025
a(n) = Sum_{k=0..floor(n/2)} (-1)^k*binomial(2*n+1-k, n-2*k). - Michael Weselcouch, Jun 17 2025
a(n) = binomial(1+2*n, n)*hypergeom([1, (1-n)/2, -n/2], [-1-2*n, 2+n], 4). - Stefano Spezia, Jun 18 2025

Extensions

Simpler definition from Alexander Adamchuk, Jul 05 2006

A079309 a(n) = C(1,1) + C(3,2) + C(5,3) + ... + C(2*n-1,n).

Original entry on oeis.org

1, 4, 14, 49, 175, 637, 2353, 8788, 33098, 125476, 478192, 1830270, 7030570, 27088870, 104647630, 405187825, 1571990935, 6109558585, 23782190485, 92705454895, 361834392115, 1413883873975, 5530599237775, 21654401079325, 84859704298201, 332818970772253
Offset: 1

Views

Author

Miklos Kristof, Feb 10 2003

Keywords

Comments

a(n) is the sum of pyramid weights of all Dyck paths of length 2n (for pyramid weight see Denise and Simion). Equivalently, a(n) is the sum of the total lengths of end branches of an ordered tree, summation being over all ordered trees with n edges. For example, the five ordered trees with 3 edges have total lengths of endbranches 3,2,3,3 and 3. - Emeric Deutsch, May 30 2003
a(n) is the number of Motzkin paths of length 2n with exactly one level segment. (A level segment is a maximal sequence of contiguous flatsteps.) Example: for n=2, the paths counted are FFFF, FFUD, UDFF, UFFD. The formula for a(n) below counts these paths by length of the level segment. - David Callan, Jul 15 2004
The inverse Catalan transform yields A024495, shifted once left. - R. J. Mathar, Jul 07 2009
From Paul Barry, Mar 29 2010: (Start)
Hankel transform is A138341.
The aerated sequence 0, 0, 1, 0, 4, 0, 14, 0, 49, ... has e.g.f. int(cosh(x-t)*Bessel_I(1,2t), t = 0..x). (End)
a(n) is the number of terms of A031443 not exceeding 4^n. - Vladimir Shevelev, Oct 01 2010
Also the number of nonempty subsets of {1..2n} with median n, bisection of A361801. The version containing n is A001700 (bisected). Replacing 2n with 2n+1 and n with n+1 gives A006134. For mean instead of median we have A212352. - Gus Wiseman, Apr 16 2023

Examples

			a(4) = C(1,1) + C(3,2) + C(5,3) + C(7,4) = 1 + 3 + 10 + 35 = 49.
G.f. = x + 4*x^2 + 14*x^3 + 49*x^4 + 175*x^5 + 637*x^6 + 2353*x^7 + ...
From _Gus Wiseman_, Apr 16 2023: (Start)
The a(1) = 1 through a(3) = 14 subsets of {1..2n} with median n:
  {1}  {2}      {3}
       {1,3}    {1,5}
       {1,2,3}  {2,4}
       {1,2,4}  {1,3,4}
                {1,3,5}
                {1,3,6}
                {2,3,4}
                {2,3,5}
                {2,3,6}
                {1,2,4,5}
                {1,2,4,6}
                {1,2,3,4,5}
                {1,2,3,4,6}
                {1,2,3,5,6}
(End)
		

Crossrefs

Equals A024718(n) - 1.
This is the even (or odd) bisection of A361801.
A007318 counts subsets by length, A327481 by mean, A013580 by median.
A359893 and A359901 count partitions by median.

Programs

  • Maple
    a := n -> add(binomial(2*j, j)/2, j=1..n): seq(a(n), n=1..24); # Zerinvary Lajos, Oct 25 2006
    a := n -> add(abs(binomial(-j, -2*j)), j=1..n): seq(a(n), n=1..24); # Zerinvary Lajos, Oct 03 2007
    f:= gfun:-rectoproc({n*a(n) +(-5*n+2)*a(n-1) +2*(2*n-1)*a(n-2)=0,a(1)=1,a(2)=4},a(n),remember):
    map(f, [$1..100]); # Robert Israel, Jun 24 2015
  • Mathematica
    Rest[CoefficientList[Series[(1/Sqrt[1-4*x]-1)/(1-x)/2, {x, 0, 20}], x]] (* Vaclav Kotesovec, Feb 13 2014 *)
    Accumulate[Table[Binomial[2n-1,n],{n,30}]] (* Harvey P. Dale, Jan 06 2021 *)
  • PARI
    {a(n) = sum(k=1, n, binomial(2*k - 1, k))}; /* Michael Somos, Feb 14 2006 */
    
  • PARI
    my(x='x+O('x^40)); Vec((1/sqrt(1-4*x)-1)/(1-x)/2) \\ Altug Alkan, Dec 24 2015

Formula

a(n) = (1/2)*(C(2, 1) + C(4, 2) + C(6, 3) + ... + C(2*n, n)) = A066796(n)/2. - Vladeta Jovovic, Feb 12 2003
G.f.: (1/sqrt(1 - 4*x) - 1)/(1 - x)/2. - Vladeta Jovovic, Feb 12 2003
Given g.f. A(x), then x * A(x - x^2) is g.f. of A024495. - Michael Somos, Feb 14 2006
a(n) = A066796(n)/2. - Zerinvary Lajos, Oct 25 2006
a(n) = Sum_{0 <= i <= j <= n} binomial(i+j, i). - Benoit Cloitre, Nov 25 2006
D-finite with recurrence n*a(n) + (-5*n+2)*a(n-1) + 2*(2*n-1)*a(n-2) = 0. - R. J. Mathar, Nov 30 2012
a(n) ~ 2^(2*n+1) / (3*sqrt(Pi*n)). - Vaclav Kotesovec, Feb 13 2014
a(n) = Sum_{k=0..n-1} A001700(k). - Doug Bell, Jun 23 2015
a(n) = -binomial(2*n+1, n)*hypergeom([1, n+3/2], [n+2], 4) - (i/sqrt(3) + 1)/2. - Peter Luschny, May 18 2018
From Gus Wiseman, Apr 18 2023: (Start)
a(n) = A024718(n) - 1.
a(n) = A231147(2n+1,n).
a(n) = A361801(2n) = A361801(2n+1). (End)
a(n) = Sum_{k=0..floor(n/2)} (-1)^k*binomial(2*n+2-k, n-2*k). - Michael Weselcouch, Jun 17 2025
a(n) = binomial(2*(1+n), n)*hypergeom([1, (1-n)/2, -n/2], [-2*(1+n), 3+n], 4). - Stefano Spezia, Jun 18 2025

Extensions

More terms from Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Feb 11 2003

A279570 Number of length n inversion sequences avoiding the patterns 110 and 120.

Original entry on oeis.org

1, 1, 2, 6, 22, 92, 423, 2091, 10950, 60120, 343453, 2029809, 12354661, 77168197, 493189283, 3217459119, 21382723456, 144518555231, 991885282987, 6904454991721, 48691257834999, 347542736059492, 2508603139285095, 18297609829743478, 134772911886028731
Offset: 0

Views

Author

Megan A. Martinez, Feb 21 2017

Keywords

Comments

A length n inversion sequence e_1e_2...e_n is a sequence of integers where 0 <= e_i <= i-1. The term a(n) counts those length n inversion sequences with no entries e_i, e_j, e_k (where i e_k and e_i > e_k. This is the same as the set of length n inversion sequences avoiding 110 and 120.

Examples

			The length 4 inversion sequences avoiding (110,120) are 0000, 0001, 0002, 0003, 0010, 0011, 0012, 0013, 0020, 0021, 0022, 0023, 0100, 0101, 0102, 0103, 0111, 0112, 0113, 0121, 0122, 0123.
		

Crossrefs

Extensions

a(10)-a(24) from Alois P. Heinz, Feb 21 2017

A279564 Number of length n inversion sequences avoiding the patterns 000 and 100.

Original entry on oeis.org

1, 1, 2, 5, 16, 60, 260, 1267, 6850, 40572, 260812, 1805646, 13377274, 105487540, 881338060, 7770957903, 72060991394, 700653026744, 7123871583656, 75561097962918, 834285471737784, 9570207406738352, 113855103776348136, 1402523725268921870, 17863056512845724036, 234910502414771617316, 3185732802058088068444, 44501675392317774477088
Offset: 0

Views

Author

Megan A. Martinez, Feb 09 2017

Keywords

Comments

A length n inversion sequence e_1e_2...e_n is a sequence of integers where 0 <= e_i <= i-1. The term a(n) counts those length n inversion sequences with no entries e_i, e_j, e_k (where i= e_j = e_k. This is the same as the set of length n inversion sequences avoiding 000 and 100.

Crossrefs

Programs

  • Maple
    b:= proc(n, i, m, s) option remember; `if`(n=0, 1, add(
          `if`(j in s, 0, b(n-1, i+1, max(m, j),
          `if`(j<=m, s union {j}, s))), j=1..i))
        end:
    a:= n-> b(n, 1, 0, {}):
    seq(a(n), n=0..15);  # Alois P. Heinz, Feb 22 2017
  • Mathematica
    b[n_, i_, m_, s_List] := b[n, i, m, s] = If[n == 0, 1, Sum[If[MemberQ[s, j], 0, b[n-1, i+1, Max[m, j], If[j <= m, s ~Union~ {j}, s]]], {j, 1, i}] ]; a[n_] := b[n, 1, 0, {}]; Table[a[n], {n, 0, 15}] (* Jean-François Alcover, Jul 10 2017, after Alois P. Heinz *)

Formula

The length 4 inversion sequences avoiding (000,100) are 0011, 0012, 0013, 0021, 0022, 0023, 0101, 0102, 0103, 0110, 0112, 0113, 0120, 0121, 0122, 0123.

Extensions

a(10)-a(23) from Alois P. Heinz, Feb 22 2017
a(24)-a(27) from Vaclav Kotesovec, Oct 08 2021

A013580 Triangle formed in same way as Pascal's triangle (A007318) except 1 is added to central element in even-numbered rows.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 4, 4, 1, 1, 5, 9, 5, 1, 1, 6, 14, 14, 6, 1, 1, 7, 20, 29, 20, 7, 1, 1, 8, 27, 49, 49, 27, 8, 1, 1, 9, 35, 76, 99, 76, 35, 9, 1, 1, 10, 44, 111, 175, 175, 111, 44, 10, 1, 1, 11, 54, 155, 286, 351, 286, 155, 54, 11, 1, 1, 12, 65, 209, 441, 637, 637, 441, 209, 65
Offset: 0

Views

Author

Martin Hecko (bigusm(AT)interramp.com)

Keywords

Comments

From Gus Wiseman, Apr 19 2023: (Start)
Appears to be the number of nonempty subsets of {1,...,n} with median k, where the median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length). For example, row n = 5 counts the following subsets:
{1} {2} {3} {4} {5}
{1,3} {1,5} {3,5}
{1,2,3} {2,4} {1,4,5}
{1,2,4} {1,3,4} {2,4,5}
{1,2,5} {1,3,5} {3,4,5}
{2,3,4}
{2,3,5}
{1,2,4,5}
{1,2,3,4,5}
Including half-steps gives A231147.
For mean instead of median we have A327481.
(End)

Examples

			Triangle begins:
   1
   1   1
   1   3   1
   1   4   4   1
   1   5   9   5   1
   1   6  14  14   6   1
   1   7  20  29  20   7   1
   1   8  27  49  49  27   8   1
   1   9  35  76  99  76  35   9   1
   1  10  44 111 175 175 111  44  10   1
   1  11  54 155 286 351 286 155  54  11   1
   1  12  65 209 441 637 637 441 209  65  12   1
		

Crossrefs

Row sums give A000975, A054106.
Central diagonal T(2n+1,n+1) appears to be A006134.
Central diagonal T(2n,n) appears to be A079309.
For partitions instead of subsets we have A359901, row sums A325347.
A000975 counts subsets with integer median.
A007318 counts subsets by length, A359893 by twice median.

Programs

  • Mathematica
    CoefficientList[CoefficientList[Series[1/(1 - (1 + y)*x)/(1 - y*x^2), {x, 0, 10}, {y, 0, 10}], x], y] // Flatten (* G. C. Greubel, Oct 10 2017 *)

Formula

G.f.: 1/(1-(1+y)*x)/(1-y*x^2). - Vladeta Jovovic, Oct 12 2003

Extensions

More terms from James Sellers

A279566 Number of length n inversion sequences avoiding the patterns 102 and 201.

Original entry on oeis.org

1, 1, 2, 6, 22, 87, 354, 1465, 6154, 26223, 113236, 494870, 2185700, 9743281, 43784838, 198156234, 902374498, 4131895035, 19012201080, 87864535600, 407664831856, 1898184887679, 8867042353912, 41543375724751, 195164372948152, 919138464708907, 4338701289961694, 20524046955770940
Offset: 0

Views

Author

Megan A. Martinez, Feb 09 2017

Keywords

Comments

A length n inversion sequence e_1e_2...e_n is a sequence of integers where 0 <= e_i <= i-1. The term a(n) counts those length n inversion sequences with no entries e_i, e_j, e_k (where i e_j < e_k and e_i <> e_k. This is the same as the set of length n inversion sequences avoiding 102 and 201.

Examples

			The length 4 inversion sequences avoiding (102, 201) are 0000, 0001, 0002, 0003, 0010, 0011, 0012, 0013, 0020, 0021, 0022, 0023, 0100, 0101, 0110, 0111, 0112, 0113, 0120, 0121, 0122, 0123
		

Crossrefs

Formula

G.f.: (-8*x^4 + 18*x^3 - 10*x^2 - 8*x + 4 + 2 * (2*x - 1) * (x^2 - 2*x + 2) * ((5*x - 1)*(x - 1))^(1/2)) / (4*x * (2*x - 1) * (x - 1) * (x - 2)^2). - Benjamin Testart, Jul 12 2024
a(n) ~ 41 * 5^(n + 3/2) / (648 * sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Nov 21 2024

Extensions

a(10)-a(11) from Alois P. Heinz, Feb 24 2017
a(12)-a(17) from Bert Dobbelaere, Dec 30 2018
a(18) and beyond from Benjamin Testart, Jul 12 2024

A047653 Constant term in expansion of (1/2) * Product_{k=-n..n} (1 + x^k).

Original entry on oeis.org

1, 2, 4, 10, 26, 76, 236, 760, 2522, 8556, 29504, 103130, 364548, 1300820, 4679472, 16952162, 61790442, 226451036, 833918840, 3084255128, 11451630044, 42669225172, 159497648600, 597950875256, 2247724108772, 8470205600640, 31991616634296, 121086752349064
Offset: 0

Views

Author

Keywords

Comments

Or, constant term in expansion of Product_{k=1..n} (x^k + 1/x^k)^2. - N. J. A. Sloane, Jul 09 2008
Or, maximal coefficient of the polynomial (1+x)^2 * (1+x^2)^2 *...* (1+x^n)^2.
a(n) = A000302(n) - A181765(n).
From Gus Wiseman, Apr 18 2023: (Start)
Also the number of subsets of {1..2n} that are empty or have mean n. The a(0) = 1 through a(3) = 10 subsets are:
{} {} {} {}
{1} {2} {3}
{1,3} {1,5}
{1,2,3} {2,4}
{1,2,6}
{1,3,5}
{2,3,4}
{1,2,3,6}
{1,2,4,5}
{1,2,3,4,5}
Also the number of subsets of {-n..n} with no 0's but with sum 0. The a(0) = 1 through a(3) = 10 subsets are:
{} {} {} {}
{-1,1} {-1,1} {-1,1}
{-2,2} {-2,2}
{-2,-1,1,2} {-3,3}
{-3,1,2}
{-2,-1,3}
{-2,-1,1,2}
{-3,-1,1,3}
{-3,-2,2,3}
{-3,-2,-1,1,2,3}
(End)

Crossrefs

Cf. A025591.
Cf. A053632; variant: A127728.
For median instead of mean we have A079309(n) + 1.
Odd bisection of A133406.
A000980 counts nonempty subsets of {1..2n-1} with mean n.
A007318 counts subsets by length, A327481 by mean.

Programs

  • Maple
    f:=n->coeff( expand( mul((x^k+1/x^k)^2,k=1..n) ),x,0);
    # second Maple program:
    b:= proc(n, i) option remember; `if`(n>i*(i+1)/2, 0,
          `if`(i=0, 1, 2*b(n, i-1)+b(n+i, i-1)+b(abs(n-i), i-1)))
        end:
    a:=n-> b(0, n):
    seq(a(n), n=0..40);  # Alois P. Heinz, Mar 10 2014
  • Mathematica
    b[n_, i_] := b[n, i] = If[n>i*(i+1)/2, 0, If[i == 0, 1, 2*b[n, i-1]+b[n+i, i-1]+b[Abs[n-i], i-1]]]; a[n_] := b[0, n]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Mar 10 2014, after Alois P. Heinz *)
    nmax = 26; d = {1}; a1 = {};
    Do[
      i = Ceiling[Length[d]/2];
      AppendTo[a1, If[i > Length[d], 0, d[[i]]]];
      d = PadLeft[d, Length[d] + 2 n] + PadRight[d, Length[d] + 2 n] +
        2 PadLeft[PadRight[d, Length[d] + n], Length[d] + 2 n];
    , {n, nmax}];
    a1 (* Ray Chandler, Mar 15 2014 *)
    Table[Length[Select[Subsets[Range[2n]],Length[#]==0||Mean[#]==n&]],{n,0,6}] (* Gus Wiseman, Apr 18 2023 *)
  • PARI
    a(n)=polcoeff(prod(k=-n,n,1+x^k),0)/2
    
  • PARI
    {a(n)=sum(k=0,n*(n+1)/2,polcoeff(prod(m=1,n,1+x^m+x*O(x^k)),k)^2)} \\ Paul D. Hanna, Nov 30 2010

Formula

Sum of squares of coefficients in Product_{k=1..n} (1+x^k):
a(n) = Sum_{k=0..n(n+1)/2} A053632(n,k)^2. - Paul D. Hanna, Nov 30 2010
a(n) = A000980(n)/2.
a(n) ~ sqrt(3) * 4^n / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Sep 11 2014
From Gus Wiseman, Apr 18 2023 (Start)
a(n) = A133406(2n+1).
a(n) = A212352(n) + 1.
a(n) = A362046(2n) + 1.
(End)

Extensions

More terms from Michael Somos, Jun 10 2000

A279565 Number of length n inversion sequences avoiding the patterns 100, 110, 120, 201, and 210.

Original entry on oeis.org

1, 1, 2, 6, 21, 81, 332, 1420, 6266, 28318, 130412, 609808, 2887582, 13818590, 66726628, 324713196, 1590853485, 7840315329, 38843186366, 193342353214, 966409013021, 4848846341569, 24412146213116, 123290812268404, 624448756434476, 3171046361310556
Offset: 0

Views

Author

Megan A. Martinez, Feb 09 2017

Keywords

Comments

A length n inversion sequence e_1e_2...e_n is a sequence of integers where 0 <= e_i <= i-1. The term a(n) counts those length n inversion sequences with no entries e_i, e_j, e_k (where i e_k. This is the same as the set of length n inversion sequences avoiding 100, 110, 120, 201, and 210.

Examples

			The length 4 inversion sequences avoiding (100, 110, 120, 201, 210) are 0000, 0001, 0002, 0003, 0010, 0011, 0012, 0013, 0020, 0021, 0022, 0023, 0101, 0102, 0103, 0111, 0112, 0113, 0121, 0122, 0123.
		

Crossrefs

Programs

  • Magma
    I:=[6, 21, 81]; [1,1,2] cat [n le 3 select I[n] else ( (n+1)*(17*n+6)*Self(n-1) +(49*n^2+11*n+22)*Self(n-2) +3*(3*n-1)*(3*n-2)*Self(n-3) )/(5*(n+2)*(n+1)) : n in [1..30]]; // G. C. Greubel, Mar 29 2019
    
  • Maple
    a:= proc(n) option remember; `if`(n<3, n!,
          ((n-1)*(17*n-28)*a(n-1) +(49*n^2-185*n+196)*a(n-2)
           +(3*(3*n-7))*(3*n-8)*a(n-3)) / (5*n*(n-1)))
        end:
    seq(a(n), n=0..30);  # Alois P. Heinz, Feb 22 2017
  • Mathematica
    a[n_] := a[n] = If[n < 3, n!, (((n - 1)*(17*n - 28)*a[n-1] + (49*n^2 - 185*n + 196)*a[n-2] + (3*(3*n - 7))*(3*n - 8)*a[n-3]) / (5*n*(n - 1)))]; Array[a, 30, 0] (* Jean-François Alcover, Nov 06 2017, after Alois P. Heinz *)
    Join[{1}, Table[(1/n)*Sum[m*Sum[Binomial[k, n-m-k]*Binomial[n+k-1, k], {k, 0, n-m}], {m, 1, n}], {n, 1, 30}]] (* G. C. Greubel, Mar 29 2019 *)
  • Maxima
    a(n):=if n=0 then 1 else sum(m*sum(binomial(k,n-m-k)*binomial(n+k-1,k),k,0,n-m),m,1,n)/n; /* Vladimir Kruchinin, Mar 26 2019 */
    
  • PARI
    my(x='x+O('x^30)); Vec(round(3/(4-4*sin(asin((27*x+11)/16)/3)))) \\ G. C. Greubel, Mar 29 2019
    
  • Sage
    [1] +[(1/n)*(sum(sum(k*binomial(j,n-k-j)*binomial(n+j-1,j) for j in (0..n-k)) for k in (1..n))) for n in (1..30)] # G. C. Greubel, Mar 29 2019

Formula

G.f.: 3/(4-4*sin(asin((27*x+11)/16)/3)). - Vladimir Kruchinin, Mar 25 2019
a(n) = (1/n)*Sum_{m=1..n} m*Sum_{k=0..n-m} C(k,n-m-k)*C(n+k-1,k), n>0, a(0)=1. - Vladimir Kruchinin, Mar 26 2019
a(n) ~ 3^(3*n + 1/2) / (2^(7/2) * sqrt(Pi) * n^(3/2) * 5^(n - 1/2)). - Vaclav Kotesovec, Oct 07 2021
Conjecture: a(n) = (v_n + v_{n+1})/2 for n > 0 with a(0) = 1 where we start with vector v of fixed length m with elements v_i = 1 and for i=1..m-2, for j=i+2..m apply v_j := Sum_{k=0..2} v_{j-k}. - Mikhail Kurkov, Sep 03 2024

Extensions

a(10)-a(25) from Alois P. Heinz, Feb 22 2017

A279561 Number of length n inversion sequences avoiding the patterns 101, 102, 201, and 210.

Original entry on oeis.org

1, 1, 2, 6, 21, 77, 287, 1079, 4082, 15522, 59280, 227240, 873886, 3370030, 13027730, 50469890, 195892565, 761615285, 2965576715, 11563073315, 45141073925, 176423482325, 690215089745, 2702831489825, 10593202603775, 41550902139551, 163099562175851
Offset: 0

Views

Author

Megan A. Martinez, Jan 17 2017

Keywords

Comments

A length n inversion sequence e_1e_2...e_n is a sequence of integers where 0 <= e_i <= i-1. The term a(n) counts those length n inversion sequences with no entries e_i, e_j, e_k (where i e_j <> e_k. This is the same as the set of length n inversion sequences avoiding 101, 102, 201, and 210.
It is conjectured that a_n also counts those length n inversion sequences with no entries e_i, e_j, e_k (where i e_k and e_i <> e_k. This is the same as the set of length n inversion sequences avoiding 021 and 120.

Examples

			The length 4 inversion sequences avoiding (101, 102, 201, 210) are 0000, 0001, 0002, 0003, 0010, 0011, 0012, 0013, 0020, 0021, 0022, 0023, 0100, 0110, 0111, 0112, 0113, 0120, 0121, 0122, 0123.
The length 4 inversion sequences avoiding (021, 120) are 0000, 0001, 0002, 0003, 0010, 0011, 0012, 0013, 0020, 0022, 0023, 0100, 0101, 0102, 0103, 0110, 0111, 0112, 0113, 0122, 0123.
		

Crossrefs

Programs

  • Maple
    a:= proc(n) option remember; `if`(n<3, 1+n*(n-1)/2,
          ((5*n^2-12*n+6)*a(n-1)-(4*n^2-10*n+6)*a(n-2))/((n-2)*n))
        end:
    seq(a(n), n=0..30);  # Alois P. Heinz, Jan 18 2017
  • Mathematica
    a[n_] := 1 + Sum[Binomial[2i, i-1], {i, 0, n-1}];
    Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Mar 28 2017 *)

Formula

a(n) = 1 + Sum_{i=1..n-1} binomial(2i, i-1).
a(n) = 1 + A057552(n-2).
G.f.: (1-4*x+sqrt(-16*x^3+20*x^2-8*x+1))/(2*(x-1)*(4*x-1)).
D-finite with recurrence: n*a(n) +(-7*n+6)*a(n-1) +2*(7*n-13)*a(n-2) +4*(-2*n+5)*a(n-3)=0. - R. J. Mathar, Feb 21 2020

A279562 Number of length n inversion sequences avoiding the patterns 100, 102, and 201.

Original entry on oeis.org

1, 1, 2, 6, 21, 78, 299, 1176, 4729, 19378, 80667, 340260, 1451277, 6248758, 27124703, 118576648, 521574769, 2306766426, 10251761219, 45759404076, 205050758165, 922104978430, 4160045001703, 18823187479504, 85400356099001, 388422301113250, 1770695668597643, 8089198184655732, 37027394471695197
Offset: 0

Views

Author

Megan A. Martinez, Feb 09 2017

Keywords

Comments

A length n inversion sequence e_1e_2...e_n is a sequence of integers where 0 <= e_i <= i-1. The term a(n) counts those length n inversion sequences with no entries e_i, e_j, e_k (where i e_j <= e_k and e_i <> e_k. This is the same as the set of length n inversion sequences avoiding 100, 102, and 201.

Examples

			The length 4 inversion sequences avoiding (100, 102, 201) are 0000, 0001, 0002, 0003, 0010, 0011, 0012, 0013, 0020, 0021, 0022, 0023, 0101, 0110, 0111, 0112, 0113, 0120, 0121, 0122, 0123.
		

Crossrefs

Formula

G.f.: (2 + x - 10*x^2 + 4*x^3 - (2-3*x)*(1 - 4*x - 4*x^2)^(1/2)) / (8*x*(1 - x)^2). - Nathan J. Britt, Jun 08 2025
a(n) ~ c * (2 + 2*sqrt(2))^n / n^(3/2), where c = 0.40413545332026258682681691461076303199449216224437... - Nathan J. Britt, Jun 08 2025

Extensions

a(10)-a(12) from Alois P. Heinz, Feb 24 2017
a(13)-a(17) from Bert Dobbelaere, Dec 30 2018
More terms from Nathan J. Britt, Jun 08 2025
Showing 1-10 of 28 results. Next