算法.md 55.6 KB
Newer Older
C
CyC2018 已提交
1
<!-- GFM-TOC -->
C
CyC2018 已提交
2 3
* [一、前言](#一前言)
* [二、算法分析](#二算法分析)
C
CyC2018 已提交
4 5 6 7
    * [数学模型](#数学模型)
    * [ThreeSum](#threesum)
    * [倍率实验](#倍率实验)
    * [注意事项](#注意事项)
C
CyC2018 已提交
8
* [三、栈和队列](#三栈和队列)
C
CyC2018 已提交
9 10
    * [](#栈)
    * [队列](#队列)
C
CyC2018 已提交
11 12 13 14 15 16 17
* [四、并查集](#四并查集)
    * [quick-find](#quick-find)
    * [quick-union](#quick-union)
    * [加权 quick-union](#加权-quick-union)
    * [路径压缩的加权 quick-union](#路径压缩的加权-quick-union)
    * [各种 union-find 算法的比较](#各种-union-find-算法的比较)
* [五、排序](#五排序)
C
CyC2018 已提交
18 19 20 21 22 23 24 25 26 27 28 29
    * [选择排序](#选择排序)
    * [冒泡排序](#冒泡排序)
    * [插入排序](#插入排序)
    * [希尔排序](#希尔排序)
    * [归并排序](#归并排序)
    * [快速排序](#快速排序)
    * [堆排序](#堆排序)
    * [桶排序](#桶排序)
    * [基数排序](#基数排序)
    * [外部排序](#外部排序)
    * [排序算法的比较](#排序算法的比较)
    * [Java 的排序算法实现](#java-的排序算法实现)
C
CyC2018 已提交
30
* [六、查找](#六查找)
C
CyC2018 已提交
31 32 33 34 35 36 37 38 39 40
    * [二分查找实现有序符号表](#二分查找实现有序符号表)
    * [二叉查找树](#二叉查找树)
    * [2-3 查找树](#2-3-查找树)
    * [红黑二叉查找树](#红黑二叉查找树)
    * [散列表](#散列表)
    * [应用](#应用)
* [参考资料](#参考资料)
<!-- GFM-TOC -->


C
CyC2018 已提交
41 42 43 44 45
# 一、前言

本文实现代码以及测试代码放在 [Algorithm](https://github.com/CyC2018/Algorithm)

# 二、算法分析
C
CyC2018 已提交
46 47 48 49 50 51 52 53 54 55 56 57

## 数学模型

### 1. 近似

N<sup>3</sup>/6-N<sup>2</sup>/2+N/3 \~ N<sup>3</sup>/6。使用 \~f(N) 来表示所有随着 N 的增大除以 f(N) 的结果趋近于 1 的函数。

### 2. 增长数量级

N<sup>3</sup>/6-N<sup>2</sup>/2+N/3 的增长数量级为 O(N<sup>3</sup>)。增长数量级将算法与它的实现隔离开来,一个算法的增长数量级为 O(N<sup>3</sup>) 与它是否用 Java 实现,是否运行于特定计算机上无关。

### 3. 内循环
C
CyC2018 已提交
58 59 60

执行最频繁的指令决定了程序执行的总时间,把这些指令称为程序的内循环。

C
CyC2018 已提交
61
### 4. 成本模型
C
CyC2018 已提交
62 63 64

使用成本模型来评估算法,例如数组的访问次数就是一种成本模型。

C
CyC2018 已提交
65
## ThreeSum
C
CyC2018 已提交
66

C
CyC2018 已提交
67
ThreeSum 用于统计一个数组中和为 0 的三元组数量。
C
CyC2018 已提交
68 69

```java
C
CyC2018 已提交
70 71 72 73 74 75 76 77 78 79 80 81
public class ThreeSum {
    public static int count(int[] nums) {
        int N = nums.length;
        int cnt = 0;
        for (int i = 0; i < N; i++)
            for (int j = i + 1; j < N; j++)
                for (int k = j + 1; k < N; k++)
                    if (nums[i] + nums[j] + nums[k] == 0)
                        cnt++;

        return cnt;
    }
C
CyC2018 已提交
82 83 84
}
```

C
CyC2018 已提交
85
该算法的内循环为 `if (nums[i] + nums[j] + nums[k] == 0)` 语句,总共执行的次数为 N(N-1)(N-2) = N<sup>3</sup>/6-N<sup>2</sup>/2+N/3,因此它的近似执行次数为 \~N<sup>3</sup>/6,增长数量级为 O(N<sup>3</sup>)。
C
CyC2018 已提交
86

C
CyC2018 已提交
87
<font size=4> **改进** </font></br>
C
CyC2018 已提交
88

C
CyC2018 已提交
89
通过将数组先排序,对两个元素求和,并用二分查找方法查找是否存在该和的相反数,如果存在,就说明存在三元组的和为 0。
C
CyC2018 已提交
90

C
CyC2018 已提交
91 92
应该注意的是,只有数组不含有相同元素才能使用这种解法,否则二分查找的结果会出错。

C
CyC2018 已提交
93
该方法可以将 ThreeSum 算法增长数量级降低为 O(N<sup>2</sup>logN)。
C
CyC2018 已提交
94 95

```java
C
CyC2018 已提交
96 97 98 99 100 101 102 103
public class ThreeSumFast {
    public static int count(int[] nums) {
        Arrays.sort(nums);
        int N = nums.length;
        int cnt = 0;
        for (int i = 0; i < N; i++)
            for (int j = i + 1; j < N; j++) {
                int target = -nums[i] - nums[j];
C
CyC2018 已提交
104
                int index = BinarySearch.search(nums, target);
C
CyC2018 已提交
105
                // 应该注意这里的下标必须大于 j,否则会重复统计。
C
CyC2018 已提交
106
                if (index > j)
C
CyC2018 已提交
107 108 109 110
                    cnt++;
            }
        return cnt;
    }
C
CyC2018 已提交
111 112
}
```
C
CyC2018 已提交
113

C
CyC2018 已提交
114 115 116
```java
public class BinarySearch {
    public static int search(int[] nums, int target) {
C
CyC2018 已提交
117 118 119 120 121 122 123 124 125 126 127 128
        int l = 0, h = nums.length - 1;
        while (l <= h) {
            int m = l + (h - l) / 2;
            if (target == nums[m])
                return m;
            else if (target > nums[m])
                l = m + 1;
            else
                h = m - 1;
        }
        return -1;
    }
C
CyC2018 已提交
129 130 131
}
```

C
CyC2018 已提交
132
## 倍率实验
C
CyC2018 已提交
133

C
CyC2018 已提交
134
如果 T(N) \~ aN<sup>b</sup>logN,那么 T(2N)/T(N) \~ 2<sup>b</sup>
C
CyC2018 已提交
135

C
CyC2018 已提交
136
例如对于暴力的 ThreeSum 算法,近似时间为 \~N<sup>3</sup>/6。进行如下实验:多次运行该算法,每次取的 N 值为前一次的两倍,统计每次执行的时间,并统计本次运行时间与前一次运行时间的比值,得到如下结果:
C
CyC2018 已提交
137

C
CyC2018 已提交
138 139 140 141 142 143 144 145
| N | Time(ms) | Ratio |
| :---: | :---: | :---: |
| 500 | 48 | / |
| 1000 | 320 | 6.7 |
| 2000 | 555 | 1.7 |
| 4000 | 4105 | 7.4 |
| 8000 | 33575 | 8.2 |
| 16000 | 268909 | 8.0 |
C
CyC2018 已提交
146

C
CyC2018 已提交
147
可以看到,T(2N)/T(N) \~ 2<sup>3</sup>,因此可以确定 T(N) \~ aN<sup>3</sup>logN。
C
CyC2018 已提交
148

C
CyC2018 已提交
149
```java
C
CyC2018 已提交
150 151 152
public class RatioTest {
    public static void main(String[] args) {
        int N = 500;
C
CyC2018 已提交
153 154 155
        int loopTimes = 7;
        double preTime = -1;
        while (loopTimes-- > 0) {
C
CyC2018 已提交
156
            int[] nums = new int[N];
C
CyC2018 已提交
157
            StopWatch.start();
C
CyC2018 已提交
158
            int cnt = ThreeSum.count(nums);
C
CyC2018 已提交
159 160 161 162 163
            System.out.println(cnt);
            double elapsedTime = StopWatch.elapsedTime();
            double ratio = preTime == -1 ? 0 : elapsedTime / preTime;
            System.out.println(N + "  " + elapsedTime + "  " + ratio);
            preTime = elapsedTime;
C
CyC2018 已提交
164 165 166
            N *= 2;
        }
    }
C
CyC2018 已提交
167 168 169
}
```

C
CyC2018 已提交
170 171 172 173 174 175 176 177 178 179 180 181 182 183 184
```java
public class StopWatch {
    private static long start;
    
    public static void start(){
        start = System.currentTimeMillis();
    }
    
    public static double elapsedTime() {
        long now = System.currentTimeMillis();
        return (now - start) / 1000.0;
    }
}
```

C
CyC2018 已提交
185
## 注意事项
C
CyC2018 已提交
186

C
CyC2018 已提交
187
### 1. 大常数
C
CyC2018 已提交
188 189 190

在求近似时,如果低级项的常数系数很大,那么近似的结果就是错误的。

C
CyC2018 已提交
191
### 2. 缓存
C
CyC2018 已提交
192 193 194

计算机系统会使用缓存技术来组织内存,访问数组相邻的元素会比访问不相邻的元素快很多。

C
CyC2018 已提交
195
### 3. 对最坏情况下的性能的保证
C
CyC2018 已提交
196 197 198

在核反应堆、心脏起搏器或者刹车控制器中的软件,最坏情况下的性能是十分重要的。

C
CyC2018 已提交
199
### 4. 随机化算法
C
CyC2018 已提交
200 201 202

通过打乱输入,去除算法对输入的依赖。

C
CyC2018 已提交
203
### 5. 均摊分析
C
CyC2018 已提交
204

C
CyC2018 已提交
205
将所有操作的总成本除于操作总数来将成本均摊。例如对一个空栈进行 N 次连续的 push() 调用需要访问数组的元素为 N+4+8+16+...+2N=5N-4(N 是向数组写入元素,其余的都是调整数组大小时进行复制需要的访问数组操作),均摊后每次操作访问数组的平均次数为常数。
C
CyC2018 已提交
206

C
CyC2018 已提交
207
# 三、栈和队列
C
CyC2018 已提交
208

C
CyC2018 已提交
209
## 栈
C
CyC2018 已提交
210

C
CyC2018 已提交
211
First-In-Last-Out
C
CyC2018 已提交
212

C
CyC2018 已提交
213 214 215 216 217 218 219 220 221 222 223 224
```java
public interface MyStack<Item> extends Iterable<Item> {
    MyStack<Item> push(Item item);

    Item pop() throws Exception;

    boolean isEmpty();

    int size();
}
```

C
CyC2018 已提交
225
### 1. 数组实现
C
CyC2018 已提交
226 227

```java
C
CyC2018 已提交
228
public class ResizingArrayStack<Item> implements MyStack<Item> {
C
CyC2018 已提交
229 230 231 232 233
    // 栈元素数组
    private Item[] a = (Item[]) new Object[1]; // 只能通过转型来创建泛型数组
    // 元素数量
    private int N = 0;

C
CyC2018 已提交
234 235
    @Override
    public MyStack<Item> push(Item item) {
C
CyC2018 已提交
236 237
        check();
        a[N++] = item;
C
CyC2018 已提交
238
        return this;
C
CyC2018 已提交
239 240
    }

C
CyC2018 已提交
241
    @Override
C
CyC2018 已提交
242 243 244
    public Item pop() throws Exception {
        if (isEmpty())
            throw new Exception("stack is empty");
C
CyC2018 已提交
245

C
CyC2018 已提交
246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268
        Item item = a[--N];
        check();
        a[N] = null; // 避免对象游离
        return item;
    }

    private void check() {
        if (N >= a.length)
            resize(2 * a.length);
        else if (N > 0 && N <= a.length / 4)
            resize(a.length / 2);
    }

    /**
     * 调整数组大小,使得栈具有伸缩性
     */
    private void resize(int size) {
        Item[] tmp = (Item[]) new Object[size];
        for (int i = 0; i < N; i++)
            tmp[i] = a[i];
        a = tmp;
    }

C
CyC2018 已提交
269
    @Override
C
CyC2018 已提交
270 271 272 273
    public boolean isEmpty() {
        return N == 0;
    }

C
CyC2018 已提交
274
    @Override
C
CyC2018 已提交
275 276 277 278 279 280 281
    public int size() {
        return N;
    }

    @Override
    public Iterator<Item> iterator() {
        // 返回逆序遍历的迭代器
C
CyC2018 已提交
282 283
        return new Iterator<Item>() {
            private int i = N;
C
CyC2018 已提交
284

C
CyC2018 已提交
285 286 287 288
            @Override
            public boolean hasNext() {
                return i > 0;
            }
C
CyC2018 已提交
289

C
CyC2018 已提交
290 291 292 293 294
            @Override
            public Item next() {
                return a[--i];
            }
        };
C
CyC2018 已提交
295
    }
C
CyC2018 已提交
296 297 298
}
```

C
CyC2018 已提交
299
### 2. 链表实现
C
CyC2018 已提交
300

C
CyC2018 已提交
301
需要使用链表的头插法来实现,因为头插法中最后压入栈的元素在链表的开头,它的 next 指针指向前一个压入栈的元素,在弹出元素使就可以通过 next 指针遍历到前一个压入栈的元素从而让这个元素称为新的栈顶元素。
C
CyC2018 已提交
302 303

```java
C
CyC2018 已提交
304
public class ListStack<Item> implements MyStack<Item> {
C
CyC2018 已提交
305 306 307 308 309 310 311 312
    private Node top = null;
    private int N = 0;

    private class Node {
        Item item;
        Node next;
    }

C
CyC2018 已提交
313 314
    @Override
    public MyStack<Item> push(Item item) {
C
CyC2018 已提交
315 316 317 318 319
        Node newTop = new Node();
        newTop.item = item;
        newTop.next = top;
        top = newTop;
        N++;
C
CyC2018 已提交
320
        return this;
C
CyC2018 已提交
321 322
    }

C
CyC2018 已提交
323
    @Override
C
CyC2018 已提交
324 325 326 327 328 329 330 331
    public Item pop() throws Exception {
        if (isEmpty())
            throw new Exception("stack is empty");
        Item item = top.item;
        top = top.next;
        N--;
        return item;
    }
C
CyC2018 已提交
332

C
CyC2018 已提交
333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351
    @Override
    public boolean isEmpty() {
        return N == 0;
    }

    @Override
    public int size() {
        return N;
    }

    @Override
    public Iterator<Item> iterator() {
        return new Iterator<Item>() {
            private Node cur = top;

            @Override
            public boolean hasNext() {
                return cur != null;
            }
C
CyC2018 已提交
352

C
CyC2018 已提交
353 354 355 356 357 358 359 360 361
            @Override
            public Item next() {
                Item item = cur.item;
                cur = cur.next;
                return item;
            }
        };
    }
}
C
CyC2018 已提交
362 363
```

C
CyC2018 已提交
364
## 队列
C
CyC2018 已提交
365

C
CyC2018 已提交
366
First-In-First-Out
C
CyC2018 已提交
367

C
CyC2018 已提交
368
下面是队列的链表实现,需要维护 first 和 last 节点指针,分别指向队首和队尾。
C
CyC2018 已提交
369

C
CyC2018 已提交
370
这里需要考虑 first 和 last 指针哪个作为链表的开头。因为出队列操作需要让队首元素的下一个元素成为队首,所以需要容易获取下一个元素,而链表的头部节点的 next 指针指向下一个元素,因此可以让 first 指针链表的开头。
C
CyC2018 已提交
371 372

```java
C
CyC2018 已提交
373 374 375 376 377 378 379 380 381 382 383 384 385
public interface MyQueue<Item> extends Iterable<Item> {
    int size();

    boolean isEmpty();

    MyQueue<Item> add(Item item);

    Item remove() throws Exception;
}
```

```java
public class ListQueue<Item> implements MyQueue<Item> {
C
CyC2018 已提交
386 387 388 389 390 391 392 393 394
    private Node first;
    private Node last;
    int N = 0;

    private class Node {
        Item item;
        Node next;
    }

C
CyC2018 已提交
395
    @Override
C
CyC2018 已提交
396 397 398 399
    public boolean isEmpty() {
        return N == 0;
    }

C
CyC2018 已提交
400
    @Override
C
CyC2018 已提交
401 402 403 404
    public int size() {
        return N;
    }

C
CyC2018 已提交
405 406
    @Override
    public MyQueue<Item> add(Item item) {
C
CyC2018 已提交
407 408 409 410 411 412 413 414 415 416 417
        Node newNode = new Node();
        newNode.item = item;
        newNode.next = null;
        if (isEmpty()) {
            last = newNode;
            first = newNode;
        } else {
            last.next = newNode;
            last = newNode;
        }
        N++;
C
CyC2018 已提交
418
        return this;
C
CyC2018 已提交
419 420
    }

C
CyC2018 已提交
421
    @Override
C
CyC2018 已提交
422 423 424 425 426 427 428 429 430 431
    public Item remove() throws Exception {
        if (isEmpty())
            throw new Exception("queue is empty");
        Node node = first;
        first = first.next;
        N--;
        if (isEmpty())
            last = null;
        return node.item;
    }
C
CyC2018 已提交
432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450

    @Override
    public Iterator<Item> iterator() {
        return new Iterator<Item>() {
            Node cur = first;

            @Override
            public boolean hasNext() {
                return cur != null;
            }

            @Override
            public Item next() {
                Item item = cur.item;
                cur = cur.next;
                return item;
            }
        };
    }
C
CyC2018 已提交
451 452 453
}
```

C
CyC2018 已提交
454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616
# 四、并查集

用于解决动态连通性问题,能动态连接两个点,并且判断两个点是否连通。

<div align="center"> <img src="../pics//9d0a637c-6a8f-4f5a-99b9-fdcfa26793ff.png" width="400"/> </div><br>

| 方法 | 描述 |
| :---: | :---: |
| UF(int N) | 构造一个大小为 N 的并查集 |
| void union(int p, int q) | 连接 p 和 q 节点 |
| int find(int p) | 查找 p 所在的连通分量 |
| boolean connected(int p, int q) | 判断 p 和 q 节点是否连通 |

```java
public abstract class UF {
    protected int[] id;

    public UF(int N) {
        id = new int[N];
        for (int i = 0; i < N; i++)
            id[i] = i;
    }

    public boolean connected(int p, int q) {
        return find(p) == find(q);
    }

    public abstract int find(int p);

    public abstract void union(int p, int q);
}
```

## quick-find

可以快速进行 find 操作,即可以快速判断两个节点是否连通。

同一连通分量的所有节点的 id 值相等。

但是 union 操作代价却很高,需要将其中一个连通分量中的所有节点 id 值都修改为另一个节点的 id 值。

<div align="center"> <img src="../pics//8f0cc500-5994-4c7a-91a9-62885d658662.png" width="350"/> </div><br>

```java
public class QuickFindUF extends UF {
    public QuickFindUF(int N) {
        super(N);
    }

    @Override
    public int find(int p) {
        return id[p];
    }

    @Override
    public void union(int p, int q) {
        int pID = find(p);
        int qID = find(q);

        if (pID == qID)
            return;

        for (int i = 0; i < id.length; i++)
            if (id[i] == pID)
                id[i] = qID;
    }
}

```

## quick-union

可以快速进行 union 操作,只需要修改一个节点的 id 值即可。

但是 find 操作开销很大,因为同一个连通分量的节点 id 值不同,id 值只是用来指向另一个节点。因此需要一直向上查找操作,直到找到最上层的节点。

<div align="center"> <img src="../pics//5d4a5181-65fb-4bf2-a9c6-899cab534b44.png" width="350"/> </div><br>

```java
public class QuickUnionUF extends UF {
    public QuickUnionUF(int N) {
        super(N);
    }

    @Override
    public int find(int p) {
        while (p != id[p])
            p = id[p];
        return p;
    }

    @Override
    public void union(int p, int q) {
        int pRoot = find(p);
        int qRoot = find(q);
        if (pRoot != qRoot)
            id[pRoot] = qRoot;
    }
}
```

这种方法可以快速进行 union 操作,但是 find 操作和树高成正比,最坏的情况下树的高度为触点的数目。

<div align="center"> <img src="../pics//bfbb11e2-d208-4efa-b97b-24cd40467cd8.png" width="150"/> </div><br>

## 加权 quick-union

为了解决 quick-union 的树通常会很高的问题,加权 quick-union 在 union 操作时会让较小的树连接较大的树上面。

理论研究证明,加权 quick-union 算法构造的树深度最多不超过 logN。

<div align="center"> <img src="../pics//a4c17d43-fa5e-4935-b74e-147e7f7e782c.png" width="200"/> </div><br>

```java
public class WeightedQuickUnionUF extends UF {

    // 保存节点的数量信息
    private int[] sz;

    public WeightedQuickUnionUF(int N) {
        super(N);
        this.sz = new int[N];
        for (int i = 0; i < N; i++)
            this.sz[i] = 1;
    }

    @Override
    public int find(int p) {
        while (p != id[p])
            p = id[p];
        return p;
    }

    @Override
    public void union(int p, int q) {
        int i = find(p);
        int j = find(q);
        if (i == j) return;
        if (sz[i] < sz[j]) {
            id[i] = j;
            sz[j] += sz[i];
        } else {
            id[j] = i;
            sz[i] += sz[j];
        }
    }
}
```

## 路径压缩的加权 quick-union

在检查节点的同时将它们直接链接到根节点,只需要在 find 中添加一个循环即可。

## 各种 union-find 算法的比较

| 算法 | union | find |
| :---: | :---: | :---: |
| quick-find | N | 1 |
| quick-union | 树高 | 树高 |
| 加权 quick-union | logN | logN |
| 路径压缩的加权 quick-union | 非常接近 1 | 非常接近 1 |

# 五、排序
C
CyC2018 已提交
617

C
CyC2018 已提交
618
待排序的元素需要实现 Java 的 Comparable 接口,该接口有 compareTo() 方法,可以用它来判断两个元素的大小关系。
C
CyC2018 已提交
619 620 621

研究排序算法的成本模型时,计算的是比较和交换的次数。

C
CyC2018 已提交
622
使用辅助函数 less() 和 swap() 来进行比较和交换的操作,使得代码的可读性和可移植性更好。
C
CyC2018 已提交
623 624

```java
C
CyC2018 已提交
625 626
private static boolean less(Comparable v, Comparable w) {
    return v.compareTo(w) < 0;
C
CyC2018 已提交
627 628
}

C
CyC2018 已提交
629 630 631 632
private static void swap(Comparable[] a, int i, int j) {
    Comparable t = a[i];
    a[i] = a[j];
    a[j] = t;
C
CyC2018 已提交
633 634 635
}
```

C
CyC2018 已提交
636
## 选择排序
C
CyC2018 已提交
637

C
CyC2018 已提交
638
选择出数组中的最小元素,将它与数组的第一个元素交换位置。再从剩下的元素中选择出最小的元素,将它与数组的第二个元素交换位置。不断进行这样的操作,直到将整个数组排序。
C
CyC2018 已提交
639

C
CyC2018 已提交
640
<div align="center"> <img src="../pics//ed7b96ac-6428-4bd5-9986-674c54c2a959.png" width="250"/> </div><br>
C
CyC2018 已提交
641 642

```java
C
CyC2018 已提交
643 644 645 646 647 648 649 650 651 652 653
public class Selection {
    public static void sort(Comparable[] a) {
        int N = a.length;
        for (int i = 0; i < N; i++) {
            int min = i;
            for (int j = i + 1; j < N; j++)
                if (less(a[j], a[min]))
                    min = j;
            swap(a, i, min);
        }
    }
C
CyC2018 已提交
654 655 656
}
```

C
CyC2018 已提交
657
选择排序需要 \~N<sup>2</sup>/2 次比较和 \~N 次交换,它的运行时间与输入无关,这个特点使得它对一个已经排序的数组也需要这么多的比较和交换操作。
C
CyC2018 已提交
658

C
CyC2018 已提交
659
## 冒泡排序
C
CyC2018 已提交
660

C
CyC2018 已提交
661
通过从左到右不断交换相邻逆序的相邻元素,在一轮的交换之后,可以让未排序的元素上浮到最右侧,是的右侧是已排序的。
C
CyC2018 已提交
662

C
CyC2018 已提交
663
在一轮交换中,如果没有发生交换,就说明数组已经是有序的,此时可以直接退出。
C
CyC2018 已提交
664 665

```java
C
CyC2018 已提交
666 667 668 669 670 671 672 673 674 675 676 677 678 679
public class Bubble {
    public static void sort(Comparable[] a) {
        int N = a.length;
        boolean hasSorted = false;
        for (int i = 0; i < N && !hasSorted; i++) {
            hasSorted = true;
            for (int j = 0; j < N - i - 1; j++) {
                if (less(a[j + 1], a[j])) {
                    hasSorted = false;
                    swap(a, j, j + 1);
                }
            }
        }
    }
C
CyC2018 已提交
680 681 682
}
```

C
CyC2018 已提交
683
## 插入排序
C
CyC2018 已提交
684 685 686

插入排序从左到右进行,每次都将当前元素插入到左侧已经排序的数组中,使得插入之后左部数组依然有序。

C
CyC2018 已提交
687
第 j 元素是通过不断向左比较并交换来实现插入过程:当第 j 元素小于第 j - 1 元素,就将它们的位置交换,然后令 j 指针向左移动一个位置,不断进行以上操作。
C
CyC2018 已提交
688

C
CyC2018 已提交
689
<div align="center"> <img src="../pics//c9a1de44-b1c0-4d13-a654-827d4ef8a723.png" width="250"/> </div><br>
C
CyC2018 已提交
690 691

```java
C
CyC2018 已提交
692 693 694 695 696 697 698
public class Insertion {
    public static void sort(Comparable[] a) {
        int N = a.length;
        for (int i = 1; i < N; i++)
            for (int j = i; j > 0 && less(a[j], a[j - 1]); j--)
                swap(a, j, j - 1);
    }
C
CyC2018 已提交
699 700 701
}
```

C
CyC2018 已提交
702
对于数组 {3, 5, 2, 4, 1},它具有以下逆序:(3, 2), (3, 1), (5, 2), (5, 4), (5, 1), (2, 1), (4, 1),插入排序每次只能交换相邻元素,令逆序数量减少 1,因此插入排序需要交换的次数为逆序数量。
C
CyC2018 已提交
703 704

插入排序的复杂度取决于数组的初始顺序,如果数组已经部分有序了,逆序较少,那么插入排序会很快。
C
CyC2018 已提交
705

C
CyC2018 已提交
706 707 708
- 平均情况下插入排序需要 \~N<sup>2</sup>/4 比较以及 \~N<sup>2</sup>/4 次交换;
- 最坏的情况下需要 \~N<sup>2</sup>/2 比较以及 \~N<sup>2</sup>/2 次交换,最坏的情况是数组是倒序的;
- 最好的情况下需要 N-1 次比较和 0 次交换,最好的情况就是数组已经有序了。
C
CyC2018 已提交
709

C
CyC2018 已提交
710
## 希尔排序
C
CyC2018 已提交
711

C
CyC2018 已提交
712
对于大规模的数组,插入排序很慢,因为它只能交换相邻的元素,每次只能将逆序数量减少 1。
C
CyC2018 已提交
713

C
CyC2018 已提交
714
希尔排序的出现就是为了改进插入排序的这种局限性,它通过交换不相邻的元素,每次可以将逆序数量减少大于 1。
C
CyC2018 已提交
715

C
CyC2018 已提交
716
希尔排序使用插入排序对间隔 h 的序列进行排序。通过不断减小 h,最后令 h=1,就可以使得整个数组是有序的。
C
CyC2018 已提交
717

C
CyC2018 已提交
718
<div align="center"> <img src="../pics//cdbe1d12-5ad9-4acb-a717-bbc822c2acf3.png" width="500"/> </div><br>
C
CyC2018 已提交
719 720

```java
C
CyC2018 已提交
721 722 723 724 725 726 727 728 729 730 731 732 733 734
public class Shell {
    public static void sort(Comparable[] a) {
        int N = a.length;
        int h = 1;
        while (h < N / 3)
            h = 3 * h + 1; // 1, 4, 13, 40, ...

        while (h >= 1) {
            for (int i = h; i < N; i++)
                for (int j = i; j >= h && less(a[j], a[j - h]); j -= h)
                    swap(a, j, j - h);
            h = h / 3;
        }
    }
C
CyC2018 已提交
735 736 737
}
```

C
CyC2018 已提交
738
希尔排序的运行时间达不到平方级别,使用递增序列 1, 4, 13, 40, ...  的希尔排序所需要的比较次数不会超过 N 的若干倍乘于递增序列的长度。后面介绍的高级排序算法只会比希尔排序快两倍左右。
C
CyC2018 已提交
739

C
CyC2018 已提交
740
## 归并排序
C
CyC2018 已提交
741 742 743

归并排序的思想是将数组分成两部分,分别进行排序,然后归并起来。

C
CyC2018 已提交
744
<div align="center"> <img src="../pics//8dfb4cc9-26da-45e7-b820-4576fa1cbb0e.png" width="350"/> </div><br>
C
CyC2018 已提交
745

C
CyC2018 已提交
746
### 1. 归并方法
C
CyC2018 已提交
747 748 749 750

归并方法将数组中两个已经排序的部分归并成一个。

```java
C
CyC2018 已提交
751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770
public class MergeSort {
    private static Comparable[] aux;

    private static void merge(Comparable[] a, int l, int m, int h) {
        int i = l, j = m + 1;

        for (int k = l; k <= h; k++)
            aux[k] = a[k]; // 将数据复制到辅助数组

        for (int k = l; k <= h; k++) {
            if (i > m)
                a[k] = aux[j++];
            else if (j > h)
                a[k] = aux[i++];
            else if (aux[i].compareTo(a[j]) <= 0)
                a[k] = aux[i++]; // 先进行这一步,保证稳定性
            else
                a[k] = aux[j++];
        }
    }
C
CyC2018 已提交
771 772 773
}
```

C
CyC2018 已提交
774
### 2. 自顶向下归并排序
C
CyC2018 已提交
775

C
CyC2018 已提交
776
<div align="center"> <img src="../pics//0c55e11c-d3ce-4cd8-b139-028aea6f40e3.png" width="450"/> </div><br>
C
CyC2018 已提交
777 778 779


```java
C
CyC2018 已提交
780 781 782
public static void sort(Comparable[] a) {
    aux = new Comparable[a.length];
    sort(a, 0, a.length - 1);
C
CyC2018 已提交
783 784
}

C
CyC2018 已提交
785 786 787 788 789 790 791
private static void sort(Comparable[] a, int l, int h) {
    if (h <= l)
        return;
    int mid = l + (h - l) / 2;
    sort(a, l, mid);
    sort(a, mid + 1, h);
    merge(a, l, mid, h);
C
CyC2018 已提交
792 793 794
}
```

C
CyC2018 已提交
795
因为每次都将问题对半分成两个子问题,而这种对半分的算法复杂度一般为 O(NlogN),因此该归并排序方法的时间复杂度也为 O(NlogN)。
C
CyC2018 已提交
796

C
CyC2018 已提交
797
### 3. 自底向上归并排序
C
CyC2018 已提交
798

C
CyC2018 已提交
799
先归并那些微型数组,然后成对归并得到的微型数组。
C
CyC2018 已提交
800 801

```java
C
CyC2018 已提交
802 803 804 805 806 807 808 809 810
 public static void sort(Comparable[] a) {
     int N = a.length;
     aux = new Comparable[N];
     for (int sz = 1; sz < N; sz += sz) {
         for (int lo = 0; lo < N - sz; lo += sz + sz) {
             merge(a, lo, lo + sz - 1, Math.min(lo + sz + sz - 1, N - 1));
         }
     }
 }
C
CyC2018 已提交
811 812
```

C
CyC2018 已提交
813
## 快速排序
C
CyC2018 已提交
814

C
CyC2018 已提交
815
### 1. 基本算法
C
CyC2018 已提交
816

C
CyC2018 已提交
817 818
- 归并排序将数组分为两个子数组分别排序,并将有序的子数组归并使得整个数组排序;
- 快速排序通过一个切分元素将数组分为两个子数组,左子数组小于等于切分元素,右子数组大于等于切分元素,将这两个子数组排序也就将整个数组排序了。
C
CyC2018 已提交
819

C
CyC2018 已提交
820
<div align="center"> <img src="../pics//ab77240d-7338-4547-9183-00215e7220ec.png" width="500"/> </div><br>
C
CyC2018 已提交
821 822

```java
C
CyC2018 已提交
823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841
public class QuickSort {
    public static void sort(Comparable[] a) {
        shuffle(a);
        sort(a, 0, a.length - 1);
    }

    private static void sort(Comparable[] a, int l, int h) {
        if (h <= l)
            return;
        int j = partition(a, l, h);
        sort(a, l, j - 1);
        sort(a, j + 1, h);
    }

    private static void shuffle(Comparable[] array) {
        List<Comparable> list = Arrays.asList(array);
        Collections.shuffle(list);
        list.toArray(array);
    }
C
CyC2018 已提交
842 843 844
}
```

C
CyC2018 已提交
845
### 2. 切分
C
CyC2018 已提交
846

C
CyC2018 已提交
847
取 a[lo] 作为切分元素,然后从数组的左端向右扫描直到找到第一个大于等于它的元素,再从数组的右端向左扫描找到第一个小于等于它的元素,交换这两个元素,并不断进行这个过程,就可以保证左指针 i 的左侧元素都不大于切分元素,右指针 j 的右侧元素都不小于切分元素。当两个指针相遇时,将切分元素 a[lo] 和 a[j] 交换位置。
C
CyC2018 已提交
848

C
CyC2018 已提交
849
<div align="center"> <img src="../pics//5aac64d3-2c7b-4f32-9e9a-1df2186f588b.png" width="300"/> </div><br>
C
CyC2018 已提交
850 851

```java
C
CyC2018 已提交
852 853 854 855 856 857 858 859 860 861 862 863
private static int partition(Comparable[] a, int l, int h) {
    int i = l, j = h + 1;
    Comparable v = a[l];
    while (true) {
        while (less(a[++i], v) && i != h) ;
        while (less(v, a[--j]) && j != l) ;
        if (i >= j)
            break;
        swap(a, i, j);
    }
    swap(a, l, j);
    return j;
C
CyC2018 已提交
864 865 866
}
```

C
CyC2018 已提交
867
### 3. 性能分析
C
CyC2018 已提交
868 869 870

快速排序是原地排序,不需要辅助数组,但是递归调用需要辅助栈。

C
CyC2018 已提交
871
快速排序最好的情况下是每次都正好能将数组对半分,这样递归调用次数才是最少的。这种情况下比较次数为 C<sub>N</sub>=2C<sub>N/2</sub>+N,复杂度为 O(NlogN)。
C
CyC2018 已提交
872

C
CyC2018 已提交
873
最坏的情况下,第一次从最小的元素切分,第二次从第二小的元素切分,如此这般。因此最坏的情况下需要比较 N<sup>2</sup>/2。为了防止数组最开始就是有序的,在进行快速排序时需要随机打乱数组。
C
CyC2018 已提交
874

C
CyC2018 已提交
875
### 4. 算法改进
C
CyC2018 已提交
876 877 878

(一)切换到插入排序

C
CyC2018 已提交
879
因为快速排序在小数组中也会递归调用自己,对于小数组,插入排序比快速排序的性能更好,因此在小数组中可以切换到插入排序。
C
CyC2018 已提交
880

C
CyC2018 已提交
881
(二)三数取中
C
CyC2018 已提交
882

C
CyC2018 已提交
883
最好的情况下是每次都能取数组的中位数作为切分元素,但是计算中位数的代价很高。人们发现取 3 个元素并将大小居中的元素作为切分元素的效果最好。
C
CyC2018 已提交
884 885 886 887 888 889 890 891

(三)三向切分

对于有大量重复元素的数组,可以将数组切分为三部分,分别对应小于、等于和大于切分元素。

三向切分快速排序对于只有若干不同主键的随机数组可以在线性时间内完成排序。

```java
C
CyC2018 已提交
892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913
public class Quick3Way {
    public static void sort(Comparable[] a) {
        sort(a, 0, a.length - 1);
    }

    private static void sort(Comparable[] a, int l, int h) {
        if (h <= l)
            return;
        int lt = l, i = l + 1, gt = h;
        Comparable v = a[l];
        while (i <= gt) {
            int cmp = a[i].compareTo(v);
            if (cmp < 0)
                swap(a, lt++, i++);
            else if (cmp > 0)
                swap(a, i, gt--);
            else
                i++;
        }
        sort(a, l, lt - 1);
        sort(a, gt + 1, h);
    }
C
CyC2018 已提交
914 915 916
}
```

C
CyC2018 已提交
917
### 5. 基于切分的快速选择算法
C
CyC2018 已提交
918

C
CyC2018 已提交
919
快速排序的 partition() 方法,会返回一个整数 j 使得 a[l..j-1] 小于等于 a[j],且 a[j+1..h] 大于等于 a[j],此时 a[j] 就是数组的第 j 大元素。
C
CyC2018 已提交
920

C
CyC2018 已提交
921
可以利用这个特性找出数组的第 k 个元素。
C
CyC2018 已提交
922 923

```java
C
CyC2018 已提交
924 925 926 927 928 929 930 931 932 933 934 935
public static Comparable select(Comparable[] a, int k) {
    int l = 0, h = a.length - 1;
    while (h > l) {
        int j = partition(a, l, h);
        if (j == k)
            return a[k];
        else if (j > k)
            h = j - 1;
        else
            l = j + 1;
    }
    return a[k];
C
CyC2018 已提交
936 937 938
}
```

C
CyC2018 已提交
939
该算法是线性级别的,因为每次正好将数组二分,那么比较的总次数为 (N+N/2+N/4+..),直到找到第 k 个元素,这个和显然小于 2N。
C
CyC2018 已提交
940

C
CyC2018 已提交
941
## 堆排序
C
CyC2018 已提交
942

C
CyC2018 已提交
943
### 1. 堆
C
CyC2018 已提交
944 945 946

堆的某个节点的值总是大于等于子节点的值,并且堆是一颗完全二叉树。

C
CyC2018 已提交
947
堆可以用数组来表示,因为堆是一种完全二叉树,而完全二叉树很容易就存储在数组中。位置 k 的节点的父节点位置为 k/2,而它的两个子节点的位置分别为 2k 和 2k+1。这里不使用数组索引为 0 的位置,是为了更清晰地描述节点的位置关系。
C
CyC2018 已提交
948

C
CyC2018 已提交
949
<div align="center"> <img src="../pics//f3080f83-6239-459b-8e9c-03b6641f7815.png" width="200"/> </div><br>
C
CyC2018 已提交
950 951

```java
C
CyC2018 已提交
952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977
public class Heap {
    private Comparable[] heap;
    private int N = 0;

    public Heap(int maxN) {
        heap = new Comparable[maxN + 1];
        N = maxN;
    }

    public boolean isEmpty() {
        return N == 0;
    }

    public int size() {
        return N;
    }

    private boolean less(int i, int j) {
        return heap[i].compareTo(heap[j]) < 0;
    }

    private void swap(int i, int j) {
        Comparable t = heap[i];
        heap[i] = heap[j];
        heap[j] = t;
    }
C
CyC2018 已提交
978 979 980
}
```

C
CyC2018 已提交
981
### 2. 上浮和下沉
C
CyC2018 已提交
982

C
CyC2018 已提交
983
在堆中,当一个节点比父节点大,那么需要交换这个两个节点。交换后还可能比它新的父节点大,因此需要不断地进行比较和交换操作,把这种操作称为上浮。
C
CyC2018 已提交
984

C
CyC2018 已提交
985
<div align="center"> <img src="../pics//33ac2b23-cb85-4e99-bc41-b7b7199fad1c.png" width="400"/> </div><br>
C
CyC2018 已提交
986 987

```java
C
CyC2018 已提交
988 989 990 991 992
private void swim(int k) {
    while (k > 1 && less(k / 2, k)) {
        swap(k / 2, k);
        k = k / 2;
    }
C
CyC2018 已提交
993 994 995
}
```

C
CyC2018 已提交
996
类似地,当一个节点比子节点来得小,也需要不断地向下进行比较和交换操作,把这种操作称为下沉。一个节点有两个子节点,应当与两个子节点中最大那么节点进行交换。
C
CyC2018 已提交
997

C
CyC2018 已提交
998
<div align="center"> <img src="../pics//72f0ff69-138d-4e54-b7ac-ebe025d978dc.png" width="400"/> </div><br>
C
CyC2018 已提交
999 1000

```java
C
CyC2018 已提交
1001 1002 1003 1004 1005 1006 1007 1008 1009 1010
private void sink(int k) {
    while (2 * k <= N) {
        int j = 2 * k;
        if (j < N && less(j, j + 1))
            j++;
        if (!less(k, j))
            break;
        swap(k, j);
        k = j;
    }
C
CyC2018 已提交
1011 1012 1013
}
```

C
CyC2018 已提交
1014
### 3. 插入元素
C
CyC2018 已提交
1015 1016 1017 1018

将新元素放到数组末尾,然后上浮到合适的位置。

```java
C
CyC2018 已提交
1019 1020 1021
public void insert(Comparable v) {
    heap[++N] = v;
    swim(N);
C
CyC2018 已提交
1022 1023 1024
}
```

C
CyC2018 已提交
1025
### 4. 删除最大元素
C
CyC2018 已提交
1026 1027 1028 1029

从数组顶端删除最大的元素,并将数组的最后一个元素放到顶端,并让这个元素下沉到合适的位置。

```java
C
CyC2018 已提交
1030 1031 1032 1033 1034 1035
public Comparable delMax() {
    Comparable max = heap[1];
    swap(1, N--);
    heap[N + 1] = null;
    sink(1);
    return max;
C
CyC2018 已提交
1036 1037 1038
}
```

C
CyC2018 已提交
1039
### 5. 堆排序
C
CyC2018 已提交
1040 1041 1042 1043 1044 1045 1046

由于堆可以很容易得到最大的元素并删除它,不断地进行这种操作可以得到一个递减序列。如果把最大元素和当前堆中数组的最后一个元素交换位置,并且不删除它,那么就可以得到一个从尾到头的递减序列,从正向来看就是一个递增序列。因此很容易使用堆来进行排序,并且堆排序是原地排序,不占用额外空间。

(一)构建堆

无序数组建立堆最直接的方法是从左到右遍历数组,然后进行上浮操作。一个更高效的方法是从右至左进行下沉操作,如果一个节点的两个节点都已经是堆有序,那么进行下沉操作可以使得这个节点为根节点的堆有序。叶子节点不需要进行下沉操作,因此可以忽略叶子节点的元素,因此只需要遍历一半的元素即可。

C
CyC2018 已提交
1047
<div align="center"> <img src="../pics//b84ba6fb-312b-4e69-8c77-fb6eb6fb38d4.png" width="300"/> </div><br>
C
CyC2018 已提交
1048 1049 1050 1051 1052

(二)交换堆顶元素与最后一个元素

交换之后需要进行下沉操作维持堆的有序状态。

C
CyC2018 已提交
1053
<div align="center"> <img src="../pics//51fb761d-8ce0-4472-92ff-2f227ac7888a.png" width="270"/> </div><br>
C
CyC2018 已提交
1054

C
CyC2018 已提交
1055
<div align="center"> <img src="../pics//1f039a45-6b91-4f31-a2c2-6c63eb8bdb56.png" width="300"/> </div><br>
C
CyC2018 已提交
1056 1057

```java
C
CyC2018 已提交
1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081
public class HeapSort {
    public static void sort(Comparable[] a) { // 数组第 0 个位置不能有元素
        int N = a.length - 1;
        for (int k = N / 2; k >= 0; k--)
            sink(a, k, N);

        while (N > 1) {
            swap(a, 1, N--);
            sink(a, 1, N);
        }
    }


    private static void sink(Comparable[] a, int k, int N) {
        while (2 * k <= N) {
            int j = 2 * k;
            if (j < N && less(a, j, j + 1))
                j++;
            if (!less(a, k, j))
                break;
            swap(a, k, j);
            k = j;
        }
    }
C
CyC2018 已提交
1082 1083 1084
}
```

C
CyC2018 已提交
1085
### 6. 分析
C
CyC2018 已提交
1086

C
CyC2018 已提交
1087
一个堆的高度为 logN,因此在堆中插入元素和删除最大元素的复杂度都为 logN。
C
CyC2018 已提交
1088

C
CyC2018 已提交
1089
对于堆排序,由于要对 N 个节点进行下沉操作,因此复杂度为 NlogN。
C
CyC2018 已提交
1090 1091 1092 1093 1094

堆排序时一种原地排序,没有利用额外的空间。

现代操作系统很少使用堆排序,因为它无法利用缓存,也就是数组元素很少和相邻的元素进行比较。

C
CyC2018 已提交
1095
## 桶排序
C
CyC2018 已提交
1096

C
CyC2018 已提交
1097
## 基数排序
C
CyC2018 已提交
1098

C
CyC2018 已提交
1099
## 外部排序
C
CyC2018 已提交
1100

C
CyC2018 已提交
1101
## 排序算法的比较
C
CyC2018 已提交
1102

C
CyC2018 已提交
1103 1104 1105 1106 1107 1108 1109 1110 1111
| 算法 | 稳定 | 原地排序 | 时间复杂度 | 空间复杂度 | 备注 |
| :---: | :---: | :---: | :---: | :---: | :---: |
| 选择排序 | no | yes | N<sup>2</sup> | 1 | |
| 插入排序 | yes | yes | N \~ N<sup>2</sup> | 1 | 时间复杂度和初始顺序有关 |
| 希尔排序 | no | yes | N 的若干倍乘于递增序列的长度 | 1 | |
| 快速排序 | no | yes | NlogN | logN | |
| 三向切分快速排序 | no | yes | N \~ NlogN | logN | 适用于有大量重复主键|
| 归并排序 | yes | no | NlogN | N | |
| 堆排序 | no | yes | NlogN | 1 | | |
C
CyC2018 已提交
1112

C
CyC2018 已提交
1113
快速排序是最快的通用排序算法,它的内循环的指令很少,而且它还能利用缓存,因为它总是顺序地访问数据。它的运行时间近似为 \~cNlogN,这里的 c 比其他线性对数级别的排序算法都要小。使用三向切分快速排序,实际应用中可能出现的某些分布的输入能够达到线性级别,而其它排序算法仍然需要线性对数时间。
C
CyC2018 已提交
1114

C
CyC2018 已提交
1115
## Java 的排序算法实现
C
CyC2018 已提交
1116

C
CyC2018 已提交
1117
Java 主要排序方法为 java.util.Arrays.sort(),对于原始数据类型使用三向切分的快速排序,对于引用类型使用归并排序。
C
CyC2018 已提交
1118

C
CyC2018 已提交
1119
# 六、查找
C
CyC2018 已提交
1120 1121 1122

符号表是一种存储键值对的数据结构,主要支持两种操作:插入一个新的键值对、根据给定键得到值。

C
CyC2018 已提交
1123
符号表分为有序和无序两种,有序符号表主要指支持 min()、max() 等根据键的大小关系来实现的操作。
C
CyC2018 已提交
1124

C
CyC2018 已提交
1125
有序符号表的键需要实现 Comparable 接口。
C
CyC2018 已提交
1126

C
CyC2018 已提交
1127
## 二分查找实现有序符号表
C
CyC2018 已提交
1128 1129 1130

使用一对平行数组,一个存储键一个存储值。

C
CyC2018 已提交
1131
rank() 方法至关重要,当键在表中时,它能够知道该键的位置;当键不在表中时,它也能知道在何处插入新键。
C
CyC2018 已提交
1132

C
CyC2018 已提交
1133
复杂度:二分查找最多需要 logN+1 次比较,使用二分查找实现的符号表的查找操作所需要的时间最多是对数级别的。但是插入操作需要移动数组元素,是线性级别的。
C
CyC2018 已提交
1134 1135

```java
C
CyC2018 已提交
1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188
public class BinarySearchST<Key extends Comparable<Key>, Value> {
    private Key[] keys;
    private Value[] values;
    private int N;

    public BinarySearchST(int capacity) {
        keys = (Key[]) new Comparable[capacity];
        values = (Value[]) new Object[capacity];
    }

    public int size() {
        return N;
    }

    public Value get(Key key) {
        int i = rank(key);
        if (i < N && keys[i].compareTo(key) == 0) {
            return values[i];
        }
        return null;
    }

    public int rank(Key key) {
        int lo = 0, hi = N - 1;
        while (lo <= hi) {
            int mid = lo + (hi - lo) / 2;
            int cmp = key.compareTo(keys[mid]);
            if (cmp == 0) return mid;
            else if (cmp < 0) hi = mid - 1;
            else lo = mid + 1;
        }
        return lo;
    }

    public void put(Key key, Value value) {
        int i = rank(key);
        if (i < N && keys[i].compareTo(key) == 0) {
            values[i] = value;
            return;
        }
        for (int j = N; j > i; j--) {
            keys[j] = keys[j - 1];
            values[j] = values[j - 1];
        }
        keys[i] = key;
        values[i] = value;
        N++;
    }

    public Key ceiling(Key key){
        int i = rank(key);
        return keys[i];
    }
C
CyC2018 已提交
1189 1190 1191
}
```

C
CyC2018 已提交
1192
## 二叉查找树
C
CyC2018 已提交
1193

C
CyC2018 已提交
1194
**二叉树**  是一个空链接,或者是一个有左右两个链接的节点,每个链接都指向一颗子二叉树。
C
CyC2018 已提交
1195

C
CyC2018 已提交
1196
<div align="center"> <img src="../pics//f9f9f993-8ece-4da7-b848-af9b438fad76.png" width="200"/> </div><br>
C
CyC2018 已提交
1197

C
CyC2018 已提交
1198
**二叉查找树** (BST)是一颗二叉树,并且每个节点的值都大于等于其左子树中的所有节点的值而小于等于右子树的所有节点的值。
C
CyC2018 已提交
1199

C
CyC2018 已提交
1200
<div align="center"> <img src="../pics//8ae4550b-f0cb-4e4d-9e2b-c550538bf230.png" width="200"/> </div><br>
C
CyC2018 已提交
1201

C
CyC2018 已提交
1202
BST 有一个重要性质,就是它的中序遍历结果递增排序。
C
CyC2018 已提交
1203

C
CyC2018 已提交
1204
<div align="center"> <img src="../pics//fbe54203-c005-48f0-8883-b05e564a3173.png" width="200"/> </div><br>
C
CyC2018 已提交
1205 1206 1207 1208

基本数据结构:

```java
C
CyC2018 已提交
1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233
public class BST<Key extends Comparable<Key>, Value> {
    private Node root;

    private class Node {
        private Key key;
        private Value val;
        private Node left, right;
        // 以该节点为根的子树中节点总数
        private int N;

        public Node(Key key, Value val, int N) {
            this.key = key;
            this.val = val;
            this.N = N;
        }
    }

    public int size() {
        return size(root);
    }

    private int size(Node x) {
        if (x == null) return 0;
        return x.N;
    }
C
CyC2018 已提交
1234 1235 1236 1237 1238
}
```

(为了方便绘图,二叉树的空链接不画出来。)

C
CyC2018 已提交
1239
### 1. get()
C
CyC2018 已提交
1240

C
CyC2018 已提交
1241 1242 1243
- 如果树是空的,则查找未命中;
- 如果被查找的键和根节点的键相等,查找命中;
- 否则递归地在子树中查找:如果被查找的键较小就在左子树中查找,较大就在右子树中查找。
C
CyC2018 已提交
1244 1245

```java
C
CyC2018 已提交
1246 1247
public Value get(Key key) {
    return get(root, key);
C
CyC2018 已提交
1248
}
C
CyC2018 已提交
1249 1250 1251 1252 1253 1254
private Value get(Node x, Key key) {
    if (x == null) return null;
    int cmp = key.compareTo(x.key);
    if (cmp == 0) return x.val;
    else if (cmp < 0) return get(x.left, key);
    else return get(x.right, key);
C
CyC2018 已提交
1255 1256 1257
}
```

C
CyC2018 已提交
1258
### 2. put()
C
CyC2018 已提交
1259 1260 1261

当插入的键不存在于树中,需要创建一个新节点,并且更新上层节点的链接使得该节点正确链接到树中。

C
CyC2018 已提交
1262
<div align="center"> <img src="../pics//107a6a2b-f15b-4cad-bced-b7fb95258c9c.png" width="200"/> </div><br>
C
CyC2018 已提交
1263 1264

```java
C
CyC2018 已提交
1265 1266
public void put(Key key, Value val) {
    root = put(root, key, val);
C
CyC2018 已提交
1267
}
C
CyC2018 已提交
1268 1269 1270 1271 1272 1273 1274 1275
private Node put(Node x, Key key, Value val) {
    if (x == null) return new Node(key, val, 1);
    int cmp = key.compareTo(x.key);
    if (cmp == 0) x.val = val;
    else if (cmp < 0) x.left = put(x.left, key, val);
    else x.right = put(x.right, key, val);
    x.N = size(x.left) + size(x.right) + 1;
    return x;
C
CyC2018 已提交
1276 1277 1278
}
```

C
CyC2018 已提交
1279
### 3. 分析
C
CyC2018 已提交
1280

C
CyC2018 已提交
1281
二叉查找树的算法运行时间取决于树的形状,而树的形状又取决于键被插入的先后顺序。最好的情况下树是完全平衡的,每条空链接和根节点的距离都为 logN。
C
CyC2018 已提交
1282

C
CyC2018 已提交
1283
<div align="center"> <img src="../pics//4d741402-344d-4b7c-be01-e57184bcad0e.png" width="200"/> </div><br>
C
CyC2018 已提交
1284

C
CyC2018 已提交
1285
在最坏的情况下,树的高度为 N。
C
CyC2018 已提交
1286

C
CyC2018 已提交
1287
<div align="center"> <img src="../pics//be7dca03-12ec-456b-8b54-b1b3161f5531.png" width="200"/> </div><br>
C
CyC2018 已提交
1288

C
CyC2018 已提交
1289
### 4. floor()
C
CyC2018 已提交
1290 1291 1292

floor(key):小于等于键的最大键

C
CyC2018 已提交
1293 1294
- 如果键小于根节点的键,那么 floor(key) 一定在左子树中;
- 如果键大于根节点的键,需要先判断右子树中是否存在 floor(key),如果存在就找到,否则根节点就是 floor(key)。
C
CyC2018 已提交
1295 1296 1297


```java
C
CyC2018 已提交
1298 1299 1300 1301
public Key floor(Key key) {
    Node x = floor(root, key);
    if (x == null) return null;
    return x.key;
C
CyC2018 已提交
1302
}
C
CyC2018 已提交
1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313
private Node floor(Node x, Key key) {
    if (x == null) return null;
    int cmp = key.compareTo(x.key);
    if (cmp == 0) return x;
    if (cmp < 0) return floor(x.left, key);
    Node t = floor(x.right, key);
    if (t != null) {
        return t;
    } else {
        return x;
    }
C
CyC2018 已提交
1314 1315 1316
}
```

C
CyC2018 已提交
1317
### 5. rank()
C
CyC2018 已提交
1318

C
CyC2018 已提交
1319
rank(key) 返回 key 的排名。
C
CyC2018 已提交
1320

C
CyC2018 已提交
1321 1322 1323
- 如果键和根节点的键相等,返回左子树的节点数;
- 如果小于,递归计算在左子树中的排名;
- 如果大于,递归计算在右子树中的排名,并加上左子树的节点数,再加上 1(根节点)。
C
CyC2018 已提交
1324 1325

```java
C
CyC2018 已提交
1326 1327
public int rank(Key key) {
    return rank(key, root);
C
CyC2018 已提交
1328
}
C
CyC2018 已提交
1329 1330 1331 1332 1333 1334
private int rank(Key key, Node x) {
    if (x == null) return 0;
    int cmp = key.compareTo(x.key);
    if (cmp == 0) return size(x.left);
    else if (cmp < 0) return rank(key, x.left);
    else return 1 + size(x.left) + rank(key, x.right);
C
CyC2018 已提交
1335 1336 1337
}
```

C
CyC2018 已提交
1338
### 6. min()
C
CyC2018 已提交
1339 1340

```java
C
CyC2018 已提交
1341 1342 1343 1344
private Node min(Node x) {
    if (x == null) return null;
    if (x.left == null) return x;
    return min(x.left);
C
CyC2018 已提交
1345 1346 1347
}
```

C
CyC2018 已提交
1348
### 7. deleteMin()
C
CyC2018 已提交
1349 1350 1351

令指向最小节点的链接指向最小节点的右子树。

C
CyC2018 已提交
1352
<div align="center"> <img src="../pics//dd15a984-e977-4644-b127-708cddb8ca99.png" width="500"/> </div><br>
C
CyC2018 已提交
1353 1354

```java
C
CyC2018 已提交
1355 1356
public void deleteMin() {
    root = deleteMin(root);
C
CyC2018 已提交
1357
}
C
CyC2018 已提交
1358 1359 1360 1361 1362
public Node deleteMin(Node x) {
    if (x.left == null) return x.right;
    x.left = deleteMin(x.left);
    x.N = size(x.left) + size(x.right) + 1;
    return x;
C
CyC2018 已提交
1363 1364 1365
}
```

C
CyC2018 已提交
1366
### 8. delete()
C
CyC2018 已提交
1367

C
CyC2018 已提交
1368 1369
- 如果待删除的节点只有一个子树,那么只需要让指向待删除节点的链接指向唯一的子树即可;
- 否则,让右子树的最小节点替换该节点。
C
CyC2018 已提交
1370

C
CyC2018 已提交
1371
<div align="center"> <img src="../pics//fa568fac-ac58-48dd-a9bb-23b3065bf2dc.png" width="400"/> </div><br>
C
CyC2018 已提交
1372 1373 1374


```java
C
CyC2018 已提交
1375 1376
public void delete(Key key) {
    root = delete(root, key);
C
CyC2018 已提交
1377
}
C
CyC2018 已提交
1378 1379 1380 1381 1382 1383 1384 1385 1386 1387 1388 1389 1390 1391 1392
private Node delete(Node x, Key key) {
    if (x == null) return null;
    int cmp = key.compareTo(x.key);
    if (cmp < 0) x.left = delete(x.left, key);
    else if (cmp > 0) x.right = delete(x.right, key);
    else {
        if (x.right == null) return x.left;
        if (x.left == null) return x.right;
        Node t = x;
        x = min(t.right);
        x.right = deleteMin(t.right);
        x.left = t.left;
    }
    x.N = size(x.left) + size(x.right) + 1;
    return x;
C
CyC2018 已提交
1393 1394 1395
}
```

C
CyC2018 已提交
1396
### 9. keys()
C
CyC2018 已提交
1397 1398 1399 1400

利用二叉查找树中序遍历的结果为递增的特点。

```java
C
CyC2018 已提交
1401 1402 1403 1404
public Iterable<Key> keys(Key lo, Key hi) {
    Queue<Key> queue = new LinkedList<>();
    keys(root, queue, lo, hi);
    return queue;
C
CyC2018 已提交
1405
}
C
CyC2018 已提交
1406 1407 1408 1409 1410 1411 1412
private void keys(Node x, Queue<Key> queue, Key lo, Key hi) {
    if (x == null) return;
    int cmpLo = lo.compareTo(x.key);
    int cmpHi = hi.compareTo(x.key);
    if (cmpLo < 0) keys(x.left, queue, lo, hi);
    if (cmpLo <= 0 && cmpHi >= 0) queue.add(x.key);
    if (cmpHi > 0) keys(x.right, queue, lo, hi);
C
CyC2018 已提交
1413 1414 1415
}
```

C
CyC2018 已提交
1416
### 10. 性能分析
C
CyC2018 已提交
1417 1418 1419

复杂度:二叉查找树所有操作在最坏的情况下所需要的时间都和树的高度成正比。

C
CyC2018 已提交
1420
## 2-3 查找树
C
CyC2018 已提交
1421

C
CyC2018 已提交
1422
<div align="center"> <img src="../pics//ff396233-1bb1-4e74-8bc2-d7c90146f5dd.png" width="250"/> </div><br>
C
CyC2018 已提交
1423

C
CyC2018 已提交
1424
2-3 查找树引入了 2- 节点和 3- 节点,目的是为了让树平衡。一颗完美平衡的 2-3 查找树的所有空链接到根节点的距离应该是相同的。
C
CyC2018 已提交
1425

C
CyC2018 已提交
1426
### 1. 插入操作
C
CyC2018 已提交
1427

C
CyC2018 已提交
1428
插入操作和 BST 的插入操作有很大区别,BST 的插入操作是先进行一次未命中的查找,然后再将节点插入到对应的空链接上。但是 2-3 查找树如果也这么做的话,那么就会破坏了平衡性。它是将新节点插入到叶子节点上。
C
CyC2018 已提交
1429 1430 1431

根据叶子节点的类型不同,有不同的处理方式。

C
CyC2018 已提交
1432
插入到 2- 节点上,那么直接将新节点和原来的节点组成 3- 节点即可。
C
CyC2018 已提交
1433

C
CyC2018 已提交
1434
<div align="center"> <img src="../pics//7f38a583-2f2e-4738-97af-510e6fb403a7.png" width="400"/> </div><br>
C
CyC2018 已提交
1435

C
CyC2018 已提交
1436
如果是插入到 3- 节点上,就会产生一个临时 4- 节点时,需要将 4- 节点分裂成 3 个 2- 节点,并将中间的 2- 节点移到上层节点中。如果上移操作继续产生临时 4- 节点则一直进行分裂上移,直到不存在临时 4- 节点。
C
CyC2018 已提交
1437

C
CyC2018 已提交
1438
<div align="center"> <img src="../pics//ef280699-da36-4b38-9735-9b048a3c7fe0.png" width="500"/> </div><br>
C
CyC2018 已提交
1439

C
CyC2018 已提交
1440
### 2. 性质
C
CyC2018 已提交
1441

C
CyC2018 已提交
1442
2-3 查找树插入操作的变换都是局部的,除了相关的节点和链接之外不必修改或者检查树的其它部分,而这些局部变换不会影响树的全局有序性和平衡性。
C
CyC2018 已提交
1443

C
CyC2018 已提交
1444
2-3 查找树的查找和插入操作复杂度和插入顺序无关,在最坏的情况下查找和插入操作访问的节点必然不超过 logN 个,含有 10 亿个节点的 2-3 查找树最多只需要访问 30 个节点就能进行任意的查找和插入操作。
C
CyC2018 已提交
1445

C
CyC2018 已提交
1446
## 红黑二叉查找树
C
CyC2018 已提交
1447

C
CyC2018 已提交
1448
2-3 查找树需要用到 2- 节点和 3- 节点,红黑树使用红链接来实现 3- 节点。指向一个节点的链接颜色如果为红色,那么这个节点和上层节点表示的是一个 3- 节点,而黑色则是普通链接。
C
CyC2018 已提交
1449

C
CyC2018 已提交
1450
<div align="center"> <img src="../pics//4f48e806-f90b-4c09-a55f-ac0cd641c047.png" width="250"/> </div><br>
C
CyC2018 已提交
1451 1452 1453

红黑树具有以下性质:

C
CyC2018 已提交
1454 1455
1. 红链接都为左链接;
2. 完美黑色平衡,即任意空链接到根节点的路径上的黑链接数量相同。
C
CyC2018 已提交
1456 1457 1458

画红黑树时可以将红链接画平。

C
CyC2018 已提交
1459
<div align="center"> <img src="../pics//3086c248-b552-499e-b101-9cffe5c2773e.png" width="300"/> </div><br>
C
CyC2018 已提交
1460 1461

```java
C
CyC2018 已提交
1462 1463 1464 1465 1466 1467 1468 1469 1470 1471 1472 1473 1474 1475 1476 1477 1478 1479 1480 1481 1482 1483 1484 1485
public class RedBlackBST<Key extends Comparable<Key>, Value> {
    private Node root;
    private static final boolean RED = true;
    private static final boolean BLACK = false;

    private class Node {
        Key key;
        Value val;
        Node left, right;
        int N;
        boolean color;

        Node(Key key, Value val, int n, boolean color) {
            this.key = key;
            this.val = val;
            N = n;
            this.color = color;
        }
    }

    private boolean isRed(Node x) {
        if (x == null) return false;
        return x.color == RED;
    }
C
CyC2018 已提交
1486 1487 1488
}
```

C
CyC2018 已提交
1489
### 1. 左旋转
C
CyC2018 已提交
1490 1491 1492

因为合法的红链接都为左链接,如果出现右链接为红链接,那么就需要进行左旋转操作。

C
CyC2018 已提交
1493
<div align="center"> <img src="../pics//9110c1a0-8a54-4145-a814-2477d0128114.png" width="450"/> </div><br>
C
CyC2018 已提交
1494 1495

```java
C
CyC2018 已提交
1496 1497 1498 1499 1500 1501 1502 1503 1504
public Node rotateLeft(Node h) {
    Node x = h.right;
    h.right = x.left;
    x.left = h;
    x.color = h.color;
    h.color = RED;
    x.N = h.N;
    h.N = 1 + size(h.left) + size(h.right);
    return x;
C
CyC2018 已提交
1505 1506 1507
}
```

C
CyC2018 已提交
1508
### 2. 右旋转
C
CyC2018 已提交
1509 1510 1511

进行右旋转是为了转换两个连续的左红链接,这会在之后的插入过程中探讨。

C
CyC2018 已提交
1512
<div align="center"> <img src="../pics//e2f0d889-2330-424c-8193-198edebecff7.png" width="450"/> </div><br>
C
CyC2018 已提交
1513 1514

```java
C
CyC2018 已提交
1515 1516 1517 1518 1519 1520 1521 1522
public Node rotateRight(Node h) {
    Node x = h.left;
    h.left = x.right;
    x.color = h.color;
    h.color = RED;
    x.N = h.N;
    h.N = 1 + size(h.left) + size(h.right);
    return x;
C
CyC2018 已提交
1523 1524 1525
}
```

C
CyC2018 已提交
1526
### 3. 颜色转换
C
CyC2018 已提交
1527

C
CyC2018 已提交
1528
一个 4- 节点在红黑树中表现为一个节点的左右子节点都是红色的。分裂 4- 节点除了需要将子节点的颜色由红变黑之外,同时需要将父节点的颜色由黑变红,从 2-3 树的角度看就是将中间节点移到上层节点。
C
CyC2018 已提交
1529

C
CyC2018 已提交
1530
<div align="center"> <img src="../pics//af4639f9-af54-4400-aaf5-4e261d96ace7.png" width="300"/> </div><br>
C
CyC2018 已提交
1531 1532

```java
C
CyC2018 已提交
1533 1534 1535 1536
void flipColors(Node h){
    h.color = RED;
    h.left.color = BLACK;
    h.right.color = BLACK;
C
CyC2018 已提交
1537 1538 1539
}
```

C
CyC2018 已提交
1540
### 4. 插入
C
CyC2018 已提交
1541 1542 1543

先将一个节点按二叉查找树的方法插入到正确位置,然后再进行如下颜色操作:

C
CyC2018 已提交
1544 1545 1546
- 如果右子节点是红色的而左子节点是黑色的,进行左旋转;
- 如果左子节点是红色的,而且左子节点的左子节点也是红色的,进行右旋转;
- 如果左右子节点均为红色的,进行颜色转换。
C
CyC2018 已提交
1547

C
CyC2018 已提交
1548
<div align="center"> <img src="../pics//08427d38-8df1-49a1-8990-e0ce5ee36ca2.png" width="400"/> </div><br>
C
CyC2018 已提交
1549 1550

```java
C
CyC2018 已提交
1551 1552 1553
public void put(Key key, Value val) {
    root = put(root, key, val);
    root.color = BLACK;
C
CyC2018 已提交
1554 1555
}

C
CyC2018 已提交
1556 1557 1558 1559 1560 1561
private Node put(Node x, Key key, Value val) {
    if (x == null) return new Node(key, val, 1, RED);
    int cmp = key.compareTo(x.key);
    if (cmp == 0) x.val = val;
    else if (cmp < 0) x.left = put(x.left, key, val);
    else x.right = put(x.right, key, val);
C
CyC2018 已提交
1562

C
CyC2018 已提交
1563 1564 1565
    if (isRed(x.right) && !isRed(x.left)) x = rotateLeft(x);
    if (isRed(x.left) && isRed(x.left.left)) x = rotateRight(x);
    if (isRed(x.left) && isRed(x.right)) flipColors(x);
C
CyC2018 已提交
1566

C
CyC2018 已提交
1567 1568
    x.N = size(x.left) + size(x.right) + 1;
    return x;
C
CyC2018 已提交
1569 1570 1571 1572 1573
}
```

可以看到该插入操作和二叉查找树的插入操作类似,只是在最后加入了旋转和颜色变换操作即可。

C
CyC2018 已提交
1574
根节点一定为黑色,因为根节点没有上层节点,也就没有上层节点的左链接指向根节点。flipColors() 有可能会使得根节点的颜色变为红色,每当根节点由红色变成黑色时树的黑链接高度加 1.
C
CyC2018 已提交
1575

C
CyC2018 已提交
1576
### 5. 分析
C
CyC2018 已提交
1577

C
CyC2018 已提交
1578
一颗大小为 N 的红黑树的高度不会超过 2logN。最坏的情况下是它所对应的 2-3 树,构成最左边的路径节点全部都是 3- 节点而其余都是 2- 节点。
C
CyC2018 已提交
1579 1580 1581

红黑树大多数的操作所需要的时间都是对数级别的。

C
CyC2018 已提交
1582
## 散列表
C
CyC2018 已提交
1583 1584 1585 1586 1587

散列表类似于数组,可以把散列表的散列值看成数组的索引值。访问散列表和访问数组元素一样快速,它可以在常数时间内实现查找和插入操作。

由于无法通过散列值知道键的大小关系,因此散列表无法实现有序性操作。

C
CyC2018 已提交
1588
### 1. 散列函数
C
CyC2018 已提交
1589

C
CyC2018 已提交
1590
对于一个大小为 M 的散列表,散列函数能够把任意键转换为 [0, M-1] 内的正整数,该正整数即为 hash 值。
C
CyC2018 已提交
1591

C
CyC2018 已提交
1592
散列表有冲突的存在,也就是两个不同的键可能有相同的 hash 值。
C
CyC2018 已提交
1593 1594 1595

散列函数应该满足以下三个条件:

C
CyC2018 已提交
1596 1597 1598
1. 一致性:相等的键应当有相等的 hash 值,两个键相等表示调用 equals() 返回的值相等。
2. 高效性:计算应当简便,有必要的话可以把 hash 值缓存起来,在调用 hash 函数时直接返回。
3. 均匀性:所有键的 hash 值应当均匀地分布到 [0, M-1] 之间,这个条件至关重要,直接影响到散列表的性能。
C
CyC2018 已提交
1599

C
CyC2018 已提交
1600
除留余数法可以将整数散列到 [0, M-1] 之间,例如一个正整数 k,计算 k%M 既可得到一个 [0, M-1] 之间的 hash 值。注意 M 必须是一个素数,否则无法利用键包含的所有信息。例如 M 为 10<sup>k</sup>,那么只能利用键的后 k 位。
C
CyC2018 已提交
1601 1602 1603

对于其它数,可以将其转换成整数的形式,然后利用除留余数法。例如对于浮点数,可以将其表示成二进制形式,然后使用二进制形式的整数值进行除留余数法。

C
CyC2018 已提交
1604
对于有多部分组合的键,每部分都需要计算 hash 值,并且最后合并时需要让每部分 hash 值都具有同等重要的地位。可以将该键看成 R 进制的整数,键中每部分都具有不同的权值。
C
CyC2018 已提交
1605 1606 1607 1608

例如,字符串的散列函数实现如下

```java
C
CyC2018 已提交
1609 1610 1611
int hash = 0;
for(int i = 0; i < s.length(); i++)
    hash = (R * hash + s.charAt(i)) % M;
C
CyC2018 已提交
1612 1613 1614 1615 1616
```

再比如,拥有多个成员的自定义类的哈希函数如下:

```java
C
CyC2018 已提交
1617
int hash = (((day * R + month) % M) * R + year) % M;
C
CyC2018 已提交
1618 1619
```

C
CyC2018 已提交
1620
R 通常取 31。
C
CyC2018 已提交
1621

C
CyC2018 已提交
1622
Java 中的 hashCode() 实现了 hash 函数,但是默认使用对象的内存地址值。在使用 hashCode() 函数时,应当结合除留余数法来使用。因为内存地址是 32 位整数,我们只需要 31 位的非负整数,因此应当屏蔽符号位之后再使用除留余数法。
C
CyC2018 已提交
1623 1624

```java
C
CyC2018 已提交
1625
int hash = (x.hashCode() & 0x7fffffff) % M;
C
CyC2018 已提交
1626 1627
```

C
CyC2018 已提交
1628
使用 Java 自带的 HashMap 等自带的哈希表实现时,只需要去实现 Key 类型的 hashCode() 函数即可。Java 规定 hashCode() 能够将键均匀分布于所有的 32 位整数,Java 中的 String、Integer 等对象的 hashCode() 都能实现这一点。以下展示了自定义类型如何实现 hashCode()。
C
CyC2018 已提交
1629 1630

```java
C
CyC2018 已提交
1631 1632 1633 1634 1635 1636 1637 1638 1639 1640 1641 1642
public class Transaction{
    private final String who;
    private final Date when;
    private final double amount;

    public int hashCode(){
        int hash = 17;
        hash = 31 * hash + who.hashCode();
        hash = 31 * hash + when.hashCode();
        hash = 31 * hash + ((Double) amount).hashCode();
        return hash;
    }
C
CyC2018 已提交
1643 1644 1645
}
```

C
CyC2018 已提交
1646
### 2. 基于拉链法的散列表
C
CyC2018 已提交
1647

C
CyC2018 已提交
1648
拉链法使用链表来存储 hash 值相同的键,从而解决冲突。此时查找需要分两步,首先查找 Key 所在的链表,然后在链表中顺序查找。
C
CyC2018 已提交
1649

C
CyC2018 已提交
1650
<div align="center"> <img src="../pics//b4252c85-6fb0-4995-9a68-a1a5925fbdb1.png" width="300"/> </div><br>
C
CyC2018 已提交
1651

C
CyC2018 已提交
1652
对于 N 个键,M 条链表 (N>M),如果哈希函数能够满足均匀性的条件,每条链表的大小趋向于 N/M,因此未命中的查找和插入操作所需要的比较次数为 \~N/M。
C
CyC2018 已提交
1653

C
CyC2018 已提交
1654
### 3. 基于线性探测法的散列表
C
CyC2018 已提交
1655

C
CyC2018 已提交
1656
线性探测法使用空位来解决冲突,当冲突发生时,向前探测一个空位来存储冲突的键。使用线性探测法,数组的大小 M 应当大于键的个数 N(M>N)。
C
CyC2018 已提交
1657

C
CyC2018 已提交
1658
<div align="center"> <img src="../pics//dbb8516d-37ba-4e2c-b26b-eefd7de21b45.png" width="400"/> </div><br>
C
CyC2018 已提交
1659 1660

```java
C
CyC2018 已提交
1661 1662 1663 1664 1665 1666 1667 1668 1669 1670 1671 1672 1673 1674 1675 1676 1677 1678 1679 1680 1681 1682 1683
public class LinearProbingHashST<Key, Value> {
    private int N;
    private int M = 16;
    private Key[] keys;
    private Value[] vals;

    public LinearProbingHashST() {
        init();
    }

    public LinearProbingHashST(int M) {
        this.M = M;
        init();
    }

    private void init() {
        keys = (Key[]) new Object[M];
        vals = (Value[]) new Object[M];
    }

    private int hash(Key key) {
        return (key.hashCode() & 0x7fffffff) % M;
    }
C
CyC2018 已提交
1684 1685 1686
}
```

C
CyC2018 已提交
1687
**(一)查找** 
C
CyC2018 已提交
1688 1689

```java
C
CyC2018 已提交
1690 1691 1692 1693 1694 1695 1696
public Value get(Key key) {
    for (int i = hash(key); keys[i] != null; i = (i + 1) % M) {
        if (keys[i].equals(key)) {
            return vals[i];
        }
    }
    return null;
C
CyC2018 已提交
1697 1698 1699
}
```

C
CyC2018 已提交
1700
**(二)插入** 
C
CyC2018 已提交
1701 1702

```java
C
CyC2018 已提交
1703 1704 1705 1706 1707 1708 1709 1710 1711 1712 1713 1714
public void put(Key key, Value val) {
    int i;
    for (i = hash(key); keys[i] != null; i = (i + 1) % M) {
        if (keys[i].equals(key)) {
            vals[i] = val;
            return;
        }
    }
    keys[i] = key;
    vals[i] = val;
    N++;
    resize();
C
CyC2018 已提交
1715 1716 1717
}
```

C
CyC2018 已提交
1718
**(三)删除** 
C
CyC2018 已提交
1719 1720 1721 1722

删除操作应当将右侧所有相邻的键值对重新插入散列表中。

```java
C
CyC2018 已提交
1723 1724 1725 1726 1727 1728 1729 1730 1731 1732 1733 1734 1735 1736 1737 1738 1739 1740 1741 1742
public void delete(Key key) {
    if (!contains(key)) return;
    int i = hash(key);
    while (!key.equals(keys[i])) {
        i = (i + 1) % M;
    }
    keys[i] = null;
    vals[i] = null;
    i = (i + 1) % M;
    while (keys[i] != null) {
        Key keyToRedo = keys[i];
        Value valToRedo = vals[i];
        keys[i] = null;
        vals[i] = null;
        N--;
        put(keyToRedo, valToRedo);
        i = (i + 1) % M;
    }
    N--;
    resize();
C
CyC2018 已提交
1743 1744 1745
}
```

C
CyC2018 已提交
1746
**(四)调整数组大小** 
C
CyC2018 已提交
1747

C
CyC2018 已提交
1748
线性探测法的成本取决于连续条目的长度,连续条目也叫聚簇。当聚簇很长时,在查找和插入时也需要进行很多次探测。例如下图中 2\~5 位置就是一个聚簇。
C
CyC2018 已提交
1749

C
CyC2018 已提交
1750
<div align="center"> <img src="../pics//386cd64f-7a9d-40e6-8c55-22b90ee2d258.png" width="400"/> </div><br>
C
CyC2018 已提交
1751

C
CyC2018 已提交
1752
α = N/M,把 α 称为利用率。理论证明,当 α 小于 1/2 时探测的预计次数只在 1.5 到 2.5 之间。
C
CyC2018 已提交
1753

C
CyC2018 已提交
1754
为了保证散列表的性能,应当调整数组的大小,使得 α 在 [1/4, 1/2] 之间。
C
CyC2018 已提交
1755 1756

```java
C
CyC2018 已提交
1757 1758 1759
private void resize() {
    if (N >= M / 2) resize(2 * M);
    else if (N <= M / 8) resize(M / 2);
C
CyC2018 已提交
1760 1761
}

C
CyC2018 已提交
1762 1763 1764 1765 1766 1767 1768 1769 1770 1771
private void resize(int cap) {
    LinearProbingHashST<Key, Value> t = new LinearProbingHashST<>(cap);
    for (int i = 0; i < M; i++) {
        if (keys[i] != null) {
            t.put(keys[i], vals[i]);
        }
    }
    keys = t.keys;
    vals = t.vals;
    M = t.M;
C
CyC2018 已提交
1772 1773 1774
}
```

C
CyC2018 已提交
1775
## 应用
C
CyC2018 已提交
1776

C
CyC2018 已提交
1777
### 1. 各种符号表实现的比较
C
CyC2018 已提交
1778

C
CyC2018 已提交
1779 1780 1781 1782 1783 1784 1785
| 算法 | 插入 | 查找 | 是否有序 |
| :---: | :---: | :---: | :---: |
| 二分查找实现的有序表 | N | logN | yes |
| 二叉查找树 | logN | logN | yes |
| 2-3 查找树 | logN | logN | yes |
| 拉链法实现的散列表 | N/M | N/M | no |
| 线性探测法试下的散列表 | 1 | 1 | no |
C
CyC2018 已提交
1786 1787 1788

应当优先考虑散列表,当需要有序性操作时使用红黑树。

C
CyC2018 已提交
1789
### 2. Java 的符号表实现
C
CyC2018 已提交
1790

C
CyC2018 已提交
1791 1792
- java.util.TreeMap:红黑树
- java.util.HashMap:拉链法的散列表
C
CyC2018 已提交
1793

C
CyC2018 已提交
1794
### 3. 集合类型
C
CyC2018 已提交
1795 1796 1797

除了符号表,集合类型也经常使用,它只有键没有值,可以用集合类型来存储一系列的键然后判断一个键是否在集合中。

C
CyC2018 已提交
1798
### 4. 稀疏向量乘法
C
CyC2018 已提交
1799

C
CyC2018 已提交
1800
当向量为稀疏向量时,可以使用符号表来存储向量中的非 0 索引和值,使得乘法运算只需要对那些非 0 元素进行即可。
C
CyC2018 已提交
1801 1802

```java
C
CyC2018 已提交
1803 1804 1805 1806 1807 1808 1809 1810 1811 1812 1813 1814 1815 1816 1817 1818 1819 1820 1821 1822 1823 1824 1825
public class SparseVector {
    private HashMap<Integer, Double> hashMap;

    public SparseVector(double[] vector) {
        hashMap = new HashMap<>();
        for (int i = 0; i < vector.length; i++) {
            if (vector[i] != 0) {
                hashMap.put(i, vector[i]);
            }
        }
    }

    public double get(int i) {
        return hashMap.getOrDefault(i, 0.0);
    }

    public double dot(SparseVector other) {
        double sum = 0;
        for (int i : hashMap.keySet()) {
            sum += this.get(i) * other.get(i);
        }
        return sum;
    }
C
CyC2018 已提交
1826 1827 1828
}
```

C
CyC2018 已提交
1829 1830 1831 1832
# 参考资料

- Sedgewick, Robert, and Kevin Wayne. _Algorithms_. Addison-Wesley Professional, 2011.