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.

A179190 Coefficient [x^n] of the Maclaurin series for 2 - sqrt(1 - 4*x - 4*x^2).

Original entry on oeis.org

1, 2, 4, 8, 24, 80, 288, 1088, 4256, 17088, 70016, 291584, 1230592, 5251584, 22623232, 98248704, 429677056, 1890700288, 8364824576, 37186449408, 166030266368, 744180244480, 3347321831424, 15104525959168, 68357598756864
Offset: 0

Views

Author

Clark Kimberling, Jul 01 2010

Keywords

Examples

			The Maclaurin series is 1 + 2*x + 4*x^2 + 8*x^3 + 24*x^4 + ...
		

Crossrefs

Programs

  • Maple
    A179190 := proc(n) if n = 0 then 1; else add( doublefactorial(2*n-2*k-3) *2^(n-k) / k! / (n-2*k)!, k=0..floor(n/2)) ; end if; end proc: # R. J. Mathar, Jul 11 2011
  • Mathematica
    Table[SeriesCoefficient[Series[2-Sqrt[1-4*t-4*t^2], {t,0,n}], n], {n, 0, 30}] (* G. C. Greubel, Jan 25 2019 *)
  • Maxima
    makelist(coeff(taylor(2-sqrt(1-4*x-4*x^2), x, 0, n), x, n), n, 0, 24); /* Bruno Berselli, Jul 04 2011 */

Formula

G.f.: 2 - sqrt(1 - 4*x - 4*x^2).
a(n) = 4*A071356(n-2), n >= 2. - R. J. Mathar, Jul 08 2010
a(n) = Sum_{k=0..floor(n/2)} (2*n - 2*k - 3)!! *2^(n-k)/(k!*(n-2k)!), n > 0. - R. J. Mathar, Jul 11 2011
a(n) ~ 2^(n - 1/4) * (1 + sqrt(2))^(n - 1/2) / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Jan 26 2019
D-finite with recurrence: n*a(n) +2*(-2*n+3)*a(n-1) +4*(-n+3)*a(n-2)=0. - R. J. Mathar, Jan 20 2020

A260774 Certain directed lattice paths.

Original entry on oeis.org

1, 6, 33, 189, 1107, 6588, 39663, 240894, 1473147, 9058554, 55954395, 346934745, 2157989445, 13459891500, 84152389833, 527224251861, 3309194474451, 20804569738218, 130987600581699, 825796890644895, 5212349717906889, 32935490120006604, 208316726580941037
Offset: 0

Views

Author

N. J. A. Sloane, Jul 30 2015

Keywords

Comments

See Dziemianczuk (2014) for precise definition.

Crossrefs

Programs

  • Maple
    b:= proc(x, y) option remember; `if`([x, y]=[0$2], 1,
          `if`(x>0, add(b(x-1, y+j), j=-1..1), 0)+
          `if`(y>0, b(x, y-1), 0)+`if`(y<0, b(x, y+1), 0))
        end:
    a:= n-> b(n, 1):
    seq(a(n), n=0..23);  # Alois P. Heinz, Sep 21 2021
  • Mathematica
    b[x_, y_] := b[x, y] = If[{x, y} == {0, 0}, 1,
         If[x > 0, Sum[b[x - 1, y + j], {j, -1, 1}], 0] +
         If[y > 0, b[x, y - 1], 0] + If[y < 0, b[x, y + 1], 0]];
    a[n_] := b[n, 1];
    Table[a[n], {n, 0, 23}] (* Jean-François Alcover, May 02 2022, after Alois P. Heinz *)

Formula

See Dziemianczuk (2014) Equation (33a) with m=1.
From Vaclav Kotesovec, Jul 15 2022: (Start)
Recurrence: (n+1)*(4*n - 3)*a(n) = 6*(4*n^2 - n - 1)*a(n-1) + 3*(n-1)*(4*n + 1)*a(n-2).
a(n) ~ (3 + 2*sqrt(3))^(n+1) / sqrt(6*Pi*n). (End)

Extensions

More terms from Lars Blomberg, Aug 01 2015

A260772 Certain directed lattice paths.

Original entry on oeis.org

1, 3, 10, 41, 190, 946, 4940, 26693, 147990, 837102, 4811860, 28027210, 165057100, 981177060, 5879570200, 35478788269, 215398416870, 1314794380374, 8064119033220, 49673222082782, 307163049317540, 1906066361809148, 11865666767361960, 74081851132379426
Offset: 0

Views

Author

N. J. A. Sloane, Jul 30 2015

Keywords

Comments

See Dziemianczuk (2014) for precise definition.

Crossrefs

Programs

  • Maple
    # A260772 satisfies a 4th-order recurrence that can be reduced
    # to a 2nd-order recurrence given in this program t:
    t := proc(n) options remember;
    if n <= 1 then
        [-1/2, 0, 1, 4][2*n+2]
      else
        (16*(n-2)*(2*n-3)*(5*n-2)*t(n-2) + (440*n^3-1056*n^2+724*n-144)*t(n-1))
           /( n*(2*n+1)*(5*n-7) )
      fi
    end:
    A260772 := proc(n)
    t(n/2) + ( (2-2*n)*t((n-1)/2)+(n+2)*t((n+1)/2) ) / (1+5*n)
    end:
    seq(A260772(i),i=0..100);
    # Mark van Hoeij, Jul 14 2022
  • Maxima
    a(n):=if n=0 then 1 else sum((-1)^j*binomial(n,j)*binomial(3*n-4*j,n-4*j+1),j,0,(n+1)/4)/n; /* Vladimir Kruchinin, Apr 04 2019 */
    
  • PARI
    a(n) = if (n==0, 1, sum(j=0, (n+1)/4, (-1)^j*binomial(n,j)*binomial(3*n-4*j, n-4*j+1))/n); \\ Michel Marcus, Apr 05 2019

Formula

G.f.: P1(x) = (2*(1-x)/3)/x - ((2*sqrt(1-5*x-2*x^2)/3)/x)*sin((Pi/6 + arccos(((20*x^3-6*x^2+15*x-2)/2)/(1-5*x-2*x^2)^(3/2))/3)). - See Dziemianczuk (2014), Proposition 11.
a(n) = (1/n)*Sum_{j=0..(n+1)/4} (-1)^j*C(n,j)*C(3*n-4*j,n-4*j+1), a(0)=1. - Vladimir Kruchinin, Apr 04 2019
n*(n+1)*(25*n^2-70*n+21)*a(n) - 30*(7*n-15)*n*a(n-1) + (-1100*n^4+5280*n^3-6424*n^2-1188*n+3816)*a(n-2) + 120*(n+2)*(n-3)*a(n-3) - 16*(n-3)*(n-4)*(25*n^2-20*n-24)*a(n-4) = 0. - Mark van Hoeij, Jul 14 2022
a(n) ~ 2^(n - 1/2) * phi^((10*n - 1)/4) / (sqrt(Pi) * 5^(1/4) * sqrt(phi^(3/2) - 2) * n^(3/2)), where phi = A001622 is the golden ratio. - Vaclav Kotesovec, Jul 15 2022

Extensions

More terms from Lars Blomberg, Aug 01 2015
Showing 1-3 of 3 results.