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