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

A256890 Triangle T(n,k) = t(n-k, k); t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 2.

Original entry on oeis.org

1, 2, 2, 4, 12, 4, 8, 52, 52, 8, 16, 196, 416, 196, 16, 32, 684, 2644, 2644, 684, 32, 64, 2276, 14680, 26440, 14680, 2276, 64, 128, 7340, 74652, 220280, 220280, 74652, 7340, 128, 256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172, 256, 512, 72076, 1637860, 10978444, 27227908, 27227908, 10978444, 1637860, 72076, 512
Offset: 0

Views

Author

Dale Gerdemann, Apr 12 2015

Keywords

Comments

Related triangles may be found by varying the function f(x). If f(x) is a linear function, it can be parameterized as f(x) = a*x + b. With different values for a and b, the following triangles are obtained:
a\b 1.......2.......3.......4.......5.......6
The row sums of these, and similarly constructed number triangles, are shown in the following table:
a\b 1.......2.......3.......4.......5.......6.......7.......8.......9
The formula can be further generalized to: t(n,m) = f(m+s)*t(n-1,m) + f(n-s)*t(n,m-1), where f(x) = a*x + b. The following table specifies triangles with nonzero values for s (given after the slash).
a\b 0 1 2 3
-2 A130595/1
-1
0
With the absolute value, f(x) = |x|, one obtains A038221/3, A038234/4, A038247/5, A038260/6, A038273/7, A038286/8, A038299/9 (with value for s after the slash).
If f(x) = A000045(x) (Fibonacci) and s = 1, the result is A010048 (Fibonomial).
In the notation of Carlitz and Scoville, this is the triangle of generalized Eulerian numbers A(r, s | alpha, beta) with alpha = beta = 2. Also the array A(2,1,4) in the notation of Hwang et al. (see page 31). - Peter Bala, Dec 27 2019

Examples

			Array, t(n, k), begins as:
   1,    2,      4,        8,        16,         32,          64, ...;
   2,   12,     52,      196,       684,       2276,        7340, ...;
   4,   52,    416,     2644,     14680,      74652,      357328, ...;
   8,  196,   2644,    26440,    220280,    1623964,    10978444, ...;
  16,  684,  14680,   220280,   2643360,   27227908,   251195000, ...;
  32, 2276,  74652,  1623964,  27227908,  381190712,  4677894984, ...;
  64, 7340, 357328, 10978444, 251195000, 4677894984, 74846319744, ...;
Triangle, T(n, k), begins as:
    1;
    2,     2;
    4,    12,      4;
    8,    52,     52,       8;
   16,   196,    416,     196,      16;
   32,   684,   2644,    2644,     684,      32;
   64,  2276,  14680,   26440,   14680,    2276,     64;
  128,  7340,  74652,  220280,  220280,   74652,   7340,   128;
  256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172,   256;
		

Crossrefs

Programs

  • Magma
    A256890:= func< n,k | (&+[(-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n: j in [0..k]]) >;
    [A256890(n,k): k in [0..n], n in [0..10]]; // G. C. Greubel, Oct 18 2022
    
  • Mathematica
    Table[Sum[(-1)^(k-j)*Binomial[j+3, j] Binomial[n+4, k-j] (j+2)^n, {j,0,k}], {n,0, 9}, {k,0,n}]//Flatten (* Michael De Vlieger, Dec 27 2019 *)
  • PARI
    t(n,m) = if ((n<0) || (m<0), 0, if ((n==0) && (m==0), 1, (m+2)*t(n-1, m) + (n+2)*t(n, m-1)));
    tabl(nn) = {for (n=0, nn, for (k=0, n, print1(t(n-k, k), ", ");); print(););} \\ Michel Marcus, Apr 14 2015
    
  • SageMath
    def A256890(n,k): return sum((-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n for j in range(k+1))
    flatten([[A256890(n,k) for k in range(n+1)] for n in range(11)]) # G. C. Greubel, Oct 18 2022

Formula

T(n,k) = t(n-k, k); t(0,0) = 1, t(n,m) = 0 if n < 0 or m < 0 else t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 2.
Sum_{k=0..n} T(n, k) = A001715(n).
T(n,k) = Sum_{j = 0..k} (-1)^(k-j)*binomial(j+3,j)*binomial(n+4,k-j)*(j+2)^n. - Peter Bala, Dec 27 2019
Modified rule of Pascal: T(0,0) = 1, T(n,k) = 0 if k < 0 or k > n else T(n,k) = f(n-k) * T(n-1,k-1) + f(k) * T(n-1,k), where f(x) = x + 2. - Georg Fischer, Nov 11 2021
From G. C. Greubel, Oct 18 2022: (Start)
T(n, n-k) = T(n, k).
T(n, 0) = A000079(n). (End)

A034910 One quarter of octo-factorial numbers.

Original entry on oeis.org

1, 12, 240, 6720, 241920, 10644480, 553512960, 33210777600, 2258332876800, 171633298636800, 14417197085491200, 1326382131865190400, 132638213186519040000, 14324927024144056320000, 1661691534800710533120000, 206049750315288106106880000
Offset: 1

Views

Author

Keywords

Comments

A034910 occurs in connection with the Vandermonde permanent of (1,3,5,7,9,...); see the Mathematica section of A203516. - Clark Kimberling, Jan 03 2012

Examples

			G.f. = x + 12*x^2 + 240*x^3 + 6720*x^4 + 241920*x^5 + 10644480*x^6 + ...
		

Crossrefs

Programs

  • Magma
    [n le 2 select 12^(n-1) else (7*n-3)*Self(n-1) +4*(n-1)*(2*n-3)*Self(n-2): n in [1..30]]; // G. C. Greubel, Oct 20 2022
    
  • Maple
    [seq((2*n)!/(n)!*2^(n-2), n=1..14)]; # Zerinvary Lajos, Sep 25 2006
  • Mathematica
    s=1;lst={s};Do[s+=n*s;AppendTo[lst, s], {n, 11, 5!, 8}];lst (* Vladimir Joseph Stephan Orlovsky, Nov 08 2008 *)
    a[ n_] := Pochhammer[ 1/2, n] 8^n / 4; (* Michael Somos, Feb 04 2015 *)
  • PARI
    {a(n) = if( n==1, 1, n>1, a(n-1) * (8*n - 4), a(n+1) / (8*n + 4))}; /* Michael Somos, Feb 04 2015 */
    
  • SageMath
    [2^(3*n-2)*rising_factorial(1/2, n) for n in range(1,40)] # G. C. Greubel, Oct 20 2022

Formula

4*a(n) = (8*n-4)(!^8) = Product_{j=1..n} (8*j-4) = 4^n*A001147(n) = 2^n*(2*n)!/n!, A001147(n) = (2*n-1)!!.
E.g.f. (-1+(1-8*x)^(-1/2))/4.
a(n) = A090802(2n-1, n). - Ross La Haye, Oct 18 2005
a(n) = ((2*n)!/n!)*2^(n-2). - Zerinvary Lajos, Sep 25 2006
G.f.: x/(1-12*x/(1-8*x/(1-20*x/(1-16*x/(1-28*x/(1-24*x/(1-36*x/(1-32*x/(1-... (continued fraction). - Philippe Deléham, Jan 07 2011
From Peter Bala, Feb 01 2015: (Start)
Recurrence equation: a(n) = (7*n - 3)*a(n-1) + 4*(n - 1)*(2*n - 3)*a(n-2).
The sequence b(n) := a(n)* Sum_{k = 0..n-1} (-1)^k/( 2^k*(2*k + 1)*binomial(2*k,k) ) beginning [1, 11, 222, 6210, 223584, ...] satisfies the same recurrence. This leads to the finite continued fraction expansion b(n)/a(n) = 1/(1 + 1/(11 + 24/(18 + 60/(25 + ... + 4*(n - 1)*(2*n - 3)/(7*n - 3) )))) for n >= 3.
Letting n tend to infinity gives the continued fraction expansion Sum_{k>=0} (-1)^k/( 2^k*(2*k + 1)*binomial(2*k,k) ) = (4/3)*log(2) = 1/(1 + 1/(11 + 24/(18 + 60/(25 + ... + 4*(n - 1)*(2*n - 3)/((7*n - 3) + ... ))))). (End)
From Peter Bala, Feb 03 2015: (Start)
This sequence satisfies several other second order recurrence equations leading to some continued fraction expansions.
1) a(n) = (9*n + 4)*a(n-1) - 4*n*(2*n - 1)*a(n-2).
This recurrence is also satisfied by the (integer) sequence c(n) := a(n)*Sum_{k = 0..n} 1/( 2^k*(2*k + 1)*binomial(2*k,k) ). From this we can obtain the continued fraction expansion Sum_{k >= 0} 1/( 2^k*(2*k + 1)*binomial(2*k,k) ) = (8/sqrt(7))*arctan(sqrt(7)/7) = (8/sqrt(7))*A195699 = 1 + 1/(12 - 24/(22 - 60/(31 - ... - 4*n*(2*n - 1)/((9*n + 4) - ... )))).
2) a(n) = (12*n + 2)*a(n-1) - 8*(2*n - 1)^2*a(n-2).
This recurrence is also satisfied by the (integer) sequence d(n) := a(n)*Sum_{k = 0..n} 1/( (2*k + 1)*2^k ). From this we can obtain the continued fraction expansion Sum_{k >= 0} 1/( (2*k + 1)*2^k ) = (1/sqrt(2))*log(3 + 2*sqrt(2)) = 1 + 2/(12 - 8*3^2/(26 - 8*5^2/(38 - ... - 8*(2*n - 1)^2/((12*n + 2) - ... )))). Cf. A002391.
3) a(n) = (4*n + 6)*a(n-1) + 8*(2*n - 1)^2*a(n-2).
This recurrence is also satisfied by the (integer) sequence e(n) := a(n)*Sum_{k = 0..n} (-1)^k/( (2*k + 1)*2^k ). From this we can obtain the continued fraction expansion Sum_{k >= 0} (-1)^k/( (2*k + 1)*2^k ) = (1/sqrt(2))*arctan(sqrt(2)/2) = 1 - 2/(12 + 8*3^2/(14 + 8*5^2/(18 + ... + 8*(2*n - 1)^2/((4*n + 6) + ... )))). Cf. A073000. (End)
a(n) = (-1)^n / (16*a(-n)) for all n in Z. - Michael Somos, Feb 04 2015
From Amiram Eldar, Jan 08 2022: (Start)
Sum_{n>=1} 1/a(n) = e^(1/8)*sqrt(2*Pi)*erf(1/(2*sqrt(2))), where erf is the error function.
Sum_{n>=1} (-1)^(n+1)/a(n) = e^(-1/8)*sqrt(2*Pi)*erfi(1/(2*sqrt(2))), where erfi is the imaginary error function. (End)

A034909 One third of octo-factorial numbers.

Original entry on oeis.org

1, 11, 209, 5643, 197505, 8492715, 433128465, 25554579435, 1712156822145, 128411761660875, 10658176217852625, 969894035824588875, 96019509546634298625, 10274087521489869952875, 1181520064971335044580625, 145326967991474210483416875, 19037832806883121573327610625
Offset: 1

Views

Author

Keywords

Crossrefs

Formula

3*a(n) = (8*n-5)(!^8) = product(8*j-5, j=1..n).
E.g.f.: (-1+(1-8*x)^(-3/8))/3.
G.f.: x/(1-11x/(1-8x/(1-19x/(1-16x/(1-27x/(1-24x/(1-35x/(1-32x/(1-... (continued fraction). - Philippe Deléham, Jan 07 2012
From Amiram Eldar, Dec 20 2022: (Start)
a(n) = A144756(n)/3.
Sum_{n>=1} 1/a(n) = 3*(e/8^5)^(1/8)*(Gamma(3/8) - Gamma(3/8, 1/8)). (End)

A034911 One fifth of octo-factorial numbers.

Original entry on oeis.org

1, 13, 273, 7917, 292929, 13181805, 698635665, 42616775565, 2940557513985, 226422928576845, 19245948929031825, 1789873250399959725, 180777198290395932225, 19704714613653156612525, 2305451609797419323665425, 288181451224677415458178125, 38328133012882096255937690625
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [n le 1 select 1 else (8*n-3)*Self(n-1): n in [1..40]]; // G. C. Greubel, Oct 20 2022
    
  • Mathematica
    With[{nn=20},CoefficientList[Series[(-1+(1-8x)^(-5/8))/5,{x,0,nn}],x] Range[0,nn]!] (* or *) FoldList[Times,Range[5,200,8]]/5 (* Harvey P. Dale, May 25 2016 *)
  • SageMath
    [8^n*rising_factorial(5/8,n)/5 for n in range(1,40)] # G. C. Greubel, Oct 20 2022

Formula

5*a(n) = (8*n-3)(!^8) = Product_{j=1..n} 8*j-3.
E.g.f.: (-1+(1-8*x)^(-5/8))/5.
G.f.: x/(1-13*x/(1-8*x/(1-21*x/(1-16*x/(1-29*x/(1-24*x/(1-37*x/(1-32*x/(1-... (continued fraction). - Philippe Deléham, Jan 07 2012
D-finite with recurrence: a(n) = (8*n-3)*a(n-1). - R. J. Mathar, Jan 28 2020
a(n) = (1/5)* 8^n * Pochhammer(n, 5/8). - G. C. Greubel, Oct 20 2022
From Amiram Eldar, Dec 20 2022: (Start)
a(n) = A147625(n+1)/5.
Sum_{n>=1} 1/a(n) = 5*(e/8^3)^(1/8)*(Gamma(5/8) - Gamma(5/8, 1/8)). (End)

A034975 One seventh of octo-factorial numbers.

Original entry on oeis.org

1, 15, 345, 10695, 417105, 19603935, 1078216425, 67927634775, 4822862069025, 381006103452975, 33147531000408825, 3149015445038838375, 324348590839000352625, 36002693583129039141375, 4284320536392355657823625, 544108708121829168543600375, 73454675596446937753386050625
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [n le 1 select 1 else (8*n-1)*Self(n-1): n in [1..40]]; // G. C. Greubel, Oct 21 2022
    
  • Mathematica
    Table[8^n*Pochhammer[7/8, n]/7, {n, 40}] (* G. C. Greubel, Oct 21 2022 *)
  • SageMath
    [8^n*rising_factorial(7/8,n)/7 for n in range(1,40)] # G. C. Greubel, Oct 21 2022

Formula

7*a(n) = (8*n-1)!^8 = Product_{j=1..n} (8*j-1) = (8*n)!/((2*n)!*2^(6*n)*3^2*5 * A045755(n)*A007696(n)*A034909(n)*A034911(n)*A034176(n)).
E.g.f.: (-1+(1-8*x)^(-7/8))/7.
G.f.: x/(1-15*x/(1-8*x/(1-23*x/(1-16*x/(1-31*x/(1-24*x/(1-39*x/(1-32*x/(1-... (continued fraction). - Philippe Deléham, Jan 07 2012
a(n) = (1/7) * 8^n * Pochhammer(n, 7/8). - G. C. Greubel, Oct 21 2022
From Amiram Eldar, Dec 20 2022: (Start)
a(n) = A049210(n)/7.
Sum_{n>=1} 1/a(n) = 7*(e/8)^(1/8)*(Gamma(7/8) - Gamma(7/8, 1/8)). (End)

A034976 One eighth of octo-factorial numbers.

Original entry on oeis.org

1, 16, 384, 12288, 491520, 23592960, 1321205760, 84557168640, 6088116142080, 487049291366400, 42860337640243200, 4114592413463347200, 427917611000188108800, 47926772432021068185600, 5751212691842528182272000, 736155224555843607330816000
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [8^(n-1)*Factorial(n): n in [1..40]]; // G. C. Greubel, Oct 20 2022
    
  • Mathematica
    Table[8^(n-1)*n!, {n,40}] (* G. C. Greubel, Oct 20 2022 *)
  • SageMath
    [8^(n-1)*factorial(n) for n in range(1,40)] # G. C. Greubel, Oct 20 2022

Formula

8*a(n) = (8*n)!^8 = Product_{j=1..n} 8*j = 8^n*n!.
E.g.f.: (-1+(1-8*x)^(-1))/8.
G.f.: x/(1-16*x/(1-8*x/(1-24*x/(1-16*x/(1-32*x/(1-24*x/(1-40*x/(1-32*x/(1-... (continued fraction). - Philippe Deléham, Jan 07 2012
From Amiram Eldar, Jan 08 2022: (Start)
Sum_{n>=1} 1/a(n) = 8*(exp(1/8)-1).
Sum_{n>=1} (-1)^(n+1)/a(n) = 8*(1-exp(-1/8)). (End)
Showing 1-6 of 6 results.