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