r=range(63) for i in r:m=[[i>0==(x*y-961)*0**m[y][x]*~-sum(sum(k[x+x%~x:x+2])for k in m[y+y%~y:y+2])for x in r]for y in r] for i in r:print"".join(" X"[_]for _ in m[i])

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