刚才在这里看到了如题所说的图像，立即想到用 Mathematica 验证一下。我选出了几个个人比较感兴趣的 k ，再用一句话便可输出所有对应 k 的图像：
kArray = {2, 3, 4, 6, 8, 10, 12, 14, 16, 18, 20, 36, 50};
For[i = 1, i <= Length[kArray], i++, Export["F:\" <> ToString[kArray[[i]]] <> ".png",
ArrayPlot[Table[Boole[Length[Divisors[x*y]] == kArray[[i]]], {x, 1, 400}, {y, 1, 400}],
PixelConstrained -> {1, 1}, Frame -> False]]];
当 k=2 时，由于只有素数才有两个约数，因此所有点都是形如 (p, 1) 或者 (1, p) 的点，其中 p 为某个素数：
图1
当 k=3 时，图像上的点大致成一条斜线，斜线上的点都是形如 (p, p) 的点，其中 p 为某个素数：
图2
当k是更大一些的奇数时，情况类似，因此不再重复了。
当 k=4 时，图像上出现了很多竖直线和水平线，这些直线的位置与素数的平方相对应：
图3
当 k=6 时，逐渐出现了一些乍看之下不好理解的元素：
图4
当 k=8 时，出现了之前不曾见过的情况：图中有一些连续的直线段：
图5
当 k=10 时：
图6
当 k=12 时：
图7
当 k=14 时：
图8
当 k=16 时，更加怪异了：
图9
当 k=18 时：
图10
当 k=20 时：
图11
当 k=36 时：
图12
当 k=50 时：
图13
研究这些图像各种有趣的模式产生的原因，无疑是一个有趣的课题。
声明：文章转自Matrix67博客，版权归原作者所有，转载仅供学习使用，不用于任何商业用途，如有侵权请联系删除，谢谢。
