OFFSET
1,1
COMMENTS
Numbers of the following forms: p[1]*p[2]^2*p[3]^2, p[1]^2*p[2]^5, p[1]*p[2]^8, p[1]^17, 2^2*p[1]*p[2]^2*p[3]^2, 2^2*p[1]^2*p[2]^5, 2^3*p[1]^2*p[2]^2, 2^3*p[1]^8, 2^4*p[1]*p[2]^2, 2^4*p[1]^5, 2^7*p[1]^2, 2^10*p[1], 2^19, where p[i] are distinct odd primes. - Robert Israel, Jun 19 2018
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000 (first 260 terms from Colin Barker)
EXAMPLE
720 is in the sequence because there are 9 solutions to x^2 - y^2 = 720, namely (x,y) = (27,3), (28,8), (29,11), (36,24), (41,31), (49,41), (63,57), (92,88), (181,179).
MAPLE
filter:= proc(n) local k;
k:= padic:-ordp(n, 2);
(k = 0 and numtheory:-tau(n)=18) or (k-1)*numtheory:-tau(n/2^k)=18
end proc:
select(filter, [$1..10^4]); # Robert Israel, Jun 19 2018
MATHEMATICA
nn = 6800;
t = Table[0, {nn}];
Do[n = x^2 - y^2; If[n <= nn, t[[n]]++], {x, nn}, {y, 0, x - 1}];
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Colin Barker, Apr 22 2015
STATUS
approved