|
EXAMPLE
|
a(6) = 26 since 26 = 1^2 + 4^2 + 3^2 + 0^2;
a(7) = 66 since 66 = 1^2 + 5^2 + 6^2 + 2^2;
a(8) = 173 since 173 = 1^2 + 6^2 + 10^2 + 6^2 + 0^2;
indicated by this array of coefficients of x^m*A(x^2)^m:
A(x)........: 1,1,1,2,5,11,26,66,173,...
x^0*A(x^2)^0: 1,__0,__0,____0,_____0,...
x^1*A(x^2)^1: __1,__1,___1,____2,____...
x^2*A(x^2)^2: ____1,__2,____3,_____6,...
x^3*A(x^2)^3: ______1,___3,____6,____...
x^4*A(x^2)^4: ________1,____4,____10,...
x^5*A(x^2)^5: ___________1,____5,____...
x^6*A(x^2)^6: ______________1,_____6,...
x^7*A(x^2)^7: _________________1,____...
x^8*A(x^2)^8: _____________________1,...
|
|
PROG
|
(PARI) {a(n)=if(n==0, 1, sum(k=0, n\2, polcoeff(x^(n-2*k)*(sum(j=0, k, a(j)*x^(2*j))+x*O(x^n))^(n-2*k), n)^2))}
|