|
Search: id:A101914
|
|
|
| A101914 |
|
G.f. satisfies: A(x) = 1/(1 + x*A(x^4)) and also the continued fraction: 1+x*A(x^5) = [1;1/x,1/x^4,1/x^16,1/x^64,...,1/x^(4^(n-1)),...]. |
|
+0
3
|
|
| 1, -1, 1, -1, 1, 0, -1, 2, -3, 3, -2, 0, 3, -6, 8, -8, 5, 1, -9, 17, -22, 20, -10, -8, 31, -51, 60, -50, 16, 38, -100, 150, -163, 119, -11, -147, 315, -432, 433, -268, -70, 522, -964, 1222, -1118, 542, 484, -1756, 2887, -3385, 2793, -879, -2176, 5678, -8472, 9186, -6672, 542, 8372, -17816, 24384, -24350, 14952
(list; graph; listen)
|
|
|
OFFSET
|
0,8
|
|
|
COMMENT
|
Sequence appears to have a rational g.f. - Ralf Stephan, May 17 2007
|
|
PROGRAM
|
(PARI) {a(n)=local(A); A=1-x; for(i=1, n\4+1, A=1/(1+x*subst(A, x, x^4)+x*O(x^n))); polcoeff(A, n, x)} (PARI) {a(n)=local(M=contfracpnqn(concat(1, vector(ceil(log(n+1)/log(4))+1, n, 1/x^(4^(n-1)))))); polcoeff(M[1, 1]/M[2, 1]+x*O(x^(5*n+1)), 5*n+1)}
|
|
CROSSREFS
|
Cf. A101912-A101913, A101915-A101918.
Sequence in context: A077869 A076585 A022906 this_sequence A099530 A074989 A123548
Adjacent sequences: A101911 A101912 A101913 this_sequence A101915 A101916 A101917
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
Paul D. Hanna (pauldhanna(AT)juno.com), Dec 20 2004
|
|
|
Search completed in 0.002 seconds
|
