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

A153641 Nonzero coefficients of the Swiss-Knife polynomials for the computation of Euler, tangent, and Bernoulli numbers (triangle read by rows).

Original entry on oeis.org

1, 1, 1, -1, 1, -3, 1, -6, 5, 1, -10, 25, 1, -15, 75, -61, 1, -21, 175, -427, 1, -28, 350, -1708, 1385, 1, -36, 630, -5124, 12465, 1, -45, 1050, -12810, 62325, -50521, 1, -55, 1650, -28182, 228525, -555731, 1, -66, 2475, -56364, 685575, -3334386, 2702765, 1
Offset: 0

Views

Author

Peter Luschny, Dec 29 2008

Keywords

Comments

In the following the expression [n odd] is 1 if n is odd, 0 otherwise.
(+) W_n(0) = E_n are the Euler (or secant) numbers A122045.
(+) W_n(1) = T_n are the signed tangent numbers, see A009006.
(+) W_{n-1}(1) n / (4^n - 2^n) = B_n gives for n > 1 the Bernoulli number A027641/A027642.
(+) W_n(-1) 2^{-n}(n+1) = G_n the Genocchi number A036968.
(+) W_n(1/2) 2^{n} are the signed generalized Euler (Springer) number, see A001586.
(+) | W_n([n odd]) | the number of alternating permutations A000111.
(+) | W_n([n odd]) / n! | for 0<=n the Euler zeta number A099612/A099617 (see Wikipedia on Bernoulli number). - Peter Luschny, Dec 29 2008
The diagonals in the full triangle (with zero coefficients) of the polynomials have the general form E(k)*binomial(n+k,k) (k>=0 fixed, n=0,1,...) where E(n) are the Euler numbers in the enumeration A122045. For k=2 we find the triangular numbers A000217 and for k=4 A154286. - Peter Luschny, Jan 06 2009
From Peter Bala, Jun 10 2009: (Start)
The Swiss-Knife polynomials W_n(x) may be expressed in terms of the Bernoulli polynomials B(n,x) as
... W_n(x) = 4^(n+1)/(2*n+2)*[B(n+1,(x+3)/4) - B(n+1,(x+1)/4)].
The Swiss-Knife polynomials are, apart from a multiplying factor, examples of generalized Bernoulli polynomials.
Let X be the Dirichlet character modulus 4 defined by X(4*n+1) = 1, X(4*n+3) = -1 and X(2*n) = 0. The generalized Bernoulli polynomials B(X;n,x), n = 1,2,..., associated with the character X are defined by means of the generating function
... t*exp(x*t)*(exp(t)-exp(3*t))/(exp(4*t)-1) = sum {n = 1..inf} B(X;n,x)*t^n/n!.
The first few values are B(X;1,x) = -1/2, B(X;2,x) = -x, B(X,3,x) = -3/2*(x^2-1) and B(X;4,x) = -2*(x^3-3*x).
In general, W_n(x) = -2/(n+1)*B(X;n+1,x).
For the theory of generalized Bernoulli polynomials associated to a periodic arithmetical function see [Cohen, Section 9.4].
The generalized Bernoulli polynomials may be used to evaluate twisted sums of k-th powers. For the present case the result is
sum{n = 0..4*N-1} X(n)*n^k = 1^k - 3^k + 5^k - 7^k + ... - (4*N-1)^k
= [B(X;k+1,4*N) - B(X;k+1,0)]/(k+1) = [W_k(0) - W_k(4*N)]/2.
For the proof apply [Cohen, Corollary 9.4.17 with m = 4 and x = 0].
The generalized Bernoulli polynomials and the Swiss-Knife polynomials are also related to infinite sums of powers through their Fourier series - see the formula section below. For a table of the coefficients of generalized Bernoulli polynomials attached to a Dirichlet character modulus 8 see A151751.
(End)
The Swiss-Knife polynomials provide a general formula for alternating sums of powers similar to the formula which are provided by the Bernoulli polynomials for non-alternating sums of powers (see the Luschny link). Sequences covered by this formula include A001057, A062393, A062392, A011934, A144129, A077221, A137501, A046092. - Peter Luschny, Jul 12 2009
The greatest common divisor of the nonzero coefficients of the decapitated Swiss-Knife polynomials is exp(Lambda(n)), where Lambda(n) is the von Mangoldt function for odd primes, symbolically:
gcd(coeffs(SKP_{n}(x) - x^n)) = A155457(n) (n>1). - Peter Luschny, Dec 16 2009
Another version is at A119879. - Philippe Deléham, Oct 26 2013

Examples

			1
x
x^2  -1
x^3  -3x
x^4  -6x^2   +5
x^5 -10x^3  +25x
x^6 -15x^4  +75x^2  -61
x^7 -21x^5 +175x^3 -427x
		

References

  • H. Cohen, Number Theory - Volume II: Analytic and Modern Tools, Graduate Texts in Mathematics. Springer-Verlag. [From Peter Bala, Jun 10 2009]

Crossrefs

W_n(k), k=0,1,...
W_0: 1, 1, 1, 1, 1, 1, ........ A000012
W_1: 0, 1, 2, 3, 4, 5, ........ A001477
W_2: -1, 0, 3, 8, 15, 24, ........ A067998
W_3: 0, -2, 2, 18, 52, 110, ........ A121670
W_4: 5, 0, -3, 32, 165, 480, ........
W_n(k), n=0,1,...
k=0: 1, 0, -1, 0, 5, 0, -61, ... A122045
k=1: 1, 1, 0, -2, 0, 16, 0, ... A155585
k=2: 1, 2, 3, 2, -3, 2, 63, ... A119880
k=3: 1, 3, 8, 18, 32, 48, 128, ... A119881
k=4: 1, 4, 15, 52, 165, 484, ........ [Peter Luschny, Jul 07 2009]

Programs

  • Maple
    w := proc(n,x) local v,k,pow,chen; pow := (a,b) -> if a = 0 and b = 0 then 1 else a^b fi; chen := proc(m) if irem(m+1,4) = 0 then RETURN(0) fi; 1/((-1)^iquo(m+1,4) *2^iquo(m,2)) end; add(add((-1)^v*binomial(k,v)*pow(v+x+1,n)*chen(k),v=0..k), k=0..n) end:
    # Coefficients with zeros:
    seq(print(seq(coeff(i!*coeff(series(exp(x*t)*sech(t),t,16),t,i),x,i-n),n=0..i)), i=0..8);
    # Recursion
    W := proc(n,z) option remember; local k,p;
    if n = 0 then 1 else p := irem(n+1,2);
    z^n - p + add(`if`(irem(k,2)=1,0,
    W(k,0)*binomial(n,k)*(power(z,n-k)-p)),k=2..n-1) fi end:
    # Peter Luschny, edited and additions Jul 07 2009, May 13 2010, Oct 24 2011
  • Mathematica
    max = 9; rows = (Reverse[ CoefficientList[ #, x]] & ) /@ CoefficientList[ Series[ Exp[x*t]*Sech[t], {t, 0, max}], t]*Range[0, max]!; par[coefs_] := (p = Partition[ coefs, 2][[All, 1]]; If[ EvenQ[ Length[ coefs]], p, Append[ p, Last[ coefs]]]); Flatten[ par /@ rows] (* Jean-François Alcover, Oct 03 2011, after g.f. *)
    sk[n_, x_] := Sum[Binomial[n, k]*EulerE[k]*x^(n-k), {k, 0, n}]; Table[CoefficientList[sk[n, x], x] // Reverse // Select[#, # =!= 0 &] &, {n, 0, 13}] // Flatten (* Jean-François Alcover, May 21 2013 *)
    Flatten@Table[Binomial[n, 2k] EulerE[2k], {n, 0, 12}, {k, 0, n/2}](* Oliver Seipel, Jan 14 2025 *)
  • Sage
    def A046978(k):
        if k % 4 == 0:
            return 0
        return (-1)**(k // 4)
    def A153641_poly(n, x):
        return expand(add(2**(-(k // 2))*A046978(k+1)*add((-1)**v*binomial(k,v)*(v+x+1)**n for v in (0..k)) for k in (0..n)))
    for n in (0..7): print(A153641_poly(n, x))  # Peter Luschny, Oct 24 2011

Formula

W_n(x) = Sum_{k=0..n}{v=0..k} (-1)^v binomial(k,v)*c_k*(x+v+1)^n where c_k = frac((-1)^(floor(k/4))/2^(floor(k/2))) [4 not div k] (Iverson notation).
From Peter Bala, Jun 10 2009: (Start)
E.g.f.: 2*exp(x*t)*(exp(t)-exp(3*t))/(1-exp(4*t))= 1 + x*t + (x^2-1)*t^2/2! + (x^3-3*x)*t^3/3! + ....
W_n(x) = 1/(2*n+2)*Sum_{k=0..n+1} 1/(k+1)*Sum_{i=0..k} (-1)^i*binomial(k,i)*((x+4*i+3)^(n+1) - (x+4*i+1)^(n+1)).
Fourier series expansion for the generalized Bernoulli polynomials:
B(X;2*n,x) = (-1)^n*(2/Pi)^(2*n)*(2*n)! * {sin(Pi*x/2)/1^(2*n) - sin(3*Pi*x/2)/3^(2*n) + sin(5*Pi*x/2)/5^(2*n) - ...}, valid for 0 <= x <= 1 when n >= 1.
B(X;2*n+1,x) = (-1)^(n+1)*(2/Pi)^(2*n+1)*(2*n+1)! * {cos(Pi*x/2)/1^(2*n+1) - cos(3*Pi*x/2)/3^(2*n+1) + cos(5*Pi*x/2)/5^(2*n+1) - ...}, valid for 0 <= x <= 1 when n >= 1 and for 0 <= x < 1 when n = 0.
(End)
E.g.f.: exp(x*t) * sech(t). - Peter Luschny, Jul 07 2009
O.g.f. as a J-fraction: z/(1-x*z+z^2/(1-x*z+4*z^2/(1-x*z+9*z^2/(1-x*z+...)))) = z + x*z^2 + (x^2-1)*z^3 + (x^3-3*x)*z^4 + .... - Peter Bala, Mar 11 2012
Conjectural o.g.f.: Sum_{n >= 0} (1/2^((n-1)/2))*cos((n+1)*Pi/4)*( Sum_{k = 0..n} (-1)^k*binomial(n,k)/(1 - (k + x)*t) ) = 1 + x*t + (x^2 - 1)*t^2 + (x^3 - 3*x)*t^3 + ... (checked up to O(t^13)), which leads to W_n(x) = Sum_{k = 0..n} 1/2^((k - 1)/2)*cos((k + 1)*Pi/4)*( Sum_{j = 0..k} (-1)^j*binomial(k, j)*(j + x)^n ). - Peter Bala, Oct 03 2016

A000464 Expansion of e.g.f. sin(x)/cos(2*x).

Original entry on oeis.org

1, 11, 361, 24611, 2873041, 512343611, 129570724921, 44110959165011, 19450718635716001, 10784052561125704811, 7342627959965776406281, 6023130568334172003579011, 5858598896811701995459355761, 6667317162352419006959182803611, 8776621742176931117228228227924441
Offset: 0

Views

Author

Keywords

Comments

From Peter Bala, Dec 22 2021: (Start)
Conjectures:
1) Taking the sequence (a(n))n>=1 modulo an integer k gives a purely periodic sequence with period dividing phi(k). For example, the sequence taken modulo 21 begins [11, 4, 20, 10, 17, 1, 11, 4, 20, 10, 17, 1, 11, 4, 20, 10, 17, 1, ...] with an apparent period of length 6, which divides phi(21) = 12.
2) For i >= 0, define a_i(n) = a(n+i). Then for each i the Gauss congruences a_i(n*p^k) == a_i(n*p^(k-1)) ( mod p^k ) hold for all prime p and positive integers n and k. If true, then for each i the expansion of exp(Sum_{n >= 1} a_i(n)*x^n/n) has integer coefficients.
3) a(m*n) == a(m)^n (mod 2^k) for k = 2*v_2(m) + 4, where v_p(i) denotes the p-adic valuation of i.
4)(i) a(2*m*n) == a(n)^(2*m) (mod 2^k) for k = v_2(m) + 4
(ii) a((2*m+1)*n) == a(n)^(2*m+1) (mod 2^k) for k = v_2(m) + 4. (End)

References

  • H. Cohen, Number Theory - Volume II: Analytic and Modern Tools, Graduate Texts in Mathematics. Springer-Verlag.
  • J. W. L. Glaisher, "On the coefficients in the expansions of cos x/ cos 2x and sin x/ cos 2x", Quart. J. Pure and Applied Math., 45 (1914), 187-222.
  • I. J. Schwatt, Intro. to Operations with Series, Chelsea, p. 278.
  • 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

Row 2 of A235606.
Cf. A064073. Bisection of A000822, A001586.

Programs

  • Maple
    a := n -> (-1)^n*2^(6*n+4)*(Zeta(0, -2*n-1, 5/8)-Zeta(0, -2*n-1, 7/8)):
    seq(a(n), n=0..12); # Peter Luschny, Oct 15 2015
  • Mathematica
    With[{nn=30},Take[CoefficientList[Series[Sin[x]/Cos[2x],{x,0,nn}],x] Range[0,nn-1]!,{2,-1,2}]] (* Harvey P. Dale, Mar 23 2012 *)
    nmax = 15; km0 = 10; d[n_, km_] := Round[(2^(4n-1/2) (2n-1)! Sum[ JacobiSymbol[2, 2k+1]/(2k+1)^(2n), {k, 0, km}])/Pi^(2n)]; dd[km_] := dd[km] = Table[d[n, km], {n, 1, nmax}]; dd[km0]; dd[km = 2*km0]; While[dd[km] != dd[km/2, km = 2*km]]; A000464 = dd[km] (* Jean-François Alcover, Feb 08 2016 *)
  • PARI
    a(n)=if(n<0, 0, n+=n+1; n!*polcoeff(sin(x+x*O(x^n))/cos(2*x+x*O(x^n)),n)) /* Michael Somos, Feb 09 2006 */

Formula

E.g.f.: Sum_{k>=0} a(k)x^(2k+1)/(2k+1)! = sin(x)/cos(2x).
a(n) = (-1)^n*L(X,-2n+1) where L(X,z) is the Dirichlet L-function L(X,z) = Sum_{k>=0} X(k)/k^z and where X(k) is the Dirichlet character Legendre(k,2) which begins 1,0,-1,0,-1,0,1,0,1,0,-1,0,-1,0,1,0,1,0,-1,0.... - Benoit Cloitre, Mar 22 2009 [This Dirichlet character is A091337. - Jianing Song, Oct 22 2023]
From Peter Bala, Mar 24 2009: (Start)
Basic hypergeometric generating function:
2*exp(-t)*Sum_{n = 0..inf} (Product_{k = 1..n} (1-exp(-16*k*t))/Product_{k = 1..n+1} (1+exp(-(16*k-8)*t))) = 1 + 11*t + 361*t^2/2! + 24611*t^3/3! + .... For other sequences with generating functions of a similar type see A000364, A002105, A002439, A079144 and A158690.
a(n) = (-1)^(n+1)*L(-2*n-1), where L(s) is the Dirichlet L-function L(s) = 1 - 1/3^s - 1/5^s + 1/7^s + - - + ... [Andrews et al., Theorem 5]. (End)
From Peter Bala, Jun 18 2009: (Start)
a(n) = (-1)^n*B_(2*n+2)(X)/(2*n+2), where B_n(X) denotes the X-Bernoulli number with X a Dirichlet character modulus 8 given by X(8*n+1) = X(8*n+7) = 1, X(8*n+3) = X(8*n+5) = -1 and X(2*n) = 0. See A161722 for the values of B_n(X).
For the theory and properties of the generalized Bernoulli numbers B_n(X) and the associated generalized Bernoulli polynomials B_n(X,x) see [Cohen, Section 9.4].
The present sequence also occurs in the evaluation of the finite sum of powers Sum_{i = 0..m-1} {(8*i+1)^n - (8*i+3)^n - (8*i+5)^n + (8*i+7)^n}, n = 1,2,... - see A151751 for details. (End)
G.f. 1/G(0) where G(k) = 1 + x - x*(4*k+3)*(4*k+4)/(1 - (4*k+4)*(4*k+5)*x/G(k+1)); (continued fraction, 2-step). - Sergei N. Gladkovskii, Aug 11 2012
G.f.: 1/E(0) where E(k) = 1 - 11*x - 32*x*k*(k+1) - 16*x^2*(k+1)^2*(4*k+3)*(4*k+5)/E(k+1) (continued fraction, 1-step). - Sergei N. Gladkovskii, Sep 17 2012
a(n) ~ (2*n+1)! * 2^(4*n+7/2) / Pi^(2*n+2). - Vaclav Kotesovec, May 03 2014
a(n) = (-1)^n*2^(6*n+4)*(Zeta(-2*n-1,5/8)-Zeta(-2*n-1,7/8)). - Peter Luschny, Oct 15 2015
From Peter Bala, May 11 2017: (Start)
G.f. A(x) = 1 + 11*x + 361*x^2 + ... = 1/(1 + x - 12*x/(1 - 20*x/(1 + x - 56*x/(1 - 72*x/(1 + x - ... - 4*n*(4*n - 1)*x/(1 - 4*n*(4*n + 1)*x/((1 + x) - ...))))))).
A(x) = 1/(1 + 9*x - 20*x/(1 - 12*x/(1 + 9*x - 72*x/(1 - 56*x/(1 + 9*x - ... - 4*n*(4*n + 1)*x/(1 - 4*n*(4*n - 1)*x/(1 + 9*x - ...))))))).
It follows that the first binomial transform of A(x) and the ninth binomial transform of A(x) have continued fractions of Stieltjes-type (S-fractions). (End)
a(n) = (-1)^(n+1)*4^(2*n+1)*E(2*n+1,1/4), where E(n,x) is the n-th Euler polynomial. Cf. A002439. - Peter Bala, Aug 13 2017
From Peter Bala, Dec 04 2021: (Start)
F(x) = exp(x)*(exp(2*x) - 1)/(exp(4*x) + 1) = x - 11*x^3/3! + 361x^5/5! - 24611*x^7/7! + ... is the e.g.f. for the sequence [1, 0, -11, 0, 361, 0, -24611, 0, ...], a signed and aerated version of this sequence.
The binomial transform exp(x)*F(x) = x + 2*x^2/2! - 8*x^3/3! - 40*x^4/4! + + - - is an e.g.f. for a signed version of A000828 (omitting the initial term). (End)
From Peter Bala, Dec 22 2021: (Start)
a(1) = 1, a(n) = (-1)^(n-1) - Sum_{k = 1..n} (-4)^k*C(2*n-1,2*k)*a(n-k).
a(n) == 1 (mod 10); a(5*n+1) == 0 mod(11);
a(n) == - 23^(n+1) (mod 108); a(n) == (7^2)*59^n (mod 144);
a(n) == 11^n (mod 240); a(n) == (11^2)*131^n (mod 360). (End)

Extensions

Better description, new reference, Aug 15 1995

A161722 Generalized Bernoulli numbers B_n(X,0), X a Dirichlet character modulus 8.

Original entry on oeis.org

0, 2, -44, 2166, -196888, 28730410, -6148123332, 1813990148894, -705775346640176, 350112935442888018, -215681051222514096220, 161537815119247080938182, -144555133640020128085896264, 152323571317104251881943249786
Offset: 0

Views

Author

Peter Bala, Jun 18 2009

Keywords

Comments

Let X be a periodic arithmetical function with period m. The generalized Bernoulli polynomials B_n(X,x) attached to X are defined by means of the generating function
(1)... t*exp(t*x)/(exp(m*t)-1) * Sum_{r = 0..m-1} X(r)*exp(r*t) = Sum_{n >= 0} B_n(X,x)*t^n/n!.
The values B_n(X,0) are generalizations of the Bernoulli numbers (case X = 1). For the theory and properties of these polynomials and numbers see [Cohen, Section 9.4]. In the present case, X is chosen to be the Dirichlet character modulus 8 given by
(2)... X(8*n+1) = X(8*n+7) = 1; X(8*n+3) = X(8*n+5) = -1; X(2*n) = 0.
The odd-indexed generalized Bernoulli numbers B_(2*n+1)(X,0) vanish. The current sequence lists the even-indexed values B_(2*n)(X,0).
The coefficients of the generalized Bernoulli polynomials B_n(X,x) are listed in A151751.

References

  • H. Cohen, Number Theory - Volume II: Analytic and Modern Tools, Graduate Texts in Mathematics. Springer-Verlag.

Crossrefs

Programs

  • Maple
    G := x*sinh(x)/cosh(2*x): ser := series(G, x, 30):
    seq((2*n)!*coeff(ser, x, 2*n), n = 0..14); # Peter Luschny, Nov 26 2020
    # After an observation of F. Chapoton in A117442:
    A161722 := proc(n) 4^n*add(binomial(2*n, k)*euler(k)*((x+1)/2)^(2*n-k), k=0..2*n);
    coeff(%, x, 1) end: seq(A161722(n), n=0..13); # Peter Luschny, Nov 26 2020
  • Mathematica
    terms = 13;
    (CoefficientList[x(Sinh[x]/Cosh[2x]) + O[x]^(2terms+3), x] Range[0, 2terms+2]!)[[ ;; ;; 2]] (* Jean-François Alcover, Nov 16 2020 *)

Formula

(1)... a(n) = (-1)^(n+1)*2*n*A000464(n-1).
The sequence of generalized Bernoulli numbers
(2)... [B_n(X,0)]n>=2 = [2,0,-44,0,2166,0,...]
has the e.g.f.
(3)... t*(exp(t)-exp(3*t)-exp(5*t)+exp(7*t))/(exp(8*t)-1),
which simplifies to
(4)... t*sinh(t)/cosh(2*t) = 2*t^2/2! - 44*t^4/4! + ....
Hence
(5)... B_(2*n)(X,0) = (-1)^(n+1)*2*n*A000464(n-1) and B_(2*n+1)(X,0) = 0.
a(n) = (-1/2)*16^n*n*euler(2*n-1, 1/4) for n >= 1 after a formula of Peter Bala in A000464. - Peter Luschny, Nov 26 2020

Extensions

Cross-reference corrected by Peter Bala, Jun 22 2009
Offset set to 0 and a(0) = 0 prepended by Peter Luschny, Nov 26 2020
Showing 1-3 of 3 results.