A230279 Number of nonnegative integer solutions to the equation x^2 - 16*y^2 = n.

1, 0, 0, 1, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 2, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 2, 0, 0, 1, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 1, 2, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 3, 0, 0, 2, 0, 0

Offset

1,9

Comments

For (x, y) to be a solution to the more general equation x^2 - d^2*y^2 = n, it can be shown that n-f^2 must be divisible by 2*f*d, where f is a divisor of n not exceeding sqrt(n). Then y = (n-f^2)/(2*f*d) and x = d*y+f.

Links

Bruno Berselli, Table of n, a(n) for n = 1..1000

Example

a(33)=2 because x^2 - 16*y^2 = 33 has two solutions: (x,y) = (17,4) and (7,1).

Prog

(PARI) a(n) = sumdiv(n, f, f^2<=n && (n-f^2)%(8*f)==0);
(Magma) d:=4; solutions:=func<i | [f: f in Divisors(i) | f le Isqrt(i) and IsZero((i-f^2) mod (2*f*d))]>; [#solutions(n): n in [1..100]]; // Bruno Berselli, Oct 15 2013

Crossrefs

Cf. A034178, A230263, A230264.

Sequence in context: A078359 A107329 A263717 * A085859 A218218 A332006

Adjacent sequences: A230276 A230277 A230278 * A230280 A230281 A230282

Keyword

nonn

Author

Colin Barker, Oct 15 2013

Status

approved