nsq.py

from PIL import Image
from fractions import gcd
from math import sqrt

w = 2048

img = Image.new("RGB", (w, w), color=(255,255,255))

for i in range(1,w):
    for j in range(1,w):
        m = i*j
        n = int(sqrt(m))
        p = n
        
        while m%p != 0:
            p += 1
        
        q = m/p

        img.putpixel((p-1,w-q-1), (0,0,0))
        img.putpixel((q-1, w-p-1), (0,0,0))

img.save("nsq.jpg")