|
|
A155046
|
|
List of pairs: first pair is (1,1); then follow (x,y) with (x+2y, x+y).
|
|
2
|
|
|
1, 1, 3, 2, 7, 5, 17, 12, 41, 29, 99, 70, 239, 169, 577, 408, 1393, 985, 3363, 2378, 8119, 5741, 19601, 13860, 47321, 33461, 114243, 80782, 275807, 195025, 665857, 470832, 1607521, 1136689, 3880899, 2744210, 9369319, 6625109, 22619537, 15994428
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,3
|
|
COMMENTS
|
List of pairs (a, b) such that (a, b*sqrt(2)) = (1 + sqrt(2))^n. In the commutative ring Z[sqrt(2)], the set { +/- (1 + sqrt(2)^n} is a multiplicative group. - Michel Lagneau, Nov 27 2015
The fractions a(2*n-1)/a(2*n) are successive convergents of the simple continued fraction of sqrt(2). - Alexander Fraebel, Sep 03 2020
|
|
LINKS
|
|
|
FORMULA
|
a(n) = 2*a(n-2) + a(n-4) for n > 4. - R. J. Mathar, Feb 19 2009
|
|
MATHEMATICA
|
LinearRecurrence[{0, 2, 0, 1}, {1, 1, 3, 2}, 40] (* Vincenzo Librandi, Mar 14 2012 *)
NestList[{#[[1]]+2#[[2]], Total[#]}&, {1, 1}, 20]//Flatten (* Harvey P. Dale, Nov 21 2020 *)
|
|
PROG
|
(Haskell)
import Data.List (transpose)
a155046 n = a155046_list !! n
a155046_list = concat $ transpose [tail a001333_list, tail a000129_list]
(PARI) Vec(x*(1+x+x^2)/(1-2*x^2-x^4) + O(x^50)) \\ Michel Marcus, Nov 28 2015
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy,tabf
|
|
AUTHOR
|
|
|
EXTENSIONS
|
First term in two pairs corrected by R. J. Mathar, Feb 19 2009
|
|
STATUS
|
approved
|
|
|
|