%I #7 May 20 2014 03:25:25
%S 29,41,53,59,61,83,89,101,107,109,113,131,137,139,149,157,173,179,181,
%T 197,211,227,229,233,241,251,257,269,277,281,283,293,307,313,317,331,
%U 337,347,349,353,373,379,389,397,401,409,419,421,433,443,449,457,461
%N Prime numbers that are the sum of three distinct positive squares.
%H Zak Seidov, <a href="/A125516/b125516.txt">Table of n, a(n) for n = 1..1000</a>
%e 29 = 2^2 + 3^2 + 4^2 = 4 + 9 + 16.
%e 89 = 2^2 + 6^2 + 7^2 = 4 + 36 + 49; also 89 = 3^2 + 4^2 + 8^2 = 9 + 16 + 64.
%e 353 = 2^2 + 5^2 + 18^2 = 4 + 25 + 324; also 353 = 4^2 + 9^2 + 16^2 = 16 + 81 + 256.
%o (PARI) {m=22;p=m^2;v=vector(m,x,x^2);w=[];for(i=1,m-2,for(j=i+1,m-1,for(k=j+1,m,if((n=v[i]+v[j]+v[k])<p&&isprime(n),w=concat(w,n)))));w=listsort(List(w),1);for(j=1,#w-1,print1(w[j],","))} /* Klaus Brockhaus, Feb 11 2007 */
%Y Cf. A088348, A126657, A126658.
%K nonn
%O 1,1
%A _Tomas Xordan_, Jan 21 2007
%E Edited, corrected and extended by _Klaus Brockhaus_, Feb 11 2007
|