|
Search: id:A105309
|
|
|
| A105309 |
|
a(n) = |b(n)|^2 = x^2 + 3*y*2 where (x,y,y,y) is the quaternion b(n) of the sequence b of quaternions defined by b(0)=1,b(1)=1, b(n) = b(n-1) + b(n-2)*(0,c,c,c) where c = 1/sqrt(3). |
|
+0
3
|
|
| 1, 1, 2, 5, 9, 20, 41, 85, 178, 369, 769, 1600, 3329, 6929, 14418, 30005, 62441, 129940, 270409, 562725, 1171042, 2436961, 5071361, 10553600, 21962241, 45703841, 95110562, 197926885, 411889609, 857150100, 1783745641, 3712008565
(list; graph; listen)
|
|
|
OFFSET
|
0,3
|
|
|
COMMENT
|
Prepending 0 and keeping the offset at 0, turns this into a divisibility sequence with g.f. x(1-x^2)/(1-x-2x^2-x^3+x^4). [From T. D. Noe (noe(AT)sspectra.com), Dec 22 2008]
Equals INVERT transform of (1, 1, 2, 0, 2, 0, 2,...) [From Gary W. Adamson (qntmpkt(AT)yahoo.com), Apr 28 2009]
|
|
LINKS
|
Eric Weisstein. "Quaternion"
|
|
FORMULA
|
a(n) = A092886(n+1)-A092886(n-1), n>0.
|
|
PROGRAM
|
(PARI) {a(n)=local(A); if(n<2, n>=0, n++; A=vector(n, i, 1); for(i=3, n, A[i]=A[i-1]+A[i-2]*I); norm(A[n]))} (Somos)
|
|
CROSSREFS
|
Sequence in context: A062092 A079117 A030137 this_sequence A097163 A117186 A155042
Adjacent sequences: A105306 A105307 A105308 this_sequence A105310 A105311 A105312
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Gerald McGarvey (gerald.mcgarvey(AT)comcast.net), Apr 25 2005
|
|
|
Search completed in 0.002 seconds
|
