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.

A166863 a(1)= 1; a(2)= 5; thereafter a(n)= a(n-1) + a(n-2) + 5.

Original entry on oeis.org

1, 5, 11, 21, 37, 63, 105, 173, 283, 461, 749, 1215, 1969, 3189, 5163, 8357, 13525, 21887, 35417, 57309, 92731, 150045, 242781, 392831, 635617, 1028453, 1664075, 2692533, 4356613, 7049151, 11405769, 18454925, 29860699, 48315629, 78176333, 126491967
Offset: 1

Views

Author

Geoff Ahiakwo, Oct 22 2009

Keywords

Examples

			a(3) = 5 + 1 + 5 = 11.
		

Crossrefs

Programs

  • Haskell
    a166863 n = a166863_list !! (n-1)
    a166863_list = 1 : zipWith (+) a166863_list (drop 3 $ map (* 2) a000045_list)
    -- Reinhard Zumkeller, Nov 17 2013
  • Mathematica
    2 * Fibonacci[Range[4,4! ]] - 5 (* Vladimir Joseph Stephan Orlovsky, Mar 19 2010 *)
    RecurrenceTable[{a[1]==1,a[2]==5,a[n]==a[n-1]+a[n-2]+5},a,{n,40}] (* or *) LinearRecurrence[{2,0,-1},{1,5,11},40] (* Harvey P. Dale, Jan 29 2021 *)

Formula

a(n) = A154691(n) - 2 = 2*A000045(n+3) - 5. - R. J. Mathar, Oct 26 2009
From R. J. Mathar, Oct 26 2009: (Start)
a(n) = 2*a(n-1) - a(n-3).
G.f: x*(1+3*x+x^2)/((x-1)* (x^2+x-1)). (End)
a(n+1) = a(n) + 2*A000045(n+2). - Reinhard Zumkeller, Nov 17 2013

Extensions

Missing value for a(29) inserted by Reinhard Zumkeller, Nov 17 2013