前幾天在上競賽的培訓課時,講到尤拉函數(Euler’s φ
function),一位協助指導的老師提供了一個用機率的「獨立事件」的觀點的說明方法,記錄如下:
定義(尤拉φ函數):
設n為自然數,φ(n)定義為小於n且與n互質的正整數的個數。
首先看一道題目:α、β、γ為正整數,試求小於n=2α×3β×5γ,且不是2,不是3,也不是5的倍數有幾個(即求φ(n)的值)?
設一個試驗為從1~30的整數中,任選一數,則樣本空間S={1, 2, ..., 30},A為2的倍數的事件,B為3的倍數的事件,C為5的倍數的事件。
設一個試驗為從1~30的整數中,任選一數,則樣本空間S={1, 2, ..., 30},A為2的倍數的事件,B為3的倍數的事件,C為5的倍數的事件。
顯然A∩B為2×3=6的倍數的事件,B∩C為3×5=15的倍數的事件,C∩A為5×2=10的倍數的事件,A∩B∩C為2×3×5=30的倍數的事件。
故P(A)=1/2,P(B)=1/3,P(C)=1/5,P(A∩B)=P(A)×P(B)=1/6,P(B∩C)=P(B)×P(C)=1/15,P(C∩A)=P(C)×P(A)=1/10,P(A∩B∩C)=P(A)×P(B)×P(C)=1/30。
所以A、B、C為獨立事件。所以若從小於n的數中選出不是2、3、5的倍數的數的事件為A’∩B’∩C’,機率為P(A’∩B’∩C’)=(1- P(A))(1-P(B))(1-P(C))=(1-1/2)(1-1/3)(1-1/5),故φ(n)=n×(1-1/2)(1-1/3)(1-1/5)
實際上,上面的題目講的是:
設p1、p2、…、pr是n所有相異的質因數。若從1~n的整數中任取一數,A1、A2、…、Ar分別是取到p1、p2、…、pr的倍數的事件,則A1’ ∩A2’ ∩…∩Ar’即為取到的整數為與n互質的數的事件。
因為p1、p2、…、pr兩兩互質且都是n的因數,所以A1、A2、…、Ar是獨立事件,當然A1’、A2’、…、Ar’
也是獨立事件,故
P(A1’ ∩A2’ ∩…∩Ar’)= P(A1’)P(A2’)…P(Ar’)=(1-1/p1) (1-1/p2)…(1-1/pr),
φ(n)=nP(A1’ ∩A2’ ∩…∩Ar’)=n(1-1/p1) (1-1/p2)…(1-1/pr)。
沒有留言:
張貼留言