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

A036968 Genocchi numbers (of first kind): expansion of 2*x/(exp(x)+1).

Original entry on oeis.org

1, -1, 0, 1, 0, -3, 0, 17, 0, -155, 0, 2073, 0, -38227, 0, 929569, 0, -28820619, 0, 1109652905, 0, -51943281731, 0, 2905151042481, 0, -191329672483963, 0, 14655626154768697, 0, -1291885088448017715, 0, 129848163681107301953
Offset: 1

Views

Author

Keywords

Comments

The sign of a(1) depends on which convention one chooses: B(n) = B_n(1) or B(n) = B_n(0) where B(n) are the Bernoulli numbers and B_n(x) the Bernoulli polynomials (see the Wikipedia article on Bernoulli numbers). The definition given is in line with B(n) = B_n(0). The convention B(n) = B_n(1) corresponds to the e.g.f. -2*x/(1+exp(-x)). - Peter Luschny, Jun 28 2013
According to Hetyei [2017], "alternation acyclic tournaments in which at least one ascent begins at each vertex, except for the largest one, are counted by the Genocchi numbers of the first kind." - Danny Rorabaugh, Apr 25 2017
Named after the Italian mathematician Angelo Genocchi (1817-1889). - Amiram Eldar, Jun 06 2021
Conjecture: For any positive integer n, -a(n+1) is the permanent of the n X n matrix M with M(j, k) = floor((2*j - k)/n), (j,k=1..n). - Zhi-Wei Sun, Sep 07 2021
A corresponding conjecture can also be made for L. Seidel's 'Genocchi numbers of second kind' A005439. - Peter Luschny, Sep 08 2021

References

  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 49.
  • 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. 73.
  • Ronald L. Graham, Donald E. Knuth and Oren Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 528.
  • Richard P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 5.8.

Crossrefs

A001469 is the main entry for this sequence. A226158 is another version.
Cf. A005439 (Genocchi numbers of second kind).

Programs

  • Maple
    a := n -> n*euler(n-1,0); # Peter Luschny, Jul 13 2009
  • Mathematica
    a[n_] := n*EulerE[n - 1, 0]; Table[a[n], {n, 1, 32}] (* Jean-François Alcover, Dec 08 2011, after Peter Luschny *)
    Range[0, 31]! CoefficientList[ Series[ 2x/(1 + Exp[x]), {x, 0, 32}], x] (* Robert G. Wilson v, Oct 26 2012 *)
    Table[(-1)^n 2 n PolyLog[1 - n, -1], {n, 1, 32}] (* Peter Luschny, Aug 17 2021 *)
  • PARI
    {a(n) = if( n<0, 0, n! * polcoeff( 2*x / (1 + exp(x + x * O(x^n))), n))}; /* Michael Somos, Jul 23 2005 */
    
  • PARI
    /* From o.g.f. (Paul D. Hanna, Aug 03 2014): */
    {a(n)=local(A=1); A=x*sum(m=0, n, m!*(-x)^m/(1-m*x)/prod(k=1,m,1 - k*x +x*O(x^n))); polcoeff(A, n)}
    for(n=1, 32, print1(a(n), ", "))
    
  • Python
    from sympy import bernoulli
    def A036968(n): return (2-(2<Chai Wah Wu, Apr 14 2023
  • Sage
    # with a(1) = -1
    [z*zeta(1-z)*(2^(z+1)-2) for z in (1..32)]  # Peter Luschny, Jun 28 2013
    
  • Sage
    def A036968_list(len):
        e, f, R, C = 4, 1, [], [1]+[0]*(len-1)
        for n in (2..len-1):
            for k in range(n, 0, -1):
                C[k] = C[k-1] / (k+1)
            C[0] = -sum(C[k] for k in (1..n))
            R.append((2-e)*f*C[0])
            f *= n; e *= 2
        return R
    print(A036968_list(34)) # Peter Luschny, Feb 22 2016
    

Formula

E.g.f.: 2*x/(exp(x)+1).
a(n) = 2*(1-2^n)*B_n (B = Bernoulli numbers). - Benoit Cloitre, Oct 26 2003
2*x/(exp(x)+1) = x + Sum_{n>=1} x^(2*n)*G_{2*n}/(2*n)!.
a(n) = Sum_{k=0..n-1} binomial(n,k) 2^k*B(k). - Peter Luschny, Apr 30 2009
From Sergei N. Gladkovskii, Dec 12 2012 to Nov 23 2013: (Start) Continued fractions:
E.g.f.: 2*x/(exp(x)+1) = x - x^2/2*G(0) where G(k) = 1 - x^2/(x^2 + 4*(2*k+1)*(2*k+3)/G(k+1)).
E.g.f.: 2/(E(0)+1) where E(k) = 1 + x/(2*k+1 - x*(2*k+1)/(x + (2*k+2)/E(k+1))).
G.f.: 2 - 1/G(0) where G(k) = 1 - x*(k+1)/(1 + x*(k+1)/(1 - x*(k+1)/(1 + x*(k+1)/G(k+1)))).
E.g.f.: 2*x/(1 + exp(x)) = 2*x-2 - 2*T(0), where T(k) = 4*k-1 + x/(2 - x/( 4*k+1 + x/(2 - x/T(k+1)))).
G.f.: 2 - Q(0)/(1-x+x^2) where Q(k) = 1 - x^4*(k+1)^4/(x^4*(k+1)^4 - (1 - x + x^2 + 2*x^2*k*(k+1))*(1 - x + x^2 + 2*x^2*(k+1)*(k+2))/Q(k+1)). (End)
a(n) = n*zeta(1-n)*(2^(n+1)-2) for n > 1. - Peter Luschny, Jun 28 2013
O.g.f.: x*Sum_{n>=0} n! * (-x)^n / (1 - n*x) / Product_{k=1..n} (1 - k*x). - Paul D. Hanna, Aug 03 2014
Sum_{n>=1} 1/a(2*n) = A321595. - Amiram Eldar, May 07 2021
a(n) = (-1)^n*2*n*PolyLog(1 - n, -1). - Peter Luschny, Aug 17 2021

A083007 a(n) = Sum_{k=0..n-1} 3^k*B(k)*C(n,k) where B(k) is the k-th Bernoulli number and C(n,k)=binomial(n,k).

Original entry on oeis.org

0, 1, -2, 1, 4, -5, -26, 49, 328, -809, -6710, 20317, 201772, -722813, -8370194, 34607305, 457941136, -2145998417, -31945440878, 167317266613, 2767413231220, -16020403322021, -291473080313162, 1848020950359841, 36679231132772824, -252778977216700025, -5435210060467425446
Offset: 0

Views

Author

Benoit Cloitre, May 31 2003

Keywords

Crossrefs

Programs

  • Maple
    A083007 := proc(n)
        3*x/(1+exp(x)+exp(2*x)) ;
        coeftayl(%,x=0,n) ;
        %*n! ;
    end proc:
    seq(A083007(n),n=0..30) ; # R. J. Mathar, Jul 13 2023
  • Mathematica
    Range[0, 15]! CoefficientList[ Series[ 3x/(1 + Exp[x] + Exp[ 2x]), {x, 0, 15}], x] (* Robert G. Wilson v, Oct 26 2012 *)
    Table[Sum[3^k BernoulliB[k]Binomial[n,k],{k,0,n-1}],{n,0,30}] (* Harvey P. Dale, May 26 2014 *)
  • PARI
    a(n)=sum(k=0,n-1,3^k*binomial(n,k)*bernfrac(k))

Formula

E.g.f.: 3x/(1+e^x+e^(2x)). - Ira M. Gessel, Jan 28 2012
From Peter Bala, Mar 01 2015: (Start)
a(2*n+1) = (-1)^(n+1)*A002111(n) for n >= 1.
a(n) = 3^n * ( B(n,1/3) - B(n,0) ), where B(n,x) denotes the n-th Bernoulli polynomial. More generally, Almkvist and Meurman show that k^n * ( B(n, 1/k) - B(n, 0) ) is an integer sequence for k = 2,3,4,..., which proves the integrality of A083008 through A083014.
a(0) = 1 and for n >= 1, a(n) = 1 - 1/(n + 1)*Sum_{k = 1..n-1} 3^(n-k)*binomial(n+1,k)*a(k) (Sury, Section 1). (End)

A083008 a(n) = Sum_{k=0..n-1} 4^k*B(k)*C(n,k) where B(k) is the k-th Bernoulli number and C(n,k) = binomial(n,k).

Original entry on oeis.org

0, 1, -3, 3, 9, -25, -99, 427, 2193, -12465, -79515, 555731, 4247577, -35135945, -313193811, 2990414715, 30461046561, -329655706465, -3777604994187, 45692713833379, 581778811909545, -7777794952988025, -108933009112011843, 1595024111042171723, 24370176181315498929
Offset: 0

Views

Author

Benoit Cloitre, May 31 2003

Keywords

Crossrefs

Programs

  • Mathematica
    Range[0, 15]! CoefficientList[ Series[ 4x/(1 + Exp[x] + Exp[ 2x] + Exp[ 3x]), {x, 0, 15}], x] (* Robert G. Wilson v, Oct 26 2012 *)
    Table[Sum[4^k*BernoulliB[k] Binomial[n, k], {k, 0, n - 1}], {n, 0, 24}] (* Michael De Vlieger, Sep 28 2016 *)
  • PARI
    a(n)=sum(k=0,n-1,4^k*binomial(n,k)*bernfrac(k))

Formula

E.g.f.: 4*x/(1+exp(x)+exp(2*x)+exp(3*x)). - Ira M. Gessel, Feb 23 2012
a(n) ~ n! * (cos(n*Pi/2)-sin(n*Pi/2)) * 2^(n+1) / Pi^n. - Vaclav Kotesovec, Mar 02 2014

Extensions

Offset changed to 0 by Seiichi Manyama, Sep 28 2016

A083009 a(n) = Sum_{k=0,n-1} 5^k*B(k)*binomial(n,k) where B(k) is the k-th Bernoulli number.

Original entry on oeis.org

0, 1, -4, 6, 16, -74, -264, 1946, 9056, -88434, -512024, 6154786, 42716496, -607884394, -4920817384, 80834386026, 747784582336, -13923204233954, -144898927180344, 3015393801263666, 34867899296006576, -801997872697905114, -10201104981227536904, 256982712667627683706
Offset: 0

Views

Author

Benoit Cloitre, May 31 2003

Keywords

Crossrefs

Programs

  • Mathematica
    Range[0, 15]! CoefficientList[ Series[ 5x/(1 + Exp[x] + Exp[ 2x] + Exp[ 3x] + Exp[ 4x]), {x, 0, 15}], x] (* Robert G. Wilson v, Oct 26 2012 *)
    Table[Sum[5^k*BernoulliB[k] Binomial[n, k], {k, 0, n - 1}], {n, 0, 23}] (* Michael De Vlieger, Sep 28 2016 *)
  • PARI
    a(n)=sum(k=0,n-1,5^k*binomial(n,k)*bernfrac(k))

Formula

E.g.f.: 5x/(1+exp(x)+exp(2x)+exp(3x)+exp(4x)). - Benoit Cloitre, Oct 26 2012 (following I. Gessel).

Extensions

Offset changed to 0 by Seiichi Manyama, Sep 28 2016

A083010 a(n) = 6^n(B_n(1/6)-B_n(0)) where B_n(x) is the n-th Bernoulli polynomial.

Original entry on oeis.org

0, 1, -5, 10, 25, -170, -575, 6370, 28225, -415826, -2294975, 41649850, 275622625, -5922729722, -45718037855, 1134081384850, 10004182986625, -281284596509858, -2791456543622015, 87722769712529770, 967282878165054625, -33597252908389628234, -407509096583935700255
Offset: 0

Views

Author

Benoit Cloitre, May 31 2003

Keywords

Crossrefs

Programs

  • Mathematica
    Range[0, 15]! CoefficientList[ Series[ 6x/(1 + Exp[x] + Exp[ 2x] + Exp[ 3x] + Exp[ 4x] + Exp[ 5x]), {x, 0, 15}], x] (* Robert G. Wilson v, Oct 26 2012 *)
  • PARI
    a(n)=sum(k=0,n-1,6^k*binomial(n,k)*bernfrac(k))
    
  • PARI
    {a(n)=if(n<1, 0, n!*polcoeff( 6*x*(exp(x+x*O(x^n))-1)/(exp(6*x +x*O(x^n))-1), n))} /* Michael Somos, Aug 02 2006 */

Formula

E.g.f.: 6x(exp(x)-1)/(exp(6x)-1). - Michael Somos, Aug 02 2006
a(n) = Sum_{k=0..n-1} 6^k*B(k)*C(n,k) where B(k) is the k-th Bernoulli number and C(n,k) = binomial(n,k).

A083011 a(n) = Sum_{k=0..n-1} 7^k*B(k)*binomial(n,k) where B(k) is the k-th Bernoulli number.

Original entry on oeis.org

0, 1, -6, 15, 36, -335, -1098, 16955, 73032, -1503963, -8075430, 204957775, 1319806188, -39666688711, -297958666242, 10337889346275, 88743928066704, -3489994294713779, -33703905982634334, 1481439997178305655, 15896303102840841780, -772269573963075710367
Offset: 0

Views

Author

Benoit Cloitre, May 31 2003

Keywords

Crossrefs

Programs

  • Mathematica
    Range[0, 15]! CoefficientList[ Series[ 7x/(1 + Exp[x] + Exp[ 2x] + Exp[ 3x] + Exp[ 4x] + Exp[ 5x] + Exp[ 6x]), {x, 0, 15}], x] (* Robert G. Wilson v, Oct 26 2012 *)
  • PARI
    a(n)=sum(k=0,n-1,7^k*binomial(n,k)*bernfrac(k))

Extensions

Offset changed to 0 by Seiichi Manyama, Sep 28 2016

A083012 a(n) = Sum_{k=0..n-1} 8^k*B(k)*binomial(n,k) where B(k) is the k-th Bernoulli number.

Original entry on oeis.org

0, 1, -7, 21, 49, -595, -1911, 39109, 165473, -4525731, -23883335, 805349237, 5097585297, -203564524787, -1503073984279, 69292329479205, 584713994953921, -30553447357629763, -290046835163027943, 16939595863125337813, 178676615255242261745
Offset: 0

Views

Author

Benoit Cloitre, May 31 2003

Keywords

Crossrefs

Programs

  • Mathematica
    Range[0, 15]! CoefficientList[ Series[ 8x/(1 + Exp[x] + Exp[ 2x] + Exp[ 3x] + Exp[ 4x] + Exp[ 5x] + Exp[ 6x] + Exp[ 7x]), {x, 0, 15}], x] (* Robert G. Wilson v, Oct 26 2012 *)
  • PARI
    a(n)=sum(k=0,n-1,8^k*binomial(n,k)*bernfrac(k))

Extensions

Offset changed to 0 by Seiichi Manyama, Sep 28 2016

A083013 a(n) = Sum_{k=0..n-1} 9^k*B(k)*binomial(n,k) where B(k) is the k-th Bernoulli number.

Original entry on oeis.org

0, 1, -8, 28, 64, -980, -3104, 81172, 339328, -11878244, -61958240, 2674671076, 16735235392, -855605816468, -6245150369696, 368601472639540, 3074742020313856, -205700802920736452, -1930357641628367072, 144338957346266943364, 1505019970814899568320
Offset: 0

Views

Author

Benoit Cloitre, May 31 2003

Keywords

Crossrefs

Programs

  • Mathematica
    Range[0, 15]! CoefficientList[ Series[ 9x/(1 + Exp[x] + Exp[ 2x] + Exp[ 3x] + Exp[ 4x] + Exp[ 5x] + Exp[ 6x] + Exp[ 7x] + Exp[ 8x]), {x, 0, 15}], x] (* Robert G. Wilson v, Oct 26 2012 *)
    Table[Sum[9^k BernoulliB[k]Binomial[n,k],{k,0,n-1}],{n,0,20}] (* Harvey P. Dale, Apr 13 2016 *)
  • PARI
    a(n)=sum(k=0,n-1,9^k*binomial(n,k)*bernfrac(k))

Extensions

Offset changed to 0 by Seiichi Manyama, Sep 28 2016
Showing 1-8 of 8 results.