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.

Previous Showing 11-16 of 16 results.

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

A160468 Triangle of polynomial coefficients related to the o.g.f.s of the RES1 polynomials.

Original entry on oeis.org

1, 1, 2, 1, 17, 26, 2, 62, 192, 60, 1, 1382, 7192, 5097, 502, 2, 21844, 171511, 217186, 55196, 2036, 2, 929569, 10262046, 20376780, 9893440, 1089330, 16356, 4, 6404582, 94582204, 271154544, 215114420, 48673180, 2567568, 16376, 1
Offset: 1

Views

Author

Johannes W. Meijer, May 24 2009

Keywords

Comments

In A160464 we defined the ES1 matrix by ES1[2*m-1,n=1] and in A094665 it was shown that the n-th term of the coefficients of matrix row ES1[1-2*m,n] for m >= 1 can be generated with the RES1(1-2*m,n) polynomials.
We define the o.g.f.s. of these polynomials by GFRES1(z,1-2*m) = sum(RES1(1-2*m,n)*z^(n-1), n=1..infinity) for m >= 1. The general expression of the o.g.f.s. is GFRES1(z,1-2*m) = (-1)*RE(z,1-2*m)/(2*p(m-1)*(z-1)^(m)). The p(m-1), m >= 1, sequence is Gould's sequence A001316.
The coefficients of the RE(z,1-2*m) polynomials lead to the triangle given above.
The E(z,n) = numer(sum((-1)^(n+1)*k^n*z^(k-1), k=1..infinity)) polynomials with n >= 1, see the Maple algorithm, lead to the Eulerian numbers A008292.
Some of our results are conjectures based on numerical evidence.

Examples

			The first few rows are:
[1]
[1]
[2, 1]
[17, 26, 2]
[62, 192, 60, 1]
The first few polynomials RE(z,m) are:
RE(z,-1) = 1
RE(z,-3) = 1
RE(z,-5) = 2+z
RE(z,-7) = 17+26*z+2*z^2
The first few GFRES1(z,m) are:
GFRES1(z,-1) = -(1/1)*(1)/(2*(z-1)^1)
GFRES1(z,-3) = -(1/2)*(1)/(2*(z-1)^2)
GFRES1(z,-5) = -(1/2)*(2+z)/(2*(z-1)^3)
GFRES1(z,-7) = -(1/4)*(17+26*z+2*z^2)/(2*(z-1)^4)
		

Crossrefs

For the Eulerian numbers E(n, k) see A008292.
The p(n) sequence equals Gould's sequence A001316.
The first right hand column of the triangle equals A048896.
The first left hand column equals A160469.
The row sums equal the absolute values of A117972.

Programs

  • Maple
    nmax := 8; mmax := nmax: T(0, x) := 1: for i from 1 to nmax do dgr := degree(T(i-1, x), x): for na from 0 to dgr do c(na) := coeff(T(i-1, x), x, na) od: T(i-1, x+1) := 0: for nb from 0 to dgr do T(i-1, x+1) := T(i-1, x+1) + c(nb)*(x+1)^nb od: for nc from 0 to dgr do ECGP(i-1, nc+1) := coeff(T(i-1, x), x, nc) od: T(i, x) := expand((2*x+1)*(x+1)*T(i-1, x+1) - 2*x^2*T(i-1, x)) od: dgr := degree (T(nmax, x), x): kmax := nmax: for k from 1 to kmax do p := k: for m from 1 to k do E(m, k) := sum((-1)^(m-q)*(q^k)*binomial(k+1, m-q), q=1..m) od: fx(p) := (-1)^(p+1) * (sum(E(r, k)*z^(k-r), r=1..k))/(z-1)^(p+1): GF(-(2*p+1)) := sort(simplify(((-1)^p* 1/2^(p+1)) * sum(ECGP(k-1, k-s)*fx(k-s), s=0..k-1)), ascending): NUMGF(-(2*p+1)) := -numer(GF(-(2*p+1))): for n from 1 to mmax+1 do A(k+1, n) := coeff(NUMGF(-(2*p+1)), z, n-1) od: od: for m from 2 to mmax do A(1, m) := 0 od: A(1, 1) := 1: FT(1) := 1: for n from 1 to nmax do for m from 1 to n do FT((n)*(n-1)/2+m+1) := A(n+1, m) end do end do: a := n-> FT(n): seq(a(n), n = 1..(nmax+1)*(nmax)/2+1);
  • Mathematica
    T[ n_, k_] := Coefficient[a[2 n]/2^IntegerExponent[(2 n)!, 2], x, n + k];
    a[0] = a[1] = 1; a[ m_] := a[m] = With[{n = m - 1}, x Sum[ a[k] a[n - k] Binomial[n, k], {k, 0, n}]]; Join[{1}, Flatten@Table[T[n, k], {n, 1, 8}, {k, 0, n - 1}]] (* Michael Somos, Apr 22 2020 *)

Extensions

Edited by Johannes W. Meijer, Sep 23 2012

A162007 Third left hand column of the EG1 triangle A162005.

Original entry on oeis.org

1, 270, 36096, 4766048, 704357760, 120536980224, 24060789342208, 5590122715250688, 1503080384197754880, 464520829174515630080, 163839204411117787938816, 65500849343294249018327040
Offset: 3

Views

Author

Johannes W. Meijer, Jun 27 2009

Keywords

Crossrefs

Third left hand column of the EG1 triangle A162005.
Other left hand columns are A000182 and A162006.
Related to A094665, A083061 and A156919.
A000079, A036289 and A100381 appear in the a(n, 3) formula.
A001789, A003472, A054849, A002409, A054851, A140325 and A140354 (scaled by 2^(m-1)) appear one by one in the a(n, m) formulas for m= 4 and higher .

Programs

  • Maple
    nmax := 14; mmax := nmax: 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 A094665(n+1, k+1) := A083061(n, k) od: od: A094665(0, 0) := 1: for n from 1 to nmax do A094665(n, 0) := 0 od: for m from 1 to mmax do A156919(0, m) := 0 end do: for n from 0 to nmax do A156919(n, 0) := 2^n end do: for n from 1 to nmax do for m from 1 to mmax do A156919(n, m) := (2*m+2)*A156919(n-1, m) + (2*n-2*m+1) * A156919(n-1, m-1) end do end do: m:=3; for n from m to nmax do a(n, m) := sum((-1)^(m-p1-1)*sum(2^(n-q-1)*binomial(n-q-1, m-p1-1) * A094665(n-1, q) * A156919(q, p1), q=1..n-m+p1), p1=0..m-1) od: seq(a(n, m), n = m..nmax);
    # Maple program edited by Johannes W. Meijer, Sep 25 2012

Formula

a(n) = sum((-1)^(m-p-1)*sum(2^(n-q-1)*binomial(n-q-1,m-p-1)*A094665(n-1,q)* A156919(q,p),q=1..n-m+p), p=0..m-1) with m = 3.

A162006 Second left hand column of the EG1 triangle A162005.

Original entry on oeis.org

1, 28, 1032, 52736, 3646208, 330545664, 38188155904, 5488365862912, 961530104709120, 201865242068910080, 50052995352723193856, 14476381898608390176768, 4831399425299156001882112
Offset: 2

Views

Author

Johannes W. Meijer, Jun 27 2009

Keywords

Crossrefs

Second left hand column of the EG1 triangle A162005.
Other left hand columns are A000182 and A162007.
Related to A094665, A083061 and A156919.
A000079 and A036289 appear in the Maple program.

Programs

  • Maple
    nmax := 14; mmax := nmax: 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 A094665(n+1, k+1) := A083061(n, k) od: od: A094665(0, 0) := 1: for n from 1 to nmax do A094665(n, 0) := 0 od: for m from 1 to mmax do A156919(0, m) := 0 end do: for n from 0 to nmax do A156919(n, 0) := 2^n end do: for n from 1 to nmax do for m from 1 to mmax do A156919(n, m) := (2*m+2)*A156919(n-1, m) + (2*n-2*m+1) * A156919(n-1, m-1) end do end do: m:=2; for n from m to nmax do a(n, m) := sum((-1)^(m-p1-1)*sum(2^(n-q-1)*binomial(n-q-1, m-p1-1) * A094665(n-1, q) * A156919(q, p1), q=1..n-m+p1), p1=0..m-1) od: seq(a(n, m), n = m..nmax);
    # Maple program edited by Johannes W. Meijer, Sep 25 2012

Formula

a(n) = sum((-1)^(m-p-1)*sum(2^(n-q-1)*binomial(n-q-1,m-p-1)*A094665(n-1,q)*A156919(q,p),q=1..n-m+p), p=0..m-1) with m = 2.

A094346 Another version of triangular array in A036970: triangle T(n,k), 0<=k<=n, read by rows; given by [0, 1, 2, 4, 6, 9, 12, 16, 20, 25, 30, ...] DELTA [1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, ...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 3, 8, 6, 0, 17, 54, 60, 24, 0, 155, 556, 762, 480, 120, 0, 2073, 8146, 12840, 10248, 4200, 720, 0, 38227, 161424, 282078, 263040, 139440, 40320, 5040, 0, 929569, 4163438, 7886580, 8240952, 5170800, 1965600, 423360, 40320
Offset: 0

Views

Author

Philippe Deléham, Jun 08 2004, Jun 13 2007

Keywords

Comments

Diagonals: A000007, A001469, A005440; A000182, A005990. Row sums: A001469.

Examples

			Triangle begins:
1;
0, 1;
0, 1, 2;
0, 3, 8, 6;
0, 17, 54, 60, 24;
0, 155, 556, 762, 480, 120;
0, 2073, 8146, 12840, 10248, 4200, 720;
0, 38227, 161424, 282078, 263040, 139440, 40320, 5040;
0, 929569, 4163438, 7886580, 8240952, 5170800, 1965600, 423360, 40320; ...
		

Crossrefs

Programs

  • Mathematica
    G[_, 1] = 1;
    G[x_, n_] := G[x, n] = (x+1)^2 G[x+1, n-1] - x^2 G[x, n-1] // Expand;
    row[0] = {1};
    row[n_] := CoefficientList[x G[x, n], x];
    Table[row[n], {n, 0, 8}] // Flatten (* Jean-François Alcover, Aug 17 2018 *)
  • PARI
    {T(n, k) = local( A = x); if( k<0 || k>n, 0, for( j = 1, n, A = x^2 * ( subst(A, x, x+1) - A)); polcoeff( A, k+1))} /* Michael Somos, Apr 10 2011 */

Formula

For n>=1, Sum_{k =1..n} T(n, k)*x^(k-1) = G(x, n), n-th Gandhi polynomial; the Gandhi polynomials are defined by G(x, n) = (x+1)^2*G(x+1, n-1) - x^2*G(x, n-1), G(x, 1) = 1. Sum_{k =0..n} T(n, k)*2^(2n-k) = A000182(n+1), tangent numbers. Sum_{k =0..n} T(n, k) = A001469(n+1), Genocchi numbers of first kind.
Sum_{k = 0..n} T(n, k)*2^(n-k) = A002105(n+1). - Philippe Deléham, Jun 10 2004
Previous Showing 11-16 of 16 results.