コンテンツにスキップ

基本的なデータ構造:配列・連結リスト・スタック・キュー

前提:計算量と O 記法:アルゴリズムの速さを入力サイズの関数で測る

生 Markdown
  • データ構造の優劣は「何ができるか」ではなく「どの操作がどれだけ速いか」で決まります。同じ列を表す配列と連結リストは、速い操作がちょうど正反対です。
  • 配列は添字によるアクセスが Θ(1)\Theta(1)、途中への挿入・削除が Θ(n)\Theta(n)。連結リストは位置を指すノードが手元にあればその直後への挿入・削除が Θ(1)\Theta(1)kk 番目へのアクセスが Θ(k)\Theta(k) です(命題 3.2命題 4.2)。
  • 容量を 2 倍ずつ増やす動的配列は、1 回の push の最悪計算量が Θ(n)\Theta(n) であっても、nn 回の合計は 3n3n 未満に収まります。償却計算量は O(1)O(1) です(定理 3.3)。
  • 「償却 O(1)O(1)」は「平均的に速い」という意味ではありません。入力の確率分布を一切仮定せず、どんな操作列に対しても合計コストが保証されます。
  • スタック(LIFO)とキュー(FIFO)は抽象データ型であり、実装に触れずに公理で規定できます。同じ公理を配列でも連結リストでも満たせます。
  • キューはリングバッファで最悪 Θ(1)\Theta(1)命題 7.2)、2 本のスタックで償却 O(1)O(1)定理 8.2)に実現できます。後者はポテンシャル法で解析します。

プログラムが扱うデータの多くは「列」です。テキストは文字の列、画像は画素の列、口座の履歴は取引の列です。列を計算機の中でどう並べるかには選択の余地があり、その選択によって処理時間が桁単位で変わります。

素朴な疑問から始めましょう。なぜ配列と連結リストという二つの表現が、どちらも半世紀以上使われ続けているのでしょうか。答えは、両者の得意な操作がちょうど正反対だからです。配列は「ii 番目を見せてほしい」に強く、「途中に割り込ませてほしい」に弱い。連結リストはその逆です。片方がすべての操作で勝っていれば、もう片方は歴史から消えていたはずです。

歴史的にも、この二つは別々の要求から生まれました。配列は最初期の高級言語 Fortran(1957 年)から言語機能として存在します。連結リストはやや遅れて、1956 年前後に Allen Newell、Cliff Shaw、Herbert Simon が定理証明プログラム Logic Theorist のために設計した情報処理言語 IPL で導入され、1958 年の LISP に受け継がれました。彼らが必要としたのは「途中に要素を割り込ませても、後ろ全体をずらさなくてよい表現」でした。

この記事では、前章 計算量と O 記法OOΘ\Theta 記法(定義 3.1[計算量と O 記法])を道具として使い、まず「何ができるか」(抽象データ型)と「どう実現するか」(データ構造)を分けます。そのうえで配列・連結リスト・スタック・キューを取り上げ、各操作の計算量を証明付きで求めます。

2. 準備:計算モデルと抽象データ型

Section titled “2. 準備:計算モデルと抽象データ型”

計算量を「証明する」には、何を 1 単位時間とみなすかを先に決めなければなりません。この記事ではワード RAM モデルを使います。これは前章の 一様コスト RAM モデル(定義 2.1)[計算量と O 記法] を、記憶領域の番地づけまで具体化したものです。

記憶領域は番地 0,1,2,0, 1, 2, \ldots で添字づけられた語の列 M[0],M[1],M[0], M[1], \ldots であり、1 語は整数か番地を 1 つ保持できます。番地 ii を指定した M[i]M[i] の読み出し・書き込み、および語同士の加減乗除と比較は、いずれも 1 単位時間で実行できるものとします。この「番地を指定した読み書きが定数時間」という仮定が、後で見る「配列の添字アクセスが Θ(1)\Theta(1)」の根拠です。

注意 2.1

ワード RAM は理想化です。実際の計算機には階層的なキャッシュがあり、1 語のアクセス時間は一定ではありません。それでもこのモデルを使うのは、アルゴリズムの本質的な差(Θ(n)\Theta(n)Θ(logn)\Theta(\log n)Θ(1)\Theta(1) か)がキャッシュの有無に影響されないからです。同じ Θ(n)\Theta(n) どうしの実測差はモデルの外にあり、注意 5.1 で扱います。

次に、データ構造を論じるための基本的な区別を導入します。

定義 2.2抽象データ型

抽象データ型(abstract data type, ADT)とは、次の 3 つの組であって、値の内部表現を一切含まないものをいう。

  1. 値の集合。
  2. 操作の名前と、その引数および返り値の型。
  3. 操作が満たすべき性質。事前条件と事後条件で書いてもよいし、操作どうしが満たす等式(公理)で書いてもよい。

抽象データ型を、具体的な記憶領域の配置と手続きによって実現したものをデータ構造という。

この区別が効くのは、同じ抽象データ型に複数のデータ構造がありうるからです。どれを選ぶかは「正しさ」では決まらず、どの操作を何回使うかと、その操作の計算量で決まります。

定義 3.1配列

長さ nn の配列とは、同じ大きさ ss(語数)の nn 個の区画を、先頭番地 bb から連続した番地に並べたものをいう。ii 番目の要素(0in10 \le i \le n-1)を A[i]A[i] と書き、その区画の先頭番地は

addr(A[i])=b+si\mathrm{addr}(A[i]) = b + s \cdot i

で与えられる。

定義に現れるこの 1 本の式が、配列の強さと弱さの両方を決めています。番地が ii の一次式で書けるので、どの要素にも計算だけで到達できます。一方、要素を 1 つ割り込ませると、それ以降の要素の ii がすべて 1 ずつずれ、番地もすべてずれます。順に確認します。

命題 3.2配列の基本操作の計算量

長さ nn の配列 AA について、ワード RAM モデルの下で次が成り立つ。

  1. 添字 ii0in10 \le i \le n-1)を指定した要素の読み出しと書き込みは Θ(1)\Theta(1) 時間で行える。
  2. 与えられた値 xxAA に現れるかどうかを判定する探索は、要素の並びに順序の仮定がない場合、最悪 Θ(n)\Theta(n) 時間を要する。先頭から順に比較する線形探索がこの計算量を達成する。
  3. 位置 ii0in0 \le i \le n)への挿入は既存要素の移動をちょうど nin - i 回、位置 ii0in10 \le i \le n-1)の要素の削除は n1in - 1 - i 回必要とし、いずれも Θ(ni)\Theta(n-i) 時間で行える。とくに最悪(i=0i = 0)は Θ(n)\Theta(n) であり、挿入位置 ii0,1,,n0, 1, \ldots, n 上一様分布に従うときの平均移動回数は n/2n/2 である。
証明(命題 3.2)

1. 定義 3.1 より addr(A[i])=b+si\mathrm{addr}(A[i]) = b + s \cdot i です。bbss は配列ごとに固定なので、この計算は乗算 1 回と加算 1 回。ワード RAM では算術演算と番地指定の読み書きがそれぞれ 1 単位時間なので、合計は定数時間、すなわち Θ(1)\Theta(1) です。nn にも ii にも依存しない点が重要です。

2. 上界は線形探索が与えます。i=0,1,,n1i = 0, 1, \ldots, n-1 の順に A[i]A[i] を読んで xx と比較し、一致すれば真を返し、最後まで一致しなければ偽を返します。各反復は 1 の意味で定数時間なので、全体で O(n)O(n) です。

下界を敵対者論法で示します。正しいアルゴリズム B\mathcal{B} が、xx を含まない入力 AA に対し、ある位置 jj を一度も読まずに「含まない」と答えたとします。AA の位置 jj だけを xx に置き換えた AA' を与えると、B\mathcal{B}jj を読まないので同じ計算をたどり、やはり「含まない」と答えます。しかし AA'xx を含むので誤りです。よって nn 箇所すべてを読む必要があり Ω(n)\Omega(n)、上界と合わせて Θ(n)\Theta(n) です。

3. 挿入を示します。挿入後の配列 AA'A[j]=A[j]A'[j] = A[j]j<ij < i)、A[i]=xA'[i] = xA[j+1]=A[j]A'[j+1] = A[j]ijn1i \le j \le n-1)を満たします。jij \ge i の要素は例外なく番地が ss だけ後ろにずれるので、少なくとも nin - i 回の書き込みが必要です。逆に、j=n1,n2,,ij = n-1, n-2, \ldots, i の順に A[j+1]A[j]A[j+1] \leftarrow A[j] と後ろから詰めれば、上書きによる損失なく nin-i 回の書き込みで済み、最後に A[i]xA[i] \leftarrow x を 1 回行えば完了します。1 回の移動は 1 の意味で定数時間なので、合計 Θ(ni)\Theta(n-i) です。削除も同様に、j=i+1,,n1j = i+1, \ldots, n-1 の順に A[j1]A[j]A[j-1] \leftarrow A[j] と前へ詰めることで n1in-1-i 回の移動、すなわち Θ(ni)\Theta(n-i) です。

最悪は i=0i = 0 で移動回数 nn、よって Θ(n)\Theta(n) です。平均は

1n+1i=0n(ni)=1n+1k=0nk=1n+1n(n+1)2=n2\frac{1}{n+1}\sum_{i=0}^{n}(n-i) = \frac{1}{n+1}\sum_{k=0}^{n}k = \frac{1}{n+1}\cdot\frac{n(n+1)}{2} = \frac{n}{2}

となり、これも Θ(n)\Theta(n) です。平均を取っても Θ(n)\Theta(n) から逃げられません。

配列の弱点はもう一つあります。長さを最初に決めなければならないことです。要素数が実行時に増えていく場合、配列そのままでは扱えません。この問題を解くのが動的配列(可変長配列)です。素朴に「1 つ増えるたびに長さ +1+1 の配列を確保して全部コピーする」と、nn 回の追加で k=1nk=Θ(n2)\sum_{k=1}^{n} k = \Theta(n^2) 回のコピーが発生します。ところが、増やし方を「+1+1」から「22 倍」に変えるだけで、合計が線形に落ちます。

定理 3.3動的配列の償却計算量

次の規則で末尾追加 push\mathrm{push} を実現するデータ構造を考える。容量 cc の配列 AA と現在の要素数 nn0nc0 \le n \le c)を保持し、初期状態を n=0n = 0c=1c = 1 とする。push(x)\mathrm{push}(x) は次のように動作する。

  • n<cn < c ならば A[n]xA[n] \leftarrow x とし、nn+1n \leftarrow n + 1 とする。
  • n=cn = c ならば、まず長さ 2c2c の配列を新たに確保して既存の cc 個の要素をそこへ移し、c2cc \leftarrow 2c としてから上の操作を行う。

このとき、空の状態から push\mathrm{push}nn 回(n1n \ge 1)行うときの要素の書き込み回数の合計は 3n3n 未満である。したがって 1 回の push\mathrm{push} あたりの償却計算量は O(1)O(1) である。ただし、個々の push\mathrm{push} の最悪計算量は Θ(n)\Theta(n) である。

証明(定理 3.3)

集約法(総コストを直接見積もり、操作回数で割る方法)で示します。

書き込みを 2 種類に分けます。A[n]xA[n] \leftarrow x という末尾への書き込みは、どの push\mathrm{push} でもちょうど 1 回起きるので合計 nn 回です。

残りは容量拡張に伴うコピーです。拡張が起きるのは push\mathrm{push} の直前に n=cn = c のときで、容量は 1,2,4,1, 2, 4, \ldots と 2 の冪をたどるので、これは直前の要素数が 2k2^kk=0,1,2,k = 0, 1, 2, \ldots)のときに限られ、コピー回数はちょうど 2k2^k です。nn 回の push\mathrm{push} の間に要素数が 2k2^k に達するのは 2k<n2^k < n を満たす kk に対してだけなので、KK をそのような最大の kk とすると、コピーの総数は

k=0K2k=2K+11<22K2n\sum_{k=0}^{K} 2^{k} = 2^{K+1} - 1 < 2 \cdot 2^{K} \le 2n

です。最後の不等号で 2K<n2^K < n を使いました。

以上より総書き込み回数は n+2n=3nn + 2n = 3n 未満です。nn 回の操作で総コストが 3n3n 未満なので、1 回あたりの償却コストは 33 未満、すなわち O(1)O(1) です。

最悪計算量については、要素数が n=c=2kn = c = 2^k の状態での push\mathrm{push}2k=n2^k = n 回のコピーを行うので、命題 3.2 の 1 より Θ(n)\Theta(n) 時間かかります。

例 3.410 回の push を数え上げる

空の状態から 10 回 push\mathrm{push} したときの書き込み回数を数えます。拡張が起きるのは直前の要素数が 1,2,4,81, 2, 4, 8 のとき、すなわち 2・3・5・9 回目の push\mathrm{push} です。

push 番号直前の (n,c)(n, c)拡張後の容量コピー回数この回の書き込み合計
1(0,1)(0, 1)01
2(1,1)(1, 1)212
3(2,2)(2, 2)423
5(4,4)(4, 4)845
9(8,8)(8, 8)1689
4, 6, 7, 8, 100各 1

総書き込み回数は 10+(1+2+4+8)=2510 + (1 + 2 + 4 + 8) = 25 で、定理の保証する上界 3×10=303 \times 10 = 30 を下回っています。9 回目の push\mathrm{push} だけで 9 回の書き込みが必要ですが、平らにならせば 1 回あたり 2.52.5 回です。

注意 3.5償却と平均は別のもの

償却計算量は確率をまったく含みません。定理 3.3 が主張しているのは「どんな push\mathrm{push} の並びに対しても、合計が 3n3n 未満」という決定論的な保証です。ハッシュ表の平均計算量(定理 5.3[探索アルゴリズム])のように「入力が一様分布なら」という仮定を置く議論とは、意味も強さも違います。この区別については 注意 6.5[探索アルゴリズム] も参照してください。

また、この定理は増やし方が定数倍であること(コピー回数が幾何級数になること)に依存しています。「毎回 +1+1」にすると kk 回目の push\mathrm{push}k1k-1 個のコピーが発生し、合計 k=1n(k1)=Θ(n2)\sum_{k=1}^{n}(k-1) = \Theta(n^2) となって償却 O(1)O(1) は崩れます。逆に、比率が 11 より大きい定数でありさえすれば 22 倍である必要はなく、1.51.5 倍でも 1.1251.125 倍でも同じ議論が通ります。

配列の弱点は「番地が添字で決まってしまう」ことでした。ならば、番地を要素自身に持たせればよい、というのが連結リストの発想です。

定義 4.1単方向連結リスト

ノードとは、値を保持するフィールド value\mathrm{value} と、別のノードへの参照または nil\mathrm{nil} を保持するフィールド next\mathrm{next} からなる区画をいう。ノード p0,p1,,pn1p_0, p_1, \ldots, p_{n-1}

pj.next=pj+1(0jn2),pn1.next=nilp_j.\mathrm{next} = p_{j+1} \quad (0 \le j \le n-2), \qquad p_{n-1}.\mathrm{next} = \mathrm{nil}

を満たすとき、この列を単方向連結リストといい、先頭ノードへの参照 head=p0\mathrm{head} = p_0 によってリスト全体を表す。空のリストは head=nil\mathrm{head} = \mathrm{nil} で表す。各ノードが prev\mathrm{prev} フィールドも持ち pj+1.prev=pjp_{j+1}.\mathrm{prev} = p_j を満たすものを双方向連結リストという。

ノードの番地には何の制約もありません。記憶領域のどこに散らばっていてもよく、つながりは next\mathrm{next} の値だけが表しています。この自由さが、そのまま計算量の特徴になります。

命題 4.2連結リストの基本操作の計算量

nn 個のノードからなる単方向連結リストについて、ワード RAM モデルの下で次が成り立つ。

  1. ノード pp への参照が与えられているとき、pp の直後への新しいノードの挿入、および pp の直後のノードの削除は、いずれも Θ(1)\Theta(1) 時間で行える。
  2. 先頭から数えて kk 番目(00 始まり、0kn10 \le k \le n-1)のノードへ到達するには、next\mathrm{next} の追跡がちょうど kk 回必要であり、Θ(k)\Theta(k) 時間かかる。最悪は Θ(n)\Theta(n) である。
  3. 単方向連結リストでは、ノード pp への参照が与えられていても、pp 自身の削除には最悪 Θ(n)\Theta(n) 時間かかる。双方向連結リストであれば Θ(1)\Theta(1) 時間で行える。
証明(命題 4.2)

1. 挿入は、新しいノード qq を用意して

q.nextp.next,p.nextqq.\mathrm{next} \leftarrow p.\mathrm{next}, \qquad p.\mathrm{next} \leftarrow q

の 2 回の書き込みで完了します。実行後、pp の直後は qqqq の直後は元の p.nextp.\mathrm{next} となり、定義 4.1 の連結条件が保たれます。削除は、p.nextnilp.\mathrm{next} \ne \mathrm{nil} のとき

p.nextp.next.nextp.\mathrm{next} \leftarrow p.\mathrm{next}.\mathrm{next}

の 1 回の書き込みで完了します。どちらも読み書きの回数が nn にも kk にも依存しないので Θ(1)\Theta(1) 時間です。命題 3.2 の 3 では位置 ii の挿入に nin-i 回の移動が必要でしたが、ここでは要素をひとつも動かしていません。

2. 上界は、head\mathrm{head} から始めて next\mathrm{next}kk 回たどる手続きが与えます。各追跡は 1 語の読み出しなので定数時間、合計 Θ(k)\Theta(k) です。

下界については、このモデルでノードの番地を知る手段が「head\mathrm{head} を読む」「既に到達したノードの next\mathrm{next} を読む」の 2 つしかないことを使います。定義 3.1 の配列と違い、kk 番目のノードの番地を kk から計算する式は存在しません。したがって到達済みノードの集合は 1 回の追跡で高々 1 個しか増えず、pkp_k に至るには少なくとも kk 回の追跡が必要です。以上より Θ(k)\Theta(k) です。

3. pp 自身の削除には、直前のノード pp' を見つけて p.nextp.nextp'.\mathrm{next} \leftarrow p.\mathrm{next} とする必要があります。単方向リストでは pp から pp' をたどれないため、head\mathrm{head} から順に next\mathrm{next}pp に一致するノードを探すしかなく、2 より最悪 Θ(n)\Theta(n) です。双方向連結リストなら p=p.prevp'= p.\mathrm{prev} が定数時間で得られるので、

p.prev.nextp.next,p.next.prevp.prevp.\mathrm{prev}.\mathrm{next} \leftarrow p.\mathrm{next}, \qquad p.\mathrm{next}.\mathrm{prev} \leftarrow p.\mathrm{prev}

の 2 回の書き込みで済み、Θ(1)\Theta(1) 時間です(両端の nil\mathrm{nil} の場合分けは 注意 4.4 の番兵で消せます)。

例 4.3挿入と削除を実際に動かす

Python で 命題 4.2 の 1 を確かめます。

class Node:
__slots__ = ("value", "next")
def __init__(self, value, next=None):
self.value = value
self.next = next
def insert_after(p: Node, x) -> Node:
"""p の直後に値 x のノードを挿入して、そのノードを返す。"""
q = Node(x, p.next) # q.next <- p.next
p.next = q # p.next <- q
return q
def delete_after(p: Node) -> None:
"""p の直後のノードを削除する。直後が無ければ何もしない。"""
if p.next is not None:
p.next = p.next.next
def to_list(head: Node) -> list:
out, cur = [], head
while cur is not None:
out.append(cur.value)
cur = cur.next
return out
head = Node(3, Node(1, Node(4)))
print(to_list(head)) # [3, 1, 4]
insert_after(head, 5)
print(to_list(head)) # [3, 5, 1, 4]
delete_after(head.next)
print(to_list(head)) # [3, 5, 4]

insert_after の本体は 2 行、delete_after は 1 行で、どちらもリストの長さに依存するループを含みません。これが Θ(1)\Theta(1) の中身です。一方 to_list は末尾まで next\mathrm{next} をたどるので Θ(n)\Theta(n) かかり、長さを知りたいだけでも全体の走査が要ります(要素数を別に保持すれば Θ(1)\Theta(1) です)。

注意 4.4番兵ノード

実装では、値を持たないダミーのノード(番兵、sentinel)を先頭に 1 つ置き、head\mathrm{head} の代わりに番兵への参照を保持する手法がよく使われます。こうすると「先頭への挿入」「先頭の削除」が「番兵の直後への挿入・削除」に統一され、命題 4.2 の 1 がそのまま適用できます。空リストの判定は「番兵の next\mathrm{next}nil\mathrm{nil} かどうか」になります。場合分けが減るぶんバグも減るので、実装するときは番兵から始めるとよいでしょう。

ここまでの結果を並べます。まず、両者のメモリ上の姿を見比べてください。

配列: 連続した番地に並ぶ。A[i] の番地は b + s·i と計算できる314159bb+sb+2sb+3sb+4sb+5sA[0]A[1]A[2]A[3]A[4]A[5]連結リスト: ノードの番地はばらばら。各ノードが値と「次への参照」を持つ3141head次への参照nil
配列と連結リストのメモリ配置。配列は番地が添字から計算でき、連結リストは参照をたどるしかない。

命題 3.2命題 4.2、および 定理 3.3 をまとめると、次の表になります。連結リストは末尾ノードへの参照も保持しているものとします。

操作動的配列単方向連結リスト
ii 番目の読み書きΘ(1)\Theta(1)Θ(i)\Theta(i)(最悪 Θ(n)\Theta(n)
値による探索(順序の仮定なし)Θ(n)\Theta(n)Θ(n)\Theta(n)
値による探索(整列済み)Θ(logn)\Theta(\log n)(二分探索)Θ(n)\Theta(n)
先頭への挿入・先頭の削除Θ(n)\Theta(n)Θ(1)\Theta(1)
末尾への追加償却 Θ(1)\Theta(1)Θ(1)\Theta(1)
末尾の削除Θ(1)\Theta(1)Θ(n)\Theta(n)(直前ノードが必要)
添字 ii で指定した位置への挿入・削除Θ(ni)\Theta(n-i)Θ(i)\Theta(i)
手元にあるノードの直後への挿入・削除Θ(ni)\Theta(n-i)Θ(1)\Theta(1)
要素以外に必要な記憶領域未使用の容量(最大で要素数と同程度)ノードごとに参照 1 個分

表の最後から 2 行目が両者の性格をもっとも鋭く表しています。「今いる場所の隣に入れる」は連結リストでは Θ(1)\Theta(1)、配列では Θ(ni)\Theta(n-i)。逆に「ii 番目に飛ぶ」は配列では Θ(1)\Theta(1)、連結リストでは Θ(i)\Theta(i) です。選択は、この 2 種類の操作をそれぞれ何回使うかで決まります。

注意 5.1定数因子とキャッシュ

Θ\Theta 記法は定数因子を隠します。実測では、隠された部分が効くことがあります。現代の CPU は主記憶を語単位ではなくキャッシュライン単位(典型的には 64 バイト)で読み込むため、配列を先頭から走査すると 1 回の読み込みで複数の要素がキャッシュに載ります。連結リストの走査では、次のノードの番地が読んでみるまで分からず、ノードが記憶領域に散らばっているほどキャッシュミスが増えます。

結果として、「途中への挿入が多いから連結リスト」という判断が、要素数が数千程度までは実測で裏切られることがあります。迷ったら、まず動的配列で書いて測ってから考えるのがよいでしょう。

ここからは抽象データ型の側に移ります。スタックは、最後に入れたものが最初に出てくる(last in, first out, LIFO)という一点だけを規定した型です。

定義 6.1スタック

スタックとは、次の操作をもつ抽象データ型である。値の型を XX、スタックの全体を S\mathcal{S} と書く。

  • emptyS\mathrm{empty} \in \mathcal{S}(空のスタック)
  • push:S×XS\mathrm{push} : \mathcal{S} \times X \to \mathcal{S}(積む)
  • pop:SS\mathrm{pop} : \mathcal{S} \to \mathcal{S}(頂上を取り除く)
  • top:SX\mathrm{top} : \mathcal{S} \to X(頂上を読む)
  • isEmpty:S{true,false}\mathrm{isEmpty} : \mathcal{S} \to \{\mathrm{true}, \mathrm{false}\}

これらは任意の SSS \in \mathcal{S}xXx \in X について次の公理を満たす。

isEmpty(empty)=true,isEmpty(push(S,x))=false,top(push(S,x))=x,pop(push(S,x))=S.\begin{aligned} &\mathrm{isEmpty}(\mathrm{empty}) = \mathrm{true}, \qquad \mathrm{isEmpty}(\mathrm{push}(S, x)) = \mathrm{false}, \\ &\mathrm{top}(\mathrm{push}(S, x)) = x, \qquad \mathrm{pop}(\mathrm{push}(S, x)) = S. \end{aligned}

top(empty)\mathrm{top}(\mathrm{empty})pop(empty)\mathrm{pop}(\mathrm{empty}) は定義されない(実装ではエラーとする)。

公理 pop(push(S,x))=S\mathrm{pop}(\mathrm{push}(S, x)) = S が LIFO そのものです。xx を積んでから取り除くと、積む前の状態にぴったり戻る、と言っています。この 4 本の等式のどこにも配列やノードは現れません。実装が何であれ、これらを満たせばスタックです。実装は 2 通りあり、どちらも既に得た結果から計算量が読めます。

動的配列による実装. 要素を A[0],,A[n1]A[0], \ldots, A[n-1] と並べ、頂上を A[n1]A[n-1] とします。push\mathrm{push} は末尾追加なので償却 Θ(1)\Theta(1)定理 3.3)、pop\mathrm{pop}nn1n \leftarrow n-1top\mathrm{top}A[n1]A[n-1] の読み出しなので最悪 Θ(1)\Theta(1) です(命題 3.2 の 1)。命題 3.2 の 3 が要求する要素の移動は、末尾での操作に限れば ni=0n-i = 0 となって発生しません。

連結リストによる実装. 頂上を先頭ノードとし、push\mathrm{push} を先頭への挿入、pop\mathrm{pop} を先頭の削除とします。番兵を置けばどちらも 命題 4.2 の 1 に帰着し、最悪 Θ(1)\Theta(1) です。償却ではなく最悪で Θ(1)\Theta(1) なので、1 回の応答時間に上限が要る場面ではこちらが有利です。代わりにノードごとに参照 1 個分の領域がかかります。

スタックの使い道を、証明のつく形で一つ見ておきます。括弧の対応判定です。

定理 6.2括弧列の判定

記号の集合を Σ={(,),[,]}\Sigma = \{\,\texttt{(}\,,\,\texttt{)}\,,\,\texttt{[}\,,\,\texttt{]}\,\} とし、対応の取れた括弧列の集合 BΣB \subseteq \Sigma^{*} を、次の 3 つの規則で生成される最小の集合として定める。

εB,wB    (w)B かつ [w]B,u,vB    uvB.\varepsilon \in B, \qquad w \in B \implies \texttt{(} w \texttt{)} \in B \ \text{かつ}\ \texttt{[} w \texttt{]} \in B, \qquad u, v \in B \implies uv \in B.

入力 wΣw \in \Sigma^{*} に対する次のアルゴリズム A\mathcal{A} を考える。空のスタック SS から始め、ww を左から 1 文字ずつ読む。

  • 読んだ文字が開き括弧なら、それを SSpush\mathrm{push} する。
  • 読んだ文字が閉じ括弧なら、SS が空であれば直ちに拒否し、そうでなければ top\mathrm{top} を見て pop\mathrm{pop} し、取り出した文字がその閉じ括弧と同種の開き括弧でなければ拒否する。

全文字を読み終えた時点で SS が空なら受理、空でなければ拒否する。このとき、A\mathcal{A}ww を受理することと wBw \in B であることは同値である。また A\mathcal{A} の実行時間は Θ(w)\Theta(|w|) である。

証明(定理 6.2)

準備(スタックの相対性). A\mathcal{A} はスタックの頂上しか参照せず、空のスタックから pop\mathrm{pop} しようとした時点で拒否します。したがって、初期スタックを σ\sigma として uu を処理する実行と、空スタックから同じ uu を処理する実行は、後者が「空からの pop\mathrm{pop}」で拒否する場合を除いて、σ\sigma より上の変化がまったく同じです。以下この事実を繰り返し使います。

(\Rightarrow) wBw \in B ならば受理する. 次の、より強い主張を BB の生成規則に関する構造帰納法で示します。

wBw \in B ならば、任意の初期スタック σ\sigma から ww を処理したとき A\mathcal{A} は拒否せず、処理後のスタックは σ\sigma に戻る。

  • w=εw = \varepsilon のとき。何も読まないので拒否せず、スタックは σ\sigma のままです。
  • w=(u)w = \texttt{(} u \texttt{)}uBu \in B)のとき。最初の (\texttt{(} を読むとスタックは σ(\sigma\texttt{(} になります。帰納法の仮定を初期スタック σ(\sigma\texttt{(}uu に適用すると、拒否せず、処理後のスタックは σ(\sigma\texttt{(} に戻ります。次に )\texttt{)} を読むと、スタックは空でなく頂上は (\texttt{(} なので pop\mathrm{pop} して種類が一致し、拒否しません。処理後のスタックは σ\sigma です。w=[u]w = \texttt{[} u \texttt{]} も同じ議論です。
  • w=uvw = uvu,vBu, v \in B)のとき。帰納法の仮定を初期スタック σ\sigmauu に適用すると、拒否せずスタックは σ\sigma に戻ります。続けて同じ仮定を vv に適用すれば、やはり拒否せずスタックは σ\sigma です。

とくに σ\sigma を空スタックとすれば、A\mathcal{A} は拒否せず、読み終えたときのスタックは空なので受理します。

(\Leftarrow) 受理するならば wBw \in B. w|w| に関する強い帰納法で示します。w=εw = \varepsilon なら生成規則の第 1 条より wBw \in B です。w1|w| \ge 1 とします。

w1w_1 が閉じ括弧だとすると、その時点でスタックは空なので A\mathcal{A} は拒否し、仮定に反します。よって w1w_1 は開き括弧です。これを cc、対応する閉じ括弧を cˉ\bar{c} と書きます。cc は 1 文字目で push\mathrm{push} され、受理時にはスタックが空なので、どこかで pop\mathrm{pop} されます。それが jj 文字目を読んだときだとします。

cc はスタックの最下段にあり jj 文字目で初めて取り除かれるので、位置 2,,j12, \ldots, j-1 の処理中はスタックに cc が残り続けます。すなわち u=w2wj1u = w_2 \cdots w_{j-1} の処理は cc より上だけで完結し、jj 文字目の直前のスタックはちょうど cc 1 個です。準備の観察より、uu を空スタックから処理しても拒否せず、処理後は空になります。u<w|u| < |w| なので、帰納法の仮定より uBu \in B です。

jj 文字目では ccpop\mathrm{pop} され、A\mathcal{A} が拒否しなかったので wj=cˉw_j = \bar{c} です。残りの v=wj+1wwv = w_{j+1} \cdots w_{|w|} の処理は空スタックから始まり、全体が受理されるので空スタックで終わります。v<w|v| < |w| より vBv \in B です。以上より w=cucˉvw = c\,u\,\bar{c}\,v で、生成規則の第 2 条から cucˉBc u \bar{c} \in B、第 3 条から wBw \in B を得ます。

計算量. 各文字につき push\mathrm{push}pop\mathrm{pop}top\mathrm{top} が高々 1 回ずつで、連結リスト実装ならこれらは最悪 Θ(1)\Theta(1) です(命題 4.2 の 1)。よって全体で Θ(w)\Theta(|w|) です。

例 6.3判定アルゴリズムを走らせる

w=([])()w = \texttt{([])()} に対する A\mathcal{A} の実行を追います。スタックは左を底として書きます。

読んだ文字動作直後のスタック
(\texttt{(}push(\texttt{(}
[\texttt{[}push([\texttt{([}
]\texttt{]}pop([\texttt{[} と一致)(\texttt{(}
)\texttt{)}pop((\texttt{(} と一致)
(\texttt{(}push(\texttt{(}
)\texttt{)}pop((\texttt{(} と一致)

読み終えてスタックが空なので受理します。実際 ([])B\texttt{([])} \in B かつ ()B\texttt{()} \in B で、生成規則の第 3 条から wBw \in B です。一方 w=(]w' = \texttt{(]} では、2 文字目で pop\mathrm{pop} した結果が (\texttt{(} となり ]\texttt{]} と種類が違うので拒否されます。w=(()w'' = \texttt{(()} では拒否は起きませんが、読み終えたときスタックに (\texttt{(} が 1 個残るので拒否されます。実装は次のとおりです。

def is_balanced(w: str) -> bool:
partner = {")": "(", "]": "["}
stack = []
for c in w:
if c in "([":
stack.append(c)
elif c in ")]":
if not stack or stack.pop() != partner[c]:
return False
else:
raise ValueError(f"unexpected symbol: {c!r}")
return not stack
print(is_balanced("([])()")) # True
print(is_balanced("(]")) # False
print(is_balanced("(()")) # False

最後の return not stack が「読み終えたときスタックが空」という受理条件に対応します。この 1 行を忘れると "(((" を受理してしまいます。

この判定は、深さ優先探索や再帰呼び出しの管理と同じ形をしています。「開いたものを、開いた順序の逆に閉じる」構造が現れる場所には、たいていスタックがあります。

キューは、最初に入れたものが最初に出てくる(first in, first out, FIFO)型です。窓口の行列と同じ規律で、待ち行列とも呼ばれます。

定義 7.1キュー

キューとは、次の操作をもつ抽象データ型である。値の型を XX、キューの全体を Q\mathcal{Q} と書く。

  • emptyQ\mathrm{empty} \in \mathcal{Q}(空のキュー)
  • enqueue:Q×XQ\mathrm{enqueue} : \mathcal{Q} \times X \to \mathcal{Q}(末尾に加える)
  • dequeue:QQ\mathrm{dequeue} : \mathcal{Q} \to \mathcal{Q}(先頭を取り除く)
  • front:QX\mathrm{front} : \mathcal{Q} \to X(先頭を読む)
  • isEmpty:Q{true,false}\mathrm{isEmpty} : \mathcal{Q} \to \{\mathrm{true}, \mathrm{false}\}

これらは任意の QQQ \in \mathcal{Q}xXx \in X について次の公理を満たす。

isEmpty(empty)=true,isEmpty(enqueue(Q,x))=false,front(enqueue(Q,x))={x(Q=empty)front(Q)(Qempty)dequeue(enqueue(Q,x))={empty(Q=empty)enqueue(dequeue(Q),x)(Qempty)\begin{aligned} &\mathrm{isEmpty}(\mathrm{empty}) = \mathrm{true}, \qquad \mathrm{isEmpty}(\mathrm{enqueue}(Q, x)) = \mathrm{false}, \\[2pt] &\mathrm{front}(\mathrm{enqueue}(Q, x)) = \begin{cases} x & (Q = \mathrm{empty}) \\ \mathrm{front}(Q) & (Q \ne \mathrm{empty}) \end{cases} \\[2pt] &\mathrm{dequeue}(\mathrm{enqueue}(Q, x)) = \begin{cases} \mathrm{empty} & (Q = \mathrm{empty}) \\ \mathrm{enqueue}(\mathrm{dequeue}(Q), x) & (Q \ne \mathrm{empty}) \end{cases} \end{aligned}

front(empty)\mathrm{front}(\mathrm{empty})dequeue(empty)\mathrm{dequeue}(\mathrm{empty}) は定義されない。

定義 6.1 と見比べてください。スタックでは pop(push(S,x))=S\mathrm{pop}(\mathrm{push}(S,x)) = S と、直前に積んだものがそのまま取れました。キューでは dequeue\mathrm{dequeue}enqueue\mathrm{enqueue} をすり抜けて奥へ潜っていきます。この「すり抜け」が FIFO の正体で、公理の右辺に再帰が現れる理由です。

素朴な配列実装はうまくいきません. 先頭を常に A[0]A[0] に置くと、enqueue\mathrm{enqueue} は末尾追加なので償却 Θ(1)\Theta(1) ですが、dequeue\mathrm{dequeue} は位置 00 の削除であり 命題 3.2 の 3 より Θ(n)\Theta(n)nn 回で Θ(n2)\Theta(n^2) です。原因は先頭を A[0]A[0] に固定したことなので、固定をやめれば解決します。

命題 7.2リングバッファ

容量 mm の配列 A[0..m1]A[0..m-1]、先頭位置 hh0hm10 \le h \le m-1)、要素数 nn0nm0 \le n \le m)を保持し、次のように操作を定める。

  • enqueue(x)\mathrm{enqueue}(x)n<mn < m のとき A[(h+n)modm]xA[(h + n) \bmod m] \leftarrow xnn+1n \leftarrow n + 1
  • dequeue()\mathrm{dequeue}()n>0n > 0 のとき戻り値を A[h]A[h] とし、h(h+1)modmh \leftarrow (h + 1) \bmod mnn1n \leftarrow n - 1
  • front()\mathrm{front}()n>0n > 0 のとき A[h]A[h] を返す。

このとき、キューの内容を先頭から順に q0,q1,,qn1q_0, q_1, \ldots, q_{n-1} とすると、常に qi=A[(h+i)modm]q_i = A[(h + i) \bmod m] が成り立ち、定義 7.1 の公理が満たされる。また 3 つの操作はいずれも最悪 Θ(1)\Theta(1) 時間である。

証明(命題 7.2)

主張の等式 qi=A[(h+i)modm]q_i = A[(h+i) \bmod m]0in10 \le i \le n-1)を不変条件として、操作回数に関する帰納法で示します。

初期状態は n=0n = 0 なので、条件は空虚に成り立ちます。

enqueue(x)\mathrm{enqueue}(x) の場合。新しい内容は q0,,qn1,xq_0, \ldots, q_{n-1}, xhh は変わりません。書き込み先 (h+n)modm(h+n) \bmod m が既存の位置 (h+i)modm(h+i) \bmod m0in10 \le i \le n-1)と重ならないことを確かめます。0i<nm10 \le i < n \le m-1 より iinn の差は mm の倍数になり得ないので、mm を法として異なります。したがって既存要素は壊れず、新しい末尾要素は番号 nn の要素として等式を満たします。

dequeue()\mathrm{dequeue}() の場合。戻り値は A[h]=q0A[h] = q_0 でキューの先頭です。新しい内容 q1,,qn1q_1, \ldots, q_{n-1} を改めて q0,,qn2q'_0, \ldots, q'_{n-2} と番号づけ、h=(h+1)modmh' = (h+1) \bmod m とすると

qi=qi+1=A[(h+i+1)modm]=A[(h+i)modm]q'_i = q_{i+1} = A[(h + i + 1) \bmod m] = A[(h' + i) \bmod m]

となって不変条件が保たれます。

front()\mathrm{front}()q0q_0 を返すことは i=0i = 0 の場合そのものです。以上より dequeue\mathrm{dequeue}front\mathrm{front} はつねに先頭要素に作用し、enqueue\mathrm{enqueue} は末尾にのみ加えるので、定義 7.1 の公理が満たされます。計算量は、各操作が定数回の加算・剰余・比較と 1 回の配列アクセス(命題 3.2 の 1 より Θ(1)\Theta(1))からなることから最悪 Θ(1)\Theta(1) です。

例 7.3容量 4 のリングバッファ

m=4m = 4、初期状態 h=0h = 0n=0n = 0 から始めます。AA の未使用区画を _ と書きます。

操作書き込み先 / 読み出し元AAhhnn
enqueue(a)\mathrm{enqueue}(a)A[(0+0)mod4]=A[0]A[(0+0) \bmod 4] = A[0]a _ _ _01
enqueue(b)\mathrm{enqueue}(b)A[1]A[1]a b _ _02
enqueue(c)\mathrm{enqueue}(c)A[2]A[2]a b c _03
dequeue()\mathrm{dequeue}()A[0]=aA[0] = aa b c _12
enqueue(d)\mathrm{enqueue}(d)A[(1+2)mod4]=A[3]A[(1+2) \bmod 4] = A[3]a b c d13
enqueue(e)\mathrm{enqueue}(e)A[(1+3)mod4]=A[0]A[(1+3) \bmod 4] = A[0]e b c d14
dequeue()\mathrm{dequeue}()A[1]=bA[1] = be b c d23

最後の状態でキューの内容は c,d,ec, d, e であり、不変条件どおり A[(2+0)mod4]=A[2]=cA[(2+0) \bmod 4] = A[2] = cA[3]=dA[3] = dA[(2+2)mod4]=A[0]=eA[(2+2) \bmod 4] = A[0] = e です。ee を書き込むときに配列の右端を越えて左端へ回り込む点が「リング」の由来で、要素は一度も移動していません。

注意 7.4満杯と空の区別、そして容量拡張

hh と「末尾の次の位置」t=(h+n)modmt = (h+n) \bmod m だけを保持する実装もありますが、その場合 n=0n = 0n=mn = m のどちらでも t=ht = h となり、空と満杯が区別できません。回避策は 2 つで、上の 命題 7.2 のように要素数 nn を明示的に持つか、区画を 1 つ常に空けておいて実質容量を m1m-1 とするかです。

満杯になったら容量 2m2m の配列を確保して先頭から詰め直します。この拡張のコストと頻度は 定理 3.3 とまったく同じ形なので、enqueue\mathrm{enqueue} は償却 Θ(1)\Theta(1) に保たれます。連結リストで先頭と末尾の両方への参照を持つ実装なら、命題 4.2 の 1 より両操作とも最悪 Θ(1)\Theta(1) です。

スタックが深さ優先探索を、キューが幅優先探索を駆動します。どちらも グラフアルゴリズム の中核で、幅優先探索が最短距離を正しく求めることはキューの FIFO 規律から従い(定理 3.2[グラフアルゴリズム])、深さ優先探索では処理中の頂点がちょうどスタックをなします(補題 4.1[グラフアルゴリズム])。また、配列を「表」として使い添字アクセスの Θ(1)\Theta(1) に全面的に依存するのが 動的計画法 です(定理 3.2[動的計画法])。

8. ポテンシャル法と 2 本のスタックによるキュー

Section titled “8. ポテンシャル法と 2 本のスタックによるキュー”

定理 3.3 では総コストを直接数え上げました(集約法)。操作が複数種類あって互いに影響し合うと、この数え上げは難しくなります。そこで使うのがポテンシャル法で、データ構造の「たまった仕事」を実数値の関数で表し、各操作でその増減を見ます。

補題 8.1ポテンシャル法

データ構造の状態の列 D0,D1,,DmD_0, D_1, \ldots, D_mD0D_0 は初期状態、DiD_iii 番目の操作の直後の状態)と、ii 番目の操作の実コスト cic_i が与えられているとする。状態の集合上の実数値関数 Φ\Phi

Φ(Di)Φ(D0)(i=0,1,,m)\Phi(D_i) \ge \Phi(D_0) \qquad (i = 0, 1, \ldots, m)

を満たすとする。償却コストを c^i:=ci+Φ(Di)Φ(Di1)\hat{c}_i := c_i + \Phi(D_i) - \Phi(D_{i-1}) で定めると

i=1mci=i=1mc^i+Φ(D0)Φ(Dm)i=1mc^i\sum_{i=1}^{m} c_i = \sum_{i=1}^{m} \hat{c}_i + \Phi(D_0) - \Phi(D_m) \le \sum_{i=1}^{m} \hat{c}_i

が成り立つ。とくに、ある定数 aa がすべての ii について c^ia\hat{c}_i \le a を満たすならば、mm 回の操作の総コストは amam 以下である。

証明(補題 8.1)

定義から

i=1mc^i=i=1mci+i=1m(Φ(Di)Φ(Di1))\sum_{i=1}^{m} \hat{c}_i = \sum_{i=1}^{m} c_i + \sum_{i=1}^{m}\bigl(\Phi(D_i) - \Phi(D_{i-1})\bigr)

です。右辺の第 2 項は望遠鏡和で、隣接する項が打ち消し合って Φ(Dm)Φ(D0)\Phi(D_m) - \Phi(D_0) になります。よって

i=1mc^i=i=1mci+Φ(Dm)Φ(D0)\sum_{i=1}^{m} \hat{c}_i = \sum_{i=1}^{m} c_i + \Phi(D_m) - \Phi(D_0)

であり、移項すれば等式を得ます。仮定より Φ(Dm)Φ(D0)\Phi(D_m) \ge \Phi(D_0) すなわち Φ(D0)Φ(Dm)0\Phi(D_0) - \Phi(D_m) \le 0 なので不等式が従います。最後の主張は c^iam\sum \hat{c}_i \le am から明らかです。

Φ\Phi は「今の状態が抱えている借金」だと思ってください。安い操作で少しずつ借金を積み、高い操作でまとめて返済します。返済の瞬間は実コストが大きくても Φ\Phi が大きく減るので、償却コストは小さくなります。

定理 8.22 本のスタックによるキュー

全操作が最悪 Θ(1)\Theta(1) 時間のスタック(例えば 命題 4.2 の 1 に基づく連結リスト実装)を 2 本用意し、in\mathrm{in}out\mathrm{out} と呼ぶ。次のように操作を定める。

  • enqueue(x)\mathrm{enqueue}(x)in\mathrm{in}xxpush\mathrm{push} する。
  • dequeue()\mathrm{dequeue}()out\mathrm{out} が空ならば、in\mathrm{in} が空になるまで「in\mathrm{in} から pop\mathrm{pop} して out\mathrm{out}push\mathrm{push}」を繰り返す。その後 out\mathrm{out} が空ならエラー、そうでなければ out\mathrm{out} から pop\mathrm{pop} して返す。
  • front()\mathrm{front}():同じ移送を行ったうえで out\mathrm{out}top\mathrm{top} を返す。

このとき、この実装は 定義 7.1 の公理を満たす。さらに、空の状態から始めた任意の mm 回の操作列の総コストは 3m3m 以下であり、1 操作あたりの償却計算量は O(1)O(1) である。ただし個々の dequeue\mathrm{dequeue} の最悪計算量は、そのときの要素数を nn として Θ(n)\Theta(n) である。

証明(定理 8.2)

正当性. キューの内容を先頭から順に q0,,qn1q_0, \ldots, q_{n-1} とし、次を不変条件とします。

ある rr0rn0 \le r \le n)が存在して、out\mathrm{out} を頂上から底へ読むと q0,q1,,qr1q_0, q_1, \ldots, q_{r-1} となり、in\mathrm{in} を底から頂上へ読むと qr,qr+1,,qn1q_r, q_{r+1}, \ldots, q_{n-1} となる。

初期状態は n=0n = 0、両方のスタックが空で、r=0r = 0 として成立します。

enqueue(x)\mathrm{enqueue}(x) では in\mathrm{in} を底から読むと qr,,qn1,xq_r, \ldots, q_{n-1}, x となり、新しいキューの内容は q0,,qn1,xq_0, \ldots, q_{n-1}, x ですから、同じ rr で不変条件が保たれます。

移送が起きるのは out\mathrm{out} が空、すなわち r=0r = 0 のときで、このとき in\mathrm{in} を底から読むと q0,,qn1q_0, \ldots, q_{n-1} です。in\mathrm{in} から pop\mathrm{pop} される順序は qn1,,q0q_{n-1}, \ldots, q_0 で、これを順に out\mathrm{out}push\mathrm{push} するので、out\mathrm{out} を底から読むと qn1,,q0q_{n-1}, \ldots, q_0、すなわち頂上からは q0,,qn1q_0, \ldots, q_{n-1} です。in\mathrm{in} は空なので r=nr = n として不変条件が成立します。順序が 2 度反転して元に戻る、というのがこの実装の仕掛けです。

out\mathrm{out} が空でないとき(r1r \ge 1)、その頂上は不変条件より q0q_0、すなわちキューの先頭です。よって front\mathrm{front}q0q_0 を返し、dequeue\mathrm{dequeue}q0q_0 を取り除きます。取り除いた後は out\mathrm{out} の頂上から q1,,qr1q_1, \ldots, q_{r-1} が並ぶので、番号を付け替えれば r=r1r' = r - 1 で不変条件が保たれます。n=0n = 0 のときは両方が空で、移送しても out\mathrm{out} は空のままなのでエラーとなり、定義 7.1dequeue(empty)\mathrm{dequeue}(\mathrm{empty}) を未定義としていることと整合します。

償却計算量. ポテンシャルを

Φ:=2(in の要素数)\Phi := 2 \cdot (\mathrm{in} \text{ の要素数})

と定めます。Φ0\Phi \ge 0 で、初期状態では Φ(D0)=0\Phi(D_0) = 0 なので、補題 8.1 の仮定 Φ(Di)Φ(D0)\Phi(D_i) \ge \Phi(D_0) が満たされます。コストは 1 回のスタック操作を 1 と数えます(仮定よりこれは最悪 Θ(1)\Theta(1) 時間です)。

  • enqueue\mathrm{enqueue}:実コストは push\mathrm{push} 1 回で c=1c = 1in\mathrm{in} の要素数が 1 増えるので ΔΦ=2\Delta\Phi = 2。よって c^=1+2=3\hat{c} = 1 + 2 = 3
  • out\mathrm{out} が空でないときの dequeue\mathrm{dequeue}:実コストは pop\mathrm{pop} 1 回で c=1c = 1in\mathrm{in} は変わらないので ΔΦ=0\Delta\Phi = 0。よって c^=1\hat{c} = 1
  • out\mathrm{out} が空で in\mathrm{in}k1k \ge 1 個あるときの dequeue\mathrm{dequeue}:移送で pop\mathrm{pop}kk 回、push\mathrm{push}kk 回、その後の pop\mathrm{pop} が 1 回なので c=2k+1c = 2k + 1in\mathrm{in} の要素数が kk から 00 になるので ΔΦ=2k\Delta\Phi = -2k。よって c^=2k+12k=1\hat{c} = 2k + 1 - 2k = 1
  • front\mathrm{front}pop\mathrm{pop} の代わりに top\mathrm{top} を読むだけなので、上の 2 つと同じ計算で c^=1\hat{c} = 1

いずれの場合も c^3\hat{c} \le 3 です。補題 8.1 より、mm 回の操作の総コストは 3m3m 以下であり、1 操作あたりの償却計算量は O(1)O(1) です。

最悪計算量. nn 個すべてが in\mathrm{in} に積まれ out\mathrm{out} が空の状態での dequeue\mathrm{dequeue}2n+12n + 1 回のスタック操作を行うので Θ(n)\Theta(n) です。償却が O(1)O(1) であることと矛盾しません。この高価な dequeue\mathrm{dequeue} が起きる前には、nn 個の要素を積むために nn 回の enqueue\mathrm{enqueue} が必要だからです。

flowchart LR
E["enqueue x は in に push する"] --> I["スタック in(底から c, d, e)"]
I -->|"out が空のとき in が空になるまで移す"| O["スタック out(頂上から c, d, e)"]
O --> D["dequeue は out の頂上を pop する"]
2 本のスタックによるキュー。in から out へ移すときに順序が反転し、最も古い要素が out の頂上に来る。

例 8.32 本のスタックを動かす

実装と、Φ\Phi の変化を追った実行例を示します。

class TwoStackQueue:
def __init__(self):
self._in = [] # 末尾が頂上
self._out = [] # 末尾が頂上
def enqueue(self, x):
self._in.append(x)
def _transfer(self):
while self._in:
self._out.append(self._in.pop())
def dequeue(self):
if not self._out:
self._transfer()
if not self._out:
raise IndexError("dequeue from empty queue")
return self._out.pop()
def front(self):
if not self._out:
self._transfer()
if not self._out:
raise IndexError("front of empty queue")
return self._out[-1]
def __len__(self):
return len(self._in) + len(self._out)
q = TwoStackQueue()
for x in "abc":
q.enqueue(x)
print(q.dequeue(), q.dequeue()) # a b
q.enqueue("d")
print(q.dequeue(), q.dequeue()) # c d

各操作の実コスト cc、ポテンシャル Φ=2in\Phi = 2\,|\mathrm{in}|、償却コスト c^=c+ΔΦ\hat{c} = c + \Delta\Phi は次のようになります。

操作直後の in\mathrm{in}(底から)直後の out\mathrm{out}(頂上から)ccΦ\Phic^\hat{c}
初期状態0
enqueue(a)\mathrm{enqueue}(a)aa123
enqueue(b)\mathrm{enqueue}(b)a,ba, b143
enqueue(c)\mathrm{enqueue}(c)a,b,ca, b, c163
dequeue()=a\mathrm{dequeue}() = ab,cb, c701
dequeue()=b\mathrm{dequeue}() = bcc101
enqueue(d)\mathrm{enqueue}(d)ddcc123
dequeue()=c\mathrm{dequeue}() = cdd121
dequeue()=d\mathrm{dequeue}() = d301

4 行目の dequeue\mathrm{dequeue} は実コスト 77pop\mathrm{pop} 3 回、push\mathrm{push} 3 回、pop\mathrm{pop} 1 回)と飛び抜けていますが、Φ\Phi66 から 00 へ落ちるので償却コストは 11 です。総実コストは 1616、操作回数は 88 で、補題 8.1 の上界 3×8=243 \times 8 = 24 を下回ります。

演習 9.1

定義 7.1 の公理だけを使って、次の等式を導いてください。途中でどの公理を使ったかを明記してください。

dequeue(enqueue(enqueue(empty,a),b))=enqueue(empty,b)\mathrm{dequeue}\bigl(\mathrm{enqueue}(\mathrm{enqueue}(\mathrm{empty}, a), b)\bigr) = \mathrm{enqueue}(\mathrm{empty}, b)

また、この結果と front\mathrm{front} の公理から、aabb の順に入れたキューから 1 つ取り出すと bb が残ることを確かめてください。

解答

Q:=enqueue(empty,a)Q := \mathrm{enqueue}(\mathrm{empty}, a) と置きます。公理 isEmpty(enqueue(Q,x))=false\mathrm{isEmpty}(\mathrm{enqueue}(Q', x)) = \mathrm{false}Q=emptyQ' = \mathrm{empty}x=ax = a に適用すると QemptyQ \ne \mathrm{empty} です。よって dequeue\mathrm{dequeue} の公理の第 2 の場合が使えて

dequeue(enqueue(Q,b))=enqueue(dequeue(Q),b)\mathrm{dequeue}(\mathrm{enqueue}(Q, b)) = \mathrm{enqueue}(\mathrm{dequeue}(Q), b)

となります。次に dequeue(Q)=dequeue(enqueue(empty,a))\mathrm{dequeue}(Q) = \mathrm{dequeue}(\mathrm{enqueue}(\mathrm{empty}, a)) に、同じ公理の第 1 の場合(内側の引数が empty\mathrm{empty})を適用すると dequeue(Q)=empty\mathrm{dequeue}(Q) = \mathrm{empty} です。代入して

dequeue(enqueue(Q,b))=enqueue(empty,b)\mathrm{dequeue}(\mathrm{enqueue}(Q, b)) = \mathrm{enqueue}(\mathrm{empty}, b)

を得ます。

さらに front\mathrm{front} の公理の第 1 の場合より front(enqueue(empty,b))=b\mathrm{front}(\mathrm{enqueue}(\mathrm{empty}, b)) = b です。先に入れた aa が先に出ていったので、確かに FIFO です。定義 6.1 の公理で同じ計算をすると pop(push(push(empty,a),b))=push(empty,a)\mathrm{pop}(\mathrm{push}(\mathrm{push}(\mathrm{empty}, a), b)) = \mathrm{push}(\mathrm{empty}, a) となり、残るのは aa です。公理の違いがそのまま挙動の違いになっています。

演習 9.2標準

単方向連結リストを、Θ(n)\Theta(n) 時間・O(1)O(1) 追加領域で反転する手続きを書いてください。新しいノードを作らず、既存のノードの next\mathrm{next} を書き換えるだけで行うこと。ループ不変条件を述べ、それを使って正当性を示してください。

解答

3 つの変数だけを使います。

def reverse(head):
prev = None
cur = head
while cur is not None:
nxt = cur.next # 先に控えないと next を上書きした瞬間に迷子になる
cur.next = prev # 向きを反転
prev = cur
cur = nxt
return prev

不変条件. 元のリストを p0,p1,,pn1p_0, p_1, \ldots, p_{n-1} とします。ループの各反復の開始時点で、ある kk0kn0 \le k \le n)について次が成り立ちます。

  • prev は、pk1,pk2,,p0p_{k-1}, p_{k-2}, \ldots, p_0 をこの順に連ねたリストの先頭(k=0k = 0 のときは nil\mathrm{nil})である。
  • curpkp_kk=nk = n のときは nil\mathrm{nil})である。
  • この 2 本のリストは互いに素で、合わせて元のノード全体をなす。

初期化. ループに入る直前は prev = Nonecur = head =p0= p_0 なので、k=0k = 0 で成立します。

維持. 反復開始時に kk で成立し、cur =pknil= p_k \ne \mathrm{nil} とします。nxtpk+1p_{k+1}(存在しなければ nil\mathrm{nil})を控えます。cur.next = prev により pkp_k の直後が pk1p_{k-1} になり、prev から始まるリストは pk,pk1,,p0p_k, p_{k-1}, \ldots, p_0 になります。この代入で失われるのは pkp_k から pk+1p_{k+1} への参照だけですが、それは nxt に控えてあります。最後に prev = curcur = nxt とすれば、k+1k+1 について不変条件が成り立ちます。

終了. ループは curnil\mathrm{nil} になったとき、すなわち k=nk = n のときに終わります。不変条件より prevpn1,,p0p_{n-1}, \ldots, p_0 を連ねたリストの先頭で、これは元のリストの反転です。

計算量. ループは各ノードにつきちょうど 1 回、合計 nn 回まわり、1 回の本体は代入 4 個で定数時間です(命題 4.2 の 1 と同じ理由です)。よって Θ(n)\Theta(n) 時間。追加で使う記憶領域は変数 3 個だけなので O(1)O(1) です。

なお nxt の退避を省くと、cur.next = prev の直後に pk+1p_{k+1} へ到達する手段が失われます。命題 4.2 の 2 の証明のとおり、ノードの番地は参照をたどってしか得られないからです。

演習 9.3

定理 3.3 の動的配列に、要素の削除 pop\mathrm{pop} と、記憶領域を返すための縮小規則を加えます。

  1. pop\mathrm{pop} の後、要素数が容量の半分以下になったら容量を半分にする」という規則を採用すると、償却 O(1)O(1) が成り立たなくなります。mm 回の操作で総コストが Θ(m2)\Theta(m^2) になる操作列を具体的に構成してください。
  2. 「要素数が容量の 1/41/4 を下回ったら容量を半分にする」という規則なら、push\mathrm{push}pop\mathrm{pop} の償却コストが定数に保たれます。連続する 2 回の容量変更の間に少なくとも何回の操作が必要かを評価して、その理由を説明してください。
解答

1. 容量 c=m0c = m_0、要素数 n=m0n = m_0(満杯)の状態から始め、push\mathrm{push}pop\mathrm{pop} を交互に繰り返します。

  • push\mathrm{push}:満杯なので容量が 2m02m_0 に拡張され、m0m_0 個がコピーされます。実行後は n=m0+1n = m_0 + 1c=2m0c = 2m_0
  • pop\mathrm{pop}:実行後 n=m0n = m_0 となり、m02m0/2=m0m_0 \le 2m_0/2 = m_0 なので縮小規則が発動し、容量が m0m_0 に戻って m0m_0 個がコピーされます。

これで最初の状態に完全に戻るので、以降同じことが繰り返されます。2 回の操作ごとに 2m02m_0 回のコピーが発生するので総コストは Θ(mm0)\Theta(m \cdot m_0) となり、m0=Θ(m)m_0 = \Theta(m) と取れば Θ(m2)\Theta(m^2) です。原因は、拡張の閾値と縮小の閾値が同じ点にあることです。

2. 閾値をずらすと、容量変更の直後に境界から遠い場所に着地します。容量変更の直後の状態を見ます。

  • 拡張の直後:容量 cc、要素数 n=c/2n = c/2
  • 縮小の直後:縮小前の容量を 2c2c、そのとき要素数は 2c/4=c/22c/4 = c/2 を下回った直後なので n=c/2n = c/211 単位の誤差を除く)。容量は cc

どちらの場合も nc/2n \approx c/2 です。次に拡張が起きるには nncc に達する必要があるので push\mathrm{push} が少なくとも c/2c/2 回、次に縮小が起きるには nnc/4c/4 を下回る必要があるので pop\mathrm{pop} が少なくとも c/4c/4 回必要です。したがって、連続する 2 回の容量変更の間には少なくとも c/4c/4 回の操作があります。

容量 cc の変更にかかるコストは Θ(c)\Theta(c) ですから、これを直前の c/4c/4 回以上の操作に割り振ると 1 回あたり定数です。形式的には、補題 8.1 のポテンシャルとして

Φ={2nc(nc/2)c/2n(n<c/2)\Phi = \begin{cases} 2n - c & (n \ge c/2) \\ c/2 - n & (n < c/2) \end{cases}

を取れば、拡張・縮小の有無による場合分けで償却コストが定数に抑えられることを確かめられます。どちらの式も容量変更の直後(n=c/2n = c/2)で 00 になり、境界へ近づくほど増える点が要点です。

  • T. H. Cormen, C. E. Leiserson, R. L. Rivest, C. Stein, Introduction to Algorithms, 4th ed., MIT Press, 2022 — 第 10 章(Elementary Data Structures)に配列・連結リスト・スタック・キューが、第 16 章(Amortized Analysis)に集約法・課金法・ポテンシャル法があります。動的配列(table doubling)の解析も第 16 章です。
  • D. E. Knuth, The Art of Computer Programming, Volume 1: Fundamental Algorithms, 3rd ed., Addison-Wesley, 1997 — §2.2(Linear Lists)。§2.2.1 がスタック・キュー・デック、§2.2.2 が逐次配置(配列)、§2.2.3 が連結配置(連結リスト)で、この記事の §3 から §7 に対応します。連結リストの歴史的経緯も §2.6 に記述があります。
  • A. V. Aho, J. E. Hopcroft, J. D. Ullman, Data Structures and Algorithms, Addison-Wesley, 1983 — 第 2 章(Basic Abstract Data Types)。抽象データ型を先に定め、実装を後から与えるという本記事の構成は、この本の流儀に沿っています。
  • R. E. Tarjan, “Amortized Computational Complexity”, SIAM Journal on Algebraic and Discrete Methods 6 (1985), 306–318. DOI: 10.1137/0606031 — 償却計算量とポテンシャル法を体系的に定式化した論文です。
  • R. Sedgewick, K. Wayne, Algorithms, 4th ed., Addison-Wesley, 2011 — §1.3(Bags, Queues, and Stacks)。可変長配列と連結リストの両方でスタック・キューを実装し、実測値を比較しています。
  • 石畑清『アルゴリズムとデータ構造』岩波書店(岩波講座 ソフトウェア科学 3)、1989 — 日本語で読める定評ある教科書です。線形リストの各種表現が詳しく扱われています。

Appendix: 主要言語の標準ライブラリとの対応

Section titled “Appendix: 主要言語の標準ライブラリとの対応”

この記事で扱ったデータ構造は、主要な言語の標準ライブラリにそのまま入っています。名前から実装が読み取れないことがあるので、対応を挙げておきます。

動的配列と連結リスト. C++ の std::vector、Java の ArrayList、Python の list が動的配列です。増加率は実装ごとに違い、Java の ArrayList はおよそ 1.51.5 倍、CPython の list はおよそ 1.1251.125 倍(容量は 0,4,8,16,25,35,46,58,72,88,0, 4, 8, 16, 25, 35, 46, 58, 72, 88, \ldots と増えます)ですが、注意 3.5 のとおり 11 より大きい定数であれば償却 O(1)O(1) は保たれます。連結リストは C++ の std::list(双方向)と std::forward_list(単方向)、Java の LinkedList(双方向)です。Python の標準ライブラリに純粋な連結リスト型がないのは、必要になる場面が少ないためで、注意 5.1 の事情も背景にあります。

スタックとキュー. C++ の std::stackstd::queue は、既存のコンテナに 定義 6.1定義 7.1 のインタフェースだけを見せるアダプタ(既定の土台は std::deque)で、抽象データ型と実装の分離がそのまま型に現れています。Java の ArrayDeque はリングバッファ(命題 7.2)で、スタックにもキューにも使えます。Python ではスタックに listappendpop)、キューに collections.deque(固定長ブロックの双方向連結リストで、両端の追加・削除が最悪 O(1)O(1))を使います。listpop(0) でキューとして使うと、命題 3.2 の 3 により 1 回あたり Θ(n)\Theta(n)nn 回で Θ(n2)\Theta(n^2) になります。実務で頻繁に見かける速度低下の原因です。

この記事の誤りを報告する ・運営: 夢現技研合同会社料金プラン利用条件特定商取引法に基づく表記

© 2026 夢現技研合同会社 ・本文の LLM への入力は自由です。コード例は MIT ライセンスです。