login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A124865
Numbers of the form p^2-q^2 with p > q prime.
6
5, 16, 21, 24, 40, 45, 48, 72, 96, 112, 117, 120, 144, 160, 165, 168, 192, 240, 264, 280, 285, 288, 312, 336, 352, 357, 360, 408, 432, 480, 504, 520, 525, 528, 552, 600, 648, 672, 720, 768, 792, 816, 832, 837, 840, 888, 912, 936, 952, 957, 960, 1008, 1032, 1080
OFFSET
1,1
COMMENTS
The only prime term is a(1) = 5.
All odd terms are of the form p^2-4.
All even terms are divisible by 8.
Numbers of the form (p^2-q^2)/8 (p, q odd primes, p>q) are listed in A124866.
Elliott & Richner call these "ans numbers". - Charles R Greathouse IV, Feb 17 2014
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
N. E. Elliott and D. Richner, An investigation of the set of ans numbers, Missouri J. of Math. Sci. 15 (2003), pp. 189-199.
Florian Luca, On the densities of some subsets of integers, Missouri Journal of Mathematical Sciences 19:3 (2007), pp. 167-170.
FORMULA
a(n) >> n log n, see Luca. - Charles R Greathouse IV, Feb 17 2014
MATHEMATICA
With[{nn=60}, Take[Union[#[[2]]^2-#[[1]]^2&/@Subsets[Prime[Range[nn]], {2}]], nn]] (* Harvey P. Dale, Aug 21 2015 *)
PROG
(PARI) is(n)=if(n%24, isprimepower(n+4)==2 || isprimepower(n+9)==2, fordiv(n/4, d, if(isprime(n/d/4+d) && isprime(n/d/4-d), return(1))); 0) \\ Charles R Greathouse IV, Feb 17 2014
CROSSREFS
Apart from a(1), a subsequence of A177713.
Cf. A045636 (numbers of the form p^2+q^2, p, q primes).
Cf. A124866 (numbers of the form (p^2-q^2)/8, p, q odd primes, p>q).
Sequence in context: A057281 A204920 A168469 * A090781 A191264 A077469
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Nov 10 2006
STATUS
approved