混淆的曼德尔布罗特函数 - 谁能解密它?

3 投票
1 回答
501 浏览
提问于 2025-04-18 15:57
_                                      =   (
                                        255,
                                      lambda
                               V       ,B,c
                             :c   and Y(V*V+B,B,  c
                               -1)if(abs(V)<6)else
               (              2+c-4*abs(V)**-0.4)/i
                 )  ;v,      x=7200,4800;C=range(v*x
                  );import  struct;P=struct.pack;M,\
            j  ='<QIIHHHH',open('M.bmp','wb').write
for X in j('BM'+P(M,v*x*3+26,26,12,v,x,1,24))or C:
            i  ,Y=_;j(P('BBB',*(lambda T:(T*80+T**9
                  *i-950*T  **99,T*70-880*T**18+701*
                 T  **9     ,T*i**(1-T**45*2)))(sum(
               [              Y(0,(A%3/3.+X%v+(X/v+
                               A/3/3.-x/2)/1j)*2.5
                             /x   -2.7,i)**2 for  \
                               A       in C
                                      [:9]])
                                        /9)
                                       )   )

我刚开始学习Python,对它能生成分形图像非常感兴趣。我自己写了一些简单的分形图,但我最近发现了一个生成曼德尔布罗特分形的脚本……它可以生成美丽的全彩图像,分辨率也可以根据你的需求调整……不过,这段代码被处理得像是曼德尔布罗特的ASCII艺术,虽然看起来很酷,但如果想轻松阅读就有点傻了。里面有一些我还没学过的Python函数,所以如果有人能把这段脚本整理成正常的Python格式,那就太好了。

正如我所说,这个艺术作品很酷,但太难读了!如果有人能帮忙做到这一点,我会非常感激。

1 个回答

3

注意:这不是一个最终的答案,而是一步一步地解释代码的过程。如果你能提供额外的步骤让事情更清楚,那就太好了,可以把它加到这个回答里。

好的,首先我们给代码加上适当的换行和缩进:

_ = (255,
    lambda V, B, c: c and Y(V*V + B, B, c-1) if(abs(V) < 6) else (2 + c - 4 * abs(V)**-0.4)/i)
v, x = 7200, 4800
C = range(v*x)
import struct
P=struct.pack
M, j = '<QIIHHHH', open('M.bmp','wb').write
for X in j('BM' + P(M, v*x*3+26, 26, 12, v, x, 1, 24)) or C:
    i, Y = _
    j(
        P('BBB',
                *(lambda T: (T * 80 + T**9 * i - 950 * T**99, T*70 - 880* T**18 + 701* T**9, T* i**(1-T**45*2)))(
                    sum([Y(0, (A%3/3. + X%v + (X/v + A/3/3. - x/2)/1j) * 2.5/x - 2.7, i)**2 for A in C[:9]]) / 9
                )
        )
    )

这样看起来好一些,但还是有点让人困惑。比如,把 P 定义为 struct.pack 的别名,还有那些一行写两个声明的地方。如果我们把这些去掉,把 lambda 定义移到循环外面,就变成这样:

import struct
i = 255
Y = lambda V, B, c: c and Y(V*V + B, B, c-1) if(abs(V) < 6) else (2 + c - 4 * abs(V)**-0.4)/i
v = 7200
x = 4800
C = range(v*x)
M = '<QIIHHHH'
color = lambda T: (T * 80 + T**9 * i - 950 * T**99, T*70 - 880* T**18 + 701 * T**9, T * i**(1-T**45*2))

f = open('M.bmp','wb')
for X in f.write('BM' + struct.pack(M, v*x*3+26, 26, 12, v, x, 1, 24)) or C:
    f.write(
        struct.pack('BBB',
            *color(sum([Y(0, (A%3/3. + X%v + (X/v + A/3/3. - x/2)/1j) * 2.5/x - 2.7, i)**2 for A in C[:9]]) / 9))
    )

现在事情开始变得清晰一些了。这个循环是在为每个像素写颜色值,而 lambda 返回的是一个包含三个值的元组,分别代表每个像素的蓝色、绿色和红色值。

import struct
i = 255
Y = lambda V, B, c: c and Y(V*V + B, B, c-1) if abs(V) < 6 else (
                    (2 + c - 4 * abs(V)**-0.4)/i)
v = 7200
x = 4800
C = range(v*x)
M = '<QIIHHHH'
color = lambda T: (T * 80 + T**9 * i - 950 * T**99,
                   T*70 - 880* T**18 + 701 * T**9,
                   T * i**(1-T**45*2))
f = open('M.bmp', 'wb')
for X in f.write('BM' + struct.pack(M, v*x*3+26, 26, 12, v, x, 1, 24)) or C:
    f.write(struct.pack('BBB',
        *color(sum(
            [Y(0, (A%3/3. + X%v + (X/v + A/3/3. - x/2)/1j) * 2.5/x - 2.7, i)**2
                for A in C[:9]]) / 9))
    )

根据你的评论,星号(*)可以把一个列表拆分成参数列表

经过很多小时的努力,这里是生成的接近100 MB的 曼德尔布罗特集合 的图像:

曼德尔布罗特集合

撰写回答