login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A096796 a(n) = n for n <= 2; for n > 2, a(n) = 2a(n-1) - a(n - floor( 1/2 + sqrt(2n) )). 3
0, 1, 2, 3, 5, 8, 13, 23, 41, 74, 135, 257, 491, 941, 1808, 3481, 6827, 13397, 26303, 51665, 101522, 199563, 395645, 784463, 1555529, 3084755, 6117845, 12134168, 24068773, 47937983, 95480321, 190176179, 378796829, 754508903, 1502899961 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

LINKS

Reinhard Zumkeller, Table of n, a(n) for n = 0..1000

MATHEMATICA

a[n_] := a[n] = If[n < 3, a[n] = n, 2a[n - 1] - a[n - Floor[1/2 + Sqrt[2n]] ]]; Table[ a[n], {n, 0, 35}] (* Robert G. Wilson v *)

PROG

(PARI) {m=35; v=vector(m+1); for(n=0, m, if(n<=2, a=n, k=n-floor(1/2+sqrt(2*n)); a=2*v[n]-v[k+1]); v[n+1]=a; print1(a, ", "))} \\ Klaus Brockhaus, Aug 20 2004

(Haskell)

a096796 n = a096796_list !! n

a096796_list = 0 : 1 : zipWith (-)

   (map (* 2) $ tail a096796_list) (map a096796 $ tail a083920_list)

-- Reinhard Zumkeller, Feb 12 2012

CROSSREFS

Cf. A005318, A096824.

Cf. A083920, A003056.

Sequence in context: A206720 A018067 A068202 * A257418 A125730 A074030

Adjacent sequences:  A096793 A096794 A096795 * A096797 A096798 A096799

KEYWORD

nonn,easy

AUTHOR

N. J. A. Sloane, Aug 17 2004

EXTENSIONS

More terms from Klaus Brockhaus and Robert G. Wilson v, Aug 20 2004

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 16 21:46 EDT 2018. Contains 313809 sequences. (Running on oeis4.)