|
PROG
|
(PARI) {a(n) = sign(n) * subst(elldivpol(ellinit([0, 17]), abs(n)), x, -2) / (if(n%2, 1, 6) * 2^(n^2\3) * 3^(n^2\4))};
(PARI) {a(n) = my(v, s); s=sign(n); n=abs(n); if(n==0, return(0)); v = concat([1, 1, -5, -23, 181, 1740, -8749], vector(max(0, n-7))); for(k=8, n, v[k] = (-23*v[k-1]*v[k-6] +348*v[k-3]*v[k-4])/v[k-7]); s*v[n]};
(PARI) {a(n) = my(v, s); s=sign(n); n=abs(n); if(n==0, return(0)); v = concat([1, 1, -5, -23], vector(max(0, n-4))); for(k=5, n, v[k] = (if(k%2, 3, 1)*v[k-1]*v[k-3] +if(k%3==2, 10, 5)*v[k-2]^2)/v[k-4]); s*v[n]};
(PARI) {a(n) = my(v, s); s=sign(n); n=abs(n); if(n==0, return(0)); v = concat([1, 1, -5, -23, 181], vector(max(0, n-5))); for(k=6, n, v[k] = (-if(k%3==1, 10, 5)*v[k-1]*v[k-4] +23*v[k-2]*v[k-3])/v[k-5]); s*v[n]};
|