跳躍表以及C++實現

原文鏈接:http://blog.csdn.net/u013011841/article/details/39158585

      問題: 2.2 若一個表從不修改,則可使用一種更爲簡單的方法來實現表的元素的查找。爲了有效地訪問第 i 個元素,向單鏈表的每個元素中添加第二個指針,使其指向表中其它元素來減少查找所需時間。

需要自己注意的是在獨立實現代碼的時候不僅僅只是需要簡單的瞭解大概的算法,而是

1)確定需要使用的數據結構,如今天的代碼中的SKNode, SkipList. 

2)如何初始化這個數據結構,這樣才能得到求解問題循環結束條件。(初始化條件很重要)

3)先自己按照算法流程,插入、查找、刪除一個一個寫好僞代碼,這樣才不至於在寫代碼的時候沒有關注點,思維混亂。可以先從最簡單也是最和核心的查找算法寫起。

4)最後是調試,碰到不對的步驟,可以添加打印,我現在對於數據結構還有太底層的調試時看不出來太多信息的,只能一步一步的試探。這也是一個軟肋,以後有空一定要再學學編程範式這種比較底層的東西,還有內存管理也很重要。

參考的博文地址如:

跳躍鏈表簡介

           二叉樹是一種常見的數據結構。它支持包括查找、插入、刪除等一系列操作。但它有一個致命的弱點,就是當數據的隨機性不夠時,會導致其樹形結構的不平衡,從而直接影響算法的效率。

              跳躍鏈表(Skip List)是1987年才誕生的一種嶄新的數據結構,它在進行查找、插入、刪除等操作時的期望時間複雜度均爲O(logn),有着近乎替代平衡樹的本領。而且最重要的一點,就是它的編程複雜度較同類的AVL樹,紅黑樹等要低得多,這使得其無論是在理解還是在推廣性上,都有着十分明顯的優勢。

              跳躍鏈表的最大優勢在於無論是查找、插入和刪除都是O(logn),不過由於跳躍鏈表的操作是基於概率形成的,那麼它操作複雜度大於O(logn)的概率爲,可以看出當n越大的時候失敗的概率越小。

              另外跳躍鏈表的實現也十分簡單,在平衡樹中是最易實現的一種結構。例如像複雜的紅黑樹,你很難在不依靠工具書的幫助下實現該算法,但是跳躍鏈表不一樣,你可以很容易在半個小時內就完成其實現。

              跳躍鏈表的空間複雜度的期望爲O(n),鏈表的層數期望爲O(logn).

如何改進普通的鏈表?

             我們先看看一個普通的鏈表


            可以看出查詢這個鏈表O(n),插入和刪除也是O(n).因此鏈表這種結構雖然節省空間,但是效率不高,那有沒有什麼辦法可以改進呢?

            我們可以增加一條鏈表做爲快速通道。這樣我們使用均勻分佈,從圖中可以看出L1層充當L0層的快速通道,底層的結點每隔固定的幾個結點出現在上面一層。


            我們這裏主要以查找操作來介紹,因爲插入和刪除操作主要的複雜度也是取決於查找,那麼兩條鏈表查找的最好的時間複雜度是多少呢?

            一次查找操作首先要在上層遍歷<=|L1|次操作,然後在下層遍歷<=(L0/L1)次操作,至多要經歷


            次操作,其中|L1|爲L1的長度,n爲L0的長度.

            那麼最好的時間複雜度,也就怎麼設置間隔距離才能使查找次數最少有

            我們對|L1|的長度求導得


            把上式代入函數,查找次數最小也就是.這意味着下層每隔個結點在上層就有一個結點作爲快速跑道。

那麼三條鏈表呢



同理那麼我們讓L2/L1=L1/L0,然後同樣列出方程,求導可得L2=,查找次數爲3*

........................................................................

第k條鏈條.....查找次數爲

我們這裏取k=logn,代入的查找次數爲2logn.

到此爲主,我們應該知道了,期望上最好的層數是logn層,而且上下層結點數比爲2,這樣查找次數常數最小,複雜度保持在O(logn)。

跳躍鏈表的結構

跳躍表由多條鏈構成(L0,L1,L2 ……,Lh),且滿足如下三個條件:

  • 每條鏈必須包含兩個特殊元素:+∞ 和 -∞(可以需要,可以不需要,我的實現是採用了-∞作爲header之後的 節點,即第一個節點,+∞作爲最後一個節點)
  • L0包含所有的元素,並且所有鏈中的元素按照升序排列。
  • 每條鏈中的元素集合必須包含於序數較小的鏈的元素集合。


結點結構源代碼

   我覺得如果不考慮空間效率是可以將數據在每層上面分開存儲,但是也可以辦到在一個n的空間存儲,跳躍表就是在這些數據之間建立links,使用links替代樹形結構,如Btreap樹堆,R-B Trees, AVL trees。 現在我給的是基於一個n空間存儲。


  1. class SKNode
  2. {
  3. public:
  4. int key;
  5. SKNode* forward[MAXLEVEL];
  6. SKNode()
  7. {
  8. key=0;
  9. for(int i =0;i<MAXLEVEL;i++)
  10. {
  11. forward[i]= NULL;
  12. }
  13. }
  14. SKNode& operator=(const SKNode* & node)
  15. {
  16. key=node->key;
  17. for(int i=0;i<MAXLEVEL;i++)
  18. {
  19. forward[i] = node->forward[i];
  20. }
  21. return *this;
  22. }
  23. };
  24. //skip list, it has a header, this header have maxlevel pointers
  25. class SkipList
  26. {
  27. public:
  28. SKNode *hdr; /* list Header */
  29. int listLevel; /* current level of list */
  30. int insert(int key);
  31. SKNode* search(int key);
  32. int deleteNode(int key);
  33. void printList();
  34. SkipList()
  35. {
  36. hdr = new SKNode;
  37. listLevel = 0;
  38. hdr->key = -INT_MAX;
  39. SKNode* end = new SKNode;
  40. SKNode* first = new SKNode;
  41. first->key=-INT_MAX;
  42. end->key=INT_MAX;
  43. for(int i =0;i<MAXLEVEL;i++)
  44. {
  45. hdr->forward[i]=first;
  46. hdr->forward[i]->forward[i] = end;
  47. }
  48. printList();
  49. }
  50. ~SkipList()
  51. {
  52. delete hdr;
  53. }
  54. };

MAXlevel可以是log(n)/log(2)

跳躍鏈表查找操作

目的:在跳躍表中查找一個元素x

   在跳躍表中查找一個元素x,按照如下幾個步驟進行:

      1、 p = hdr; 從最上層的鏈(Lh)的開頭開始,進入到各層

      2、 for i(listLevel-1, 0) , 從最上層的鏈(Lh)的開頭開始, 如L3到L0

      3、假設當前位置爲p,p初始值爲hdr, p在i層的下一個節點爲q = p->forward[i], 它向右指向的節點爲q(p與q不一定相鄰)。將x和q.key做比較

          (1)if  x>q.key 在i層繼續向右走,

               then    p = p->forward[i].

          (2) else x<=q.key,  

             then   i--, 即將當前指針直接下移到下一層。

      for循環結束,i變爲了0, 即在最後一層, p = p->foward[0]

    4、判斷p,

         (1) if p!= NULL && p->key = x

                 return p

         (2) return NULL

如我們查找29, 先從入口到l3,29>11,p = 11, 判斷知道往下走L2,29<30, 繼續往下L1,29<30,L1,29>15, 前進,29<=29, 結束, p在15的位置,所以最後指針調整到p = p->forward[0], 判斷p值。

(今天因爲我從listLevel開始,浪費了好多調試時間。。。。)


元素53的查找路徑

    下面是C++實現代碼:

  1. SKNode* SkipList::search(int key)
  2. {
  3. SKNode* current = new SKNode;
  4. current = hdr;
  5. int i = listLevel-1;
  6. for(;i>=0;i--)
  7. {
  8. while(current->forward[i]->key != INT_MAX && key>current->forward[i]->key)//key大於下一個數據的值。轉到本層下一個元素
  9. {
  10. current = current->forward[i];
  11. }
  12. //否則i--,轉到下一層
  13. }
  14. current = current->forward[0];
  15. if(current!= NULL && current->key == key)
  16. {
  17. cout<<"find"<<key<<endl;
  18. return current;
  19. }
  20. return NULL;
  21. }


跳躍鏈表插入操作

目的:向跳躍表中插入一個元素x

     首先明確,向跳躍表中插入一個元素,相當於在表中插入一列從S0中某一位置出發向上的連續一段元素。有兩個參數需要確定,即插入列的位置以及它的“高度”。

     1)關於在L0層的插入的位置,我們先利用跳躍表的查找功能, x<=q.key,所以x的位置一定是在x之後,q之前。 同樣可以推論在L1,L2 ,L listLevel層的位置是在循環中生成他們這層最後的位置,就是在search的while之後記錄一個這個位置爲S[i]。最後需要在所有的S[i]之後重連數據之間的鏈接。 (但是需要注意的是爲了,如果我們不想加入重複數據,需要判斷p->forward[0]的值,如果相等,就是找到了,不需要再插入。當然如果我們不介意重複數據,也可以不加這個判斷。)

     2)需要插入的高度,決定在L1-listLevel的哪些位置加。 而插入列的“高度”較前者來說顯得更加重要,也更加難以確定。由於它的不確定性,使得不同的決策可能會導致截然不同的算法效率。爲了使插入數據之後,保持該數據結構進行各種操作均爲O(logn)複雜度的性質,我們引入隨機化算法(Randomized Algorithms)。


僞代碼如下:

Skip List Insertion(x)

    p = hdr;

    newlevel = getlevel(); 

    s[listLevel] = hdr(注意需要全部初始化爲hdr,爲了newlevel增長了,但是增長的層次s[i]卻沒有數據,沒有初始化,應該從頭結點開始)

    for i(listLevel-1, 0) 

          while(q!=+∞ && x>key)

                    then    p = p->forward[i].

                   //else x<=q.key,  

                   //then   i--

           s[i]=p; (s[i]爲i層探索的最後一個節點,最後需要在這之後插入x)

   last = p->forward[0] //判斷是否相等,。。。。

//插入數據,重連鏈表

if(newlevel>level) level = newlevel

for  i(newlevel-1 - 0)

    node->forward[i]= s[i]->forward[i]

    s[i]->forward[i] = node


     我定義一個隨機決策模塊,它的大致內容如下,但是這個代碼不能保證完全隨機,其實每次的運行結果都是一樣的:

 

  1. int getInsertLevel()
  2. {
  3. int upcount = 0;
  4. for(int i=0;i<MAXLEVEL;i++)
  5. {
  6. int num = rand()%10;
  7. if(num<5)
  8. {
  9. upcount++;
  10. }
  11. }
  12. return upcount;
  13. }



如插入43,查找路徑如下。

                  43的下一個數接到的40下一個數45。

                  40的下一個數接到43


紫色的箭頭表示更新過的指針

  1. int SkipList::insert(int key)
  2. {
  3. int level = getInsertLevel();
  4. SKNode* node = new SKNode;
  5. node->key=key;
  6. SKNode *s[MAXLEVEL];
  7. SKNode* current = new SKNode;
  8. SKNode* last = new SKNode;
  9. for(int i =0;i<MAXLEVEL;i++)
  10. {
  11. s[i]=hdr->forward[i];//initiation
  12. }
  13. current = last = hdr;
  14. cout<<"hdr"<<hdr->key<<endl;
  15. int i = listLevel-1;
  16. for(;i>=0;i--)
  17. {
  18. while(current->forward[i]->key != INT_MAX && key>current->forward[i]->key)//key大於下一個數據的值。轉到本層下一個元素
  19. {
  20. current = current->forward[i];
  21. }
  22. s[i] = current;//保存每一層位置上的最後指針的前驅
  23. }
  24. last=current->forward[0];
  25. if(last != NULL && last->key == key)
  26. {
  27. cout<<"inset key:"<<key<<"already existed"<<endl;
  28. return 0;
  29. }
  30. if(level>listLevel)//更新層數
  31. {
  32. listLevel = level;
  33. }
  34. for(int k = 0; k <listLevel;k++)
  35. {
  36. node->forward[k]=s[k]->forward[k];
  37. s[k]->forward[k]=node;
  38. }
  39. if(level>listLevel)
  40. {
  41. listLevel = level;
  42. }
  43. return 1;
  44. }

跳躍鏈表的刪除

目的:從跳躍表中刪除一個元素x

   刪除鏈表和插入幾乎一模一樣的,只是在最後重接鏈表不同:

   在跳躍表中查找到這個元素的位置,如果未找到,則退出 

  否則將該元素所在整列從表中刪除 

   將多餘的“空鏈”刪除 

Skip List Deletion(x)

    p = hdr;

    //newlevel = getlevel(); 

    s[listLevel] = hdr(注意需要全部初始化爲hdr,爲了newlevel增長了,但是增長的層次s[i]卻沒有數據,沒有初始化,應該從頭結點開始)

    fori(listLevel-1, 0) 

          while(q!=+∞ && x>key)

                    then    p = p->forward[i].

                   //else x<=q.key,  

                   //then   i--

           s[i]=p; (s[i]爲i層探索的最後一個節點,最後需要在這之後插入x)

   last = p->forward[0] //判斷是否相等,。。。。

if(last->key != x)

 return

//刪除數據,重連鏈表

for  i(listlevel-1 - 0)

   s[i]->forward[i]=s[i]->forward[i]->forward[i]; 

這段代碼如下:

 

  1. intSkipList::deleteNode(int key)
  2. {
  3. SKNode *s[MAXLEVEL];
  4. SKNode* current = new SKNode;
  5. SKNode* last = new SKNode;
  6. for(inti =0;i<MAXLEVEL;i++)
  7. {
  8. s[i]=hdr->forward[i];//initiation
  9. }
  10. current = last = hdr;
  11. for(inti = listLevel-1;i>=0;i--)
  12. {
  13. while(current->forward[i]->key != INT_MAX && key>current->forward[i]->key)//key大於下一個數據的值。轉到本層下一個元素
  14. {
  15. current = current->forward[i];
  16. }
  17. s[i] = current;//保存每一層位置上的最後指針的前驅
  18. }
  19. last=current->forward[0];
  20. if(last->key != key)
  21. {
  22. cout<<"delete key:"<<key<<"does not existed"<<endl;
  23. return 0;
  24. }
  25. for(inti = 0; i<listLevel;i++)
  26. {
  27. s[i]->forward[i]=s[i]->forward[i]->forward[i];
  28. }
  29. return 1;
  30. }

整個C++程序如下

  1. #include <iostream>
  2. #include <vector>
  3. using namespace std;
  4. #define MAXLEVEL 4 //最多2 power n=16個數
  5. /*skip list node,they are keys and pointers*/
  6. classSKNode
  7. {
  8. public:
  9. int key;
  10. SKNode* forward[MAXLEVEL];
  11. SKNode()
  12. {
  13. key=0;
  14. for(inti =0;i<MAXLEVEL;i++)
  15. {
  16. forward[i]= NULL;
  17. }
  18. }
  19. SKNode& operator=(constSKNode* & node)
  20. {
  21. key=node->key;
  22. for(inti=0;i<MAXLEVEL;i++)
  23. {
  24. forward[i] = node->forward[i];
  25. }
  26. return *this;
  27. }
  28. };
  29. //skip list, it has a header, this header have maxlevel pointers
  30. classSkipList
  31. {
  32. public:
  33. SKNode *hdr; /* list Header */
  34. intlistLevel; /* current level of list */
  35. int insert(int key);
  36. SKNode* search(int key);
  37. intdeleteNode(int key);
  38. voidprintList();
  39. SkipList()
  40. {
  41. hdr = new SKNode;
  42. listLevel = 0;
  43. hdr->key = -INT_MAX;
  44. SKNode* end = new SKNode;
  45. SKNode* first = new SKNode;
  46. first->key=-INT_MAX;
  47. end->key=INT_MAX;
  48. for(inti =0;i<MAXLEVEL;i++)
  49. {
  50. hdr->forward[i]=first;
  51. hdr->forward[i]->forward[i] = end;
  52. }
  53. printList();
  54. }
  55. ~SkipList()
  56. {
  57. deletehdr;
  58. }
  59. };
  60. intgetInsertLevel()
  61. {
  62. intupcount = 0;
  63. for(inti=0;i<MAXLEVEL;i++)
  64. {
  65. intnum = rand()%10;
  66. if(num<5)
  67. {
  68. upcount++;
  69. }
  70. }
  71. returnupcount;
  72. }
  73. SKNode* SkipList::search(int key)
  74. {
  75. SKNode* current = new SKNode;
  76. current = hdr;
  77. inti = listLevel-1;
  78. for(;i>=0;i--)
  79. {
  80. while(current->forward[i]->key != INT_MAX && key>current->forward[i]->key)//key大於下一個數據的值。轉到本層下一個元素
  81. {
  82. current = current->forward[i];
  83. }
  84. //否則i--,轉到下一層
  85. }
  86. current = current->forward[0];
  87. if(current!= NULL && current->key == key)
  88. {
  89. cout<<"find"<<key<<endl;
  90. return current;
  91. }
  92. return NULL;
  93. }
  94. intSkipList::insert(int key)
  95. {
  96. int level = getInsertLevel();
  97. SKNode* node = new SKNode;
  98. node->key=key;
  99. SKNode *s[MAXLEVEL];
  100. SKNode* current = new SKNode;
  101. SKNode* last = new SKNode;
  102. for(inti =0;i<MAXLEVEL;i++)
  103. {
  104. s[i]=hdr->forward[i];//initiation
  105. }
  106. current = last = hdr;
  107. cout<<"hdr"<<hdr->key<<endl;
  108. inti = listLevel-1;
  109. for(;i>=0;i--)
  110. {
  111. while(current->forward[i]->key != INT_MAX && key>current->forward[i]->key)//key大於下一個數據的值。轉到本層下一個元素
  112. {
  113. current = current->forward[i];
  114. }
  115. s[i] = current;//保存每一層位置上的最後指針的前驅
  116. }
  117. last=current->forward[0];
  118. if(last != NULL && last->key == key)
  119. {
  120. cout<<"inset key:"<<key<<"already existed"<<endl;
  121. return 0;
  122. }
  123. if(level>listLevel)//更新層數
  124. {
  125. listLevel = level;
  126. }
  127. for(int k = 0; k <listLevel;k++)
  128. {
  129. node->forward[k]=s[k]->forward[k];
  130. s[k]->forward[k]=node;
  131. }
  132. if(level>listLevel)
  133. {
  134. listLevel = level;
  135. }
  136. return 1;
  137. }
  138. intSkipList::deleteNode(int key)
  139. {
  140. SKNode *s[MAXLEVEL];
  141. SKNode* current = new SKNode;
  142. SKNode* last = new SKNode;
  143. for(inti =0;i<MAXLEVEL;i++)
  144. {
  145. s[i]=hdr->forward[i];//initiation
  146. }
  147. current = last = hdr;
  148. for(inti = listLevel-1;i>=0;i--)
  149. {
  150. while(current->forward[i]->key != INT_MAX && key>current->forward[i]->key)//key大於下一個數據的值。轉到本層下一個元素
  151. {
  152. current = current->forward[i];
  153. }
  154. s[i] = current;//保存每一層位置上的最後指針的前驅
  155. }
  156. last=current->forward[0];
  157. if(last->key != key)
  158. {
  159. cout<<"delete key:"<<key<<"does not existed"<<endl;
  160. return 0;
  161. }
  162. for(inti = 0; i<listLevel;i++)
  163. {
  164. s[i]->forward[i]=s[i]->forward[i]->forward[i];
  165. }
  166. return 1;
  167. }
  168. voidSkipList::printList()
  169. {
  170. SKNode* current = hdr;
  171. for(inti = listLevel -1;i>=0;i--)
  172. {
  173. current = hdr->forward[i];
  174. cout<<"level "<<i<<"................................"<<endl;
  175. while(current->forward[i] != NULL)//key大於下一個數據的值。轉到本層下一個元素
  176. {
  177. cout<<" "<<current->key;
  178. current = current->forward[i];
  179. }
  180. cout<<" "<<current->key<<endl;
  181. }
  182. }
  183. int main()
  184. {
  185. SkipListsk;
  186. constint n = 7;
  187. intnum[n]={30,15,45,37,11,53,17};
  188. cout<<"test insert............."<<endl;
  189. for(inti = 0;i<n;i++)
  190. {
  191. sk.insert(num[i]);
  192. }
  193. sk.printList();
  194. cout<<"test search............"<<endl;
  195. sk.search(17);
  196. cout<<"test delete................."<<endl;
  197. sk.deleteNode(30);
  198. sk.printList();
  199. system("pause");
  200. return 0;
  201. }



 

跳躍鏈表的搜索時間複雜度爲O(logn)

定理:n個元素的跳躍鏈表的每一次搜索的時間複雜度有很高的概率爲O(logn).

高概率:事件E以很高的概率發生意味着對於a>=1,存在一個合適的常數使得事件E發生的概率Pr{E}>=1-O(1/n^a).

其中a是任意選擇的一個數,不同的a影響搜索時間複雜度的常數,即a*O(logn),這個在後面介紹.

我們要證跳躍鏈表的時間複雜度,不能只是證明一次搜索的複雜度爲,是要證明全部的搜索都是O(logn),因爲這是基於概率的算法,如果光一次有效率並沒有多大作用.

我們定義時間Ei爲某一次搜索失敗的概率,那麼假設k次搜素,我們先假定失敗的概率爲O(1/n^a),其中至少有一次失敗的概率爲



可以估算出k次有一次失敗的概率爲1/n^(a-c),那麼我們只要讓a>=c+1或者a取無窮大,就可以證明每一次搜索都具有高概率成功。



跳躍鏈表的層數有高概率爲O(logn)

            類似上面的方法,對於n個元素,如果有一個層數超過O(logn)就算失敗。那麼對於某一個元素超過clogn層,即失敗的概率爲Pr{E}.那麼對於一次搜索失敗的概率爲

令a=c-1,則只要a>=1時,就有高概率的可能使得層數爲O(logn)


跳躍鏈表單次查找複雜度大於O(logn)的概率

每完成一次查找,都肯定要從最頂層移動到最下面一層,這每改變一次層數是由概率選擇時候的p處於扔硬幣中的正面決定的。.既然上面知道層數高概率爲clogn層,那麼扔正面的次數爲clogn-1次.

我們假設扔了clogn個正面,超過10*clogn次是反面,則有


因爲9-log(10e)中的9是線性增長,要遠遠大於log(10e)中的對數增長,因此超過10clogn的概率隨着10的增長變得越來越.所以全部的操作都在10logn以內,我們使用k替代10作爲常數,即查找次數爲klogn,爲O(logn).


發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章