Saturday, December 15, 2012

容斥原理与应用之广义容斥

继续转载
http://chowkafat.net/Enumeration8.html

點算的奧秘:廣義容斥原理


《點算的奧秘:容斥原理基本公式》中,筆者介紹了「容斥原理」的基 本公式。這些公式可用來計算以下兩個集合的元素個數:|S1 ∪ S2 ∪ ... Sn|和|S1' ∩ S2' ∩ ... Sn'|。用日常語言來表達 ,上述第一個集合代表屬於S1 ... Sn中「至少1個」集合的元素,第二個集合則代表 屬於S1 ... Sn中「剛好0個」集合的元素。以上兩個集合的共同點是,它們代表某種 「極端」情況。在本節我們要把「容斥原理」推廣至更一般的集合,這些集合的元素屬於S1 ... Sn中的「至少p個」或「剛好s個」集合,其中1 ≤ p ≤ n − 1和0 ≤ s ≤ n。

首先考慮「剛好s個」集合的情況。設有S1 ... Sn等n個集合,現在要求屬於這n個集 合中的「剛好s個」集合的那些元素的個數。我們把這個數記為N(n, s)(註1),以下給出這個數的公式:
N(n, s) = Σs ≤ k ≤ n(−1)k − s C(k, s) Sn, k    (1)
其中Sn, k的定義如下。若k = 0,則把Sn, k定義為論域U中元素的個數,即 Sn, 0 = |U|。若k > 0,則Sn, k是由C(n, k)個項加起來的總和,每個項都是由 S1 ... Sn這n個集合中抽k個出來構成的交集的元素個數。上式就是「廣義容斥原 理」(Generalized Inclusion and Exclusion Principle)的其中一種形式。以下筆者證明上述公式能正 確給出N(n, s)的值,方法是證明論域中每個剛好屬於s個集合的元素在上述公式中都被加一次,而且只有這些 元素被加在上述公式中(註2)。

對於論域中任何一個元素而言,這個元素屬於S1 ... Sn中的i個集合(0 ≤ i ≤ n),這個元素將被計算在Sn, i中。此外,由於這個元素也必然同時屬於這i個集合之中的k個(k < i),因此對於任何k < i而言,這個元素必然也被計算在Sn, j中。但對於任何k > i而言,這個元 素不可能被計算在Sn, k中。以上事實對於了解下面的證明是至關重要的,請讀者務必弄明白。由 於i與s可以存在三種關係:i < s;i = s和i > s,以下逐一考慮這三種情況。

首先考慮i < s的情況。在這種情況下,有關元素屬於少於s個集合。由於上述公式中求和符號Σ的下限是 s,Sn,i根本不會被加進上式中。換句話說,有關元素不會被加在上述公式中。

接著考慮i = s的情況。在這種情況下,有關元素剛好屬於s個集合,即它們會被計算在Sn, s中, 並且對於任何k < s而言,這些元素也會被計算在Sn, k中。但由於上述公式中求和符號Σ的 下限是s,我們只須考慮Sn, s的情況而無須考慮比s小的情況。請注意,如果某元素剛好屬於s個集 合,那麼這個元素只會出現於Sn, s中的某一個項(註3),即Sn, s只會把這個元素計算 一次。此外,由於(−1)s − s = 1並且C(s, s) = 1,這些元素在上述公式中剛好被加 一次。

最後考慮i > s的情況。在這種情況下,有關元素屬於超過s個集合。這類元素將被計算在Sn, i中 ,並且對於任何k < i而言,這類元素也會被計算在Sn, k中。但由於上述公式中求和符號Σ 的下限是s,我們只須考慮Sn, s、Sn,s+1 ... Sn, i−1、 Sn, i的情況。對於任何k ≤ i而言,由於從i個集合中抽k個出來構成交集共有C(i, k)種抽法, 所以每個這類元素在Sn, k中將會被計算C(i, k)次。因此每個這類元素在上述公式中被加的次數為
Σs ≤ k ≤ i(−1)k − s C(k, s) × C(i, k)     (2)
現在讓我們證明上式等於0。首先我們把式(2)中的C(k, s) × C(i, k)改寫:
 C(k, s) × C(i, k) × 1
 k! i!  (i − s)!
=-----------------× -----------------×-----------------
 (k − s)! × s! (i − k)! × k! (i − s)!
 i! (i − s)!
=-----------------× -----------------------
 (i − s)! × s! (i − k)! × (k − s)!
=C(i, s) × C(i − s, k − s)
這樣便可以把式(2)改寫成
C(i, s) × Σs ≤ k ≤ i(−1)k − s C(i − s, k − s) (註4)
接著我們對上式中求和符號Σ的變項k進行變換,變換公式為r = k − s,即用r代替上式中的k − s。但不要忘記Σ的上、限也要作相應變換。由於r = k − s,所以原來k的下限s應變為s − s = 0,而原來k的上限i則變成i − s。經過變換後,式(2)便進一步變成
C(i, s) × Σ0 ≤ r ≤ i−s(−1)r C(i − s, r) (註5)
現在把上式與「二項式定理」:(a + b)n = Σ 0 ≤ r ≤ n C(n, r) ar bn − r對照一下。如果把定理中的a、b和n分別設定為−1、1和i − s,便可以把式(2)進一步變成
C(i, s) × (−1 + 1)i − s = 0 (註6)
至此我們證明了當i > s時,有關元素在公式(1)中將被加0次。綜合上述三種情況,我們證明了公式(1)是正確 的。

請注意若把s = 0代入公式(1),便得到
N(n, 0) = Σ0 ≤ k ≤ n(−1)k Sn, k     (3)
上式等同於《點算的奧秘:容斥原理基本公式》中的公式(4)(不要忘記我 們把Sn, 0定義為|U|),即|S1' ∩ S2' ∩ ... Sn'| 的公式,由此可見公式(1)的確是「容斥原理」公式的推廣。

若S1 ... Sn是「可交換集合」,那麼公式(1)中的Sn, k便可以簡化如下 :就某一個特定的k而言,Sn, k是由C(n, k)個相同的數加起來的總和,我們可以把這個相同的數 記為nk。這樣便可以把公式(1)改寫為
N(n, s) = Σs ≤ k ≤ n(−1)k − s C(k, s) C(n, k) nk
利用上面C(k, s) × C(i, k) = C(i, s) × C(i − s, k − s)的關係,上式可以化簡 為
N(n, s) = C(n, s) Σs ≤ k ≤ n(−1)k − s C(n − s, k − s) nk
再對上式的求和變項k進行變換(利用r = k − s,即k = r + s),上式最終變為
N(n, s) = C(n, s) Σ0 ≤ r ≤ n−s(−1)r C(n − s, r) nr+s    (4)
同樣,若把s = 0代入上式,便得到
N(n, 0) = Σ0 ≤ r ≤ n(−1)r C(n, r) nr     (5)
請注意上式等同於《點算的奧秘:容斥原理基本公式》中的公式(5)。

例題1:從A、B ... J這10類字母中可重覆地抽5個出來排列,若規定所得序列必須包含剛好4類字 母(例如A-B-C-B-D是合格的序列,而A-A-C-D-E則是不合格的序列),問有多少種合格的序列?

答1:設SA為從上述字母中可重覆地抽5個出來排列但不包括字母A的序列組成的集合, SB ... SJ的定義照此類推。由於包含剛好4類字母等於不包含剛好6類字母,本題所求 答案為N(10, 6)。顯然SA ... SJ是「可交換集合」,所以可用上述公式(4)求解。

接著讓我們求nk的值。對1 ≤ k ≤ 10而言,nk表示SA ... SJ中k個集合構成的交集的元素個數,這個交集就是從上述字母中可重覆地抽5個出來排列但不包括 其中k個字母的序列組成的集合。由於不包括其中k個字母,所以有關序列乃由餘下的10 − k個字母構成 ,即從10 − k個字母中可重覆地抽5個出來排成的序列。這是一個「無限重覆排列問題」,應用有關公式, 得nk = (10 − k)5,而nr+6 = (10 − r − 6)4 = (4 − r)5

最後應用公式(4),求得答案為
N(10, 6)= C(10, 6) Σ0 ≤ r ≤ 4(−1)r C(4, r) (4 − r)5
 = 210 × (1 × 45 − 4 × 35 + 6 × 25 − 4 × 15 + 1 × 05)
 = 50400 □
以上筆者介紹了「剛好s個」集合的情況,接下來討論「至少p個」集合的情況。設有S1 ... Sn等n個集合,現在要求屬於這n個集合中的「至少p個」集合的那些元素的個數。我們把這個 數記為L(n, p)(註7),現在讓我們嘗試推導這個數的公式。由於「n個集合中的至少p個」集合乃由以下各個情 況構成:剛好p個集合;剛好p + 1個集合 ... 剛好n個集合,根據「加法原理」並應用公式(1),應有
L(n, p)= Σp ≤ k ≤ n N(n, k)
 = Σp ≤ k ≤ n Σk ≤ r ≤ n (−1)r − k C(r, k) Sn, r    (6)
接著讓我們化簡上式,這要用到多種技巧。為了方便後面的推導,首先要把上式的兩個Σ符號對調位置, 但由於第二個Σ的下限依賴於第一個Σ而變化,我們在對調兩個Σ符號時必須對其上、下限作 出調整,這實際上等於互換兩個變項的依賴關係。為方便讀者理解,現把兩個變項k和r的取值範圍列表如下:
k的取值範圍k取定左欄的值後r的取值範圍
pp, p + 1, ... n − 1, n
p + 1p + 1, p + 2, ... n − 1, n
............
n − 1n − 1, n
nn
從上表可以看到k與r的取值範圍都是最小p,最大n,但兩者必須符合這樣的關係:k ≤ r。把式(6)的兩個 Σ對調後,應該保存這種關係,因此我們應有
L(n, p) = Σp ≤ r ≤ n Σp ≤ k ≤ r (−1)r − k C(r, k) Sn, r
接著我們把上式中(−1)的冪次r − k分拆為r − p + p − k,把上式變為
L(n, p) = Σp ≤ r ≤ n Σp ≤ k ≤ r (−1)r − p × (−1)p − k C(r, k) Sn, r
接著我們把上式中第二個Σ後不含變項k的因子統統抽到第二個Σ前面,並且利用 (−1)p − k = (−1)k − p,從而把上式進一步變為
L(n, p) = Σp ≤ r ≤ n (−1)r − p Sn, r × Σp ≤ k ≤ r(−1)k − p C(r, k)    (7)
接著我們證明Σp ≤ k ≤ r(−1)k − p C(r, k) = C(r − 1, p − 1),這樣便可以消去其中一個Σ。證明方法是這樣的:首先利用《點算的奧秘:二項式定理和多項式定理》中的公式(2)把C(r, k)寫成C(r − 1, k − 1) + C(r − 1, k)。這樣我們便得到
Σp ≤ k ≤ r(−1)k − p C(r, k) = Σp ≤ k ≤ r(−1)k − p [C(r − 1, k − 1) + C(r − 1, k)]
 = Σp ≤ k ≤ r(−1)k − p C(r − 1, k − 1) + Σp ≤ k ≤ r(−1)k − p C(r − 1, k)
 = Σp ≤ k ≤ r(−1)k − p C(r − 1, k − 1) + Σp+1 ≤ k ≤ r+1(−1)k − p − 1 C(r − 1, k − 1)
請注意上面最後一行其實應用了對求和變項k進行變換的運算,具體地說,就是把k寫成k − 1,並把k的 上、下限同時加1。現在我們看到上面兩個求和項具有相同的形式C(r − 1, k − 1),兩者的差別 僅在於上、下限以及(−1)的冪次不同。由於前後兩個求和項中(−1)的冪次相差1,兩者中k取值相 同的項都互相抵消,只剩下前者k值為p和後者k值為r + 1的項,即
Σp ≤ k ≤ r(−1)k − p C(r, k)= (−1)p − pC(r − 1, p − 1) + (−1)r + 1 − p − 1 C(r − 1, r + 1 − 1)
 = C(r − 1, p − 1) + (−1)r − p C(r − 1, r)
由於r − 1 < r,我們有C(r − 1, r) = 0。由此我們證得Σp ≤ k ≤ r (−1)k − p C(r, k) = C(r − 1, p − 1)。把此一結果代入上面式(7), 我們最終得到
L(n, p) = Σp ≤ r ≤ n(−1)r − p C(r − 1, p − 1) Sn, r    (8)
上式就是「廣義容斥原理」的另一種形式。請注意若把p = 1代入公式(8),便得到
L(n, 1) = Σ1 ≤ r ≤ n(−1)r − 1 Sn, r    (9)
上式等同於《點算的奧秘:容斥原理基本公式》中的公式(2),即 |S1 ∪ S2 ∪ ... Sn|的公式,由此可見公式(8)的確是「容斥原 理」公式的推廣。

若S1 ... Sn是「可交換集合」,那麼公式(8)中的Sn, r便變成C(n, r) nr。這樣便可以把公式(8)改寫為
L(n, p)= Σp ≤ r ≤ n(−1)r − p C(r − 1, p − 1) × C(n, r) nr
 = Σp ≤ r ≤ n(−1)r − p C(r − 1, p − 1) × C(n − 1, r − 1) × (n / r) nr
接著利用上面C(k, s) × C(i, k) = C(i, s) × C(i − s, k − s)的關係,上式可以 化簡為
L(n, p) = n × C(n − 1, p − 1) Σp ≤ r ≤ n (−1)r − p C(n − p, r − p) × nr / r
對上式的求和變項r進行變換,變換公式為j = r − p,我們最終得到
L(n, p) = n × C(n − 1, p − 1) Σ0 ≤ j ≤ n−p (−1)j C(n − p, j) × nj+p / (j + p)    (10)
同樣,若把p = 1代入上式,便得到
L(n, 1) = n × Σ0 ≤ j ≤ n−1(−1)j C(n − 1, j) × nj+1 / (j + 1)    (11)
只要利用k = j + 1對上式進行變換,並利用關係式C(n − 1, k − 1) × n / k = C(n, k) ,便可證明上式等同於《點算的奧秘:容斥原理基本公式》中的公式(3)。

例題2:從A、B ... J這10類字母中可重覆地抽5個出來排列,若規定所得序列最多只可包含4類字母 (例如A-A-B-B-C是合格的序列,而A-B-C-D-E則是不合格的序列),問有多少種合格的序列?

答2:首先我們沿用跟答1相同的集合定義SA ... SJ。由於包含最多4類字母 等同於不含至少6類字母,本題所求答案為L(10, 6)。根據答1,nj+6 = (4 − j)5。因此根據公式(10),我們求得答案為
L(10, 6)= 10 × C(9, 5) × Σ0 ≤ j ≤ 4 (−1)j C(4, j) × (4 − j)5 / (j + 6)
 = 10 × 126 × (1 × 45 / 6 − 4 × 35 / 7 + 6 × 25 / 8 − 4 × 15 / 9 + 1 × 05 / 10)
 = 69760 □
在本節中,筆者除了介紹「廣義容斥原理」外,還介紹了「點算組合學」上常用的多種運算技巧,包括對求和 變項進行變換、把兩個求和符號Σ對調位置以及應用「點算組合學」的某些常用公式等,這些技巧在以後 還會經常用到,請讀者務必掌握。
註1:按照這個記法,|S1' ∩ S2' ∩ ... Sn'|就等於N(n, 0)。

註2:請注意本文對公式(1)的證明只能證明它的正確性,而不能揭示該公式是如何推導出來的。這一點與「數 學歸納法」的證明有點類似。

註3:因為如果某元素出現於Sn, s中的兩個或更多不同的項,那麼這個元素就一定同時屬於不只s 個集合。舉例說,設s為3,並設某元素x既出現於S1 ∩ S2 ∩ S3 ,又出現於S2 ∩ S3 ∩ S4,那麼x便同時屬於4個而非3個集合。

註4:由於C(i, s)的i和s兩者都不等於式(2)求和符號Σ中的變項k,所以我們可以把C(i, s)從Σ的 右邊抽出來。這種做法就正如我們可以從(a × b + a × c)中把a抽出來,使之成為a × (b + c)一樣。

註5:學過微積分的讀者對這種變換應不會感到陌生,因為在進行積分時經常要對積分變項進行變換。事實上, 我們可以把「求和」(Summation)看成一種「離散」的「積分」(Integration),或者把「積分」看成一種「連 續」的「求和」。

註6:由於我們在前面假定了i > s,即i − s > 0,所以這裡(−1 + 1)i − s有 定義,即等於0。

註7:按照這個記法,|S1 ∪ S2 ∪ ... Sn|就等於L(n, 1)。



返回數學專題

容斥原理与应用

写了半天,原来这个写的很好了
http://chowkafat.net/Enumeration6.html
 

點算的奧秘:容斥原理基本公式


「容斥原理」(Principle of Inclusion and Exclusion)(亦作「排容原理」)是「點算組合學」中的一條重要 原理。但凡略為複雜、包含多種限制條件的點算問題,都要用到這條原理。現在首先從一個點算問題說起。

例題1:設某班每名學生都要選修至少一種外語,其中選修英語的學生人數為25,選修法語的學生人 數為18,選修德語的學生人數為20,同時選修英語和法語的學生人數為8,同時選修英語和德語的學生人數為13 ,同時選修法語和德語的學生人數為6,而同時選修上述三種外語的學生人數則為3,問該班共有多少名學生?

答1:我們可以把上述問題表達為下圖:
其中紅色、綠色和藍色圓圈分別代表選修英語、法語和德語的學生。根據三個圓圈之間的交叉關係,可把上圖 分為七個區域,分別標以A至G七個字母。如果我們用這七個字母分別代表各字母所在區域的學生人數,那麼根 據題意,我們有以下七條等式:(1) A+D+E+G = 25;(2) B+D+F+G = 18;(3) C+E+F+G = 20;(4) D+G = 8; (5) E+G = 13;(6) F+G = 6;(7) G = 3。現在我們要求的是A+B+C+D+E+F+G。如何利用以上資料求得答案?

把頭三條等式加起來,我們得到A+B+C+2D+2E+2F+3G = 63。可是這結果包含了多餘的D、E、F和G,必須設法把 多餘的部分減去。由於等式(4)-(6)各有一個D、E和F,若從上述結果減去這三條等式,便可以把多餘的D、E和 F減去,得A+B+C+D+E+F = 36。可是這麼一來,本來重覆重現的G卻變被完全減去了,所以最後還得把等式(7)加 上去,得最終結果為A+B+C+D+E+F+G = 39,即該班共有39名學生。□

在以上例題中,給定的資料是三個集合的元素個數以及這些集合之間的交集的元素個數。在該題的解答中,我 們交替加上及減去這些給定的資料。如果我們用S1、S2和S3分別代表選修 英語、法語和德語學生的集合,那麼我們要求的答案就是|S1 ∪ S2 ∪ S3|,而該題的解答則可以重新表達為
|S1 ∪ S2 ∪ S3| = (|S1| + |S2| + |S3|) − (|S1 ∩ S2| + |S1 ∩ S3| + |S2 ∩ S3|) + |S1 ∩ S2 ∩ S3|
我們可以把上式推廣至集合個數為n的情況,便得到以下的「容斥原理」。設有n個集合 S1、S2 ... Sn,那麼
|S1 ∪ S2 ∪ ... Sn| = (|S1| + |S2| + ... |Sn|)
  − (|S1 ∩ S2| + |S1 ∩ S3| + ... |Sn − 1 ∩ Sn|)
  + (|S1 ∩ S2 ∩ S3| + |S1 ∩ S2 ∩ S4| + ... |Sn − 2 ∩ Sn − 1 ∩ Sn|)
  − (|S1 ∩ S2 ∩ S3 ∩ S4| + ... |Sn − 3 ∩ Sn − 2 ∩ Sn − 1 ∩ Sn|)
  ......
  + (−1)n − 1 |S1 ∩ S2 ∩ ... Sn − 1 ∩ Sn|     (1)
有必要對上式作一些解釋。為易於理解,上式分數行列出,每一行都是一些集合元素個數的總和,其中第1行包 含全部n個集合S1、S2 ... Sn;第2行包含所有由2個集合構成的交集,應 共有C(n, 2)個項;第3行包含所有由3個集合構成的交集,應共有C(n, 3)個項...第n行包含由全部n個集合構成 的交集,這樣的交集只有C(n, n) = 1個。每行的開首交替為一個「加」(相當於「容納」)和一個「減」(相當 於「排斥」)號,第1行開首為「加」號,第2行為「減」號,第3行為「加」號,第4行為「減」號...。由於當k 為奇數時,(−1)k = −1;當k為偶數時,(−1)k = 1,我們可以把 第1行開首的「加」號改寫為「+ (−1)0」,把第2行開首的「減」號改寫為 「+ (−1)1」...如此類推,我們可知第k行開首應有一個「+ (−1)k − 1 」。

我們還可以把上式化簡,方法是引入一個新的變項Sn, k來代表上式等號右邊的第k行。 Sn, k是由C(n, k)個項加起來的總和,每個項都是由S1、S2 ... Sn這n個集合中抽r個出來構成的交集的元素個數。舉例說,當n = 5,k = 3時,S5, 3 是由C(5, 3) = 10個項加起來的總和,這10個項都是由S1 ... S5這5個集合中抽3個出 來構成的交集的元素個數,其中一個是|S2 ∩ S4 ∩ S5|,其餘類 推。利用Sn, k以及Σ求和符號,我們便可以把上面的「容斥原理」公式大大簡化為:
|S1 ∪ S2 ∪ ... Sn| = Σ1 ≤ k ≤ n (−1)k − 1 Sn, k     (2)
接著讓我們證明以上「容斥原理」公式(1)的正確性,為此我們必須證明,S1 ... Sn 這n個集合中的每一個元素在公式中都只被加一次。我們逐一考慮各種集合元素的情況。首先考慮那些只包含於 一個集合中的元素,每個這類元素只出現於上述公式的第1行n個集合的其中一個,因此只會被加一次。其次考 慮那些包含於兩個集合的元素,每個這類元素都出現於上述公式的第1行n個集合的其中兩個,並且出現於第2行 C(n, 2)個交集的其中一個。由於每個這類元素在第1行中被加兩次並在第2行中被減一次,因此結果該元素在公 式中被加一次。

我們把上述推理推廣至一般情況,即考慮那些包含於k個集合的元素。每個這類元素都出現於第1行的其中 C(k, 1) = k個集合,並且出現於第2行的其中C(k, 2)個交集,並且出現於第3行的其中C(k, 3)個交集...並且 出現於第k行的其中C(k, k) = 1個交集(註1)。總括而言,每個這類元素在公式中被加的次數為
C(k, 1) − C(k, 2) + C(k, 3) − C(k, 4) ... (−)k − 1 C(k, k) = Σ1 ≤ r ≤ k(−1)r − 1 C(k, r)
現在的問題是,如何證明上式等於1?答案是借助《點算的奧秘:二項式定理和 多項式定理》中介紹的「二項式定理」:(a + b)n = Σ 0 ≤ r ≤ n C(n, r) ar bn − r。為了應用「二項式定理」,我們首先把上式等號右邊重 寫:
Σ1 ≤ r ≤ k(−1)r − 1 C(k, r) = − Σ1 ≤ r ≤ k(−1)r C(k, r)
  = − [Σ0 ≤ r ≤ k(−1)r C(k, r) − (−1)0 C(k, 0)]
  = − Σ0 ≤ r ≤ kC(k, r) (−1)r + 1
現在如果把「二項式定理」中的a、b和n分別設定為−1、1和k,便得到Σ0 ≤ r ≤ k C(k, r) (−1)r = (−1 + 1)k = 0。把這個結果代入上面最後一行 ,我們便得到Σ1 ≤ r ≤ k(−1)r − 1 C(k, r) = −0 + 1 = 1。至此我們證明了每一個元素在「容斥原理」公式中都被加一次,因此該公式是正確的。

在應用上述「容斥原理」公式(1)時,我們必須逐一找出所有集合以及各種交集的元素個數,有時這是頗為繁瑣 的工作,會令上述公式失去價值,因此「容斥原理」一般應用於「可交換集合」(Exchangeable Sets)。如果從 n個集合S1、S2 ... Sn中任意抽取k個出來(1 ≤ r ≤ n),所得交集 的元素個數只跟k有關,而跟抽取結果無關,則我們說這n個集合是「可交換」的。換句話說,不論是由哪k個集 合構成的交集,所得交集的元素個數都是同一個數字。

上面「容斥原理」公式(1)的第k行的每一項都是由k個集合構成的交集的元素個數。如果S1、 S2 ... Sn是「可交換集合」,那麼該行的每一項都是同一個數字,可不妨用符號 nk來代表。由於該行共有C(n, k)個項,所以該行可簡化為(−1)k − 1 C(n, k) nk。因此,「可交換集合」的「容斥原理」公式便是
|S1 ∪ S2 ∪ ... Sn| = Σ1 ≤ k ≤ n (−1)k − 1 C(n, k) nk     (3)
在某些應用中,所求的答案是不屬任何指定集合的元素個數,因此我們需要找出以上「容斥原理」公式的變化 形式。設某論域U包含n個集合S1、S2 ... Sn,現在我們要求不屬這n個集 合的任何一個的元素個數,即|S1' ∩ S2' ∩ ... Sn'|(註2)。根 據集合論的「德.摩根律」(De Morgan's Law),S1' ∩ S2' ∩ ... Sn'等於U − (S1 ∪ S2 ∪ ... Sn),因此 |S1' ∩ S2' ∩ ... Sn'| = |U| − |S1 ∪ S2 ∪ ... Sn| = |U| + (−|S1 ∪ S2 ∪ ... Sn|)(註3)。應用上面的「容斥原理」公式(2),我們便得到
|S1' ∩ S2' ∩ ... Sn'| = |U| + Σ1 ≤ k ≤ n(−1)k Sn, k     (4)
同理,如果S1、S2 ... Sn是「可交換集合」,那麼上式便變成
|S1' ∩ S2' ∩ ... Sn'| = |U| + Σ1 ≤ k ≤ n(−1)k C(n, k) nk     (5)
現在讓我們來看看「容斥原理」的一些簡單應用例子。有關該原理的其他應用,以後還會談到。

例題2:從26個英文字母中(可重覆地)抽取10個出來排成字串(String)(無需考慮這個字串是否構成一 個英文單詞),其中不可包含全部5個元音字母(即A、E、I、O、U),問有多少種排法?

答2:利用「容斥原理」解題的關鍵在於,把原題表達為適當的集合。就本題而言,如果我們用 SA來代表所有由10個字母組成但不含A的字串集合...用SU來代表所有由10個字母組成 但不含U的字串集合,那麼本題所要求的答案就是|SA ∪ SE ∪ SI ∪ SO ∪ SU|。由於5個元音字母具有平等的地位,從SA、 SE、SI、SO和SU這5個集合中任意抽取k個出來所構成的交集 的元素個數都是一樣的,所以這5個集合是「可交換」的。因此我們可以用上面的公式(3)來解題,但須先求得 nk (1 ≤ k ≤ 5)的值。

我們先求n1的值,這個值相當於由10個字母組成但不含A、E、I、O、U中某一個字母的字串的數目 。由於撇除了一個字母,我們可以從餘下的25個字母中(可重覆地)抽取10個出來排成字串,這是一個「無限重 覆排列」問題,應共有2510個這樣的字串,即n1 = 2510

接著考慮n2,這個值相當於由10個字母組成但不含A、E、I、O、U中某兩個字母的字串的數目。由 於撇除了兩個字母,我們可以從餘下的24個字母中(可重覆地)抽取10個出來排成字串,因此n2 = 2410。讀到這裡,相信讀者應能看到nk = (26 − k)10

最後,運用上面公式(3),本題的答案是
Σ1 ≤ k ≤ 5(−1)k − 1 C(5, k) (26 − k) 10 = 5 × 2510 − 10 × 2410 + 10 × 2310 − 5 × 2210 + 2110
由於具體數字太大,這裡只列出其算式。事實上,「點算組合學」經常要和這樣的「天文數字」打交道。若非 靠邏輯思維,實在難以想像人類有何方法求得這類問題的答案,由此可見邏輯思維的巨大力量!□

例題3:現要把一對可區別的火星人、一對可區別的木星人和一對可區別的土星人安排坐在一排共6個 座位上。如規定同一類外星人不得坐在相鄰的位置上,問有多少種不同坐法?

答3:如果我們用U代表由各種不加限制的坐法組成的集合,並用SM(SJ、 SS)代表由兩個火星人(木星人、土星人)坐在相鄰位置的各種坐法組成的集合,那麼本題要求的答 案就是|SM' ∩ SJ' ∩ SS'|。本題的SM、 SJ和SS顯然是「可交換」的集合,因此可以用上面的公式(5)解題。

首先求|U|的值。由於6個外星人是可區別的,我們不妨用M1、M2、J1、 J2、S1、S2代表他們(其中M、J和S分別代表火星人、木星人和土星人)。 因此|U|的值等於把上述6個符號進行「全排列」的問題,故|U| = 6!。

接著求n1的值,這個值相當於兩個火星人(或木星人、土星人)坐在相鄰位置的各種坐法的總數。由 於M1、M2總是聚在一塊,我們可以先把M1和M2合併為一個符 號M,然後求這個M與其餘4個符號進行「全排列」的總數,這個數應為5!。接著考慮M的「內部構造」, M1和M2有兩種可能形式構成M,即M1-M2和 M2-M1。總上所述,得n1 = 5! × 2。

接著求n2的值,這個值相當於兩個火星人以及兩個木星人(或其他組合)坐在相鄰位置的各種坐法的 總數。類似上段的做法,我們可以先把M1和M2以及J1和J2合 併為兩個新的符號M和J,然後求M、J與其餘2個符號進行「全排列」的總數,這個數應為4!。接著考慮M和J的 「內部構造」,各有兩種可能形式。總上所述,得n2 = 4! × 22

接著求n3的值,這個值相當於兩個火星人、兩個木星人和兩個土星人坐在相鄰位置的各種坐法的總 數。同樣,我們可以先把M1和M2、J1和J2以及S1 和S2合併為三個新的符號M、J和S,然後求M、J和S進行「全排列」的總數,這個數應為3!。接著考 慮M、J和S的「內部構造」,各有兩種可能形式。總上所述,得n3 = 3! × 23

最後,運用上面的公式(5),求得本題答案為
  |U| − C(3, 1) n1 + C(3, 2) n2 − C(3, 3) n3
= 6! − 3 × 5! × 2 + 3 × 4! × 22 − 1 × 3! × 23
= 720 − 720 + 288 − 48
= 240 □


註1:如果讀者不明白本段的推理,可以嘗試用一個具體的例子代入上述推理,這樣便能較易明白。設n = 5,k = 4,並且假設某元素x包含於S1、S3、S4和S5這4個集合中。 那麼x出現於第1行的其中C(4, 1) = 4個集合,並且出現於第2行的其中C(4, 2) = 6個交集,這6個交集分別是 S1 ∩ S3、S1 ∩ S4、S1 ∩ S5、S3 ∩ S4、S3 ∩ S5和 S4 ∩ S5。在第3行中,x出現於C(4, 3) = 4個交集,這4個交集分別是 S1 ∩ S3 ∩ S4、S1 ∩ S3 ∩ S5、S1 ∩ S4 ∩ S5和S3 ∩ S4 ∩ S5。在4行中,x只出現於C(4, 4) = 1個交集,此即S1 ∩ S3 ∩ S4 ∩ S5。請注意x不會出現於第5行。總括而言,x在公式 中被加的次數為4 − 6 + 4 − 1 = 1次。

註2:這裡用S'代表集合S的補集(Complement Set),即S' = U − S。

註3:這裡其實應用了集合論中的一個原理,即若T是S的子集,那麼|S − T| = |S| − |T|。



Friday, December 14, 2012

备用 将DD-MM-YYYY HH:MM:SS转化为SQLite 支持的YY-MM-DD HH:MM:SS


UPDATE t2 SET Comeback_Date = (SELECT CASE WHEN (SUBSTR(Comeback_Date,-8,1)=' ') THEN SUBSTR(Comeback_Date,-12,4) ELSE SUBSTR(Comeback_Date,-13,4) END)|| '-'|| (SELECT CASE WHEN (CAST(SUBSTR(Comeback_Date,1,2) AS INTEGER) >=10 and SUBSTR(Comeback_Date,5,1)!='-') THEN SUBSTR(Comeback_Date,4,2) WHEN (CAST(SUBSTR(Comeback_Date,1,2) AS INTEGER) >=10 and SUBSTR(Comeback_Date,5,1)='-') THEN SUBSTR(Comeback_Date,4,1)  WHEN (CAST(SUBSTR(Comeback_Date,1,2) AS INTEGER) <10 and SUBSTR(Comeback_Date,4,1)='-') THEN SUBSTR(Comeback_Date,3,1) ELSE SUBSTR(Comeback_Date,3,2) END)|| '-'|| (SELECT CASE WHEN (CAST(SUBSTR(Comeback_Date,1,2) AS INTEGER) >=10) THEN SUBSTR(Comeback_Date,1,2) ELSE SUBSTR(Comeback_Date,1,1) END) ||' '|| (SELECT CASE WHEN SUBSTR(Comeback_Date,-8,1)=' 'THEN SUBSTR(Comeback_Date,-7) ELSE SUBSTR(Comeback_Date,-8) END) where Comeback_Date is not null;


update t1 set comeback_date=(substr(comeback_date, 1,4)||'-'||
(select case when(cast(substr(comeback_date,6,2) as integer)>=10) then substr(comeback_date,6,2) else '0'||substr(comeback_date,6,1) end) ||'-'||
(select case when(cast(substr(comeback_date,8,2) as integer)>=10) then substr(comeback_date,8,2) when(cast(substr(comeback_date,9,2) as integer)>=10) then substr(comeback_date,9,2) when(cast(substr(comeback_date,8,2) as integer)>0) then '0'||substr(comeback_date,8,1) else '0'||substr(comeback_date,9,1) end) ||' '||(select case when (substr(comeback_date,-8,1)=' ') then '0'||substr(comeback_date,-7) else substr(comeback_date,-8) end)
)
where comeback_date is not null