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.

A001871 Expansion of 1/(1 - 3*x + x^2)^2.

Original entry on oeis.org

1, 6, 25, 90, 300, 954, 2939, 8850, 26195, 76500, 221016, 632916, 1799125, 5082270, 14279725, 39935214, 111228804, 308681550, 853904015, 2355364650, 6480104231, 17786356776, 48715278000, 133167004200, 363372003625, 989900286774
Offset: 0

Views

Author

Keywords

Comments

Convolution of A001906(n), n >= 1 (even-indexed Fibonacci numbers) with itself.
A001787 and this sequence arise in counting ordered trees of height at most k where only the rightmost branch at the root actually achieves this height and the count is by the number of edges, with k = 3 for A001787 and k = 4 for this sequence.
Gives the number of 3412-avoiding permutations containing exactly one subsequence of type 321. - Dan Daly (ddaly(AT)du.edu), Apr 24 2008

References

  • 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

Partial sums of A001870 (one half of odd-indexed A001629(n), n >= 2, Fibonacci convolution).

Programs

  • Magma
    I:=[1, 6, 25, 90]; [n le 4 select I[n] else 6*Self(n-1)-11*Self(n-2)+6*Self(n-3)-Self(n-4): n in [1..30]]; // Vincenzo Librandi, Jun 10 2012
    
  • Maple
    f:= gfun:-rectoproc({a(n)=6*a(n-1)-11*a(n-2)+6*a(n-3)-a(n-4),
    a(0)=1,a(1)=6,a(2)=25,a(3)=90},a(n),remember):
    map(f, [$0..50]); # Robert Israel, May 05 2017
    # alternative
    A001871 := proc(n)
        option remember ;
        if n <= 3 then
            op(n+1,[1,6,25,90]) ;
        else
            6*procname(n-1)-11*procname(n-2)+6*procname(n-3)-procname(n-4) ;
        end if;
    end proc:
    seq(A001871(n),n=0..10) ; # R. J. Mathar, Dec 16 2024
  • Mathematica
    CoefficientList[Series[1/(1-3x+x^2)^2,{x,0,40}],x] (* Vincenzo Librandi, Jun 10 2012 *)
  • PARI
    a(n)=((4*n+2)*fibonacci(2*n)+(7*n+5)*fibonacci(2*n+1))/5
    
  • PARI
    Vec(1/(1-3*x+x^2)^2 + O(x^100)) \\ Altug Alkan, Oct 31 2015

Formula

a(n) = (2*(2*n+1)*F(2*(n+1))+3*(n+1)*F(2*n+1))/5 with F(n) = A000045 (Fibonacci numbers).
a(n) = -a(-4-n) = ((4*n+2)*F(2*n) + (7*n+5)*F(2*n+1))/5 with F(n) = A000045 (Fibonacci numbers).
a(n) = (2*a(n-1) + (n+1)*F(2n+4))/3, where F(n) = A000045 (Fibonacci numbers). - Emeric Deutsch, Oct 08 2002
G.f.: 1/(1 - 3*x + x^2)^2. - Simon Plouffe in his 1992 dissertation
a(n) = (Sum_{k=0..n} S(k, 3)*S(n-k, 3)), where S(n, x) = U(n, x/2) is the n-th Chebyshev polynomial of the 2nd kind, A049310. - Paul Barry, Nov 14 2003
a(n) = Sum_{k=1..n+1} F(2k)*F(2(n-k+2)), where F(k) is the k-th Fibonacci number. - Dan Daly (ddaly(AT)du.edu), Apr 24 2008
a(n) = 6*a(n-1) - 11*a(n-2) + 6*a(n-3) - a(n-4). - Vincenzo Librandi, Mar 14 2011
a(n) = 2*A001870(n) - A238846(n). - Philippe Deléham, Mar 06 2014
a(n) ~ (7 + 3*sqrt(5))*n*cos(n*arccos(3/2))/5. - Stefano Spezia, Mar 29 2022
From Peter Bala, Nov 05 2024: (Start)
a(n) = Sum_{k = 0..n} (n + 2*k + 1)*binomial(n+k, 2*k).
a(n) = (n+1) * hypergeom([-n, n+1, (n+3)/2], [1/2, (n+1)/2], -1/4).
Second-order recurrence: n*a(n) = 3*(n + 1)*a(n-1) - (n + 2)*a(n-2) with a(0) = 1 and a(1) = 6. (End)
E.g.f.: exp(3*x/2)*(5*(5 + 18*x)*cosh(sqrt(5)*x/2) + sqrt(5)*(9 + 40*x)*sinh(sqrt(5)*x/2))/25. - Stefano Spezia, Mar 04 2025

Extensions

Additional comments from Wolfdieter Lang, Apr 07 2000

A001870 Expansion of (1-x)/(1 - 3*x + x^2)^2.

Original entry on oeis.org

1, 5, 19, 65, 210, 654, 1985, 5911, 17345, 50305, 144516, 411900, 1166209, 3283145, 9197455, 25655489, 71293590, 197452746, 545222465, 1501460635, 4124739581, 11306252545, 30928921224, 84451726200, 230204999425
Offset: 0

Views

Author

Keywords

Comments

a(n) = ((n+1)*F(2*n+3)+(2*n+3)*F(2*(n+1)))/5 with F(n)=A000045(n) (Fibonacci numbers). One half of odd-indexed A001629(n), n >= 2 (Fibonacci convolution).
Convolution of F(2n+1) (A001519) and F(2n+2) (A001906(n+1)). - Graeme McRae, Jun 07 2006
Number of reentrant corners along the lower contours of all directed column-convex polyominoes of area n+3 (a reentrant corner along the lower contour is a vertical step that is followed by a horizontal step). a(n) = Sum_{k=0..ceiling((n+1)/2)} k*A121466(n+3,k). - Emeric Deutsch, Aug 02 2006
From Wolfdieter Lang, Jan 02 2012: (Start)
a(n) = A024458(2*n), n >= 1 (bisection, even arguments).
a(n) is also the odd part of the bisection of the half-convolution of the sequence A000045(n+1), n >= 0, with itself. See a comment on A201204 for the definition of the half-convolution of a sequence with itself. There one also finds the rule for the o.g.f. which in this case is Chato(x)/2 with the o.g.f. Chato(x) = 2*(1-x)/(1-3*x+x^2)^2 of A001629(2*n+3), n >= 0.
(End)

References

  • 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

a(n) = A060921(n+1, 1)/2.
Partial sums of A030267. First differences of A001871.
Cf. A121466.
Cf. A023610.

Programs

  • GAP
    F:=Fibonacci;; List([0..30], n-> ((n+1)*F(2*n+3)+(2*n+3)*F(2*(n+1)))/5); # G. C. Greubel, Jul 15 2019
  • Haskell
    a001870 n = a001870_list !! n
    a001870_list = uncurry c $ splitAt 1 $ tail a000045_list where
       c us vs'@(v:vs) = (sum $ zipWith (*) us vs') : c (v:us) vs
    -- Reinhard Zumkeller, Oct 31 2013
    
  • Magma
    I:=[1, 5, 19, 65]; [n le 4 select I[n] else 6*Self(n-1) -11*Self(n-2)+6*Self(n-3)-Self(n-4): n in [1..30]]; // Vincenzo Librandi, Jun 10 2012
    
  • Maple
    A001870:=-(-1+z)/(z**2-3*z+1)**2; # Simon Plouffe in his 1992 dissertation.
  • Mathematica
    CoefficientList[Series[(1-x)/(1-3*x+x^2)^2,{x,0,40}],x] (* Vincenzo Librandi, Jun 10 2012 *)
    LinearRecurrence[{6,-11,6,-1},{1,5,19,65},30] (* Harvey P. Dale, Aug 17 2013 *)
    With[{F=Fibonacci}, Table[((n+1)*F[2*n+3]+(2*n+3)*F[2*n+2])/5, {n,0,30}]] (* G. C. Greubel, Jul 15 2019 *)
  • PARI
    Vec((1-x)/(1-3*x+x^2)^2+O(x^30)) \\ Charles R Greathouse IV, Sep 23 2012
    
  • Sage
    f=fibonacci; [((n+1)*f(2*n+3)+(2*n+3)*f(2*n+2))/5 for n in (0..30)] # G. C. Greubel, Jul 15 2019
    

Formula

a(n) = Sum_{k=1..n+1} k*binomial(n+k+1, 2k). - Emeric Deutsch, Jun 11 2003
a(n) = 6*a(n-1) - 11*a(n-2) + 6*a(n-3) - a(n-4). - Vincenzo Librandi, Jun 10 2012
a(n) = (A238846(n) + A001871(n))/2. - Philippe Deléham, Mar 06 2014
a(n) = ((2*n-1)*Fibonacci(2*n) - n*Fibonacci(2*n-1))/5 [Czabarka et al.]. - N. J. A. Sloane, Sep 18 2018
E.g.f.: exp(3*x/2)*(5*(5 + 11*x)*cosh(sqrt(5)*x/2) + sqrt(5)*(13 + 25*x)*sinh(sqrt(5)*x/2))/25. - Stefano Spezia, Mar 04 2025

Extensions

More terms from Christian G. Bower

A238941 Triangle T(n,k), read by rows given by (1, 1, 1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, -1, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 5, 8, 4, 1, 13, 21, 13, 6, 1, 34, 55, 40, 25, 7, 1, 89, 144, 120, 90, 33, 9, 1, 233, 377, 354, 300, 132, 51, 10, 1, 610, 987, 1031, 954, 483, 234, 62, 12, 1, 1597, 2584, 2972, 2939, 1671, 951, 308, 86, 13, 1, 4181, 6765, 8495, 8850, 5561, 3573, 1345, 480, 100, 15, 1
Offset: 0

Views

Author

Philippe Deléham, Mar 07 2014

Keywords

Comments

Row sums are A025192(n).

Examples

			Triangle begins:
1;
1, 1;
2, 3, 1;
5, 8, 4, 1;
13, 21, 13, 6, 1;
34, 55, 40, 25, 7, 1;
89, 144, 120, 90, 33, 9, 1;
233, 377, 354, 300, 132, 51, 10, 1;
		

Crossrefs

Cf. Columns: A001519, A001906, A238846, A001871.
Cf. Diagonals: A000012, A032766.

Programs

  • Mathematica
    nmax=10; Column[CoefficientList[Series[CoefficientList[Series[(1 - 2*x + x*y)/(1 - 3*x + x^2 - x^2*y^2), {x, 0, nmax }], x], {y, 0, nmax}], y]] (* Indranil Ghosh, Mar 14 2017 *)

Formula

G.f. for the column k: x^k*(1-2*x)^A059841(k)/(1-3*x+x^2)^A008619(k).
G.f.: (1-2*x+x*y)/(1-3*x+x^2-x^2*y^2).
T(n,k) = 3*T(n-1,k) + T(n-2,k-2) - T(n-2,k), T(0,0) = T(1,0) = T(1,1) = 1, T(n,k) = 0 if k<0 or if k>n.
Sum_{k = 0..n} T(n,k)*x^k = A000007(n), A001519(n), A025192(n), A030195(n+1) for x = -1, 0, 1, 2 respectively.
Sum_{k = 0..n} T(n,k)*3^k = A015525(n) + A015525(n+1).

Extensions

Data section corrected and extended by Indranil Ghosh, Mar 14 2017
Showing 1-3 of 3 results.