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.

A004171 a(n) = 2^(2n+1).

Original entry on oeis.org

2, 8, 32, 128, 512, 2048, 8192, 32768, 131072, 524288, 2097152, 8388608, 33554432, 134217728, 536870912, 2147483648, 8589934592, 34359738368, 137438953472, 549755813888, 2199023255552, 8796093022208, 35184372088832, 140737488355328, 562949953421312
Offset: 0

Views

Author

Keywords

Comments

Same as Pisot sequences E(2, 8), L(2, 8), P(2, 8), T(2, 8). See A008776 for definitions of Pisot sequences.
In the Chebyshev polynomial of degree 2n, a(n) is the coefficient of x^2n. - Benoit Cloitre, Mar 13 2002
1/2 - 1/8 + 1/32 - 1/128 + ... = 2/5. - Gary W. Adamson, Mar 03 2009
From Adi Dani, May 15 2011: (Start)
Number of ways of placing an even number of indistinguishable objects in n + 1 distinguishable boxes with at most 3 objects in box.
Number of compositions of even natural numbers into n + 1 parts less than or equal to 3 (0 is counted as part). (End)
Also the number of maximal cliques in the (n+1)-Sierpinski tetrahedron graph for n > 0. - Eric W. Weisstein, Dec 01 2017
Assuming the Collatz conjecture is true, any starting number eventually leads to a power of 2. A number in this sequence can never be the first power of 2 in a Collatz sequence except of course for the Collatz sequence starting with that number. For example, except for 8, 4, 2, 1, any Collatz sequence that includes 8 must also include 16 (e.g., 5, 16, 8, 4, 2, 1). - Alonso del Arte, Oct 01 2019
First differences of A020988, and thus the "wavelengths" of the local maxima in A020986. See the Brillhart and Morton link, pp. 855-856. - John Keith, Mar 04 2021

Examples

			G.f. = 2 + 8*x + 32*x^2 + 128*x^3 + 512*x^4 + 2048*x^5 + 8192*x^6 + 32768*x^7 + ...
From _Adi Dani_, May 15 2011: (Start)
a(1) = 8 because all compositions of even natural numbers into 2 parts less than or equal to 3 are:
  for 0: (0, 0)
  for 2: (0, 2), (2, 0), (1, 1)
  for 4: (1, 3), (3, 1), (2, 2)
  for 6: (3, 3).
a(2) = 32 because all compositions of even natural numbers into 3 parts less than or equal to 3 are:
  for 0: (0, 0, 0)
  for 2: (0, 0, 2), (0, 2, 0), (2, 0, 0), (0, 1, 1), (1, 0, 1) , (1, 1, 0)
  for 4: (0, 1, 3), (0, 3, 1), (1, 0, 3), (1, 3, 0), (3, 0, 1), (3, 1, 0), (0, 2, 2), (2, 0, 2), (2, 2, 0), (1, 1, 2), (1, 2, 1), (2, 1, 1)
  for 6: (0, 3, 3), (3, 0, 3), (3, 3, 0), (1, 2, 3), (1, 3, 2), (2, 1, 3), (2, 3, 1), (3, 1, 2), (3, 2, 1), (2, 2, 2)
  for 8: (2, 3, 3), (3, 2, 3), (3, 3, 2).
(End)
		

References

  • Adi Dani, Quasicompositions of natural numbers, Proceedings of III congress of mathematicians of Macedonia, Struga Macedonia 29 IX -2 X 2005 pages 225-238.

Crossrefs

Absolute value of A009117. Essentially the same as A081294.
Cf. A132020, A164632. Equals A000980(n) + 2*A181765(n). Cf. A013776.

Programs

Formula

a(n) = 2*4^n.
a(n) = 4*a(n-1).
1 = 1/2 + Sum_{n >= 1} 3/a(n) = 3/6 + 3/8 + 3/32 + 3/128 + 3/512 + 3/2048 + ...; with partial sums: 1/2, 31/32, 127/128, 511/512, 2047/2048, ... - Gary W. Adamson, Jun 16 2003
From Philippe Deléham, Nov 23 2008: (Start)
a(n) = 2*A000302(n).
G.f.: 2/(1-4*x). (End)
a(n) = A081294(n+1) = A028403(n+1) - A000079(n+1) for n >= 1. a(n-1) = A028403(n) - A000079(n). - Jaroslav Krizek, Jul 27 2009
E.g.f.: 2*exp(4*x). - Ilya Gutkovskiy, Nov 01 2016
a(n) = A002063(n)/3 - A000302(n). - Zhandos Mambetaliyev, Nov 19 2016
a(n) = Sum_{k = 0..2*n} (-1)^(k+n)*binomial(4*n + 2, 2*k + 1); a(2*n) = Sum_{k = 0..2*n} binomial(4*n + 2, 2*k + 1) = A013776(n). - Peter Bala, Nov 25 2016
Product_{n>=0} (1 - 1/a(n)) = A132020. - Amiram Eldar, May 08 2023

A081294 Expansion of (1-2*x)/(1-4*x).

Original entry on oeis.org

1, 2, 8, 32, 128, 512, 2048, 8192, 32768, 131072, 524288, 2097152, 8388608, 33554432, 134217728, 536870912, 2147483648, 8589934592, 34359738368, 137438953472, 549755813888, 2199023255552, 8796093022208, 35184372088832
Offset: 0

Views

Author

Paul Barry, Mar 17 2003

Keywords

Comments

Binomial transform of A046717. Second binomial transform of A000302 (with interpolated zeros). Partial sums are A007583.
Counts closed walks of length 2n at a vertex of the cyclic graph on 4 nodes C_4. With interpolated zeros, counts closed walks of length n at a vertex of the cyclic graph on 4 nodes C_4. - Paul Barry, Mar 10 2004
In general, Sum_{k=0..n} Sum_{j=0..n} C(2*(n-k), j)*C(2*k, j)*r^j has expansion (1 - (r+1)*x)/(1 - (r+3)*x - (r-1)*(r+3)*x^2 + (r-1)^3*x^3). - Paul Barry, Jun 04 2005 [corrected by Jason Yuen, Jan 20 2025]
a(n) is the number of binary strings of length 2n with an even number of 0's (and hence an even number of 1's). - Toby Gottfried, Mar 22 2010
Number of compositions of n where there are 2 sorts of part 1, 4 sorts of part 2, 8 sorts of part 3, ..., 2^k sorts of part k. - Joerg Arndt, Aug 04 2014
a(n) is also the number of permutations simultaneously avoiding 231 and 321 in the classical sense which can be realized as labels on an increasing strict binary tree with 2n-1 nodes. See A245904 for more information on increasing strict binary trees. - Manda Riehl Aug 07 2014
INVERT transform of powers of 2 (A000079). - Alois P. Heinz, Feb 11 2021
a(n) is the number of elements in an n-interval of the binomial poset of even-sized subsets of positive integers, cf. Stanley reference and second formula by Paul Barry. Each multichain 0 = x_0 <= x_1 <= x_2 = 1 in such an n-interval corresponds to a closed walk described above by Paul Barry. More generally, each multichain 0 = x_0 <= x_1 <= ... <= x_k = 1 corresponds to a closed walk of length 2n on the k-dimensional hypercube, cf. A054879, A092812, A121822. - Geoffrey Critzer, Apr 21 2023

Examples

			G.f. = 1 + 2*x + 8*x^2 + 32*x^3 + 128*x^4 + 512*x^5 + 2048*x^6 + 8192*x^7 + ...
		

References

  • Richard P. Stanley, Enumerative Combinatorics, Vol 1, second edition, Example 3.18.3-f, page 323.

Crossrefs

Row sums of triangle A136158.
Cf. A000079, A081295, A009117, A016742, A054879, A092812, A121822. Essentially the same as A004171.

Programs

  • Magma
    [(4^n+0^n)/2: n in [0..30]]; // Vincenzo Librandi, Jul 26 2011
    
  • Magma
    R:=PowerSeriesRing(Rationals(), 25); Coefficients(R!( (1-2*x)/(1-4*x))); // Marius A. Burtea, Jan 20 2020
    
  • Maple
    a:= n-> 2^max(0, (2*n-1)):
    seq(a(n), n=0..30);  # Alois P. Heinz, Jul 20 2017
  • Mathematica
    CoefficientList[Series[(1-2x)/(1-4x),{x,0,40}],x] (* or *)
    Join[{1}, NestList[4 # &, 2, 40]] (* Harvey P. Dale, Apr 22 2011 *)
  • PARI
    a(n)=1<Charles R Greathouse IV, Jul 25 2011
    
  • PARI
    x='x+O('x^100); Vec((1-2*x)/(1-4*x)) \\ Altug Alkan, Dec 21 2015

Formula

G.f.: (1-2*x)/(1-4*x).
a(n) = 4*a(n-1) n > 1, with a(0)=1, a(1)=2.
a(n) = (4^n+0^n)/2 (i.e., 1 followed by 4^n/2, n > 0).
E.g.f.: exp(2*x)*cosh(2*x) = (exp(4*x)+exp(0))/2. - Paul Barry, May 10 2003
a(n) = Sum_{k=0..n} C(2*n, 2*k). - Paul Barry, May 20 2003
a(n) = A001045(2*n+1) - A001045(2*n-1) + 0^n/2. - Paul Barry, Mar 10 2004
a(n) = 2^n*A011782(n); a(n) = gcd(A011782(2n), A011782(2n+1)). - Paul Barry, Jan 12 2005
a(n) = Sum_{k=0..n} Sum_{j=0..n} C(2*(n-k), j)*C(2*k, j). - Paul Barry, Jun 04 2005
a(n) = Sum_{k=0..n} A038763(n,k). - Philippe Deléham, Sep 22 2006
a(n) = Integral_{x=0..4} p(n,x)^2/(Pi*sqrt(x(4-x))) dx, where p(n,x) is the sequence of orthogonal polynomials defined by C(2*n,n): p(n,x) = (2*x-4)*p(n-1,x) - 4*p(n-2,x), with p(0,x)=1, p(1,x)=-2+x. - Paul Barry, Mar 01 2007
a(n) = ((2+sqrt(4))^n + (2-sqrt(4))^n)/2. - Al Hakanson (hawkuu(AT)gmail.com), Nov 22 2008
a(n) = A000079(n) * A011782(n). - Philippe Deléham, Dec 01 2008
a(n) = A004171(n-1) = A028403(n) - A000079(n) for n >= 1. - Jaroslav Krizek, Jul 27 2009
a(n) = Sum_{k=0..n} A201730(n,k)*3^k. - Philippe Deléham, Dec 06 2011
a(n) = Sum_{k=0..n} A134309(n,k)*2^k = Sum_{k=0..n} A055372(n,k). - Philippe Deléham, Feb 04 2012
G.f.: Q(0), where Q(k) = 1 - 2*x/(1 - 2/(2 - 1/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, Apr 29 2013
E.g.f.: 1/2 + exp(4*x)/2 = (Q(0)+1)/2, where Q(k) = 1 + 4*x/(2*k+1 - 2*x*(2*k+1)/(2*x + (k+1)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, Apr 29 2013
a(n) = ceiling( 2^(2n-1) ). - Wesley Ivan Hurt, Jun 30 2013
G.f.: 1 + 2*x/(1 + x)*( 1 + 5*x/(1 + 4*x)*( 1 + 8*x/(1 + 7*x)*( 1 + 11*x/(1 + 10*x)*( 1 + ... )))). - Peter Bala, May 27 2017
Sum_{n>=0} 1/a(n) = 5/3. - Amiram Eldar, Aug 18 2022
Sum_{n>=0} a(n)*x^n/A000680(n) = E(x)^2 where E(x) = Sum_{n>=0} x^n/A000680(n). - Geoffrey Critzer, Apr 21 2023

A002432 Denominators of zeta(2*n)/Pi^(2*n).

Original entry on oeis.org

2, 6, 90, 945, 9450, 93555, 638512875, 18243225, 325641566250, 38979295480125, 1531329465290625, 13447856940643125, 201919571963756521875, 11094481976030578125, 564653660170076273671875, 5660878804669082674070015625, 62490220571022341207266406250
Offset: 0

Views

Author

Keywords

Comments

Also denominators in expansion of Psi(x).
For n >= 1 a(n) is always divisible by 3 (by the von Staudt-Clausen theorem, see A002445).
A comment due to G. Campbell: The original approach taken by Euler was to derive the infinite product for sin(Pi*x)/(Pi*x) equal to (1 - x^2/1^2) (1 - x^2/2^2)(1 - x^2/3^2) ... treating sin(Pi*x)/(Pi*x) as if it were a polynomial. Differentiating the logarithm of both sides and equating coefficients gives all of the zeta function values for 2, 4, 6, 8, .... - M. F. Hasler, Mar 29 2015
Note that 2n+1 divides a(n) for every n. If 2n+1 > 9 is composite, then (2n+1)^2 divides a(n). If 2n+1 is prime, then (2n+1)^2 does not divide a(n). My theorem: for n > 4, (2n+1)^2 divides a(n) if and only if the number 2n+1 is composite. - Thomas Ordowski, Nov 07 2022

Examples

			(zeta(2n)/Pi^(2n), n = 0, 1, 2, 3, ...) = (-1/2, 1/6, 1/90, 1/945, 1/9450, 1/93555, 691/638512875, 2/18243225, 3617/325641566250, ...), i.e.: zeta(0) = -1/2, zeta(2) = Pi^2/6, zeta(4) = Pi^4/90, zeta(6) = Pi^6/945, zeta(8) = Pi^8/9450, zeta(10) = Pi^10/93555, zeta(12) = Pi^12*691/638512875, ...
In Maple, series(Psi(x),x,20) gives -1*x^(-1) + (-gamma) + 1/6*Pi^2*x + (-Zeta(3))*x^2 + 1/90*Pi^4*x^3 + (-Zeta(5))*x^4 + 1/945*Pi^6*x^5 + (-Zeta(7))*x^6 + 1/9450*Pi^8*x^7 + (-Zeta(9))*x^8 + 1/93555*Pi^10*x^9 + ...
a(5) = 93555 = 10!/(2^9 * B(10)) = 3628800/(512*5/66). - _Frank Ellermann_, Apr 03 2020
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 88.
  • A. Fletcher, J. C. P. Miller, L. Rosenhead and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 84.
  • 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

Cf. A046988 (numerators), A006003.

Programs

  • Maple
    seq(denom(Zeta(2*n)/Pi^(2*n)),n=0..24); # Martin Renner, Sep 07 2016
    A002432List := proc(len) series(-x*cot(x)/2, x, 2*len+1):
    seq(denom(coeff(%, x, n)), n=0..2*len-1, 2) end:
    A002432List(17); # Peter Luschny, Jun 07 2020
  • Mathematica
    Table[Denominator[Zeta[2 n]/Pi^(2 n)], {n, 0, 30}] (* Artur Jasinski, Mar 11 2010 *)
    Denominator[Zeta[2*Range[0, 20]]] (* Harvey P. Dale, Sep 05 2013 *)
  • PARI
    a(n)=numerator(bestappr(Pi^(2*n)/zeta(2*n))) \\ Requires sufficient realprecision. The standard value of 38 digits yields erroneous values for n>9. \p99 is more than enough to get the 3 lines of displayed data. - M. F. Hasler, Mar 29 2015
    
  • PARI
    a002432(n) = denominator(polcoeff((1-x*cotan(x))/2,n*2))
    default(seriesprecision, 33); for(i=0,16,print1(a002432(i),",")) \\ Chris Boyd, Dec 21 2015

Formula

Sum_{n>=1} 2/(n^2 + 1) = Pi*coth(Pi)-1. 2*Sum_{k>=1} (-1)^(k + 1)/n^(2*k) = 2/(n^2+1). - Shmuel Spiegel (shmualm(AT)hotmail.com), Aug 13 2001
zeta(2n)/(2i * ( log(1-i)-log(1+i) ))^(2n) = zeta(2n)/(-i*log(-1))^(2n). - Eric Desbiaux, Dec 12 2008
zeta(2n) = Sum_{k >= 1} k^(-2n) = (-1)^(n-1)*B_{2n}*2^(2n-1)*Pi^(2n)/(2n)!.
a(n) = -A046988(n)*A010050(n)*A002445(n)/(A009117(n)*A000367(n))
a(n) = sqrt(denominator(Sum_{i>=1} A000005(i)/i^2n)). - Enrique Pérez Herrero, Jan 19 2012
Sum_{k >= 1} zeta(2k)*x^(2k) = (1-Pi*x*cot(Pi*x))/2. - Chris Boyd, Dec 21 2015
a(n) = denominator([x^(2*n)] -x*cot(x)/2). - Peter Luschny, Jun 07 2020

Extensions

Formula and link from Henry Bottomley, Mar 10 2000
Formula corrected by Bjoern Boettcher, May 15 2003
Corrected and edited by M. F. Hasler, Mar 29 2015
a(0) = 2 prepended by Peter Luschny, Jun 07 2020

A086872 Triangle T(n, k) read by rows; given by [1, 2, 3, 4, 5, 6, ..] DELTA [1, 4, 9, 16, 25, 36, ...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 1, 1, 3, 8, 5, 15, 75, 121, 61, 105, 840, 2478, 3128, 1385, 945, 11025, 51030, 115350, 124921, 50521, 10395, 166320, 1105335, 3859680, 7365633, 7158128, 2702765
Offset: 0

Views

Author

Philippe Deléham, Aug 20 2003, Aug 17 2007

Keywords

Examples

			Triangle begins:
1;
1, 1;
3, 8, 5;
15, 75, 121, 61;
105, 840, 2478, 3128, 1385;
945, 11025, 51030, 115350, 124921, 50521;
10395, 166320, 1105335, 3859680, 7365633, 7158128, 2702765 ; ...
		

Crossrefs

Cf. A000182 (row sums), A000364 (first diagonal), A001147 (first column), A084938, A261065 (2nd column).

Formula

Sum( k>=0, T(n, k)*(-1)^k ) = 0; if n>0.
Sum( k>=0, T(n, k)*(-1/2)^k ) = (1/2)^n.
Sum_{k, 0<=k<=n}T(n,k)*x^(n-k) = (-1)^n*A121822(n), (-1)^n*A092812(n), (-1)^n*A054879(n), A009117(n), A033999(n), A000007(n), A000364(n), A000182(n+1) for x = -6, -5, -4, -3, -2, -1, 0, 1 respectively .

A206022 Riordan array (1, x*exp(arcsinh(-2*x))).

Original entry on oeis.org

1, 0, 1, 0, -2, 1, 0, 2, -4, 1, 0, 0, 8, -6, 1, 0, -2, -8, 18, -8, 1, 0, 0, 0, -32, 32, -10, 1, 0, 4, 8, 30, -80, 50, -12, 1, 0, 0, 0, 0, 128, -160, 72, -14, 1, 0, -10, -16, -28, -112, 350, -280, 98, -16, 1, 0, 0, 0
Offset: 0

Views

Author

Philippe Deléham, Feb 02 2012

Keywords

Comments

Riordan array (1, x*(sqrt(1+4x^2)-2x)); inverse of Riordan array (1, x/sqrt(1-4x)), see A205813.
The g.f. for row sums (1,1,-1,-1,3,1,-9,1,27,13,-81,67,243,...) is (1+2*x^2+x*sqrt(1+4*x^2))/(1+3*x^2).
Triangle T(n,k), read by rows, given by (0, -2, 1, -1, 1, -1, 1, -1, 1, -1, ...) DELTA (1, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938.

Examples

			Triangle begins:
  1
  0,   1
  0,  -2,   1
  0,   2,  -4,   1
  0,   0,   8,  -6,    1,
  0,  -2,  -8,  18,   -8,    1
  0,   0,   0, -32,   32,  -10,     1
  0,   4,   8,  30,  -80,   50,   -12,    1
  0,   0,   0,   0,  128, -160,    72,  -14,    1
  0, -10, -16, -28, -112,  350,  -280,   98,  -16,   1
  0,   0,   0,   0,    0, -512,   768, -448,  128, -18,   1
  0,  28,  40,  54,   96,  420, -1512, 1470, -672, 162, -20, 1
		

Crossrefs

Cf. A104624 (column k=1).

Formula

T(n,n) = 1, T(n+1,n) = -2n = -A005843(n), T(n+2,n) = 2*n^2 = A001105(n), T(n+3,n) = -A130809(n+1), T(2n,n) = A009117(n), T(2n+3,1) = (-1)^n*2*A000108(n).
T(n,k) = T(n-2,k-2) - 4*T(n-2,k-1), for k >= 2.
Showing 1-5 of 5 results.