--- comments: true --- # 8.1 ヒープ ヒープ(heap)は、特定の条件を満たす完全二分木であり、主に次の 2 種類に分けられます。 - 最小ヒープ(min heap):任意のノードの値 $\leq$ その子ノードの値。 - 最大ヒープ(max heap):任意のノードの値 $\geq$ その子ノードの値。 { class="animation-figure" }
図 8-1 最小ヒープと最大ヒープ
ヒープは完全二分木の特殊な例であり、次の性質を持ちます。 - 最下層のノードは左から順に埋められ、ほかの層のノードはすべて埋まっています。 - 二分木の根ノードを「ヒープ頂点」、最下層で最も右にあるノードを「ヒープ底」と呼びます。 - 最大ヒープ(最小ヒープ)では、ヒープ頂点の要素(根ノード)の値が最大(最小)です。 ## 8.1.1 ヒープの基本操作 ここで注意したいのは、多くのプログラミング言語が提供しているのは優先度付きキュー(priority queue)であり、これは優先度順に並ぶキューとして定義される抽象データ構造だということです。 実際には、**ヒープは通常、優先度付きキューの実装に用いられ、最大ヒープは要素が大きい順に取り出される優先度付きキューに相当します**。利用の観点では、「優先度付きキュー」と「ヒープ」は等価なデータ構造とみなせます。そのため、本書では両者を特に区別せず、まとめて「ヒープ」と呼びます。 ヒープの基本操作を以下の表に示します。メソッド名はプログラミング言語によって異なります。表 8-1 ヒープの操作効率
図 8-2 ヒープの表現と格納
インデックスの対応式は関数にまとめておくと、後続で使いやすくなります: === "Python" ```python title="my_heap.py" def left(self, i: int) -> int: """左子ノードのインデックスを取得""" return 2 * i + 1 def right(self, i: int) -> int: """右子ノードのインデックスを取得""" return 2 * i + 2 def parent(self, i: int) -> int: """親ノードのインデックスを取得""" return (i - 1) // 2 # 切り捨て除算 ``` === "C++" ```cpp title="my_heap.cpp" /* 左子ノードのインデックスを取得 */ int left(int i) { return 2 * i + 1; } /* 右子ノードのインデックスを取得 */ int right(int i) { return 2 * i + 2; } /* 親ノードのインデックスを取得 */ int parent(int i) { return (i - 1) / 2; // 切り捨て除算 } ``` === "Java" ```java title="my_heap.java" /* 左子ノードのインデックスを取得 */ int left(int i) { return 2 * i + 1; } /* 右子ノードのインデックスを取得 */ int right(int i) { return 2 * i + 2; } /* 親ノードのインデックスを取得 */ int parent(int i) { return (i - 1) / 2; // 切り捨て除算 } ``` === "C#" ```csharp title="my_heap.cs" /* 左子ノードのインデックスを取得 */ int Left(int i) { return 2 * i + 1; } /* 右子ノードのインデックスを取得 */ int Right(int i) { return 2 * i + 2; } /* 親ノードのインデックスを取得 */ int Parent(int i) { return (i - 1) / 2; // 切り捨て除算 } ``` === "Go" ```go title="my_heap.go" /* 左子ノードのインデックスを取得 */ func (h *maxHeap) left(i int) int { return 2*i + 1 } /* 右子ノードのインデックスを取得 */ func (h *maxHeap) right(i int) int { return 2*i + 2 } /* 親ノードのインデックスを取得 */ func (h *maxHeap) parent(i int) int { // 切り捨て除算 return (i - 1) / 2 } ``` === "Swift" ```swift title="my_heap.swift" /* 左子ノードのインデックスを取得 */ func left(i: Int) -> Int { 2 * i + 1 } /* 右子ノードのインデックスを取得 */ func right(i: Int) -> Int { 2 * i + 2 } /* 親ノードのインデックスを取得 */ func parent(i: Int) -> Int { (i - 1) / 2 // 切り捨て除算 } ``` === "JS" ```javascript title="my_heap.js" /* 左子ノードのインデックスを取得 */ #left(i) { return 2 * i + 1; } /* 右子ノードのインデックスを取得 */ #right(i) { return 2 * i + 2; } /* 親ノードのインデックスを取得 */ #parent(i) { return Math.floor((i - 1) / 2); // 切り捨て除算 } ``` === "TS" ```typescript title="my_heap.ts" /* 左子ノードのインデックスを取得 */ left(i: number): number { return 2 * i + 1; } /* 右子ノードのインデックスを取得 */ right(i: number): number { return 2 * i + 2; } /* 親ノードのインデックスを取得 */ parent(i: number): number { return Math.floor((i - 1) / 2); // 切り捨て除算 } ``` === "Dart" ```dart title="my_heap.dart" /* 左子ノードのインデックスを取得 */ int _left(int i) { return 2 * i + 1; } /* 右子ノードのインデックスを取得 */ int _right(int i) { return 2 * i + 2; } /* 親ノードのインデックスを取得 */ int _parent(int i) { return (i - 1) ~/ 2; // 切り捨て除算 } ``` === "Rust" ```rust title="my_heap.rs" /* 左子ノードのインデックスを取得 */ fn left(i: usize) -> usize { 2 * i + 1 } /* 右子ノードのインデックスを取得 */ fn right(i: usize) -> usize { 2 * i + 2 } /* 親ノードのインデックスを取得 */ fn parent(i: usize) -> usize { (i - 1) / 2 // 切り捨て除算 } ``` === "C" ```c title="my_heap.c" /* 左子ノードのインデックスを取得 */ int left(MaxHeap *maxHeap, int i) { return 2 * i + 1; } /* 右子ノードのインデックスを取得 */ int right(MaxHeap *maxHeap, int i) { return 2 * i + 2; } /* 親ノードのインデックスを取得 */ int parent(MaxHeap *maxHeap, int i) { return (i - 1) / 2; // 切り捨て } ``` === "Kotlin" ```kotlin title="my_heap.kt" /* 左子ノードのインデックスを取得 */ fun left(i: Int): Int { return 2 * i + 1 } /* 右子ノードのインデックスを取得 */ fun right(i: Int): Int { return 2 * i + 2 } /* 親ノードのインデックスを取得 */ fun parent(i: Int): Int { return (i - 1) / 2 // 切り捨て除算 } ``` === "Ruby" ```ruby title="my_heap.rb" ### 左子ノードのインデックスを取得 ### def left(i) 2 * i + 1 end ### 右子ノードのインデックスを取得 ### def right(i) 2 * i + 2 end ### 親ノードのインデックスを取得 ### def parent(i) (i - 1) / 2 # 切り捨て除算 end ``` ### 2. ヒープ頂点の要素にアクセス ヒープ頂点の要素は二分木の根ノード、すなわちリストの先頭要素です: === "Python" ```python title="my_heap.py" def peek(self) -> int: """ヒープ先頭要素にアクセス""" return self.max_heap[0] ``` === "C++" ```cpp title="my_heap.cpp" /* ヒープ先頭要素にアクセス */ int peek() { return maxHeap[0]; } ``` === "Java" ```java title="my_heap.java" /* ヒープ先頭要素にアクセス */ int peek() { return maxHeap.get(0); } ``` === "C#" ```csharp title="my_heap.cs" /* ヒープ先頭要素にアクセス */ int Peek() { return maxHeap[0]; } ``` === "Go" ```go title="my_heap.go" /* ヒープ先頭要素にアクセス */ func (h *maxHeap) peek() any { return h.data[0] } ``` === "Swift" ```swift title="my_heap.swift" /* ヒープ先頭要素にアクセス */ func peek() -> Int { maxHeap[0] } ``` === "JS" ```javascript title="my_heap.js" /* ヒープ先頭要素にアクセス */ peek() { return this.#maxHeap[0]; } ``` === "TS" ```typescript title="my_heap.ts" /* ヒープ先頭要素にアクセス */ peek(): number { return this.maxHeap[0]; } ``` === "Dart" ```dart title="my_heap.dart" /* ヒープ先頭要素にアクセス */ int peek() { return _maxHeap[0]; } ``` === "Rust" ```rust title="my_heap.rs" /* ヒープ先頭要素にアクセス */ fn peek(&self) -> Option図 8-3 要素をヒープに追加する手順
ノード総数を $n$ とすると、木の高さは $O(\log n)$ です。したがって、ヒープ化操作のループ回数は高々 $O(\log n)$ であり、**要素をヒープに追加する操作の時間計算量は $O(\log n)$** です。コードは以下のとおりです: === "Python" ```python title="my_heap.py" def push(self, val: int): """要素をヒープに追加""" # ノードを追加 self.max_heap.append(val) # 下から上へヒープ化 self.sift_up(self.size() - 1) def sift_up(self, i: int): """ノード i から始めて、下から上へヒープ化""" while True: # ノード i の親ノードを取得 p = self.parent(i) # 「根ノードを越えた」または「ノードの修復が不要」になったらヒープ化を終了 if p < 0 or self.max_heap[i] <= self.max_heap[p]: break # 2 つのノードを交換 self.swap(i, p) # ループで下から上へヒープ化 i = p ``` === "C++" ```cpp title="my_heap.cpp" /* 要素をヒープに追加 */ void push(int val) { // ノードを追加 maxHeap.push_back(val); // 下から上へヒープ化 siftUp(size() - 1); } /* ノード i から始めて、下から上へヒープ化 */ void siftUp(int i) { while (true) { // ノード i の親ノードを取得 int p = parent(i); // 「根ノードを越えた」または「ノードの修復が不要」になったらヒープ化を終了 if (p < 0 || maxHeap[i] <= maxHeap[p]) break; // 2 つのノードを交換 swap(maxHeap[i], maxHeap[p]); // ループで下から上へヒープ化 i = p; } } ``` === "Java" ```java title="my_heap.java" /* 要素をヒープに追加 */ void push(int val) { // ノードを追加 maxHeap.add(val); // 下から上へヒープ化 siftUp(size() - 1); } /* ノード i から始めて、下から上へヒープ化 */ void siftUp(int i) { while (true) { // ノード i の親ノードを取得 int p = parent(i); // 「根ノードを越えた」または「ノードの修復が不要」になったらヒープ化を終了 if (p < 0 || maxHeap.get(i) <= maxHeap.get(p)) break; // 2 つのノードを交換 swap(i, p); // ループで下から上へヒープ化 i = p; } } ``` === "C#" ```csharp title="my_heap.cs" /* 要素をヒープに追加 */ void Push(int val) { // ノードを追加 maxHeap.Add(val); // 下から上へヒープ化 SiftUp(Size() - 1); } /* ノード i から始めて、下から上へヒープ化 */ void SiftUp(int i) { while (true) { // ノード i の親ノードを取得 int p = Parent(i); // 「根ノードを越えた」または「ノードの修復が不要」な場合は、ヒープ化を終了する if (p < 0 || maxHeap[i] <= maxHeap[p]) break; // 2 つのノードを交換 Swap(i, p); // ループで下から上へヒープ化 i = p; } } ``` === "Go" ```go title="my_heap.go" /* 要素をヒープに追加 */ func (h *maxHeap) push(val any) { // ノードを追加 h.data = append(h.data, val) // 下から上へヒープ化 h.siftUp(len(h.data) - 1) } /* ノード i から始めて、下から上へヒープ化 */ func (h *maxHeap) siftUp(i int) { for true { // ノード i の親ノードを取得 p := h.parent(i) // 「根ノードを越えた」または「ノードの修復が不要」になったらヒープ化を終了 if p < 0 || h.data[i].(int) <= h.data[p].(int) { break } // 2 つのノードを交換 h.swap(i, p) // ループで下から上へヒープ化 i = p } } ``` === "Swift" ```swift title="my_heap.swift" /* 要素をヒープに追加 */ func push(val: Int) { // ノードを追加 maxHeap.append(val) // 下から上へヒープ化 siftUp(i: size() - 1) } /* ノード i から始めて、下から上へヒープ化 */ func siftUp(i: Int) { var i = i while true { // ノード i の親ノードを取得 let p = parent(i: i) // 「根ノードを越えた」または「ノードの修復が不要」になったらヒープ化を終了 if p < 0 || maxHeap[i] <= maxHeap[p] { break } // 2 つのノードを交換 swap(i: i, j: p) // ループで下から上へヒープ化 i = p } } ``` === "JS" ```javascript title="my_heap.js" /* 要素をヒープに追加 */ push(val) { // ノードを追加 this.#maxHeap.push(val); // 下から上へヒープ化 this.#siftUp(this.size() - 1); } /* ノード i から始めて、下から上へヒープ化 */ #siftUp(i) { while (true) { // ノード i の親ノードを取得 const p = this.#parent(i); // 「根ノードを越えた」または「ノードの修復が不要」になったらヒープ化を終了 if (p < 0 || this.#maxHeap[i] <= this.#maxHeap[p]) break; // 2 つのノードを交換 this.#swap(i, p); // ループで下から上へヒープ化 i = p; } } ``` === "TS" ```typescript title="my_heap.ts" /* 要素をヒープに追加 */ push(val: number): void { // ノードを追加 this.maxHeap.push(val); // 下から上へヒープ化 this.siftUp(this.size() - 1); } /* ノード i から始めて、下から上へヒープ化 */ siftUp(i: number): void { while (true) { // ノード i の親ノードを取得 const p = this.parent(i); // 「根ノードを越えた」または「ノードの修復が不要」になったらヒープ化を終了 if (p < 0 || this.maxHeap[i] <= this.maxHeap[p]) break; // 2 つのノードを交換 this.swap(i, p); // ループで下から上へヒープ化 i = p; } } ``` === "Dart" ```dart title="my_heap.dart" /* 要素をヒープに追加 */ void push(int val) { // ノードを追加 _maxHeap.add(val); // 下から上へヒープ化 siftUp(size() - 1); } /* ノード i から始めて、下から上へヒープ化 */ void siftUp(int i) { while (true) { // ノード i の親ノードを取得 int p = _parent(i); // 「根ノードを越えた」または「ノードの修復が不要」になったらヒープ化を終了 if (p < 0 || _maxHeap[i] <= _maxHeap[p]) { break; } // 2 つのノードを交換 _swap(i, p); // ループで下から上へヒープ化 i = p; } } ``` === "Rust" ```rust title="my_heap.rs" /* 要素をヒープに追加 */ fn push(&mut self, val: i32) { // ノードを追加 self.max_heap.push(val); // 下から上へヒープ化 self.sift_up(self.size() - 1); } /* ノード i から始めて、下から上へヒープ化 */ fn sift_up(&mut self, mut i: usize) { loop { // ノード i はすでにヒープの先頭ノードなので、ヒープ化を終了する if i == 0 { break; } // ノード i の親ノードを取得 let p = Self::parent(i); // 「ノードの修復が不要」になったら、ヒープ化を終了 if self.max_heap[i] <= self.max_heap[p] { break; } // 2 つのノードを交換 self.swap(i, p); // ループで下から上へヒープ化 i = p; } } ``` === "C" ```c title="my_heap.c" /* 要素をヒープに追加 */ void push(MaxHeap *maxHeap, int val) { // 通常は、これほど多くのノードを追加すべきではない if (maxHeap->size == MAX_SIZE) { printf("heap is full!"); return; } // ノードを追加 maxHeap->data[maxHeap->size] = val; maxHeap->size++; // 下から上へヒープ化 siftUp(maxHeap, maxHeap->size - 1); } /* ノード i から始めて、下から上へヒープ化 */ void siftUp(MaxHeap *maxHeap, int i) { while (true) { // ノード i の親ノードを取得 int p = parent(maxHeap, i); // 「根ノードを越えた」または「ノードの修復が不要」になったらヒープ化を終了 if (p < 0 || maxHeap->data[i] <= maxHeap->data[p]) { break; } // 2 つのノードを交換 swap(maxHeap, i, p); // ループで下から上へヒープ化 i = p; } } ``` === "Kotlin" ```kotlin title="my_heap.kt" /* 要素をヒープに追加 */ fun push(_val: Int) { // ノードを追加 maxHeap.add(_val) // 下から上へヒープ化 siftUp(size() - 1) } /* ノード i から始めて、下から上へヒープ化 */ fun siftUp(it: Int) { // Kotlin の関数引数は不変のため、一時変数を作成する var i = it while (true) { // ノード i の親ノードを取得 val p = parent(i) // 「根ノードを越えた」または「ノードの修復が不要」になったらヒープ化を終了 if (p < 0 || maxHeap[i] <= maxHeap[p]) break // 2 つのノードを交換 swap(i, p) // ループで下から上へヒープ化 i = p } } ``` === "Ruby" ```ruby title="my_heap.rb" ### 要素をヒープに挿入 ### def push(val) # ノードを追加 @max_heap << val # 下から上へヒープ化 sift_up(size - 1) end ### ノード i から下から上へヒープ化 ### def sift_up(i) loop do # ノード i の親ノードを取得 p = parent(i) # 「根ノードを越えた」または「ノードの修復が不要」になったらヒープ化を終了 break if p < 0 || @max_heap[i] <= @max_heap[p] # 2 つのノードを交換 swap(i, p) # ループで下から上へヒープ化 i = p end end ``` ??? pythontutor "コードの可視化" ### 4. ヒープ頂点の要素を取り出す ヒープ頂点の要素は二分木の根ノード、すなわちリストの先頭要素です。もし先頭要素をそのまま削除すると、二分木内のすべてのノードのインデックスが変化してしまい、その後のヒープ化による修復が困難になります。要素インデックスの変動をできるだけ小さくするため、次の手順を取ります。 1. ヒープ頂点の要素とヒープ底の要素を交換する(根ノードと最も右の葉ノードを交換する)。 2. 交換後、ヒープ底をリストから削除する(すでに交換済みであるため、実際に削除されるのは元のヒープ頂点の要素であることに注意)。 3. 根ノードから開始し、**上から下へヒープ化**を行う。 次の図のように、**「上から下へのヒープ化」の方向は「下から上へのヒープ化」と逆**です。根ノードの値を 2 つの子ノードと比較し、最大の子ノードと根ノードを交換します。その後、この操作を繰り返し、葉ノードを越えるか交換不要のノードに達した時点で終了します。 === "<1>" { class="animation-figure" } === "<2>" { class="animation-figure" } === "<3>" { class="animation-figure" } === "<4>" { class="animation-figure" } === "<5>" { class="animation-figure" } === "<6>" { class="animation-figure" } === "<7>" { class="animation-figure" } === "<8>" { class="animation-figure" } === "<9>" { class="animation-figure" } === "<10>" { class="animation-figure" }図 8-4 ヒープ頂点の要素を取り出す手順
要素をヒープに追加する操作と同様に、ヒープ頂点の要素を取り出す操作の時間計算量も $O(\log n)$ です。コードは以下のとおりです: === "Python" ```python title="my_heap.py" def pop(self) -> int: """要素をヒープから取り出す""" # 空判定の処理 if self.is_empty(): raise IndexError("ヒープが空です") # 根ノードと最も右の葉ノードを交換(先頭要素と末尾要素を交換) self.swap(0, self.size() - 1) # ノードを削除 val = self.max_heap.pop() # 上から下へヒープ化 self.sift_down(0) # ヒープ先頭要素を返す return val def sift_down(self, i: int): """ノード i から始めて、上から下へヒープ化""" while True: # ノード i, l, r のうち値が最大のノードを ma とする l, r, ma = self.left(i), self.right(i), i if l < self.size() and self.max_heap[l] > self.max_heap[ma]: ma = l if r < self.size() and self.max_heap[r] > self.max_heap[ma]: ma = r # ノード i が最大、またはインデックス l, r が範囲外なら、ヒープ化は不要なので抜ける if ma == i: break # 2 つのノードを交換 self.swap(i, ma) # ループで上から下へヒープ化 i = ma ``` === "C++" ```cpp title="my_heap.cpp" /* 要素をヒープから取り出す */ void pop() { // 空判定の処理 if (isEmpty()) { throw out_of_range("ヒープが空です"); } // 根ノードと最も右の葉ノードを交換(先頭要素と末尾要素を交換) swap(maxHeap[0], maxHeap[size() - 1]); // ノードを削除 maxHeap.pop_back(); // 上から下へヒープ化 siftDown(0); } /* ノード i から始めて、上から下へヒープ化 */ void siftDown(int i) { while (true) { // ノード i, l, r のうち値が最大のノードを ma とする int l = left(i), r = right(i), ma = i; if (l < size() && maxHeap[l] > maxHeap[ma]) ma = l; if (r < size() && maxHeap[r] > maxHeap[ma]) ma = r; // ノード i が最大、またはインデックス l, r が範囲外なら、ヒープ化は不要なので抜ける if (ma == i) break; swap(maxHeap[i], maxHeap[ma]); // ループで上から下へヒープ化 i = ma; } } ``` === "Java" ```java title="my_heap.java" /* 要素をヒープから取り出す */ int pop() { // 空判定の処理 if (isEmpty()) throw new IndexOutOfBoundsException(); // 根ノードと最も右の葉ノードを交換(先頭要素と末尾要素を交換) swap(0, size() - 1); // ノードを削除 int val = maxHeap.remove(size() - 1); // 上から下へヒープ化 siftDown(0); // ヒープ先頭要素を返す return val; } /* ノード i から始めて、上から下へヒープ化 */ void siftDown(int i) { while (true) { // ノード i, l, r のうち値が最大のノードを ma とする int l = left(i), r = right(i), ma = i; if (l < size() && maxHeap.get(l) > maxHeap.get(ma)) ma = l; if (r < size() && maxHeap.get(r) > maxHeap.get(ma)) ma = r; // ノード i が最大、またはインデックス l, r が範囲外なら、ヒープ化は不要なので抜ける if (ma == i) break; // 2 つのノードを交換 swap(i, ma); // ループで上から下へヒープ化 i = ma; } } ``` === "C#" ```csharp title="my_heap.cs" /* 要素をヒープから取り出す */ int Pop() { // 空判定の処理 if (IsEmpty()) throw new IndexOutOfRangeException(); // 根ノードと最も右の葉ノードを交換(先頭要素と末尾要素を交換) Swap(0, Size() - 1); // ノードを削除 int val = maxHeap.Last(); maxHeap.RemoveAt(Size() - 1); // 上から下へヒープ化 SiftDown(0); // ヒープ先頭要素を返す return val; } /* ノード i から始めて、上から下へヒープ化 */ void SiftDown(int i) { while (true) { // ノード i, l, r のうち値が最大のノードを ma とする int l = Left(i), r = Right(i), ma = i; if (l < Size() && maxHeap[l] > maxHeap[ma]) ma = l; if (r < Size() && maxHeap[r] > maxHeap[ma]) ma = r; // 「ノード i が最大」または「葉ノードを越えた」場合は、ヒープ化を終了する if (ma == i) break; // 2 つのノードを交換 Swap(i, ma); // ループで上から下へヒープ化 i = ma; } } ``` === "Go" ```go title="my_heap.go" /* 要素をヒープから取り出す */ func (h *maxHeap) pop() any { // 空判定の処理 if h.isEmpty() { fmt.Println("error") return nil } // 根ノードと最も右の葉ノードを交換(先頭要素と末尾要素を交換) h.swap(0, h.size()-1) // ノードを削除 val := h.data[len(h.data)-1] h.data = h.data[:len(h.data)-1] // 上から下へヒープ化 h.siftDown(0) // ヒープ先頭要素を返す return val } /* ノード i から始めて、上から下へヒープ化 */ func (h *maxHeap) siftDown(i int) { for true { // ノード i, l, r のうち値が最大のノードを max とする l, r, max := h.left(i), h.right(i), i if l < h.size() && h.data[l].(int) > h.data[max].(int) { max = l } if r < h.size() && h.data[r].(int) > h.data[max].(int) { max = r } // ノード i が最大、またはインデックス l, r が範囲外なら、ヒープ化は不要なので抜ける if max == i { break } // 2 つのノードを交換 h.swap(i, max) // ループで上から下へヒープ化 i = max } } ``` === "Swift" ```swift title="my_heap.swift" /* 要素をヒープから取り出す */ func pop() -> Int { // 空判定の処理 if isEmpty() { fatalError("ヒープが空です") } // 根ノードと最も右の葉ノードを交換(先頭要素と末尾要素を交換) swap(i: 0, j: size() - 1) // ノードを削除 let val = maxHeap.remove(at: size() - 1) // 上から下へヒープ化 siftDown(i: 0) // ヒープ先頭要素を返す return val } /* ノード i から始めて、上から下へヒープ化 */ func siftDown(i: Int) { var i = i while true { // ノード i, l, r のうち値が最大のノードを ma とする let l = left(i: i) let r = right(i: i) var ma = i if l < size(), maxHeap[l] > maxHeap[ma] { ma = l } if r < size(), maxHeap[r] > maxHeap[ma] { ma = r } // ノード i が最大、またはインデックス l, r が範囲外なら、ヒープ化は不要なので抜ける if ma == i { break } // 2 つのノードを交換 swap(i: i, j: ma) // ループで上から下へヒープ化 i = ma } } ``` === "JS" ```javascript title="my_heap.js" /* 要素をヒープから取り出す */ pop() { // 空判定の処理 if (this.isEmpty()) throw new Error('ヒープが空です'); // 根ノードと最も右の葉ノードを交換(先頭要素と末尾要素を交換) this.#swap(0, this.size() - 1); // ノードを削除 const val = this.#maxHeap.pop(); // 上から下へヒープ化 this.#siftDown(0); // ヒープ先頭要素を返す return val; } /* ノード i から始めて、上から下へヒープ化 */ #siftDown(i) { while (true) { // ノード i, l, r のうち値が最大のノードを ma とする const l = this.#left(i), r = this.#right(i); let ma = i; if (l < this.size() && this.#maxHeap[l] > this.#maxHeap[ma]) ma = l; if (r < this.size() && this.#maxHeap[r] > this.#maxHeap[ma]) ma = r; // ノード i が最大、またはインデックス l, r が範囲外なら、ヒープ化は不要なので抜ける if (ma === i) break; // 2 つのノードを交換 this.#swap(i, ma); // ループで上から下へヒープ化 i = ma; } } ``` === "TS" ```typescript title="my_heap.ts" /* 要素をヒープから取り出す */ pop(): number { // 空判定の処理 if (this.isEmpty()) throw new RangeError('Heap is empty.'); // 根ノードと最も右の葉ノードを交換(先頭要素と末尾要素を交換) this.swap(0, this.size() - 1); // ノードを削除 const val = this.maxHeap.pop(); // 上から下へヒープ化 this.siftDown(0); // ヒープ先頭要素を返す return val; } /* ノード i から始めて、上から下へヒープ化 */ siftDown(i: number): void { while (true) { // ノード i, l, r のうち値が最大のノードを ma とする const l = this.left(i), r = this.right(i); let ma = i; if (l < this.size() && this.maxHeap[l] > this.maxHeap[ma]) ma = l; if (r < this.size() && this.maxHeap[r] > this.maxHeap[ma]) ma = r; // ノード i が最大、またはインデックス l, r が範囲外なら、ヒープ化は不要なので抜ける if (ma === i) break; // 2 つのノードを交換 this.swap(i, ma); // ループで上から下へヒープ化 i = ma; } } ``` === "Dart" ```dart title="my_heap.dart" /* 要素をヒープから取り出す */ int pop() { // 空判定の処理 if (isEmpty()) throw Exception('ヒープが空です'); // 根ノードと最も右の葉ノードを交換(先頭要素と末尾要素を交換) _swap(0, size() - 1); // ノードを削除 int val = _maxHeap.removeLast(); // 上から下へヒープ化 siftDown(0); // ヒープ先頭要素を返す return val; } /* ノード i から始めて、上から下へヒープ化 */ void siftDown(int i) { while (true) { // ノード i, l, r のうち値が最大のノードを ma とする int l = _left(i); int r = _right(i); int ma = i; if (l < size() && _maxHeap[l] > _maxHeap[ma]) ma = l; if (r < size() && _maxHeap[r] > _maxHeap[ma]) ma = r; // ノード i が最大、またはインデックス l, r が範囲外なら、ヒープ化は不要なので抜ける if (ma == i) break; // 2 つのノードを交換 _swap(i, ma); // ループで上から下へヒープ化 i = ma; } } ``` === "Rust" ```rust title="my_heap.rs" /* 要素をヒープから取り出す */ fn pop(&mut self) -> i32 { // 空判定の処理 if self.is_empty() { panic!("index out of bounds"); } // 根ノードと最も右の葉ノードを交換(先頭要素と末尾要素を交換) self.swap(0, self.size() - 1); // ノードを削除 let val = self.max_heap.pop().unwrap(); // 上から下へヒープ化 self.sift_down(0); // ヒープ先頭要素を返す val } /* ノード i から始めて、上から下へヒープ化 */ fn sift_down(&mut self, mut i: usize) { loop { // ノード i, l, r のうち値が最大のノードを ma とする let (l, r, mut ma) = (Self::left(i), Self::right(i), i); if l < self.size() && self.max_heap[l] > self.max_heap[ma] { ma = l; } if r < self.size() && self.max_heap[r] > self.max_heap[ma] { ma = r; } // ノード i が最大、またはインデックス l, r が範囲外なら、ヒープ化は不要なので抜ける if ma == i { break; } // 2 つのノードを交換 self.swap(i, ma); // ループで上から下へヒープ化 i = ma; } } ``` === "C" ```c title="my_heap.c" /* 要素をヒープから取り出す */ int pop(MaxHeap *maxHeap) { // 空判定の処理 if (isEmpty(maxHeap)) { printf("heap is empty!"); return INT_MAX; } // 根ノードと最も右の葉ノードを交換(先頭要素と末尾要素を交換) swap(maxHeap, 0, size(maxHeap) - 1); // ノードを削除 int val = maxHeap->data[maxHeap->size - 1]; maxHeap->size--; // 上から下へヒープ化 siftDown(maxHeap, 0); // ヒープ先頭要素を返す return val; } /* ノード i から始めて、上から下へヒープ化 */ void siftDown(MaxHeap *maxHeap, int i) { while (true) { // ノード i, l, r のうち値が最大のノードを max とする int l = left(maxHeap, i); int r = right(maxHeap, i); int max = i; if (l < size(maxHeap) && maxHeap->data[l] > maxHeap->data[max]) { max = l; } if (r < size(maxHeap) && maxHeap->data[r] > maxHeap->data[max]) { max = r; } // ノード i が最大、またはインデックス l, r が範囲外なら、ヒープ化は不要なので抜ける if (max == i) { break; } // 2 つのノードを交換 swap(maxHeap, i, max); // ループで上から下へヒープ化 i = max; } } ``` === "Kotlin" ```kotlin title="my_heap.kt" /* 要素をヒープから取り出す */ fun pop(): Int { // 空判定の処理 if (isEmpty()) throw IndexOutOfBoundsException() // 根ノードと最も右の葉ノードを交換(先頭要素と末尾要素を交換) swap(0, size() - 1) // ノードを削除 val _val = maxHeap.removeAt(size() - 1) // 上から下へヒープ化 siftDown(0) // ヒープ先頭要素を返す return _val } /* ノード i から始めて、上から下へヒープ化 */ fun siftDown(it: Int) { // Kotlin の関数引数は不変のため、一時変数を作成する var i = it while (true) { // ノード i, l, r のうち値が最大のノードを ma とする val l = left(i) val r = right(i) var ma = i if (l < size() && maxHeap[l] > maxHeap[ma]) ma = l if (r < size() && maxHeap[r] > maxHeap[ma]) ma = r // ノード i が最大、またはインデックス l, r が範囲外なら、ヒープ化は不要なので抜ける if (ma == i) break // 2 つのノードを交換 swap(i, ma) // ループで上から下へヒープ化 i = ma } } ``` === "Ruby" ```ruby title="my_heap.rb" ### 要素をヒープから取り出す ### def pop # 空判定の処理 raise IndexError, "ヒープが空です" if is_empty? # 根ノードと最も右の葉ノードを交換(先頭要素と末尾要素を交換) swap(0, size - 1) # ノードを削除 val = @max_heap.pop # 上から下へヒープ化 sift_down(0) # ヒープ先頭要素を返す val end ### ノード i から上から下へヒープ化 ### def sift_down(i) loop do # ノード i, l, r のうち値が最大のノードを ma とする l, r, ma = left(i), right(i), i ma = l if l < size && @max_heap[l] > @max_heap[ma] ma = r if r < size && @max_heap[r] > @max_heap[ma] # ノード i が最大、またはインデックス l, r が範囲外なら、ヒープ化は不要なので抜ける break if ma == i # 2 つのノードを交換 swap(i, ma) # ループで上から下へヒープ化 i = ma end end ``` ??? pythontutor "コードの可視化" ## 8.1.3 ヒープの代表的な応用 - **優先度付きキュー**:ヒープは、優先度付きキューを実装するための代表的なデータ構造です。キューへの追加と取り出しの時間計算量はいずれも $O(\log n)$ で、ヒープ構築は $O(n)$ であり、これらの操作はいずれも非常に効率的です。 - **ヒープソート**:与えられたデータ群からヒープを構築し、要素の取り出しを繰り返すことで整列済みデータを得られます。ただし、通常はより洗練された方法でヒープソートを実装します。詳しくは「ヒープソート」の章を参照してください。 - **最大の $k$ 個の要素を取得**:これは古典的なアルゴリズム問題であると同時に、典型的な応用でもあります。たとえば、人気上位 10 件のニュースをホットトピックとして選んだり、売上上位 10 件の商品を選んだりする場面です。