|
PROGRAM
|
(PARI) {a(n)=local(A); if(n<0, 0, A=sum(k=1, sqrtint(n), 2*x^k^2, 1+x*O(x^n)); polcoeff( 8*A^7 -7*A^3*subst(A, x, -x)^4, n))} /* Michael Somos Oct 24 2006 */
(PARI) {a(n)= if(n<1, n==0, 2*qfrep([2, -1, 0, 0, 0, 0, 0; -1, 2, -1, 0, 0, 0, 0; 0, -1, 2, -1, 0, 0, 0; 0, 0, -1, 2, -1, 0, -1; 0, 0, 0, -1, 2, -1, 0; 0, 0, 0, 0, -1, 2, 0; 0, 0, 0, -1, 0, 0, 2], n, 1)[n])} /* Michael Somos Jun 11 2007 */
|