|
Search: id:A108165
|
|
|
| A108165 |
|
Let beta = positive real root of x^3-x^2-x-1; b(n) = 1 + ceiling((n-1)*beta); c(n)=b(n)-b(n-1); a(n)=a(n-1)+c(n). |
|
+0 1
|
|
| 2, 5, 9, 12, 15, 19, 22, 26, 29, 32, 36, 39, 42, 46, 49, 53, 56, 59, 63, 66, 70, 73, 76, 80, 83, 86, 90, 93, 97, 100, 103, 107, 110, 114, 117, 120, 124, 127, 130, 134, 137, 141, 144, 147, 151, 154, 157, 161, 164, 168, 171, 174, 178, 181, 185, 188, 191, 195, 198, 201
(list; graph; listen)
|
|
|
OFFSET
|
0,1
|
|
|
COMMENT
|
Tribonacci version of A001950 using beta, the positive real root of x^3-x^2-x-1, as the constant.
|
|
MATHEMATICA
|
NSolve[x^3 - x^2 - x - 1 = 0, x] beta = 1.8392867552141612 a[n_] = 1 + Ceiling[(n - 1)*beta^2] (* A007066-like*) aa = Table[a[n], {n, 1, 100}] (*A076662-like*) b = Table[a[n] - a[n - 1], {n, 2, Length[aa]}] F[1] = 2; F[n_] := F[n] = F[n - 1] + b[[n]] (* A000195-like*) c = Table[F[n], {n, 1, Length[b] - 1}]
|
|
CROSSREFS
|
Cf. A007066, A076662, A001950.
Sequence in context: A050904 A057471 A099434 this_sequence A063957 A047385 A086814
Adjacent sequences: A108162 A108163 A108164 this_sequence A108166 A108167 A108168
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Roger L. Bagula (rlbagulatftn(AT)yahoo.com), Jun 13 2005
|
|
EXTENSIONS
|
Partially edited by N. J. A. Sloane (njas(AT)research.att.com), May 04 2007
|
|
|
Search completed in 0.002 seconds
|