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

A160469 The left hand column of the triangle A160468.

Original entry on oeis.org

1, 1, 2, 17, 62, 1382, 21844, 929569, 6404582, 443861162, 18888466084, 1936767361654, 58870668456604, 8374643517010684, 689005380505609448, 129848163681107301953, 1736640792209901647222, 418781231495293038913922
Offset: 1

Views

Author

Johannes W. Meijer, May 24 2009

Keywords

Comments

Resembles A002430, the numerators of the Taylor series for tan(x). The first difference occurs at a(12). (Its resemblance to this sequence led to the conjecture A160469(n) = A002430(n)*A089170(n-1).)

Crossrefs

Equals the first left hand column of A160468.
Equals A002430(n)*A089170(n-1).
Equals (A002430(n)/A036279(n))*(A117972(n)/A000265(n)).
Equals A048896(n-1)*A002425(n).
Cf. A156769 (which resembles the denominators of the Taylor series for tan(x)).

Formula

a(n) = A002430(n)*A089170(n-1) with A002430 (n) = numer((-1)^(n-1)*2^(2*n)*(2^(2*n)-1)* bernoulli(2*n)/(2*n)!) and A089170 (n-1) = numer(2*bernoulli(2*n)* (4^n-1)/(2*n))/ numer((4^n-1)*bernoulli(2*n)/(2*n)!) for n = 1, 2, 3, ....

A048896 a(n) = 2^(A000120(n+1) - 1), n >= 0.

Original entry on oeis.org

1, 1, 2, 1, 2, 2, 4, 1, 2, 2, 4, 2, 4, 4, 8, 1, 2, 2, 4, 2, 4, 4, 8, 2, 4, 4, 8, 4, 8, 8, 16, 1, 2, 2, 4, 2, 4, 4, 8, 2, 4, 4, 8, 4, 8, 8, 16, 2, 4, 4, 8, 4, 8, 8, 16, 4, 8, 8, 16, 8, 16, 16, 32, 1, 2, 2, 4, 2, 4, 4, 8, 2, 4, 4, 8, 4, 8, 8, 16, 2, 4, 4, 8, 4, 8, 8, 16, 4, 8, 8, 16, 8, 16, 16, 32, 2, 4, 4
Offset: 0

Views

Author

Keywords

Comments

a(n) = 2^A048881 = 2^{maximal power of 2 dividing the n-th Catalan number (A000108)}. [Comment corrected by N. J. A. Sloane, Apr 30 2018]
Row sums of triangle A128937. - Philippe Deléham, May 02 2007
a(n) = sum of (n+1)-th row terms of triangle A167364. - Gary W. Adamson, Nov 01 2009
a(n), n >= 1: Numerators of Maclaurin series for 1 - ((sin x)/x)^2, A117972(n), n >= 2: Denominators of Maclaurin series for 1 - ((sin x)/x)^2, the correlation function in Montgomery's pair correlation conjecture. - Daniel Forgues, Oct 16 2011
For n > 0: a(n) = A007954(A007931(n)). - Reinhard Zumkeller, Oct 26 2012
a(n) = A261363(2*(n+1), n+1). - Reinhard Zumkeller, Aug 16 2015
From Gus Wiseman, Oct 30 2022: (Start)
Also the number of coarsenings of the (n+1)-th composition in standard order. The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions. See link for sequences related to standard compositions. For example, the a(10) = 4 coarsenings of (2,1,1) are: (2,1,1), (2,2), (3,1), (4).
Also the number of times n+1 appears in A357134. For example, 11 appears at positions 11, 20, 33, and 1024, so a(10) = 4.
(End)

Examples

			From _Omar E. Pol_, Jul 21 2009: (Start)
If written as a triangle:
  1;
  1,2;
  1,2,2,4;
  1,2,2,4,2,4,4,8;
  1,2,2,4,2,4,4,8,2,4,4,8,4,8,8,16;
  1,2,2,4,2,4,4,8,2,4,4,8,4,8,8,16,2,4,4,8,4,8,8,16,4,8,8,16,8,16,16,32;
  ...,
the first half-rows converge to Gould's sequence A001316.
(End)
		

Crossrefs

This is Guy Steele's sequence GS(3, 5) (see A135416).
Equals first right hand column of triangle A160468.
Equals A160469(n+1)/A002425(n+1).
Standard compositions are listed by A066099.
The opposite version (counting refinements) is A080100.
The version for Heinz numbers of partitions is A317141.

Programs

  • Haskell
    a048896 n = a048896_list !! n
    a048896_list = f [1] where f (x:xs) = x : f (xs ++ [x,2*x])
    -- Reinhard Zumkeller, Mar 07 2011
    
  • Haskell
    import Data.List (transpose)
    a048896 = a000079 . a000120
    a048896_list = 1 : concat (transpose
       [zipWith (-) (map (* 2) a048896_list) a048896_list,
        map (* 2) a048896_list])
    -- Reinhard Zumkeller, Jun 16 2013
    
  • Magma
    [Numerator(2^n / Factorial(n+1)): n in [0..100]]; // Vincenzo Librandi, Apr 12 2014
  • Maple
    a := n -> 2^(add(i,i=convert(n+1,base,2))-1): seq(a(n), n=0..97); # Peter Luschny, May 01 2009
  • Mathematica
    NestList[Flatten[#1 /. a_Integer -> {a, 2 a}] &, {1}, 4] // Flatten (* Robert G. Wilson v, Aug 01 2012 *)
    Table[Numerator[2^n / (n + 1)!], {n, 0, 200}] (* Vincenzo Librandi, Apr 12 2014 *)
    Denominator[Table[BernoulliB[2*n] / (Zeta[2*n]/Pi^(2*n)), {n, 1, 100}]] (* Terry D. Grant, May 29 2017 *)
    Table[Denominator[((2 n)!/2^(2 n + 1)) (-1)^n], {n, 1, 100}]/4 (* Terry D. Grant, May 29 2017 *)
    2^IntegerExponent[CatalanNumber[Range[0,100]],2] (* Harvey P. Dale, Apr 30 2018 *)
  • PARI
    a(n)=if(n<1,1,if(n%2,a(n/2-1/2),2*a(n-1)))
    
  • PARI
    a(n) = 1 << (hammingweight(n+1)-1); \\ Kevin Ryde, Feb 19 2022
    

Formula

a(n) = 2^A048881(n).
a(n) = 2^k if 2^k divides A000108(n) but 2^(k+1) does not divide A000108(n).
It appears that a(n) = Sum_{k=0..n} binomial(2*(n+1), k) mod 2. - Christopher Lenard (c.lenard(AT)bendigo.latrobe.edu.au), Aug 20 2001
a(0) = 1; a(2*n) = 2*a(2*n-1); a(2*n+1) = a(n).
a(n) = (1/2) * A001316(n+1). - Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Mar 26 2004
It appears that a(n) = Sum_{k=0..2n} floor(binomial(2n+2, k+1)/2)(-1)^k = 2^n - Sum_{k=0..n+1} floor(binomial(n+1, k)/2). - Paul Barry, Dec 24 2004
a(n) = Sum_{k=0..n} (T(n,k) mod 2) where T = A039598, A053121, A052179, A124575, A126075, A126093. - Philippe Deléham, May 02 2007
a(n) = numerator(b(n)), where sin(x)^2/x = Sum_{n>0} b(n)*(-1)^n x^(2*n-1). - Vladimir Kruchinin, Feb 06 2013
a((2*n+1)*2^p-1) = A001316(n), p >= 0 and n >= 0. - Johannes W. Meijer, Feb 12 2013
a(n) = numerator(2^n / (n+1)!). - Vincenzo Librandi, Apr 12 2014
a(2n) = (2n+1)!/(n!n!)/A001803(n). - Richard Turk, Aug 23 2017
a(2n-1) = (2n-1)!/(n!(n-1)!)/A001790(n). - Richard Turk, Aug 23 2017

Extensions

New definition from N. J. A. Sloane, Mar 01 2008

A160464 The Eta triangle.

Original entry on oeis.org

-1, -11, 2, -114, 29, -2, -3963, 1156, -122, 4, -104745, 32863, -4206, 222, -4, -3926745, 1287813, -184279, 12198, -366, 4, -198491580, 67029582, -10317484, 781981, -30132, 562, -4
Offset: 2

Views

Author

Johannes W. Meijer, May 24 2009

Keywords

Comments

The ES1 matrix coefficients are defined by ES1[2*m-1,n] = 2^(2*m-1) * int(y^(2*m-1)/(cosh(y))^(2*n),y=0..infinity)/(2*m-1)! for m = 1, 2, 3, .. and n = 1, 2, 3 .. .
This definition leads to ES1[2*m-1,n=1] = 2*eta(2*m-1) and the recurrence relation ES1[2*m-1,n] = ((2*n-2)/(2*n-1))*(ES1[2*m-1,n-1] - ES1[2*m-3,n-1]/(n-1)^2) which we used to extend our definition of the ES1 matrix coefficients to m = 0, -1, -2, .. . We discovered that ES1[ -1,n] = 0.5 for n = 1, 2, .. . As usual eta(m) = (1-2^(1-m))*zeta(m) with eta(m) the Dirichlet eta function and zeta(m) the Riemann zeta function.
The coefficients in the columns of the ES1 matrix, for m = 1, 2, 3, .. , and n = 2, 3, 4 .. , can be generated with the polynomials GF(z,n) for which we found the following general expression GF(z;n) = ((-1)^(n-1)*r(n)*CFN1(z,n)*GF(z;n=1) + ETA(z,n))/p(n).
The CFN1(z,n) polynomials depend on the central factorial numbers A008955.
The ETA(z,n) are the Eta polynomials which lead to the Eta triangle.
The zero patterns of the Eta polynomials resemble a UFO. These patterns resemble those of the Zeta, Beta and Lambda polynomials, see A160474, A160480 and A160487.
The first Maple algorithm generates the coefficients of the Eta triangle. The second Maple algorithm generates the ES1[2*m-1,n] coefficients for m= 0, -1, -2, -3, .. .
The M(n) sequence, see the second Maple algorithm, leads to Gould's sequence A001316 and a sequence that resembles the denominators of the Taylor series for tan(x), A156769(n).
Some of our results are conjectures based on numerical evidence, see especially A160466.

Examples

			The first few rows of the triangle ETA(n,m) with n=2,3,.. and m=1,2,... are
  [ -1]
  [ -11, 2]
  [ -114, 29, -2]
  [ -3963, 1156, -122, 4].
The first few ETA(z,n) polynomials are
  ETA(z,n=2) = -1;
  ETA(z,n=3) = -11+2*z^2;
  ETA(z,n=4) = -114 + 29*z^2 - 2*z^4.
The first few CFN1(z,n) polynomials are
  CFN1(z,n=2) = (z^2-1);
  CFN1(z,n=3) = (z^4 - 5*z^2 + 4);
  CFN1(z,n=4) = (z^6 - 14*z^4 + 49*z^2 - 36).
The first few generating functions GF(z;n) are:
  GF(z;n=2) = ((-1)*2*(z^2 - 1)*GF(z;n=1) + (- 1))/3;
  GF(z;n=3) = (4*(z^4 - 5*z^2+4) *GF(z;n=1) + (-11 + 2*z^2))/30;
  GF(z;n=4) = ((-1)*4*(z^6 - 14*z^4 + 49*z^2 - 36)*GF(z;n=1) + (-114 + 29*z^2 - 2*z^4))/315.
		

References

  • Mohammad K. Azarian, Problem 1218, Pi Mu Epsilon Journal, Vol. 13, No. 2, Spring 2010, p. 116. Solution published in Vol. 13, No. 3, Fall 2010, pp. 183-185.

Crossrefs

The r(n) sequence equals A062383 (n>=1).
The p(n) sequence equals A160473(n) (n>=2).
The GCS(n) sequence equals the Geometric Connell sequence A049039(n).
The M(n-1) sequence equals A001316(n-1)/A156769(n) (n>=1).
The q(n) sequence leads to A081729 and the 'gossip sequence' A007456.
The first right hand column equals A053644 (n>=1).
The first left hand column equals A160465.
The row sums equal A160466.
The CFN1(z, n) and the cfn1(n, k) lead to A008955.
Cf. A094665 and A160468.
Cf. the Zeta, Beta and Lambda triangles A160474, A160480 and A160487.
Cf. A162440 (EG1 matrix).

Programs

  • Maple
    nmax:=8; c(2 ):= -1/3: for n from 3 to nmax do c(n) := (2*n-2)*c(n-1)/(2*n-1)-1/((n-1)*(2*n-1)) end do: for n from 2 to nmax do GCS(n-1) := ln(1/(2^(-(2*(n-1)-1-floor(ln(n-1)/ ln(2))))))/ln(2); p(n) := 2^(-GCS(n-1))*(2*n-1)!; ETA(n, 1) := p(n)*c(n); ETA(n, n) := 0 end do: mmax:=nmax: for m from 2 to mmax do for n from m+1 to nmax do q(n) := (1+(-1)^(n-3)*(floor(ln(n-1)/ln(2)) - floor(ln(n-2)/ln(2)))): ETA(n, m) := q(n)*((-1)*ETA(n-1, m-1)+(n-1)^2*ETA(n-1, m)) end do end do: seq(seq(ETA(n,m), m=1..n-1), n=2..nmax);
    # End first program.
    nmax1:=20; m:=1; ES1row:=1-2*m; with (combinat): cfn1 := proc(n, k): sum((-1)^j*stirling1(n+1, n+1-k+j) * stirling1(n+1, n+1-k-j), j=-k..k) end proc: mmax1:=nmax1: for m1 from 1 to mmax1 do M(m1-1) := 2^(2*m1-2)/((2*m1-1)!); ES1[-2*m1+1,1] := 2*(1-2^(1-(1-2*m1)))*(-bernoulli(2*m1)/(2*m1)) od: for n from 2 to nmax1 do for m1 from 1 to mmax1-n+1 do ES1[1-2*m1, n] := (-1)^(n-1)*M(n-1)*sum((-1)^(k+1)*cfn1(n-1,k-1)* ES1[2*k-2*n-2*m1+1, 1], k=1..n) od: od: seq(ES1[1-2*m, n], n=1..nmax1-m+1);
    # End second program.

Formula

We discovered an interesting relation between the Eta triangle coefficients ETA(n,m) = q(n)*((-1)*ETA(n-1,m-1)+(n-1)^2*ETA(n-1,m)), for n = 3, 4, ... and m = 2, 3, ... , with
q(n) = 1 + (-1)^(n-3)*(floor(log(n-1)/log(2)) - floor(log(n-2)/log(2))) for n = 3, 4, ....
See A160465 for ETA(n,m=1) and furthermore ETA(n,n) = 0 for n = 2, 3, ....
The generating functions GF(z;n) of the coefficients in the matrix columns are defined by
GF(z;n) = sum_{m>=1} ES1[2*m-1,n] * z^(2*m-2), with n = 1, 2, 3, .... This leads to
GF(z;n=1) = (2*log(2) - Psi(z) - Psi(-z) + Psi(1/2*z) + Psi(-1/2*z)); Psi(z) is the digamma-function.
GF(z;n) = ((2*n-2)/(2*n-1)-2*z^2/((n-1)*(2*n-1)))*GF(z;n-1)-1/((n-1)*(2*n-1)).
We found for GF(z;n), for n = 2, 3, ..., the following general expression:
GF(z;n) = ((-1)^(n-1)*r(n)*CFN1(z,n)*GF(z;n=1) + ETA(z,n) )/p(n) with
r(n) = 2^floor(log(n-1)/log(2)+1) and
p(n) = 2^(-GCS(n))*(2*n-1)! with
GCS(n) = log(1/(2^(-(2*(n-1)-1-floor(log(n-1)/ log(2))))))/log(2).

A094665 Another version of triangular array in A083061: triangle T(n,k), 0<=k<=n, read by rows; given by [0, 1, 3, 6, 10, 15, 21, 28, ...] DELTA [1, 2, 3, 4, 5, 6, 7, 8, ...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 0, 1, 0, 1, 3, 0, 4, 15, 15, 0, 34, 147, 210, 105, 0, 496, 2370, 4095, 3150, 945, 0, 11056, 56958, 111705, 107415, 51975, 10395, 0, 349504, 1911000, 4114110, 4579575, 2837835, 945945, 135135, 0, 14873104, 85389132, 197722980, 244909665, 178378200, 77567490, 18918900, 2027025
Offset: 0

Views

Author

Philippe Deléham, Jun 07 2004, Jun 12 2007

Keywords

Comments

Define polynomials P(n,x) = x(2x+1)P(n-1,x+1) - 2x^2P(n-1,x), P(0,x) = 1. Sequence gives triangle read by rows, defined by P(n,x) = Sum_{k = 0..n} T(n,k)*x^k. - Philippe Deléham, Jun 20 2004
From Johannes W. Meijer, May 24 2009: (Start)
In A160464 we defined the coefficients of the ES1 matrix by ES1[2*m-1,n=1] = 2*eta(2*m-1) and the recurrence relation ES1[2*m-1,n] = ((2*n-2)/(2*n-1))*(ES1[2*m-1,n-1] - ES1[2*m-3,n-1]/(n-1)^2) for m the positive and negative integers and n >= 1. As usual eta(m) = (1-2^(1-m))*zeta(m) with eta(m) the Dirichlet eta function and zeta(m) the Riemann zeta function. It is well-known that ES1[1-2*m,n=1] = (4^m-1)*(-bernoulli(2*m))/m for m >= 1. and together with the recurrence relation this leads to ES1[-1,n] = 0.5 for n >= 1.
We discovered that the n-th term of the row coefficients ES1[1-2*m,n] for m >= 1, can be generated with the rather simple polynomials RES1(1-2*m,n) = (-1)^(m+1)*ECGP(1-2*m, n)/2^m. This discovery was enabled by the recurrence relation for the RES1(1-2*m,n) which we derived from the recurrence relation for the ES1[2*m-1,n] coefficients and the fact that RES1(-1,n) = 0.5. The coefficients of the ECGP(1-2*m,n) polynomials led to this triangle and subsequently to triangle A083061. (End)
From David Callan, Jan 03 2011: (Start)
T(n,k) is the number of increasing 0-2 trees (A002105) on 2n edges in which the minimal path from the root has length k.
Proof. The number a(n,k) of such trees satisfies the recurrence a(0,0)=1, a(1,1)=1 and, counting by size of the subtree rooted at the smaller child of the root,
a(n,k) = Sum_{j=1..n-1} C(2n-1,j)*a(j,k-1)*a(n-1-j)
for 2<=k<=n, where a(n) = Sum_{k>=0} a(n,k) is the reduced tangent number A002105 (indexed from 0). The recurrence translates into the differential equation
F_x(x,y) = y*F(x,y)*G(x)
for the GF F(x,y) = Sum_{n,k>=0} a(n,k)x^(2n)/(2n)!*y^k, where G(x):=Sum_{n>=0} a(n)x^(2n+1)/(2n+1)! is known to be sqrt(2)*tan(x/sqrt(2)). The differential equation has solution F(x,y) = sec(x/sqrt(2))^(2y). (End)

Examples

			Triangle begins:
.1;
.0, 1;
.0, 1, 3;
.0, 4, 15, 15;
.0, 34, 147, 210, 105;
.0, 496, 2370, 4095, 3150, 945;
.0, 11056, 56958, 111705, 107415, 51975, 10395;
.0, 349504, 1911000, 4114110, 4579575, 2837835, 945945, 135135;
From _Johannes W. Meijer_, May 24 2009: (Start)
The first few ECGP(1-2*m,n) polynomials are: ECGP(-1,n) = 1; ECGP(-3,n) = n; ECGP(-5,n) = n + 3*n^2; ECGP(-7,n) = 4*n + 15*n^2+ 15*n^3 .
The first few RES1(1-2*m,n) are: RES1(-1,n) = (1/2)*(1); RES1(-3,n) = (-1/4)*(n); RES1(-5,n) = (1/8)*(n+3*n^2); RES1(-7,n) = (-1/16)*(4*n+15*n^2+15*n^3).
(End)
		

Crossrefs

From Johannes W. Meijer, May 24 2009 and Jun 27 2009: (Start)
A001147, A001880, A160470, A160471 and A160472 are the first five right hand columns.
Appears in A162005, A162006 and A162007.
(End)

Programs

  • Maple
    nmax:=7; imax := nmax: T1(0, x) := 1: T1(0, x+1) := 1: for i from 1 to imax do T1(i, x) := expand((2*x+1) * (x+1) * T1(i-1, x+1) - 2 * x^2 * T1(i-1, x)): dx:=degree(T1(i, x)): for k from 0 to dx do c(k) := coeff(T1(i, x), x, k) od: T1(i, x+1) := sum(c(j1)*(x+1)^(j1), j1=0..dx) od: for i from 0 to imax do for j from 0 to i do A083061(i, j) := coeff(T1(i, x), x, j) od: od: for n from 0 to nmax do for k from 0 to n do T(n+1, k+1) := A083061(n, k) od: od: T(0, 0):=1: for n from 1 to nmax do T(n, 0):=0 od: seq(seq(T(n, k), k=0..n), n=0..nmax);
    # Johannes W. Meijer, Jun 27 2009, revised Sep 23 2012
  • Mathematica
    nmax = 8;
    T[n_, k_] := SeriesCoefficient[Sec[x/Sqrt[2]]^(2y), {x, 0, 2n}, {y, 0, k}]* (2n)!;
    Table[T[n, k], {n, 0, nmax}, {k, 0, n}] // Flatten (* Jean-François Alcover, Aug 10 2018 *)

Formula

Sum_{k = 0..n} T(n, k) = A002105(n+1).
Sum_{k = 0..n} T(n, k)*2^(n-k) = A000364(n); Euler numbers.
Sum_{k = 0..n} T(n, k)*(-2)^(n-k) = 1.
RES1(1-2*m,n) = n^2*RES1(3-2*m,n)-n*(2*n+1)*RES1(3-2*m,n+1)/2 for m >= 2, with RES1(-1,n) = 0.5 for n >= 1. - Johannes W. Meijer, May 24 2009
G.f.: Sum_{n,k>=0} T(n,k)x^n/n!*y^k = sec(x/sqrt(2))^(2y).

Extensions

Term corrected by Johannes W. Meijer, Sep 23 2012

A083061 Triangle of coefficients of a companion polynomial to the Gandhi polynomial.

Original entry on oeis.org

1, 1, 3, 4, 15, 15, 34, 147, 210, 105, 496, 2370, 4095, 3150, 945, 11056, 56958, 111705, 107415, 51975, 10395, 349504, 1911000, 4114110, 4579575, 2837835, 945945, 135135, 14873104, 85389132, 197722980, 244909665, 178378200, 77567490
Offset: 0

Views

Author

Hans J. H. Tuenter, Apr 19 2003

Keywords

Comments

This polynomial arises in the setting of a symmetric Bernoulli random walk and occurs in an expression for the even moments of the absolute distance from the origin after an even number of timesteps. The Gandhi polynomial, sequence A036970, occurs in an expression for the odd moments.
When formatted as a square array, first row is A002105, first column is A001147, second column is A001880.
Another version of the triangle T(n,k), 0<=k<=n, read by rows; given by [0, 1, 3, 6, 10, 15, 21, 28, ...] DELTA [1, 2, 3, 4, 5, 6, 7, 8, 9, ...] = 1; 0, 1; 0, 1, 3; 0, 4, 15, 15; 0, 34, 147, 210, 105; ... where DELTA is the operator defined in A084938. - Philippe Deléham, Jun 07 2004
In A160464 we defined the coefficients of the ES1 matrix. Our discovery that the n-th term of the row coefficients ES1[1-2*m,n] for m>=1, can be generated with rather simple polynomials led to triangle A094665 and subsequently to this one. - Johannes W. Meijer, May 24 2009
Related to polynomials defined in A160485 by a shift of +-1/2 and scaling by a power of 2. - Richard P. Brent, Jul 15 2014

Examples

			Triangle starts (with an additional first column 1,0,0,...):
[1]
[0,      1]
[0,      1,       3]
[0,      4,      15,      15]
[0,     34,     147,     210,     105]
[0,    496,    2370,    4095,    3150,     945]
[0,  11056,   56958,  111705,  107415,   51975,  10395]
[0, 349504, 1911000, 4114110, 4579575, 2837835, 945945, 135135]
		

Crossrefs

From Johannes W. Meijer, May 24 2009 and Jun 27 2009: (Start)
A002105 equals the row sums (n>=2) and the first left hand column (n>=1).
A001147, A001880, A160470, A160471 and A160472 are the first five right hand columns.
Appears in A162005, A162006 and A162007.
(End)

Programs

  • Maple
    imax := 6;
    T1(0, x) := 1:
    T1(0, x+1) := 1:
    for i from 1 to imax do
        T1(i, x) := expand((2*x+1) * (x+1) * T1(i-1, x+1) - 2*x^2*T1(i-1, x)):
        dx := degree(T1(i, x)):
        for k from 0 to dx do
            c(k) := coeff(T1(i, x), x, k)
        od:
        T1(i, x+1) := sum(c(j1)*(x+1)^(j1), j1 = 0..dx):
    od:
    for i from 0 to imax do
        for j from 0 to i do
            a(i, j) := coeff(T1(i, x), x, j)
        od:
    od:
    seq(seq(a(i, j), j = 0..i), i = 0..imax);
    # Johannes W. Meijer, Jun 27 2009, revised Sep 23 2012
  • Mathematica
    b[0, 0] = 1;
    b[n_, k_] := b[n, k] = Sum[2^j*(Binomial[k + j, 1 + j] + Binomial[k + j + 1, 1 + j])*b[n - 1, k - 1 + j], {j, Max[0, 1 - k], n - k}];
    a[0, 0] = 1;
    a[n_, k_] := b[n, k]/2^(n - k);
    Table[a[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jun 19 2018, after Philippe Deléham *)
  • Sage
    # uses[fr2_row from A088874]
    A083061_row = lambda n: [(-1)^(n-k)*m*2^(-n+k) for k,m in enumerate(fr2_row(n))]
    for n in (0..7): print(A083061_row(n)) # Peter Luschny, Sep 19 2017

Formula

Let T(i, x)=(2x+1)(x+1)T(i-1, x+1)-2x^2T(i-1, x), T(0, x)=1; so that T(1, x)=1+3x; T(2, x)=4+15x+15x^2; T(3, x)=34+147x+210x^2+105x^3, etc. Then the (i, j)-th entry in the table is the coefficient of x^j in T(i, x).
a(n, k)*2^(n-k) = A085734(n, k). - Philippe Deléham, Feb 27 2005

A117972 Numerator of zeta'(-2n), n >= 0.

Original entry on oeis.org

1, -1, 3, -45, 315, -14175, 467775, -42567525, 638512875, -97692469875, 9280784638125, -2143861251406875, 147926426347074375, -48076088562799171875, 9086380738369043484375, -3952575621190533915703125
Offset: 0

Views

Author

Eric W. Weisstein, Apr 06 2006

Keywords

Comments

In A160464 the coefficients of the ES1 matrix are defined. This matrix led to the discovery that the successive differences of the ES1[1-2*m,n] coefficients for m = 1, 2, 3, ..., are equal to the values of zeta'(-2n), see also A094665 and A160468. - Johannes W. Meijer, May 24 2009
A048896(n), n >= 1: Numerators of Maclaurin series for 1 - ((sin x)/x)^2,
a(n), n >= 2: Denominators of Maclaurin series for 1 - ((sin x)/x)^2, the correlation function in Montgomery's pair correlation conjecture. - Daniel Forgues, Oct 16 2011
From Andrey Zabolotskiy, Sep 23 2021: (Start)
zeta'(-2n), which is mentioned in the Name, is irrational. For n > 0, a(n) is the numerator of the rational fraction g(n) = Pi^(2n)*zeta'(-2n)/zeta(2n+1). The denominator is 4*A048896(n-1). g(n) = f(n) for n > 0, where f(n) is given in the Formula section. Also, f(n) = Bernoulli(2n)/z(n)/4 (see Formula section) for all n.
For n = 0, zeta'(0) = -log(2Pi)/2, g(0) can be set to 0 because of the infinite denominator. However, a(0) is set to 1 because it is the numerator of f(0).
It seems that -4*f(n)*alpha_n = A000182(n), where alpha_n = A191657(n, p(n)) / A191658(n, p(n)) [where p(n) = A000041(n)] is the n-th "elementary coefficient" from the paper by Izaurieta et al. (End)

Examples

			-1/4, 3/4, -45/8, 315/4, -14175/8, 467775/8, -42567525/16, ...
-zeta(3)/(4*Pi^2), (3*zeta(5))/(4*Pi^4), (-45*zeta(7))/(8*Pi^6), (315*zeta(9))/(4*Pi^8), (-14175*zeta(11))/(8*Pi^10), ...
		

Crossrefs

From Johannes W. Meijer, May 24 2009: (Start)
Absolute values equal row sums of A160468. (End)

Programs

  • Maple
    # Without rational arithmetic
    a := n -> (-1)^n*(2*n)!*2^(add(i,i=convert(n,base,2))-2*n);
    # Peter Luschny, May 02 2009
  • Mathematica
    Table[Numerator[(2 n)!/2^(2 n + 1) (-1)^n], {n, 0, 30}]
  • Maxima
    L:taylor(1/x*sin(sqrt(x))^2,x,0,15); makelist(denom(coeff(L,x,n))*(-1)^(n+1),n,0,15); /* Vladimir Kruchinin, May 30 2011 */

Formula

a(n) = numerator(f(n)) where f(n) = (2*n)!/2^(2*n + 1)(-1)^n, from the Mathematica code.
From Terry D. Grant, May 28 2017: (Start)
|a(n)| = A049606(2n).
a(n) = -numerator(Bernoulli(2n)/z(n)) where Bernoulli(2n) = A000367(n) / A002445(n) and z(n) = A046988(n) / A002432(n) for n > 0. (End) [Corrected by Andrey Zabolotskiy, Sep 23 2021]

Extensions

First term added, offset changed and edited by Johannes W. Meijer, May 15 2009

A117973 a(n) = 2^(wt(n)+1), where wt() = A000120().

Original entry on oeis.org

2, 4, 4, 8, 4, 8, 8, 16, 4, 8, 8, 16, 8, 16, 16, 32, 4, 8, 8, 16, 8, 16, 16, 32, 8, 16, 16, 32, 16, 32, 32, 64, 4, 8, 8, 16, 8, 16, 16, 32, 8, 16, 16, 32, 16, 32, 32, 64, 8, 16, 16, 32, 16, 32, 32, 64, 16, 32, 32, 64, 32, 64, 64, 128, 4, 8, 8, 16, 8, 16, 16, 32, 8, 16, 16, 32, 16, 32
Offset: 0

Views

Author

Eric W. Weisstein, Apr 06 2006

Keywords

Comments

Denominator of Zeta'(-2n).
If Gould's sequence A001316 is written as a triangle, this is what the rows converge to. In other words, let S_0 = [2], and construct S_{n+1} by following S_n with 2*S_n. Then this is S_{oo}. - N. J. A. Sloane, May 30 2009
In A160464 the coefficients of the ES1 matrix are defined. This matrix led to the discovery that the successive differences of the ES1[1-2*m,n] coefficients for m = 1, 2, 3, ..., are equal to the values of Zeta'(-2n), see also A094665 and A160468. - Johannes W. Meijer, May 24 2009

Examples

			-zeta(3)/(4*Pi^2), (3*zeta(5))/(4*Pi^4), (-45*zeta(7))/(8*Pi^6), (315*zeta(9))/(4*Pi^8), (-14175*zeta(11))/(8*Pi^10), ...
		

Crossrefs

Programs

  • Maple
    S := [2]; S := [op(S), op(2*S)]; # repeat ad infinitum! - N. J. A. Sloane, May 30 2009
    a := n -> 2^(add(i,i=convert(n,base,2))+1); # Peter Luschny, May 02 2009
  • Mathematica
    Denominator[(2*n)!/2^(2*n + 1)]

Formula

For n>=0, a(n) = 2 * A001316(n). - N. J. A. Sloane, May 30 2009
For n>0, a(n) = 4 * A048896(n). - Peter Luschny, May 02 2009
a(0) = 2; for n>0, write n = 2^i + j where 0 <= j < 2^i; then a(n) = 2*a(j).
a((2*n+1)*2^p-1) = 2^(p+1) * A001316(n), p >= 0. - Johannes W. Meijer, Jan 28 2013

Extensions

Entry revised by N. J. A. Sloane, May 30 2009
Showing 1-7 of 7 results.