|
Search: id:A096796
|
|
|
| A096796 |
|
a(n) = n for n <= 2; for n > 2, a(n) = 2a(n-1) - a(n - floor( 1/2 + sqrt(2n) )). |
|
+0
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; listen)
|
|
|
OFFSET
|
0,3
|
|
|
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}] (from Robert G. Wilson v)
|
|
PROGRAM
|
(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
|
|
CROSSREFS
|
Cf. A005318, A096824.
Sequence in context: A052968 A018067 A068202 this_sequence A125730 A074030 A024318
Adjacent sequences: A096793 A096794 A096795 this_sequence A096797 A096798 A096799
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
njas, Aug 17 2004
|
|
EXTENSIONS
|
More terms from Klaus Brockhaus (klaus-brockhaus(AT)t-online.de) and Robert G. Wilson v (rgwv(AT)rgwv.com), Aug 20 2004
|
|
|
Search completed in 0.002 seconds
|
