%I #26 Jul 13 2013 12:04:29
%S 8,13,18,29,34,50,53,58,74,98,125,130,146,170,173,178,194,218,242,290,
%T 293,298,314,365,370,386,458,482,530,533,538,554,698,722,818,845,850,
%U 866,962,965,970,986,1058,1082,1202,1250,1322,1370,1373,1378,1394,1418
%N Numbers of the form p^2 + q^2, with p and q prime, in exactly one way.
%C A045698(a(n)) = 1. - _Reinhard Zumkeller_, Jul 29 2012
%H Reinhard Zumkeller, <a href="/A214723/b214723.txt">Table of n, a(n) for n = 1..1000</a>
%H <a href="/index/Su#ssq">Index entries for sequences related to sums of squares</a>
%e a(1) = 8 = 2^2 + 2^2, and no other p^2 + q^2 sums to 8.
%e Both 7^2 + 17^2 and 13^2 + 13^2 sum to 338, so 338 is not in this sequence.
%t nn = 2000; ps = Prime[Range[PrimePi[Sqrt[nn]]]]; t = Flatten[Table[ps[[i]]^2 + ps[[j]]^2, {i, Length[ps]}, {j, i, Length[ps]}]]; t = Select[t, # <= nn &]; Sort[Transpose[Select[Tally[t], #[[2]] == 1 &]][[1]]] (* _T. D. Noe_, Jul 26 2012 *)
%o (Haskell)
%o import Data.List (elemIndices)
%o a214723 n = a214723_list !! (n-1)
%o a214723_list = elemIndices 1 a045698_list
%o -- _Reinhard Zumkeller_, Jul 29 2012
%Y Subsequence of A045636 (numbers of the form p^2 + q^2).
%K nonn
%O 1,1
%A _J. Stauduhar_, Jul 26 2012
|