python 数学公式_python – 使用指数有效地计算数学公式

让我们首先修复原始代码并定义一个函数来计算小n的F(n).我们还将打印出F的前几个值.下面的所有代码都是针对Python 3的;如果您使用的是Python 2,则需要进行一些小的更改,例如将str.maketrans替换为string.maketrans并使用xrange替换.

swap_ac = str.maketrans({ord('a'): 'c', ord('c'): 'a'})

def F(n):

s = ''

for _ in range(n):

s = s + 'a' + s[::-1].translate(swap_ac)

return s

for n in range(7):

print("F({}) = {!r}".format(n, F(n)))

这给出了以下输出：

F(0) = ''

F(1) = 'a'

F(2) = 'aac'

F(3) = 'aacaacc'

F(4) = 'aacaaccaaaccacc'

F(5) = 'aacaaccaaaccaccaaacaacccaaccacc'

F(6) = 'aacaaccaaaccaccaaacaacccaaccaccaaacaaccaaaccacccaacaacccaaccacc'

此时有几点意见：

> F(n)是长度为2 ** n-1的字符串.这意味着F(n)增长很快.计算F(50)已经需要一些严肃的硬件：即使我们每位存储一个字符,我们也需要超过100太字节来存储完整的字符串. F(200)具有比太阳系中估计的原子更多的特征.所以直接计算F(10 ** 2017)的想法是可笑的：我们需要一种不同的方法.

>通过构造,每个F(n)是F(n 1)的前缀.所以我们真正拥有的是一个定义良好的无限字符串,其中每个F(n)只给出了该无限字符串的前2 ** n-1个字符,我们正在寻找计算其第k个字符.并且出于任何实际目的,F(10 ** 2017)也可能是无限的字符串：例如,当我们进行计算时,我们不需要检查k< 2 **(10 ** 2017)-1,因为超过这个的k甚至不能用这个宇宙中的正常二进制表示法表示.

幸运的是,字符串的结构非常简单,直接计算第k个字符非常简单.当我们在偶数和奇数位置查看角色时,主要线索就出现了：

>>> F(6)[::2]

'acacacacacacacacacacacacacacacac'

>>> F(6)[1::2]

'aacaaccaaaccaccaaacaacccaaccacc'

偶数位置的字符只是在a和c之间交替(并且根据结构可以直接证明这是真的).因此,如果我们的k是偶数,我们可以简单地看看k / 2是奇数还是偶数来确定我们是否会得到a或c.

奇怪的位置怎么样？好吧F(6)[1 :: 2]应该看起来有些熟悉：它只是F(5)：

>>> F(6)[1::2] == F(5)

True

同样,直接证明(例如,通过归纳)这不仅仅是巧合,而且对于所有非负n来说F(n 1)[1 :: 2] == F(n).

我们现在有一种有效的方法来计算无限字符串中的第k个字符：如果k是偶数,我们只看k / 2的奇偶校验.如果k是奇数,那么我们知道位置k处的字符等于位置(k-1)/ 2处的字符.所以这是计算该角色的第一个解决方案：

def char_at_pos(k):

"""

Return the character at position k of the string F(n), for any

n satisfying 2**n-1 > k.

"""

while k % 2 == 1:

k //= 2

return 'ac'[k//2%2]

并检查这是否正确：

>>> ''.join(char_at_pos(i) for i in range(2**6-1))

'aacaaccaaaccaccaaacaacccaaccaccaaacaaccaaaccacccaacaacccaaccacc'

>>> ''.join(char_at_pos(i) for i in range(2**6-1)) == F(6)

True

但我们可以做得更好.我们有效地盯着k的二进制表示,删除所有尾随的’1’和下一个’0′,然后简单地查看下一位以确定我们是否有’a’或’c’.识别尾随1可以通过位操作技巧来完成.这为我们提供了以下半混淆的无循环解决方案,我将其留给您放松：

def char_at_pos2(k):

"""

Return the character at position k of the string F(n), for any

n satisfying 2**n-1 > k.

"""

return 'ac'[k//(1+(k+1^k))%2]

再次,让我们检查一下：

>>> F(20) == ''.join(char_at_pos2(i) for i in range(2**20-1))

True

最后评论：这是一个非常着名且经过充分研究的序列：它被称为龙曲线序列,或者说是regular paper-folding sequence,并且在整数序列的在线百科全书中是sequence A014577.某些Google搜索可能会为您提供许多其他方法来计算其元素.另见this codegolf question.