|
%I #7 Dec 17 2022 06:51:08
%S 0,1,-1,1,0,0,1,0,1,0,1,1,1,1,2,1,2,3,2,4,3,3,5,4,4,5,4,5,4,5,5,5,5,6,
%T 5,6,7,6,8,7,8,6,8,7,7,8,7,8,7,8,8,8,8,9,8,9,10,9,11,10,11,9,11,10,10,
%U 11,10,11,10,11,11,11,11,12,11,12,13,12,14,13,14,12,14,13,13,14,13,14,13
%N a(0) = 0, a(n) = A029826(n+1) for n <= 20, otherwise a(n) = a(n - 1 - n % 20) + A029826(2 + n % 20 ).
%t p[x_] = (x^(10) + x^9 - x^7 - x^6 - x^5 - x^4 - x^3 + x + 1);
%t q[x_] = Expand[x^10*p[1/x]];
%t a = Table[SeriesCoefficient[Series[1/q[x], {x, 0, 50}], n], {n, 0, 20}];
%t b[0] := 0;
%t b[n_] := b[n] = If[n <= 20, a[[n]], b[n - 1 - Mod[n, 20]] + a[[1 + Mod[n, 20]]]];
%t Table[b[n], {n, 0, 100}]
%Y Cf. A029826.
%K sign,base,less
%O 0,15
%A _Roger L. Bagula_, Feb 06 2010
|
