login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A064551 Ado [Simone Caramel]'s Fibonacci function: define the Fibonacci sequence by f(0) = 1, f(1) = 1, f(n) = f(n-1)+f(n-2); then a(0) = 1, a(n) = a(n-1) + 2*(f(n)-n), n > 0. 3
1, 1, 1, 1, 3, 9, 23, 51, 103, 195, 353, 619, 1061, 1789, 2981, 4925, 8087, 13221, 21547, 35039, 56891, 92271, 149541, 242231, 392233, 634969, 1027753, 1663321, 2691723, 4355745, 7048223, 11404779, 18453871, 29859579, 48314441, 78175075 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,5

COMMENTS

A Pickover sequence with properties analogous to the primes.

REFERENCES

Ado [Simone Caramel], Postings in egroups and newsgroups.

LINKS

T. D. Noe, Table of n, a(n) for n = 0..500

Index entries for linear recurrences with constant coefficients, signature (4,-5,1,2,-1).

FORMULA

G.f.: (1 - 3x + 2x^2 + x^3 + x^4)/((x-1)^3 (x^2 + x - 1)). - T. D. Noe, Oct 12 2007

a(n) = 4*a(n-1) - 5*a(n-2) + a(n-3) + 2*a(n-4) - a(n-5). - T. D. Noe, Oct 12 2007

a(n) = (1/5)*2^(-n)*(-15*2^n + (10-4*sqrt(5))*(1-sqrt(5))^n + (1+sqrt(5))^n*(10+4*sqrt(5))) - n - n^2. - Jean-François Alcover, May 28 2013

a(n) = a(n-1) - 2 * A065220(n), n > 0. - Reinhard Zumkeller, Sep 13 2013

a(n) = 2*F(n+3) - n^2 - n - 3 = 1 + 2*Sum_{k=1..n} F(k+1) - k = 1 + 2*Sum_{k=1..n} A001924(k-3), F=A000045. - Ehren Metcalfe, Dec 27 2018

MAPLE

a := proc(n) option remember: f := n->fibonacci(n+1): if n=0 then RETURN(1) fi: a(n-1)+2*(f(n)-n): end: for n from 1 to 60 do printf(`%d, `, a(n)) od:

MATHEMATICA

a[0] = f[0] = f[1] = 1; f[n_] := f[n] = f[n - 1] + f[n - 2]; a[n_] := a[n] = a[n - 1] + 2*(f[n] - n); Table[ a[n], {n, 0, 40} ]

LinearRecurrence[{4, -5, 1, 2, -1}, {1, 1, 1, 1, 3}, 50] (* Harvey P. Dale, Sep 27 2011 *)

PROG

(ARIBAS): function a064551(maxarg: integer); var n, r, rm, q, qm1, qm2: integer; begin qm2 := 0; qm1 := 0; rm := 0; for n := 0 to maxarg do if n < 2 then q := 1; else q := qm1 + qm2; end; qm2 := qm1; qm1 := q; if n = 0 then r := 1; else r := rm + 2*(q - n); end; rm := r; write(r, " "); end; end; a064551(35).

(Haskell)

a064551 n = a064551_list !! n

a064551_list = 1 : zipWith (+) a064551_list

(map (* 2) $ zipWith (-) (drop 2 a000045_list) [1..])

-- Reinhard Zumkeller, Sep 13 2013

CROSSREFS

Cf. A064550, A064552, A000045.

Sequence in context: A230119 A146693 A227259 * A005783 A263330 A146440

Adjacent sequences: A064548 A064549 A064550 * A064552 A064553 A064554

KEYWORD

nonn,nice,easy

AUTHOR

Roger L. Bagula, Oct 08 2001

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 30 08:53 EST 2023. Contains 359942 sequences. (Running on oeis4.)