|
%I #6 Nov 21 2013 12:48:22
%S 17,34,68,85,136,153,170,221,272,289,306,340,425,442,493,544,578,612,
%T 629,680,697,765,833,850,884,901,986,1037,1088,1105,1156,1224,1241,
%U 1258,1360,1377,1394,1445,1513,1530,1649,1666,1700,1717,1768
%N Sums of two squares and divisible by 17.
%e Sums not divisible by 17 are indicated by asterisks:
%e ....|...1....4....9...16....25....36....49....64....81...100...121
%e ------------------------------------------------------------------
%e 1...|...*....*....*...17.....*.....*.....*.....*.....*.....*.....*
%e 4...|...*....*....*....*.....*.....*.....*....68....85.....*.....*
%e 9...|...*....*....*....*....34.....*.....*.....*.....*.....*.....*
%e 16..|..17....*....*....*.....*.....*.....*.....*.....*.....*.....*
%e 25..|...*....*...34....*.....*.....*.....*.....*.....*.....*.....*
%e 36..|...*....*....*....*.....*.....*....85.....*.....*...136.....*
%e 49..|...*....*....*....*.....*....85.....*.....*.....*.....*...170
%e 64..|...*...68....*....*.....*.....*.....*.....*.....*.....*.....*
%e 81..|...*...85....*....*.....*.....*.....*.....*.....*.....*.....*
%e 100.|...*....*....*....*.....*...136.....*.....*.....*.....*...221
%e 121.|...*....*....*....*.....*.....*...170.....*.....*...221.....*
%t ssq[upto_]:=Module[{max=Floor[Sqrt[upto]],seq},seq=Select[Union[ Total/@Tuples[ Range[max]^2,{2}]],Divisible[#,17]&];Select[seq,#<=upto&]]; ssq[1800] (* _Harvey P. Dale_, Mar 24 2011 *)
%K nonn
%O 1,1
%A Jun Mizuki (suzuki32(AT)sanken.osaka-u.ac.jp), Oct 26 2004
|
