分治策略 - 典型實例 - 選擇問題


選擇問題最常見的問題有:

1.1選最大

選擇算法
統一描述:設L是n個算法的集合,從L中選出第k小的元素,1<=k<=n,當L中元素按從小到大排好序后,排在第k個位置的數,就是第k小的數。
下面介紹 順序比較法
算法Findmax
輸入:n個數的數組L
輸出:max,k

max <- L[1]; k <- 1
for i <- 2 to n do      //for循環執行n-1次
    if max < L[i]
    then max <- L[i]
         k <- i
return max, k

算法Findmax第二行,for循環執行n-1次,所以 \(W(n)=n-1\)
這個算法是選最大問題在時間上最優的算法(對於選最小問題,只需對算法稍加改動,就可以得到順序比較的Findmin算法)

1.2同時選最大和最小的算法

設計思想:先選最大,然后把最大的從L中刪除,接着選最小。
算法:(利用Findmax和Findmin)
輸入:n個數的數組L
輸出:max,min

if n=1 then return L[1]作為max和min
else Findmax
    從L中刪除max
    Findmin

算法執行的比較次數:\(W(n)=n-1+n-2=2n-3\)

分組比賽的方法
基本思想:首先將L中的元素兩兩一組,分成\(\lfloor n/2 \rfloor\)組(當n是奇數時有一個元素輪空)。每組中的2個元素進行比較,得到組內較大和較少數,把至多\(\lfloor n/2 \rfloor + 1\)(當n為奇數時,把被輪空的元素加進來)個小組中較大的元素放在一起,運行Findmax,得到L中的最大元素,同理得到L中的最小元素。
算法FindMaxMin

將n個元素兩兩一組分成n/2(下取整)組
每組比較,得到n/2(下取整)個較小和n/2(下取整)個較大的數      //比較n/2(下取整)次
在n/2(下取整)個(n為奇數時,是n/2(下取整)+1)較小中找最小min
在n/2(下取整)個(n為奇數時,是n/2(下取整)+1)較大中找最大max

行3和行4都執行\(\lceil n/2 \rceil -1\)
所以\(W(n)=\lfloor n/2 \rfloor +2 \lceil n/2 \rceil-2=n+\lceil n/2 \rceil-2=\lceil 3n/2 \rceil-2\)
此算法效率更高,是所有同時找最大和最小算法中事件復雜度最低的算法。

1.3找第二大

2次調用Findmax算法
\(W(n)=n-1+n-2=2n-3\)

錦標賽算法
把數組中的元素兩兩一組,划分為\(\lfloor n/2 \rfloor\)組(n為奇數時1個元素輪空),每組組內兩個元素比大小,大的進入下一輪(n為奇數時,輪空的元素也進入下一輪)。
所以下一輪有\(\lceil n/2 \rceil\)個元素。繼續每組組內比大小,然后大的進入下一輪,直至找出max。篩掉\(n-1\)個元素,比較\(n-1\)次。

找第二大,不可能再用Findmax了,如果用Findmax,就又比較n-2次了,和上面的算法兩次調用Findmax一樣。
所以,我們可以利用找max時比較所產生的記錄幫我們減少比較次數。在比賽前為每個元素設定一個指針,指向一個鏈表,把比較后比它小的元素都記錄進它的鏈表中,找出max后,在max的鏈表中用Findmax找出整個數組第二大的元素。
算法FindSecond
輸入:n個數的數組L
輸出:second

k <- n
將k個元素兩兩一組,分成k/2(下取整)組
每組的兩個數比較,找到較大的
將被淘汰的較小的數在淘汰它的數所指向的鏈表中做記錄
if k為奇數 then k <- k/2(下取整)+1
           else k <- k/2(下取整)
if k>1 then goto 2
max <- 剩下的一個數
second <- max的鏈表中的最大

此時,第一輪兩兩比較的比較次數是\(n-1\),但還不知道max的鏈表中有多少個元素,
所以接下來求解max所淘汰掉的元素個數(這部分的工作量)
設本輪參與比較的有\(t\)個元素,經過分組淘汰后進入下一輪的元素數至多是\(\lceil t/2 \rceil\),下下一輪就是\(\lceil \lceil t/2 \rceil /2 \rceil = \lceil t/2^2 \rceil\)
假設k輪淘汰后只剩max,則\(\lceil n/2^k \rceil = 1\).
\(n=2^d\),那么\(k=d=logn=\lceil logn \rceil\)
所以max進行了\(d\)次比較,\(W(n)=n-1+\lceil logn \rceil\)
對於找第二大的問題,算法Findmax是時間復雜度最低的算法。

上面,我們只討論了一些特例情況,基本上使用順序比較或分組比較的方法,沒有明確的分治算法的特征。下面考慮一般性的選第k小問題的算法,用到分治策略。

2選第k小

輸入:數組S,S的長度n,正整數k,1<=k<=n
輸出:第k小的數
選第k小,可以使用排序算法,從小到大排好序后,選擇第k個即可,最好的排序算法的時間復雜度是\(O(nlogn)\)
下面考慮性能更好的分治算法,就是\(O(n)\)時間的算法,方便起見,假設數組S中元素彼此不等。
以S中的某個元素\(m^*\)作為划分標准,將S划分為兩個子數組S1和S2,把這個數組中比\(m^*\)小的都放入\(S_1\)的數組中,數組\(S_1\)的元素個數是\(|S_1|\)個;把這個數組中比\(m^*\)大的都放入\(S_2\)的數組中,數組\(S_2\)的元素個數是\(|S_2|\)個。

  • \(k<|S_1|\),則原問題歸納為在數組\(S_1\)中找第\(k\)小的子問題。
  • \(k=|S_1|+1\),則\(m^*\)就是要找的第\(k\)小元素。
  • \(k>|S_1|+1\),則原問題歸納為在數組\(S_2\)中找第\(n-|S_1|-1\)小的子問題。
    算法的關鍵是如何確定這個划分\(S\)的標准\(m^*\),它要具有以下 特征
  1. 尋找\(m^*\)的時間代價不能高於\(O(nlogn)\),如果直接尋找\(m^*\),時間應是\(O(n)\)。設選擇算法的時間復雜度為\(T(n)\),遞歸調用這個算法在\(S\)上的一個真子集M上尋找\(m^*\),應該使用\(T(cn)\)時間(\(c<1\),反映\(M\)的規模小於\(S\)
  2. 通過\(m^*\)划分的兩個子問題的大小分別記作\(|S_1|\)\(|S_2|\),每次遞歸調用時,子問題規模與原問題規模\(n\)的比都不超過\(d\)\(d<1\)),調用時間\(T(dn)\),並且應保證\(c+d<1\),否則方程\(T(n)=T(cn)+T(dn)+O(n)\)的解不會達到\(O(n)\)
    下面的分治算法采用遞歸調用的方法尋找\(m^*\),先將\(S\)分組,5個元素一組,共分成\(\lceil n/5 \rceil\)個組。在每組中取中位數,把這\(\lceil n/5 \rceil\)個中位數放入集合\(M\)中,然后使用選擇算法選出集合\(M\)中的中位數\(m^*\)。這次遞歸調用的子問題規模是原問題規模的\(1/5\),\(1/5\)即為上文的特征(1)的\(c\)
    算法Select(S,k)
    輸入:數組S,S的長度n,正整數k,1<=k<=n
    輸出:第k小的數
將S划分成5個一組,共n/5(上取整)個組
每組中找一個中位數,把這些中位數放到集合M中
m* <- Select(M,|M|/2)                  //選M中的中位數m*,將S中的數組划分成A,B,C,D四個集合
把A和D中的每個元素與m*比較,小的構成S1,大的構成S2
S1 <- S1並C ; S2 <- S2並B
if k=|S1|+1 then 輸出m*
else if k<=|S1|
    then Select(S1,k)
    else Select(S2,k-|S1|-1)

image
左下方的集合\(C\)中所以元素全部小於\(m^*\),右上角的集合B中的所有元素全部大於\(m^*\)。僅僅對於\(A\)\(D\)中的元素,我們不能確定它們是否大於或小於\(m^*\),所以在算法的行4加以比較,完成\(S1\)\(S2\)的划分。
下面分析時間復雜度
假設n是5的倍數,且\(n/5是奇數\),即\(n/5=2r+1\)
所以\(|A|=|D|=2r\)\(|B|=|C|=3r+2\)\(n=10r+5\)
若A,D的元素都小於\(m^*\),那么它們都加入至S1中,且下一步算法又在這個子問題上遞歸調用,這對應了歸約后子問題的規模的上界,也正好是時間復雜度最壞的情況。類似的,如果A,D的元素都大於\(m^*\),也會出現類似的情況。
以前者為例時,子問題的大小是
\(|A|+|C|+|D|=7r+2=7\frac{n-5}{10}+2=7\frac{n}{10}-1.5<\frac{7n}{10}\)
上式表明子問題規模最大不超過原問題的\(7/10\),這個參數\(7/10\)就是前文特征(2)中的\(d\)
所以最壞情況下的時間復雜度的遞推式:
\(W(n)<=W(\frac{n}{5})+W(\frac{7n}{10})+cn\)
image
image

\[\begin{aligned} W(n) & <= tn+0.9tn+0.9^2 tn+... \\ & = tn(1+0.9+0.9^2+...) \\ & = O(n) \end{aligned} \]

所以從上面的分析可以看出,這樣找的\(m^*\)完全滿足算法要求,兩次遞歸調用的參數分別是\(c=0.2\)\(d=0.7\)\(c+d=0.9<1\)
所以Select算法可以把時間復雜度降至線性時間。
分組時也可以選\(3\)個或\(7\)個元素,元素數的改變可能會改變\(m^*\)的特征,從而使得針對某些分組方法得到的\(c+d\)的值不再小於\(1\),這就會增加運算時間。

求解選擇問題的時間復雜度最低的算法就是O(n)時間的算法!


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM