|
|
A115637
|
|
A divide-and-conquer sequence.
|
|
3
|
|
|
1, 0, 5, 4, 1, 0, 21, 20, 17, 16, 5, 4, 1, 0, 85, 84, 81, 80, 69, 68, 65, 64, 21, 20, 17, 16, 5, 4, 1, 0, 341, 340, 337, 336, 325, 324, 321, 320, 277, 276, 273, 272, 261, 260, 257, 256, 85, 84, 81, 80, 69, 68, 65, 64, 21, 20, 17, 16, 5, 4, 1, 0, 1365, 1364, 1361, 1360, 1349
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (1/(1-x))*Sum_{k>=0} 4^k*x^(2^(k+1)-2)/(1+x^(2^k)); the g.f. G(x) satisfies G(x) - 4(1+x)*x^2*G(x^2) = 1/(1-x^2).
|
|
PROG
|
(PARI)
up_to = 1024;
A115633array(n, k) = (((-1)^n)*if(n==k, 1, if((k+k+2)==n, -4, if((k+1)==n, -(1+(-1)^k)/2, 0))));
A115637list(up_to) = { my(mA115633=matrix(up_to, up_to, n, k, A115633array(n-1, k-1)), mA115636 = matsolve(mA115633, matid(up_to)), v = vector(up_to)); for(n=1, up_to, v[n] = vecsum(mA115636[n, ])); (v); };
v115637 = A115637list(up_to+1);
|
|
CROSSREFS
|
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|