亚洲免费在线-亚洲免费在线播放-亚洲免费在线观看-亚洲免费在线观看视频-亚洲免费在线看-亚洲免费在线视频

用Python實現(xiàn)斐波那契(Fibonacci)函數(shù)

系統(tǒng) 1550 0

Fibonacci斐波那契數(shù)列,很簡單,就是一個遞歸嘛,學(xué)任何編程語言可能都會做一下這個。

最近在玩Python,在粗略的看了一下Learning Python和Core Python之后,偶然發(fā)現(xiàn)網(wǎng)上有個帖子Python程序員的進化寫的很有意思。于是打算仿照一篇,那篇帖子用了十余種方法完成一個階乘函數(shù),我在這里會用九種不同的風(fēng)格寫出一個Fibonacci函數(shù)。

要求很簡單,輸入n,輸出第n個Fibonacci數(shù),n為正整數(shù)

下面是這九種不同的風(fēng)格:

1)第一次寫程序的Python程序員:

            
def fib(n):
  return nth fibonacci number

          

說明:
第一次寫程序的人往往遵循人類語言的語法而不是編程語言的語法,就拿我一個編程很猛的哥們來說,他寫的第一個判斷閏年的程序,里面直接是這么寫的:如果year是閏年,輸出year是閏年,否則year不是閏年。

2)剛學(xué)Python不久的的C程序員:

            
def fib(n):#{
 if n<=2 :
  return 1;
 else:
  return fib(n-1)+fib(n-2);
#}
          

說明:
在剛接觸Python時,用縮進而非大括號的方式來劃分程序塊這種方式我是很不適應(yīng)的,而且每個語句后面沒有結(jié)束符,所以每次寫完一個Python函數(shù)之后干的第一件事一般就是一邊注釋大括號,一邊添加漏掉的冒號。

3)懶散的Python程序員:

            
def fib(n):
  return 1 and n<=2 or fib(n-1)+fib(n-2)
          

說明:
看了Learning Python之后,才知道Python沒有三元操作符?,不過鑒于Python里bool值比較特殊(有點像C,非零即真,非空即真),再加上Python的邏輯語句也是支持短路求值(Short-Circuit Evaluation)的,這就可以寫出一個仿?語句出來。

4)更懶的Python程序員:

            
 fib=lambda n:1 if n<=2 else fib(n-1)+fib(n-2)

          

說明:
lambda關(guān)鍵字我曾在C#和Scheme里面用過,Python里面的lambda比C#里簡便,并很像Scheme里的用法,所以很快就適應(yīng)了。在用Python Shell聲明一些小函數(shù)時經(jīng)常用這種寫法。

5)剛學(xué)完數(shù)據(jù)結(jié)構(gòu)的Python程序員:

            
def fib(n):
 x,y=0,1
 while(n):
  x,y,n=y,x+y,n-1
 return x
          

說明:
前面的Fibonacci函數(shù)都是樹形遞歸的實現(xiàn),哪怕是學(xué)一點算法就應(yīng)該知道這種遞歸的低效了。在這里從樹形遞歸改為對應(yīng)的迭代可以把效率提升不少。
Python的元組賦值特性是我很喜歡的一個東東,這玩意可以把代碼簡化不少。舉個例子,以前的tmp=a;a=b;b=tmp;可以直接用一句a,b=b,a實現(xiàn),既簡潔又明了。

6)正在修SICP課程的Python程序員:

            
def fib(n):
  def fib_iter(n,x,y):
   if n==0 : return x
   else : return fib_iter(n-1,y,x+y)

  return fib_iter(n,0,1)
          

說明:
在這里我使用了Scheme語言中很常見的尾遞歸(Tail-recursion)寫法。Scheme里面沒有迭代,但可以用不變量和尾遞歸來模擬迭代,從而實現(xiàn)相同的效果。不過我還不清楚Python有沒有對尾遞歸做相應(yīng)的優(yōu)化,回頭查一查。
PS:看過SICP的同學(xué),一眼就能看出,這個程序其實就是SICP第一章里的一個例子。

7)好耍小聰明的Python程序員:

            
fib=lambda n,x=0,y=1:x if not n else f(n-1,y,x+y)
          

說明:
基本的邏輯和上面的例子一樣,都是尾遞歸寫法。主要的區(qū)別就是利用了Python提供的默認參數(shù)和三元操作符,從而把代碼簡化至一行。至于默認參數(shù),學(xué)過C++的同學(xué)都知道這玩意,至于C#4.0也引入了這東東。

8)剛修完線性代數(shù)的Python程序員:

            
def fib(n):
 def m1(a,b):
  m=[[],[]]
  m[0].append(a[0][0]*b[0][0]+a[0][1]*b[1][0])
  m[0].append(a[0][0]*b[0][1]+a[0][1]*b[1][1])
  m[1].append(a[1][0]*b[0][0]+a[1][1]*b[1][0])
  m[1].append(a[1][0]*b[1][0]+a[1][1]*b[1][1])
  return m
 def m2(a,b):
  m=[]
  m.append(a[0][0]*b[0][0]+a[0][1]*b[1][0])
  m.append(a[1][0]*b[0][0]+a[1][1]*b[1][0])
  return m
 return m2(reduce(m1,[[[0,1],[1,1]] for i in range(n)]),[[0],[1]])[0]
          

說明:
這段代碼就不像之前的代碼那樣清晰了,所以先介紹下原理(需要一點線性代數(shù)知識):
首先看一下之前的迭代版本的Fibonacci函數(shù),很容易可以發(fā)現(xiàn)存在一個變換:y->x, x+y->y。換一個角度,就是[x,y]->[y,x+y]。
在這里,我聲明一個二元向量[x,y]T,它通過一個變換得到[y,x+y]T,可以很容易得到變換矩陣是[[1,0],[1,1]],也就是說:[[1,0],[1,1]]*[x,y]T=[y,x+y]T
令二元矩陣A=[[1,0],[1,1]],二元向量x=[0,1]T,容易知道Ax的結(jié)果就是下一個Fibonacci數(shù)值,即:
Ax=[fib(1),fib(2)]T
亦有:
Ax=[fib(2),fib(3)]T
………………
以此類推,可以得到:

            
A?x=[fib(n),fib(n-1)]T

          

也就是說可以通過對二元向量[0,1]T進行n次A變換,從而得到[fib(n),fib(n+1)]T,從而得到fib(n)。

在這里我定義了一個二元矩陣的相乘函數(shù)m1,以及一個在二元向量上的變換m2,然后利用reduce操作完成一個連乘操作得到A?x,最后得到fib(n)。

9)準備參加ACM比賽的Python程序員:

            
 
def fib(n):
 lhm=[[0,1],[1,1]]
 rhm=[[0],[1]]
 em=[[1,0],[0,1]]
 #multiply two matrixes
 def matrix_mul(lhm,rhm):
  #initialize an empty matrix filled with zero
  result=[[0 for i in range(len(rhm[0]))] for j in range(len(rhm))]
  #multiply loop
  for i in range(len(lhm)):
   for j in range(len(rhm[0])):
    for k in range(len(rhm)):
     result[i][j]+=lhm[i][k]*rhm[k][j]
  return result
 
 def matrix_square(mat):
  return matrix_mul(mat,mat)
 #quick transform
 def fib_iter(mat,n):
  if not n:
   return em
  elif(n%2):
   return matrix_mul(mat,fib_iter(mat,n-1))
  else:
   return matrix_square(fib_iter(mat,n/2))
 return matrix_mul(fib_iter(lhm,n),rhm)[0][0]


          

說明:

看過上一個fib函數(shù)就比較容易理解這一個版本了,這個版本同樣采用了二元變換的方式求fib(n)。不過區(qū)別在于這個版本的復(fù)雜度是lgn,而上一個版本則是線性的。

這個版本的不同之處在于,它定義了一個矩陣的快速求冪操作fib_iter,原理很簡單,可以類比自然數(shù)的快速求冪方法,所以這里就不多說了。

PS:雖然說是ACM版本,不過說實話我從來沒參加過那玩意,畢竟自己算法太水了,那玩意又太高端……只能在這里YY一下鳥~

python中,最基本的那種遞歸(如下fib1)效率太低了,只要n數(shù)字大了運算時間就會很長;而通過將計算的指保存到一個dict中,后面計算時直接拿來使用,這種方式成為備忘(memo),如下面的fib2函數(shù)所示,則會發(fā)現(xiàn)效率大大提高。

在n=10以內(nèi)時,fib1和fab2運行時間都很短看不出差異,但當(dāng)n=40時,就太明顯了,fib1運行花了35秒,fab2運行只花費了0.00001秒。
n=40時,輸出如下:

            
jay@jay-linux:~/workspace/python.git/py2014$ python fibonacci.py 
2014-10-16 16:28:35.176396
fib1(40)=102334155
2014-10-16 16:29:10.479953
fib2(40)=102334155
2014-10-16 16:29:10.480035

          

這兩個計算Fibonacci數(shù)列的函數(shù),如下:https://github.com/smilejay/python/blob/master/py2014/fibonacci.py

            
import datetime

def fib1(n):
  if n == 0:
    return 0
  elif n == 1:
    return 1
  else:
    return fib1(n - 1) + fib1(n - 2)
 
known = {0: 0, 1: 1}
 
def fib2(n):
  if n in known:
    return known[n]
 
  res = fib2(n - 1) + fib2(n - 2)
  known[n] = res
  return res

if __name__ == '__main__':
  n = 40
  print(datetime.datetime.now())
  print('fib1(%d)=%d' % (n, fib1(n)))
  print(datetime.datetime.now())
  print('fib2(%d)=%d' % (n, fib2(n)))
  print(datetime.datetime.now())


          

后記:

由于剛學(xué)習(xí)Python沒多久,所以對其各種特性的掌握還不夠熟練。與其說是我在用Python寫程序,倒不如說我是在用C,C++,C#或是Scheme來寫程序。至于傳說中的Pythonic way,我現(xiàn)在還沒有什么體會,畢竟還沒用Python寫過什么真正的程序。
Learning Python和Core Python都是不錯的Python入門書籍,前者更適合沒有編程基礎(chǔ)的人閱讀。
Python是最好的初學(xué)編程入門語言,沒有之一。所以它可以取代Scheme成為MIT的計算機編程入門語言。


更多文章、技術(shù)交流、商務(wù)合作、聯(lián)系博主

微信掃碼或搜索:z360901061

微信掃一掃加我為好友

QQ號聯(lián)系: 360901061

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧,狠狠點擊下面給點支持吧,站長非常感激您!手機微信長按不能支付解決辦法:請將微信支付二維碼保存到相冊,切換到微信,然后點擊微信右上角掃一掃功能,選擇支付二維碼完成支付。

【本文對您有幫助就好】

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描上面二維碼支持博主2元、5元、10元、自定義金額等您想捐的金額吧,站長會非常 感謝您的哦!!!

發(fā)表我的評論
最新評論 總共0條評論
主站蜘蛛池模板: 欧美成人性毛片免费版 | 久久伊人免费视频 | 2020久久精品永久免费 | 精品资源在线 | 久久亚洲精品中文字幕二区 | 狠狠激情五月综合婷婷俺 | 这里只有精品视频在线 | 国产福利视频一区美女 | 亚洲午夜一区二区三区 | 久久青草免费97线频观 | 亚洲一区欧洲一区 | 久久嫩模| 亚洲欧美久久 | 中文字幕一区二区日产乱码 | 97色精品视频在线观看免费 | 国内精品视频九九九九 | 一级有奶水毛片免费看 | 狠狠色噜噜狠狠狠狠 | 色综合久久中文字幕 | 亚洲两性| 国产小网站 | 亚洲精品香蕉一区二区在线观看 | 欧美永久免费 | 青春草国产成人精品久久 | 草久久久 | 成人a大片高清在线观看 | 成年人色视频 | 色噜噜狠狠色综合久 | 国产精品日本一区二区在线看 | 国产成人在线播放视频 | 狠狠久久综合伊人不卡 | 国产福利在线视频 | 欧美日韩顶级毛片www免费看 | 夜夜夜夜夜夜夜工噜噜噜 | 婷婷久久综合 | 又爽又黄又无遮挡的视频在线观看 | 国内夫妇精品对白在线播放 | 99久久精品无码一区二区毛片 | 国产精品亚洲综合一区 | 久久99爱视频| 四虎影院网站 |