关于谢尔宾斯基地毯的讲解

        和谢尔宾斯基三角形一样，谢尔宾斯基地毯也是数学家谢尔宾斯基提出的一个分形图形，谢尔宾斯基地毯和谢尔宾斯基三角形基本类似，不同之处在于谢尔宾斯基地毯采用的是正方形进行分形构造，而谢尔宾斯基三角形采用的等边三角形进行分形构造。谢尔宾斯基地毯和它本身的一部分完全相似，减掉一块会破坏自相似性。

​

 ​

理解了这些后，运行以下代码： 

# 舞台自动重置,kpen和p可直接使用.
p.width(3)
p.speed(10)
import turtle
n = int(input())
# 画好最外围的正方形，背景设为蓝色
base = int(200 / (3 ** n))  # 设置每个小正方形长度base
length = base * (3 ** n)  # 最外围的正方形长度
p.pensize(0.1)
turtle.tracer(False)  # 画图过程太长我们跳过不看
p.pencolor('blue')
p.fillcolor('blue')
p.begin_fill()  # 填充蓝色
for i in range(1, 5):
    p.forward(length)
    p.left(90)
p.hide()
p.end_fill()
def carpet(N):  # 对于点a（x，y），我们令a点位小正方形的左下角顶点
    def check(n, x, y):  # 判断对于坐标（x，y），判断该处的小正方形是否需要被挖掉，需要则返回False
        if n <= base:  # 结束条件
            return True
        n2 = n // 3
        if n2 <= x < n2 * 2 and n2 <= y < n2 * 2:  # 小正方形处于大正方形中央应该满足的条件
            return False
        return check(n2, x % n2, y % n2)  # 不处于中央的点，递归至下一层的小正方形
        # 取余的操作可以理解为将零点移至小正方形的左下角顶点

    for y in range(0, N, base):
        for x in range(0, N, base):
            if check(N, x, y):
                continue
            else:            # 将（x，y）处小正方形填为空白
                p.penup()
                p.setx(x)
                p.sety(y)
                p.pendown()
                p.pensize(0.1)
                p.pencolor('white')
                p.fillcolor('white')
                p.begin_fill()
                for i in range(1, 5):
                    p.forward(base)
                    p.left(90)
                p.end_fill()
                p.penup()


carpet(length)

p.draw()

运行如下： 

​