|
| |
|
|
A096824
|
|
a(n) = n for n <= 2; for n > 2, a(n) = 2a(n-1) - a(n - floor( 1/2 + sqrt(2(n-1)) )).
|
|
3
| |
|
|
0, 1, 2, 3, 4, 6, 9, 14, 24, 42, 75, 136, 258, 492, 942, 1809, 3482, 6828, 13398, 26304, 51666, 101523, 199564, 395646, 784464, 1555530, 3084756, 6117846, 12134169, 24068774, 47937984, 95480322, 190176180, 378796830, 754508904, 1502899962
(list; graph; refs; listen; history; 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[2(n - 1)]]]]; Table[ a[n], {n, 0, 35}] (from Robert G. Wilson v Aug 20 2004)
|
|
|
PROG
| (PARI) {m=36; v=vector(m+1); for(n=0, m, if(n<=2, a=n, k=n-floor(1/2+sqrt(2*(n-1))); a=2*v[n]-v[k+1]); v[n+1]=a; print1(a, ", "))} - Klaus Brockhaus, Aug 20 2004
(Haskell)
a096824 n = a096824_list !! n
a096824_list = 0 : 1 : 2 : zipWith (-)
(map (* 2) $ drop 2 a096824_list) (map a096824 $ tail a122797_list)
-- Reinhard Zumkeller, Feb 12 2012
|
|
|
CROSSREFS
| Cf. A005318, A096796.
Cf. A122797, A003056.
Sequence in context: A001115 A173278 A173289 * A192488 A089797 A081237
Adjacent sequences: A096821 A096822 A096823 * A096825 A096826 A096827
|
|
|
KEYWORD
| nonn,easy,changed
|
|
|
AUTHOR
| N. J. A. Sloane (njas(AT)research.att.com), 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
|
| |
|
|
