r=range(-23,24) for y in r: s="" for x in r:p=x*x+y*y or 1;s+=" [ ]"[p**.5%16in[3,7,11]or 10103**~-p%p==1::2] print s

Note that non-ascii characters in the above source code will be escaped (such as \x9f).

download

return top