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

A001469 Genocchi numbers (of first kind); unsigned coefficients give expansion of x*tan(x/2).

Original entry on oeis.org

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

Views

Author

Keywords

Comments

The Genocchi numbers satisfy Seidel's recurrence: for n>1, 0 = Sum_{j=0..[n/2]} C(n,2j)*a(n-j). - Ralf Stephan, Apr 17 2004
The (n+1)st Genocchi number is the number of Dumont permutations of the first kind on 2n letters. In a Dumont permutation of the first kind, each even integer must be followed by a smaller integer and each odd integer is either followed by a larger integer or is the last element. - Ralf Stephan, Apr 26 2004
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

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 49.
  • L. Euler, Institutionum Calculi Differentialis, volume 2 (1755), para. 181.
  • 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.
  • A. Genocchi, Intorno all'espressione generale de'numeri Bernulliani, Ann. Sci. Mat. Fis., 3 (1852), 395-405.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 528.
  • 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).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 5.8.

Crossrefs

a(n) = -A065547(n, 1) and A065547(n+1, 2) for n >= 1.

Programs

  • Magma
    [2*(1 - 4^n) * Bernoulli(2*n): n in [1..25]]; // Vincenzo Librandi, Oct 15 2018
    
  • Maple
    A001469 := proc(n::integer) (2*n)!*coeftayl( 2*x/(exp(x)+1), x=0,2*n) end proc:
    for n from 1 to 20 do print(A001469(n)) od : # R. J. Mathar, Jun 22 2006
  • Mathematica
    a[n_] := 2*(1-4^n)*BernoulliB[2n]; Table[a[n], {n, 17}] (* Jean-François Alcover, Nov 24 2011 *)
    a[n_] := 2*n*EulerE[2*n-1, 0]; Table[a[n], {n, 17}] (* Jean-François Alcover, Jul 02 2013 *)
    Table[4 n PolyLog[1 - 2 n, -1], {n, 1, 19}] (* Peter Luschny, Aug 17 2021 *)
  • PARI
    a(n)=if(n<1,0,n*=2; 2*(1-2^n)*bernfrac(n))
    
  • PARI
    {a(n)=polcoeff(sum(m=0, n, m!^2*(-x)^(m+1)/prod(k=1, m, 1-k^2*x+x*O(x^n))), n)} /* Paul D. Hanna, Jul 21 2011 */
    
  • Python
    from sympy import bernoulli
    def A001469(n): return (2-(2<<(m:=n<<1)))*bernoulli(m) # Chai Wah Wu, Apr 14 2023
  • Sage
    # Algorithm of L. Seidel (1877)
    # n -> [a(1), ..., a(n)] for n >= 1.
    def A001469_list(n) :
        D = [0]*(n+2); D[1] = -1
        R = []; b = False
        for i in(0..2*n-1) :
            h = i//2 + 1
            if b :
                for k in range(h-1, 0, -1) : D[k] -= D[k+1]
            else :
                for k in range(1, h+1, 1) :  D[k] -= D[k-1]
            b = not b
            if not b : R.append(D[h])
        return R
    A001469_list(17) # Peter Luschny, Jun 29 2012
    

Formula

a(n) = 2*(1-4^n)*B_{2n} (B = Bernoulli numbers).
x*tan(x/2) = Sum_{n>=1} x^(2*n)*abs(a(n))/(2*n)! = x^2/2 + x^4/24 + x^6/240 + 17*x^8/40320 + 31*x^10/725760 + O(x^11).
E.g.f.: 2*x/(1 + exp(x)) = x + Sum_{n>=1} a(2*n)*x^(2*n)/(2*n)! = -x^2/2! + x^4/4! - 3 x^6/6! + 17 x^8/8! + ...
O.g.f.: Sum_{n>=0} n!^2*(-x)^(n+1) / Product_{k=1..n} (1-k^2*x). - Paul D. Hanna, Jul 21 2011
a(n) = Sum_{k=0..2n-1} 2^k*B(k)*binomial(2*n,k) where B(k) is the k-th Bernoulli number. - Benoit Cloitre, May 31 2003
abs(a(n)) = Sum_{k=0..2n} (-1)^(n-k+1)*Stirling2(2n, k)*A059371(k). - Vladeta Jovovic, Feb 07 2004
G.f.: -x/(1+x/(1+2x/(1+4x/(1+6x/(1+9x/(1+12x/(1+16x/(1+20x/(1+25x/(1+...(continued fraction). - Philippe Deléham, Nov 22 2011
E.g.f.: E(x) = 2*x/(exp(x)+1) = x*(1-(x^3+2*x^2)/(2*G(0)-x^3-2*x^2)); G(k) = 8*k^3 + (12+4*x)*k^2 + (4+6*x+2*x^2)*k + x^3 + 2*x^2 + 2*x - 2*(x^2)*(k+1)*(2*k+1)*(x+2*k)*(x+2*k+4)/G(k+1); (continued fraction, Euler's kind, 1-step). - Sergei N. Gladkovskii, Jan 18 2012
a(n) = (-1)^n*(2*n)!*Pi^(-2*n)*4*(1-4^(-n))*Li{2*n}(1). - Peter Luschny, Jun 29 2012
Asymptotic: abs(a(n)) ~ 8*Pi*(2^(2*n)-1)*(n/(Pi*exp(1)))^(2*n+1/2)*exp(1/2+(1/24)/n-(1/2880)/n^3+(1/40320)/n^5+...). - Peter Luschny, Jul 24 2013
G.f.: x/(T(0)-x) -1, where T(k) = 2*x*k^2 + 4*x*k + 2*x - 1 - x*(-1+x+2*x*k+x*k^2)*(k+2)^2/T(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 17 2013
G.f.: -1 + x/(T(0)+x), where T(k) = 1 + (k+1)*(k+2)*x/(1+x*(k+2)^2/T(k+1)); (continued fraction). - Sergei N. Gladkovskii, Nov 17 2013
a(n) = 4*n*PolyLog(1 - 2*n, -1). - Peter Luschny, Aug 17 2021

A000795 Salié numbers: expansion of cosh x / cos x = Sum_{n >= 0} a(n)*x^(2n)/(2n)!.

Original entry on oeis.org

1, 2, 12, 152, 3472, 126752, 6781632, 500231552, 48656756992, 6034272215552, 929327412759552, 174008703107274752, 38928735228629389312, 10255194381004799025152, 3142142941901073853366272, 1107912434323301224813002752, 445427836895850552387642130432
Offset: 0

Views

Author

Keywords

Comments

Named after the German mathematician Hans Salié (1902-1978). - Amiram Eldar, Jun 10 2021

Examples

			cosh x / cos x = Sum_{n>=0} a(n)*x^(2n)/(2n)! = 1 + x^2 + (1/2)*x^4 + (19/90)*x^6 + (31/360)*x^8 + (3961/113400)*x^10 + ...
G.f. = 1 + 2*x + 12*x^2 + 252*x^3 + 3472*x^4 + 126752*x^5 + 6781632*x^6 + ...
		

References

  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 86, Problem 32.
  • Hans Salié, Arithmetische Eigenschaften der Koeffizienten einer speziellen Hurwitzschen Potenzreihe, Wiss. Z. Karl-Marx-Univ. Leipzig Math.-Natur. Reihe 12 (1963), pp. 617-618.
  • 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

A005647(n) = a(n)/2^n.

Programs

  • Maple
    A000795 := proc(n)
            (2*n)!*coeftayl( cosh(x)/cos(x),x=0,2*n) ;
    end proc: # R. J. Mathar, Oct 20 2011
  • Mathematica
    max = 16; se = Series[ Cosh[x] / Cos[x], {x, 0, 2*max} ]; a[n_] := SeriesCoefficient[ se, 2*n ]*(2*n)!; Table[ a[n], {n, 0, max} ] (* Jean-François Alcover, Apr 02 2012 *)
    With[{nn=40},Take[CoefficientList[Series[Cosh[x]/Cos[x],{x,0,nn}],x] Range[ 0,nn]!,{1,-1,2}]] (* Harvey P. Dale, May 11 2012 *)
    a[ n_] := If[ n < 0, 0, With[ {m = 2 n}, m! SeriesCoefficient[ Cosh[ x] / Cos[ x], {x, 0, m}]]]; (* Michael Somos, Aug 15 2015 *)
  • Sage
    # Generalized algorithm of L. Seidel (1877)
    def A000795_list(n) :
        R = []; A = {-1:0, 0:0}
        k = 0; e = 1
        for i in range(n) :
            Am = 1 if e == 1 else 0
            A[k + e] = 0
            e = -e
            for j in (0..i) :
                Am += A[k]
                A[k] = Am
                k += e
            if e == -1 : R.append(A[-i//2])
        return R
    A000795_list(10) # Peter Luschny, Jun 02 2012

Formula

a(n) = Sum_{k=0..n} binomial(2n, 2k)*A000364(n-k). - Philippe Deléham, Dec 16 2003
a(n) = Sum_{k>=0} (-1)^(n+k)*2^(2n-k)*A065547(n, k). - Philippe Deléham, Feb 26 2004
a(n) = Sum_{k>=0} A086646(n, k). - Philippe Deléham, Feb 26 2004
G.f.: 1 / (1 - (1^2+1)*x / (1 - 2^2*x / (1 - (3^2+1)*x / (1 - 4^2*x / (1 - (5^2+1)*x / (1 - 6^2*x / ...)))))). - Michael Somos, May 12 2012
G.f.: Q(0)/(1-2*x), where Q(k) = 1 - 8*x^2*(2*k^2+2*k+1)*(k+1)^2/( 8*x^2*(2*k^2+2*k+1)*(k+1)^2 - (1 - 8*x*k^2 - 4*x*k -2*x)*(1 - 8*x*k^2 - 20*x*k -14*x)/Q(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Oct 22 2013
a(n) ~ (2*n)! * 2^(2*n+2) * cosh(Pi/2) / Pi^(2*n+1). - Vaclav Kotesovec, Mar 08 2014
a(n) = 1 - Sum_{k=1..n} (-1)^k * binomial(2*n,2*k) * a(n-k). - Ilya Gutkovskiy, Mar 09 2022

A006846 Hammersley's polynomial p_n(1).

Original entry on oeis.org

1, 1, 2, 7, 41, 376, 5033, 92821, 2257166, 69981919, 2694447797, 126128146156, 7054258103921, 464584757637001, 35586641825705882, 3136942184333040727, 315295985573234822561, 35843594275585750890976, 4575961401477587844760793, 651880406652100451820206941
Offset: 0

Views

Author

Keywords

Comments

Equals column 0 of triangle A104027. Also equals column 0 of triangle A104030 (offset 1). Both A104027 and A104030 involve the trinomial coefficients. - Paul D. Hanna, Mar 06 2005

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Julia
    function A006846list(len::Int)  # Algorithm of L. Seidel (1877)
        R = Array{BigInt}(len)
        A = fill(BigInt(0), len+1); A[1] = 1
        for n in 1:len
            for k in n:-1:2 A[k] += A[k+1] end
            for k in 2:1:n A[k] += A[k-1] end
            R[n] = A[n]
        end
        return R
    end
    println(A006846list(20)) # Peter Luschny, Jan 02 2018
  • Maple
    A006846 := proc(n)
        option remember ;
        if n =0 then
            return 1;
        else
            add(binomial(2*n,2*m)*procname(m)/(-4)^(n-m),m=0..n-1) ;
            (3/4)^n-% ;
        end if
    end proc:
    seq(A006846(n),n=0..20) ; # R. J. Mathar, Jan 10 2018
  • Mathematica
    h[n_, x_] := Sum[c[k] x^k, {k, 0, n}]; eq[n_] := SolveAlways[h[n, x*(x-1)] == EulerE[2*n, x], x]; a[n_] := Sum[(-1)^(n+k)*c[k], {k, 0, n}] /. eq[n] // First; Table[a[n], {n, 0, 15}]   (* Jean-François Alcover, Oct 02 2013, after Philippe Deléham *)
  • PARI
    {a(n)=local(X=x+x*O(x^(2*n))); round((2*n)!*polcoeff(cosh(sqrt(3)*X/2)/cos(X/2),2*n))} \\ Paul D. Hanna
    

Formula

a(n) = Sum_{k>=0} (-1)^(n+k)*A065547(n, k) = Sum_{k>=0} A085707(n, k). - Philippe Deléham, Feb 26 2004
E.g.f.: cosh(sqrt(3)*x/2)/cos(x/2) = Sum_{n>=0} a(n)*x^(2n)/(2n)!. - Paul D. Hanna, Feb 27 2005
a(n) = (-1)^n*A104027(n, 0). a(n+1) = (-1)^(n+1)*A104030(n, 0). - Paul D. Hanna, Mar 06 2005
G.f.: 1/(1-x/(1-x/(1-3x/(1-4x/(1-7x/(1-.../(1-ceiling((n+1)^2/4)*x/(1-... (continued fraction). - Paul Barry, Feb 24 2010
a(n) ~ 4*cosh(sqrt(3)*Pi/2) * (2*n)! / Pi^(2*n+1). - Vaclav Kotesovec, Jun 07 2021

A098435 Triangle of Salie numbers T(n,k) for negative n,k, n < k.

Original entry on oeis.org

1, -1, 1, 2, -3, 1, -8, 13, -6, 1, 56, -92, 45, -10, 1, -608, 1000, -493, 115, -15, 1, 9440, -15528, 7662, -1799, 245, -21, 1, -198272, 326144, -160944, 37817, -5180, 462, -28, 1, 5410688, -8900224, 4392080, -1032088, 141465, -12684, 798, -36, 1
Offset: 1

Views

Author

Ralf Stephan, Sep 08 2004

Keywords

Comments

Inverse matrix of A054142. - Paul Barry, Jan 21 2005
Essentially the same as the triangle giving by [0,-1,-1,-4,-4,-9,-9,-16,-16,-25,...] DELTA[1,0,1,0,1,0,1,0,1,0,...] = 1; 0,1; 0,-1,1; 0,2,-3,1; 0,-8,13,-6,1; 0,56,-92,45,-10,1; ... where DELTA is the operator defined in A084938. - Philippe Deléham, Aug 30 2006

Examples

			   1;
  -1,   1;
   2,  -3,   1;
  -8,  13,  -6,   1;
  56, -92,  45, -10,  1;
		

Crossrefs

T(-1, k) = (-1)^k*A005439(k-1). Row sums are zero.

Programs

  • Mathematica
    rows = 9; A054142 = Table[ PadRight[ Table[ Binomial[2*n-k, k], {k, 0, n}], rows], {n, 0, rows-1}]; inv = Inverse[A054142]; Table[ Take[inv[[n]], n], {n, 1, rows}] // Flatten (* Jean-François Alcover, Oct 02 2013, after Paul Barry *)

Formula

See A065547 for formulas.

A098157 Triangle T(n,k) with diagonals T(n,n-k) = binomial(n+1,2k).

Original entry on oeis.org

1, 1, 1, 0, 3, 1, 0, 1, 6, 1, 0, 0, 5, 10, 1, 0, 0, 1, 15, 15, 1, 0, 0, 0, 7, 35, 21, 1, 0, 0, 0, 1, 28, 70, 28, 1, 0, 0, 0, 0, 9, 84, 126, 36, 1, 0, 0, 0, 0, 1, 45, 210, 210, 45, 1, 0, 0, 0, 0, 0, 11, 165, 462, 330, 55, 1, 0, 0, 0, 0, 0, 1, 66, 495, 924, 495, 66, 1, 0, 0, 0, 0, 0, 0, 13, 286, 1287, 1716, 715, 78, 1
Offset: 0

Views

Author

Paul Barry, Aug 29 2004

Keywords

Comments

Row sums are A000079. Diagonal sums are A062200. Inverse is A065547, less the first column.
Number of permutations of length n avoiding simultaneously the patterns 123 and 132 with k descents. A descent in a permutation a(1)a(2)...a(n) is position i such that a(i)>a(i+1). - Tian Han, Nov 16 2023

Examples

			Rows begin:
 {1},
 {1,1},
 {0,3,1},
 {0,1,6,1},
 {0,0,5,10,1},
 {0,0,1,15,15,1},
  ...
		

Crossrefs

Programs

  • Mathematica
    Table[Binomial[n+1, 2(n-k)],{n,0,11},{k,0,n}]//Flatten (* Stefano Spezia, Nov 16 2023 *)

Formula

T(n, k) = binomial(n+1, 2(n-k)) with 0 <= k <= n.
G.f.: (1 + x - q*x)/(1 - 2*q*x - q*x^2 + q^2*x^2). - Tian Han, Nov 16 2023
Previous Showing 11-15 of 15 results.