Bài giảng Toán rời rạc 2

pdf 124 trang phuongnguyen 6210
Bạn đang xem 20 trang mẫu của tài liệu "Bài giảng Toán rời rạc 2", để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên

Tài liệu đính kèm:

  • pdfbai_giang_toan_roi_rac_2.pdf

Nội dung text: Bài giảng Toán rời rạc 2

  1. HỌC VIỆN CÔNG NGHỆ BƯU CHÍNH VIỄN THÔNG  KHOA CÔNG NGHỆ THÔNG TIN 1 BÀI GIẢNG TOÁN RỜI RẠC 2 PTIT Hà Nội 2013
  2. LỜI GIỚI THIỆU Toán rời rạc là một lĩnh vực nghiên cứu và xử lý các đối tượng rời rạc dùng để đếm các đối tượng, và nghiên cứu mối quan hệ giữa các tập rời rạc. Một trong những yếu tố làm Toán rời rạc trở nên quan trọng là việc lưu trữ, xử lý thông tin trong các hệ thống máy tính về bản chất là rời rạc. Chính vì lý do đó, Toán học rời rạc là một môn học bắt buộc mang tính chất kinh điển của các ngành Công nghệ thông tin và Điện tử Viễn thông. Tài liệu hướng dẫn môn học Toán học rời rạc được xây dựng được xây dựng dựa trên cơ sở kinh nghiệm giảng dạy môn học và kế thừa từ giáo trình [1, 2]. Tài liệu được trình bày thành hai phần. Trong đó, phần I trình bày những kiến thức cơ bản về lý thuyết tổ hợp thông qua việc giải quyết bốn bài toán cơ bản đó là: Bài toán đếm, Bài toán tồn tại, Bài toán liệt kê và Bài toán tối ưu. Phần II trình bày những kiến thức cơ bản về Lý thuyết đồ thị: khái niệm, định nghĩa, các thuật toán trên đồ thị, đồ thị Euler, đồ thị Hamilton. Một số bài toán có ứng dụng thực tiễn quan trọng khác của lý thuyết đồ thị cũng được chú trọng giải quyết đó là Bài toán tô màu đồ thị, Bài toán tìm đường đi ngắn nhất và Bài toán luồng cực đại trong mạng. Trong mỗi phần của tài liệu, chúng tôi cố gắng trình bày ngắn gọn trực tiếp vào bản chất của vấn đề, đồng thời cài đặt hầu hết các thuật toán bằng ngôn ngữ lập trình C nhằm đạt được hai mục tiêu chính cho người học: Nâng cao tư duy toán học trong phân tích, thiết kế thuật toán và rèn luyện kỹ năng lập trình với những thuật toán phức tạp. Mặc dù đã rất cẩn trọng trong quá trPTITình biên soạn, tuy nhiên tài liệu không tránh khỏi những thiếu sót và hạn chế. Chúng tôi rất mong được sự góp ý quí báu của tất cả đọc giả và các bạn đồng nghiệp. Hà nội, tháng 11 năm 2013 2
  3. MỤC LỤC CHƯƠNG 1. MỘT SỐ KHÁI NIỆM CƠ BẢN CỦA ĐỒ THỊ 7 1.1. Định nghĩa và khái niệm 7 1.2. Một số thuật ngữ cơ bản trên đồ thị vô hướng 10 1.2.1. Bậc của đỉnh 10 1.2.2. Đường đi, chu trình, đồ thị liên thông 11 1.3. Một số thuật ngữ cơ bản trên đồ thị có hướng 13 1.3.1. Bán bậc của đỉnh 13 1.3.2. Đồ thị có hướng liên thông mạnh, liên thông yếu 13 1.4. Một số dạng đồ thị đặc biệt 15 1.5. Những điểm cần ghi nhớ 16 CHƯƠNG II. BIỂU DIỄN ĐỒ THỊ TRÊN MÁY TÍNH 17 2.1.Biểu diễn đồ thị bằng ma trận kề 17 2.1.1. Ma trận kề của đồ thị vô hướng 17 2.1.2. Ma trận kề của đồ thị có hướng 18 2.1.3. Ma trận trọng số 19 2.1.4. Qui ước khuôn dạng lưu trữ ma trận kề 20 2.2. Biểu diễn đồ thị bằng danh sách cạnh (cung ) 20 2.2.1. Biểu diễn đồ thị vô hướng bằng danh sách cạnh 20 2.2.2. Biểu diễn đồ thị có hướng bằng danh sách cạnh 21 2.2.3. Biểu diễn đồ thị trọng số bằng danh sách cạnh 22 2.2.4. Qui ước khuôn dạng lưu trữ danh sách cạnh 22 2.2.5. Cấu trúc dữ liệu biểu diễn danh sách cạnh 23 2.3. Biểu diễn đồ thị bằng danh sách kề 24 2.3.1. Biểu diễn danh sách kề dựa vào mảng 25 2.3.2. Biểu diễn danh sách kề bằng danh sách liên kết 25 2.3.3. Qui ước khuôn dạng lưu trữ danh sách kề: 26 2.4. Những điểm cần ghi nhớPTIT 26 BÀI TẬP 27 CHƯƠNG 3. TÌM KIẾM TRÊN ĐỒ THỊ 31 3.1. Thuật toán tìm kiếm theo chiều sâu (Depth First Search) 31 3.1.1.Biểu diễn thuật toán DFS(u) 31 3.1.2. Độ phức tạp thuật toán 32 3.1.3. Kiểm nghiệm thuật toán 33 3.1.4. Cài đặt thuật toán 35 3.2. Thuật toán tìm kiếm theo chiều rộng (Breadth First Search) 37 3.2.1. Biểu diễn thuật toán 37 3.2.2. Độ phức tạp thuật toán 38 3
  4. 3.2.3. Kiểm nghiệm thuật toán 38 3.2.4. Cài đặt thuật toán 39 3.3. Ứng dụng của thuật toán DFS và BFS 41 3.3.1. Xác định thành phần liên thông của đồ thị 41 a) Đặt bài toán 41 b) Mô tả thuật toán 41 c) Kiểm nghiệm thuật toán 42 d) Cài đặt thuật toán 43 3.3.2. Tìm đường đi giữa các đỉnh trên đồ thị 44 a) Đặt bài toán 44 b) Mô tả thuật toán 44 c) Kiểm nghiệm thuật toán 46 d) Cài đặt thuật toán 47 3.3.3. Tính liên thông mạnh trên đồ thị có hướng 49 a) Đặt bài toán 49 b) Mô tả thuật toán 49 c) Kiểm nghiệm thuật toán 49 d) Cài đặt thuật toán 51 3.3.4. Duyệt các đỉnh trụ 53 a) Đặt bài toán 53 b) Mô tả thuật toán 53 c) Kiểm nghiệm thuật toán 53 d) Cài đặt thuật toán 54 3.3.5. Duyệt các cạnh cầu 56 a) Đặt bài toán 56 b) Mô tả thuật toán 56 c) Kiểm nghiệm thuật toán 57 d) Cài đặt thuật toán 58 3.4. Một số bài toán quan trọng khác 61 2.4.1. Duyệt các thành phần liên thông mạnh của đồ thị 61 2.4.2. Bài toán định chiều đồ thị 61 3.5. Một số điểm cần ghi nhớ 62 BÀI TẬP PTIT 63 CHƯƠNG 4. ĐỒ THỊ EULER, ĐỒ THỊ HAMIL TON 67 4.1. Đồ thị Euler, đồ thị nửa Euler 67 4.2. Thuật toán tìm chu trình Euler 67 4.2.1. Chứng minh đồ thị là Euler 68 4.2.2. Biểu diễn thuật toán tìm chu trình Euler 69 4.2.3. Kiểm nghiệm thuật toán 70 4.2.4. Cài đặt thuật toán 70 4.3. Thuật toán tìm đường đi Euler 72 4.3.1. Chứng minh đồ thị là nửa Euler 72 4.3.2. Thuật toán tìm đường đi Euler 74 4
  5. 4.3.3. Kiểm nghiệm thuật toán 74 4.3.4. Cài đặt thuật toán 76 4.4. Đồ thị Hamilton 77 4.4.1. Thuật toán tìm tất cả các chu trình Hamilton 78 4.4.2. Kiểm nghiệm thuật toán 79 4.4.3. Cài đặt thuật toán 79 4.4.3. Cài đặt thuật toán 81 4.5. Những điểm cần ghi nhớ 82 BÀI TẬP 83 CHƯƠNG 5. CÂY KHUNG CỦA ĐỒ THỊ 86 5.1. Cây và một số tính chất cơ bản 86 5.2. Xây dựng cây khung của đồ thị dựa vào thuật toán DFS 87 5.2.1. Mô tả thuật toán 87 5.2.2. Kiểm nghiệm thuật toán 88 5.2.3. Cài đặt thuật toán 89 5.3. Xây dựng cây khung của đồ thị dựa vào thuật toán BFS 90 5.3.1. Cài đặt thuật toán 91 5.3.2. Kiểm nghiệm thuật toán 91 5.3.3. Cài đặt thuật toán 92 5.4. Bài toán xây dựng cây khung có độ dài nhỏ nhất 94 5.4.1. Đặt bài toán 94 5.4.2. Thuật toán Kruskal 95 a) Mô tả thuật toán 95 b) Kiểm nghiệm thuật toán 96 c) Cài đặt thuật toán 97 5.4.2. Thuật toán Prim 99 a) Mô tả thuật toán 100 b) Kiểm nghiệm thuật toán 100 c) Cài đặt thuật toán 101 5.5. Những nội dung cần ghi nhớ 103 BÀI TẬP PTIT 104 CHƯƠNG 6. BÀI TOÁN TÌM ĐƯỜNG ĐI NGẮN NHẤT 106 6.1. Phát biểu bài toán 106 6.2. Thuật toán Dijkstra 106 6.2.1. Mô tả thuật toán 107 6.2.2. Kiểm nghiệm thuật toán 107 6.2.3. Cài đặt thuật toán 109 6.3.Thuật toán Bellman-Ford 111 6.3.1. Mô tả thuật toán 111 6.3.2. Kiểm nghiệm thuật toán 112 6.3.3. Cài đặt thuật toán 114 5
  6. 6.4.Thuật toán Floy 116 6.4.1. Mô tả thuật toán 116 6.4.2. Cài đặt thuật toán 117 6.5. Những nội dung cần ghi nhớ 119 BÀI TẬP 120 PTIT 6
  7. CHƯƠNG 1. MỘT SỐ KHÁI NIỆM CƠ BẢN CỦA ĐỒ THỊ Nội dung chính của chương này đề cập đến những khái niệm cơ bản nhất của đồ thị, bao gồm:  Định nghĩa và ví dụ.  Phân loại đồ thị vô hướng, đồ thị có hướng, đơn đồ thị, đa đồ thị.  Khái niệm về bậc và bán bậc của đỉnh.  Khái niệm về đường đi, chu trình và tính liên thông của đồ thị.  Bài tập. Bạn đọc có thể tìm thấy những kiến thức sâu hơn và rộng hơn trong các tài liệu [1], [2], [3]. 1.1. Định nghĩa và khái niệm Đồ thị (Graph) là một cấu trúc dữ liệu rời rạc bao gồm các đỉnh và các cạnh nối các cặp đỉnh này. Chúng ta phân biệt đồ thị thông qua kiểu và số lượng cạnh và hướng của mỗi cạnh nối giữa các cặp đỉnh của đồ thị. Để minh chứng cho các loại đồ thị, chúng ta xem xét một số ví dụ về các loại mạng máy tính bao gồm: mỗi máy tính là một đỉnh, mỗi cạnh là những kênh điện thoại được nối giữa hai máy tính với nhau. Hình 1.1, là sơ đồ của mạng máy tính loại 1. San Francisco Detroit Chicago New York DenverPTIT Los Angeles Washington Hình 1.1. Đơn đồ thị vô hướng. Trong mạng máy tính này, mỗi máy tính là một đỉnh của đồ thị, mỗi cạnh vô hướng biểu diễn các đỉnh nối hai đỉnh phân biệt, không có hai cặp đỉnh nào nối cùng một cặp đỉnh. Mạng loại này có thể biểu diễn bằng một đơn đồ thị vô hướng. Định nghĩa 1. Đơn đồ thị vô hướng G = bao gồm V là tập các đỉnh, E là tập các cặp không có thứ tự gồm hai phần tử khác nhau của V gọi là các cạnh. 7
  8. Trong trường hợp giữa hai máy tính nào đó thường xuyên truyền tải nhiều thông tin, người ta nối hai máy tính bởi nhiều kênh thoại khác nhau. Mạng máy tính đa kênh thoại có thể được biểu diễn như Hình 1.2. San Francisco Detroit Chicago New York Denver Los Angeles Washington Hình 1.2. Đa đồ thị vô hướng. Trên Hình 1.2, giữa hai máy tính có thể được nối với nhau bởi nhiều hơn một kênh thoại. Với mạng loại này, chúng ta không thể dùng đơn đồ thị vô hướng để biểu diễn. Đồ thị loại này là đa đồ thị vô hướng. Định nghĩa 2. Đa đồ thị vô hướng G = bao gồm V là tập các đỉnh, E là họ các cặp không có thứ tự gồm hai phần tử khác nhau của V gọi là tập các cạnh. e1 E, e2 E được gọi là cạnh bội nếu chúng cùng tương ứng với một cặp đỉnh. Rõ ràng, mọi đơn đồ thị đều là đa đồ thị, nhưng không phải đa đồ thị nào cũng là đơn đồ thị vì giữa hai đỉnh có thể có nhiều hơn một cạnh nối giữa chúng với nhau. Trong nhiều trường hợp, có máy tính có thể nối nhiều kênh thoại với chính nó. Với loại mạng này, ta không thể dùng đa đồ thị để biểu diễn mà phải dùng giả đồ thị vô hướng. Giả đồ thị vô hướng được mô tả như trong Hình 1.3. Định nghĩa 3. Giả đồ thị vô hướng G = bao gồm V là tập đỉnh, E là họ các cặp không có thứ tự gồm hai phần tử (hai phần tử không nhất thiết phải khác nhau) trong V được gọi là các cạnh. Cạnh e được gọi là khuyên nếu có dạng e =(u, u), trong đó u là đỉnh nào đó thuộc V. San Francisco PTIT Detroit Chicago New York Denver Los Angeles Washington Hình 1.3. Giả đồ thị vô hướng. Trong nhiều mạng, các kênh thoại nối giữa hai máy tính có thể chỉ được phép truyền tin theo một chiều. Chẳng hạn máy tính đặt tại San Francisco được phép truy nhập tới máy tính đặt tại Los Angeles, nhưng máy tính đặt tại Los Angeles không được phép 8
  9. truy nhập ngược lại San Francisco. Hoặc máy tính đặt tại Denver có thể truy nhập được tới máy tính đặt tại Chicago và ngược lại máy tính đặt tại Chicago cũng có thể truy nhập ngược lại máy tính tại Denver. Để mô tả mạng loại này, chúng ta dùng khái niệm đơn đồ thị có hướng. Đơn đồ thị có hướng được mô tả như trong Hình 1.4. San Francisco Detroit Chicago New York Denver Los Angeles Washington Hình 1.4. Đơn đồ thị có hướng. Định nghĩa 4. Đơn đồ thị có hướng G = bao gồm V là tập các đỉnh, E là tập các cặp có thứ tự gồm hai phần tử của V gọi là các cung. Đồ thị có hướng trong Hình 1.4 không chứa các cạnh bội. Nên đối với các mạng đa kênh thoại một chiều, đồ thị có hướng không thể mô tả được mà ta dùng khái niệm đa đồ thị có hướng. Mạng có dạng đa đồ thị có hướng được mô tả như trong Hình 1.5. San Francisco Detroit Chicago New York Denver Los Angeles Washington Hình 5.5. Đa đồ thị có hướng. Định nghĩa 5. Đa đồ thị có hướng G = bao gồm V là tập đỉnh, E là cặp có thứ tự gồm hai phần tử của V đPTITược gọi là các cung. Hai cung e1, e2 tương ứng với cùng một cặp đỉnh được gọi là cung lặp. Từ những dạng khác nhau của đồ thị kể trên, chúng ta thấy sự khác nhau giữa các loại đồ thị được phân biệt thông qua các cạnh của đồ thị có thứ tự hay không có thứ tự, các cạnh bội, khuyên có được dùng hay không. Ta có thể tổng kết các loại đồ thị thông qua Bảng 1. Bảng 1. Phân biệt các loại đồ thị Loại đồ thị Cạnh Có cạnh bội Có khuyên 1. Đơn đồ thị vô hướng Vô hướng Không Không 2. Đa đồ thị vô hướng Vô hướng Có Không 3. Giả đồ thị vô hướng Vô hướng Có Có 9
  10. 4. Đơn đồ thị có hướng Có hướng Không Không 5. Đa đồ thị có hướng Có hướng Có Có 1.2. Một số thuật ngữ cơ bản trên đồ thị vô hướng Cho đồ thị vô hướng G = , trong đó V là tập đỉnh, E là tập cạnh. Ta bắt đầu làm quen với một số khái niệm cơ bản dưới đây. 1.2.1. Bậc của đỉnh Định nghĩa 1. Hai đỉnh u và v của đồ thị vô hướng G = được gọi là kề nhau nếu (u,v) là cạnh thuộc đồ thị G. Nếu e =(u, v) là cạnh của đồ thị G thì ta nói cạnh này liên thuộc với hai đỉnh u và v, hoặc ta nói cạnh e nối đỉnh u với đỉnh v, đồng thời các đỉnh u và v sẽ được gọi là đỉnh đầu của cạnh (u,v). Định nghĩa 2. Ta gọi bậc của đỉnh v trong đồ thị vô hướng là số cạnh liên thuộc với nó và ký hiệu là deg(v). b c d a f e g Hình 1.6 Đồ thị vô hướng G. Ví dụ 1. Xét đồ thị trong Hình 1.6, ta có: deg(a) = 2, deg(b) =deg(c) = deg(f) = 4; deg(e) = 3, deg(d) = 1, deg(g)=0. Đỉnh có bậc 0 được gọi làPTIT đỉnh cô lập. Đỉnh bậc 1 được gọi là đỉnh treo. Vì vậy : Đỉnh g là đỉnh cô lập của đồ thị Đỉnh d là đỉnh treo của đồ thị. Định lý 1. Giả sử G = là đồ thị vô hướng với m cạnh. Khi đó 2m deg(v) . v V Chứng minh. Rõ ràng mỗi cạnh e=(u,v) bất kỳ, được tính một lần trong deg(u) và một lần trong deg(v). Từ đó suy ra số tổng tất cả các bậc bằng hai lần số cạnh. Hệ quả. Trong đồ thị vô hướng G= , số các đỉnh bậc lẻ là một số chẵn. Chứng minh. Gọi O là tập các đỉnh bậc chẵn và V là tập các đỉnh bậc lẻ. Từ định lý 1 ta suy ra: 10
  11. 2m deg(v) deg(v) deg(v) v V v O v U Do deg(v) là chẵn với v là đỉnh trong O nên tổng thứ hai trong vế phải cũng là một số chẵn. 1.2.2. Đường đi, chu trình, đồ thị liên thông Định nghĩa 1. Đường đi độ dài n từ đỉnh u đến đỉnh v trên đồ thị vô hướng G= là dãy x0, x1, . . ., xn-1, xn , trong đó n là số nguyên dương, x0=u, xn=v, (xi, xi+1) E, i =0, 1, 2, . . ., n-1. Đường đi như trên còn có thể biểu diễn thành dãy các cạnh (x0, x1), (x1,x2) , . . ., (xn-1, xn). Đỉnh u là đỉnh đầu, đỉnh v là đỉnh cuối của đường đi. Đường đi có đỉnh đầu trùng với đỉnh cuối (u=v) được gọi là chu trình. Đường đi hay chu trình được gọi là đơn nếu như không có cạnh nào lặp lại. Ví dụ 1. Tìm các đường đi, chu trình trong đồ thị vô hướng như trong Hình 1.7. a, d, c, f, e là đường đi đơn độ dài 4. d, e, c, a không là đường đi vì (e,c) không phải là cạnh của đồ thị. Dãy b, c, f, e, b là chu trình độ dài 4. Đường đi a, b, e, d, a, b có độ dài 5 không phải là đường đi đơn vì cạnh (a,b) có mặt hai lần. a b c d e f Hình 1.7. Đường đi trên đồ thị. Định nghĩa 2. Đồ thị vô hướng được gọi là liên thông nếu luôn tìm được đường đi giữa hai đỉnh bất kỳ của nó. PTIT Trong trường hợp đồ thị G= không liên thông, ta có thể phân rã G thành một số đồ thị con liên thông mà chúng đôi một không có đỉnh chung. Mỗi đồ thị con như vậy được gọi là một thành phần liên thông của G. Như vậy, đồ thị liên thông khi và chỉ khi số thành phần liên thông của nó là 1. Đối với đồ thị vô hướng, đường đi từ đỉnh u đến đỉnh v cũng giống như đường đi từ đỉnh v đến đỉnh u. Chính vì vậy, nếu tồn tại đỉnh u V sao cho u có đường đi đến tất cả các đỉnh còn lại của đồ thị thì ta kết luận được đồ thị là liên thông. 11
  12. Ví dụ 2. Tìm các thành phần liên thông của đồ thị Hình 1.8 dưới đây. Số thành phần liên thông của G là 3. Thành phần liên thông thứ nhất gồm các đỉnh 1, 2, 3, 4, 6, 7. Thành phần liên thông thứ hai gồm các đỉnh 5, 8, 9, 10. Thành phần liên thông thứ ba gồm các đỉnh 11, 12, 13. Định nghĩa 3. Cạnh e E được gọi là cầu nếu loại bỏ e làm tăng thành phần liên thông của đồ thị. Đỉnh u V được gọi là đỉnh trụ nếu loại bỏ u cùng với các cạnh nối với u làm tăng thành phần liên thông của đồ thị. Ví dụ 3. Tìm các cạnh cầu và đỉnh trụ của đồ thị Hình 1.8. 2 6 8 7 1 4 3 5 10 11 9 13 12 Hình 1.8. Đồ thị vô hướng G Lời giải. Cạnh (5, 9) là cầu vì nếu loại bỏ (5, 9) thì số thành phần liên thông của đồ thị tăng từ 3 lên 4. Cạnh (5, 10) là cầu vì nếu loại bỏ (5, 10) thì số thành phần liên thông của đồ thị tăng từ 3 lên 4. Cạnh (6, 7) là cầu vì nếu loại bỏ (6, 7) thì số thành phần liên thông của đồ thị tăng từ 3 lên 4. Cạnh (8, 10) là cầu vì nếu loại bỏ (8, 10) thì số thành phần liên thông của đồ thị tăng từ 3 lên 4. Các cạnh còn lại PTITkhông là cầu vì nếu loại bỏ cạnh không làm tăng thành phần liên thông của đồ thị. Đỉnh 5 là đỉnh trụ vì nếu loại bỏ đỉnh 5 cùng với các cạnh nối với đỉnh 5 số thành phần liên thông của đồ thị tăng từ 3 lên 4. Đỉnh 6 là đỉnh trụ vì nếu loại bỏ đỉnh 6 cùng với các cạnh nối với đỉnh 6 số thành phần liên thông của đồ thị tăng từ 3 lên 4. Đỉnh 10 là đỉnh trụ vì nếu loại bỏ đỉnh 10 cùng với các cạnh nối với đỉnh 10 số thành phần liên thông của đồ thị tăng từ 3 lên 4. Các đỉnh còn lại không là trụ vì nếu loại bỏ đỉnh cùng với các cạnh nối với đỉnh không làm tăng thành phần liên thông của đồ thị. 12
  13. 1.3. Một số thuật ngữ cơ bản trên đồ thị có hướng Cho đồ thị có hướng G = , trong đó V là tập đỉnh, E là tập cạnh. Ta bắt đầu làm quen với một số khái niệm cơ bản dưới đây. 1.3.1. Bán bậc của đỉnh Định nghĩa 1. Nếu e=(u,v) là cung của đồ thị có hướng G thì ta nói hai đỉnh u và v là kề nhau, và nói cung (u, v) nối đỉnh u với đỉnh v, hoặc nói cung này đi ra khỏi đỉnh u và đi vào đỉnh v. Đỉnh u được gọi là đỉnh đầu, đỉnh v được gọi là đỉnh cuối của cung (u,v). Định nghĩa 2. Ta gọi bán bậc ra của đỉnh v trên đồ thị có hướng là số cung của đồ thị đi ra khỏi v và ký hiệu là deg+(v). Ta gọi bán bậc vào của đỉnh v trên đồ thị có hướng là số cung của đồ thị đi vào v và ký hiệu là deg-(v). a b c e d Hình 1.9. Đồ thị có hướng G. Ví dụ 2. Xét đồ thị có hướng trong Hình 1.10, ta có deg+(a) = 2, deg+(b) = 2, deg+(c) = 0, deg+(d) = 1, deg+(e) = 1. deg-(a) = 1, deg-(b) = 1, deg-(c) = 2, deg-(d) = 2, deg-(e) = 1. Do mỗi cung (u,v) được tính một lần trong bán bậc vào của đỉnh v và một lần trong bán bậc ra của đỉnh u nên ta có: Định lý 1. Giả sử G = là đồ thị có hướng. Khi đó deg (v) deg (v) | E | . PTIT v V v V Rất nhiều tính chất của đồ thị có hướng không phụ thuộc vào hướng trên các cung của nó. Vì vậy, trong nhiều trường hợp, ta bỏ qua các hướng trên cung của đồ thị. Đồ thị vô hướng nhận được bằng cách bỏ qua hướng trên các cung được gọi là đồ thị vô hướng tương ứng với đồ thị có hướng đã cho. 1.3.2. Đồ thị có hướng liên thông mạnh, liên thông yếu Khái niệm đường đi và chu trình trên đồ thị có hướng được định nghĩa hoàn toàn tương tự, chỉ có điều khác biệt duy nhất là ta phải chú ý tới các cung của đồ thị. 13
  14. Định nghĩa 1. Đường đi độ dài n từ đỉnh u đến đỉnh v trong đồ thị có hướng G= là dãy x0, x1, . . ., xn , trong đó, n là số nguyên dương, u = x0, v = xn, (xi, xi+1) E. Đường đi như trên có thể biểu diễn thành dãy các cung : (x0, x1), (x1, x2), . . ., (xn-1, xn). Đỉnh u được gọi là đỉnh đầu, đỉnh v được gọi là đỉnh cuối của đường đi. Đường đi có đỉnh đầu trùng với đỉnh cuối (u=v) được gọi là một chu trình. Đường đi hay chu trình được gọi là đơn nếu như không có hai cạnh nào lặp lại. Đối với đồ thị vô hướng, đường đi từ đỉnh u đến đỉnh v cũng giống như đường đi từ đỉnh v đến đỉnh u. Đối với đồ thị có hướng, đường đi từ đỉnh u đến đỉnh v có thể không phải là đường đi từ v đến u. Chính vì vậy, đồ thị vô hướng đưa ra hai khái niệm liên thông mạnh và liên thông yếu như sau. Định nghĩa 2. Đồ thị có hướng G= được gọi là liên thông mạnh nếu giữa hai đỉnh bất kỳ u V, v V đều có đường đi từ u đến v. Như vậy, để chứng tỏ một đồ thị có hướng liên thông mạnh ta cần chứng tỏ mọi cặp đỉnh của đồ thị đều có đường đi đến nhau. Điều này hoàn toàn khác biệt với tính liên thông của đồ thị vô hướng. Định nghĩa 3. Ta gọi đồ thị vô hướng tương ứng với đồ thị có hướng G= là đồ thị tạo bởi G và bỏ hướng của các cạnh trong G. Khi đó, đồ thị có hướng G= được gọi là liên thông yếu nếu đồ thị vô hướng tương ứng với nó là liên thông. Ví dụ 1. Hình 1.10: Đồ thị G1 là liên thông mạnh, đồ thị G2 là liên thông yếu. a b c a b c e d e d PTITG1 G2 Hình 1.10. Đồ thị có hướng liên thông mạnh, liên thông yếu Định nghĩa 4. Đồ thị vô hướng G= được gọi là định chiều được nếu ta có thể biến đổi các cạnh trong G thành các cung tương ứng để nhận được một đồ thị có hướng liên thông mạnh. Định lý 1. Đồ thị vô hướng G= định chiều được khi và chỉ khi các cạnh của nó không phải là cầu. Bạn đọc có thể tìm hiểu phần chứng minh định lý trong các tài liệu [1, 2, 3]. 14
  15. 1.4. Một số dạng đồ thị đặc biệt Dưới đây là một số dang đơn đồ thị vô hướng đặc biệt có nhiều ứng dụng khác nhau của thực tế. Đồ thị đầy đủ. Đồ thị đầy đủ n đỉnh, ký hiệu là Kn, là đơn đồ thị vô hướng mà giữa hai đỉnh bất kỳ của nó đều có cạnh nối. Ví dụ đồ thị K3, K4, K5 trong Hình 1.11. 2 1 1 2 1 3 2 3 4 3 5 4 K3 K4 K5 Hình 1.11. Đồ thị K3, K4, K5. Đồ thị vòng. Đồ thị vòng Cn (n 3) có các cạnh (1,2), (2,3), ,(n-1,n), (n,1). Ví dụ đồ thị C3, C4, C5 trong Hình 1.12. 2 1 1 2 1 3 2 3 4 3 5 4 C3 C4 C5 Hình 1.12. Đồ thị C3, C4, C5. Đồ thị bánh xe. Đồ thị bánh xe Wn thu được bằng cách bổ sung một đỉnh nối với tất cả các đỉnh của Cn. Ví dụ đồ thị W3, W4, W5 trong Hình 1.13. PTIT 2 1 1 2 51 1 6 3 4 2 3 4 3 5 4 Hình 1.13. Đồ thị C3, C4, C5. Đồ thị hai phía. Đồ thị G = được gọi là đồ thị hai phía nếu tập đỉnh V của nó có thể phân hoạch thành hai tập X và Y sao cho mỗi cạnh của đồ thị chỉ có dạng (x, y), trong đó x X và y Y. Ví dụ đồ thị K2,3, K33, K3,5 trong Hình 1.14. 15
  16. 8 3 1 6 7 1 1 4 2 5 2 6 2 3 4 5 5 3 4 Hình 1.13. Đồ thị K2,3, K3,3, K3,5. 1.5. Những điểm cần ghi nhớ  Nắm vững và phân biệt rõ các loại đồ thị: đơn đồ thị, đa đồ thị, đồ thị vô hướng, đồ thị có hướng, đồ thị trọng số.  Nắm vững những khái niệm cơ bản trên đồ thị vô hướng.  Nắm vững những khái niệm cơ bản trên đồ thị có hướng.về đồ thị.  Nắm vững các khái niệm đường đi, chu trình, liên thông, liên thông mạnh, liên thông yếu.  Nắm vững các loại đồ thị : đồ thị đầy đủ, đồ thị vòng, đồ thị bánh xe, đồ thị hai phía PTIT 16
  17. CHƯƠNG II. BIỂU DIỄN ĐỒ THỊ TRÊN MÁY TÍNH Để lưu trữ đồ thị và thực hiện các thuật toán khác nhau, ta cần phải biểu diễn đồ thị trên máy tính, đồng thời sử dụng những cấu trúc dữ liệu thích hợp để mô tả đồ thị. Việc chọn cấu trúc dữ liệu nào để biểu diễn đồ thị có tác động rất lớn đến hiệu quả thuật toán. Vì vậy, lựa chọn cấu trúc dữ liệu thích hợp biểu diễn đồ thị sẽ phụ thuộc vào từng bài toán cụ thể. Nội dung chính của chương bao gồm:  Biểu diễn đồ thị bằng ma trận kề.  Biểu diễn đồ thị bằng danh sách cạnh.  Biểu diễn đồ thị bằng danh sách kề.  Biểu diễn đồ thị bằng ma trận liên thuộc.  Bài tập Chương 2. Bạn đọc có thể tìm thấy những kiến thức sâu hơn và rộng hơn trong các tài liệu [1], [2], [3]. 2.1.Biểu diễn đồ thị bằng ma trận kề Cấu trúc dữ liệu phổ dụng nhất để biểu diễn đồ thị là biểu diễn đồ thị bằng ma trận. Về lý thuyết, người ta đã chứng minh được mỗi ma trận vuông (0,1) cấp n đều đẳng cấu với một đơn đồ thị vô hướng hoặc có hướng. Mục này, chúng ta sẽ xem xét phương pháp biểu diễn các loại đồ thị khác nhau bằng ma trận kề. 2.1.1. Ma trận kề của đồ thị vô hướng Xét đồ thị đơn vô hướng G = , với tập đỉnh V = {1, 2, . . ., n}, tập cạnh E = {e1, e2, , em}. Ta gọi ma trận kề của đồ thị G là ma trận có các phần tử hoặc bằng 0 hoặc bằng 1 theo qui định như sau: A = { aij: aij = 1 nếu (i, j)PTIT E, aij = 0 nếu (i,j) E; i, j =1, 2, . . ., n}. Ví dụ 1. Biểu diễn đồ thị trong Hình 2.1 dưới đây bằng ma trận kề. 2 5 0 1 0 0 0 0 1 0 1 1 1 0 1 6 1 1 0 1 0 0 0 1 1 0 1 1 0 1 0 1 0 1 3 4 0 0 1 1 0 0 Hình 2.1. Ma trận kề biểu diễn đồ thị vô hướng. 17
  18. Tính chất ma trận kề đối với đồ thị vô hướng: n n a) Tổng các phần tử của ma trận bằng hai lần số cạnh : aij 2m (m là số i 1 j 1 cạnh của đồ thị. n b) Tổng các phần tử của hàng u là bậc của đỉnh u: deg(u) auj . Ví dụ với ma j 1 trận kề biểu diễn đồ thị Hình 2.1, tổng các phần tử của hàng 1 là bậc của đỉnh 1, vì vậy deg(1)=2; tổng các phần tử của hàng 2 là bậc của đỉnh 2, vì vậy deg(2)=3. n c) Tổng các phần tử của cột u là bậc của đỉnh u: deg(u) a ju . Ví dụ với ma j 1 trận kề biểu diễn đồ thị Hình 2.1, tổng các phần tử của cột 1 là bậc của đỉnh 1, vì vậy deg(1)=2; tổng các phần tử của cột 2 là bậc của đỉnh 2, vì vậy deg(2)=3. p d) Nếu ký hiệu aij ,i, j 1,2, ,n là các phần tử của ma trận. Khi đó, p p A = A.A. . . A (p lần); aij ,i, j 1,2, ,n , cho ta số đường đi khác nhau từ đỉnh i đến đỉnh j qua p-1 đỉnh trung gian. 2.1.2. Ma trận kề của đồ thị có hướng Ma trận kề của đồ thị có hướng cũng được định nghĩa hoàn toàn tương tự, chúng ta chỉ cần lưu ý tới hướng của cạnh. Ma trận kề của đồ thị có hướng là không đối xứng. Ví dụ 2. Tìm ma trận kề của đồ thị có hướng trong Hình 2.2. 2 5 0 1 0 0 0 0 0 0 1 1 1 0 1 0 0 0 0 0 1 PTIT6 0 0 1 0 1 0 0 0 0 0 0 1 3 4 0 0 0 1 0 0 Hình 2.2. Ma trận kề của đồ thị có hướng. Tính chất của ma trận kề của đồ thị có hướng: n n a) Tổng các phần tử của ma trận bằng số cạnh : aij m (m là số cạnh của đồ i 1 j 1 thị. 18
  19. n b) Tổng các phần tử của hàng u là bán đỉnh bậc ra của đỉnh u: deg (u) auj . j 1 Ví dụ với ma trận kề biểu diễn đồ thị Hình 2.2, tổng các phần tử của hàng 1 là bán đỉnh bậc a của đỉnh 1, vì vậy deg+(1)=1; tổng các phần tử của hàng 2 là bán đỉnh bậc ra của đỉnh 3, vì vậy deg+(2)=3. n c) Tổng các phần tử của cột u là bán đỉnh bậc vào của đỉnh u: deg (u) a ju . j 1 Ví dụ với ma trận kề biểu diễn đồ thị Hình 2.2, tổng các phần tử cột 1 là bán đỉnh bậc vào của đỉnh 1, vì vậy deg-(1)=1; tổng các phần tử của cột 2 là bán đỉnh bậc vào của đỉnh 2, vì vậy deg-(2)=1. p p d) Nếu ký hiệu aij ,i, j 1,2, ,n là các phần tử của ma trận. Khi đó, A = A.A. . . p A (p lần); aij ,i, j 1,2, ,n , cho ta số đường đi khác nhau từ đỉnh i đến đỉnh j qua p-1 đỉnh trung gian. 2.1.3. Ma trận trọng số Trong rất nhiều ứng dụng khác nhau của lý thuyết đồ thị, mỗi cạnh e =(u,v) của nó được gán bởi một số c(e) = c(u,v) gọi là trọng số của cạnh e. Đồ thị trong trường hợp như vậy gọi là đồ thị trọng số. Trong trường hợp đó, ma trận kề của đồ thị được thay bởi ma trận trọng số c= c[i,j], i, j= 1, 2, . . ., n. c[i,j] = c(i,j) nếu (i, j) E, c[i,j] =  nếu (i, j) E. Trong đó,  nhận các giá trị: 0, , - tuỳ theo từng tình huống cụ thể của thuật toán. Ví dụ 3. Ma trận kề của đồ thị có trọng số trong Hình 2.3. 2 3 5 5 4 5 8 6 3 6 5 1 8 6 2 7 5 2 3 PTIT 4 3 7 4 3 Hình 2.3. Ma trận kề của đồ thị có hướng. Ưu điểm của ma trận kề: Đơn giản dễ cài đặt trên máy tính bằng cách sử dụng một mảng hai chiều để biểu diễn ma trận kề; Dễ dàng kiểm tra được hai đỉnh u, v có kề với nhau hay không bằng đúng một phép so sánh (a[u][v] 0?);và chúng ta chỉ mất đúng một phép so sánh. 19
  20. Nhược điểm của ma trận kề: Lãng phí bộ nhớ: bất kể số cạnh nhiều hay ít ta cần n2 đơn vị bộ nhớ để biểu diễn; Không thể biểu diễn được với các đồ thị có số đỉnh lớn (ví dụ triệu đỉnh); Để xem xét đỉnh đỉnh u có những đỉnh kề nào cần mất n phép so sánh kể cả đỉnh u là đỉnh cô lập hoặc đỉnh treo. 2.1.4. Qui ước khuôn dạng lưu trữ ma trận kề Để thuận tiện cho những nội dung kế tiếp, ta qui ước khuôn dạng dữ liệu biểu diễn đồ thị dưới dạng ma trận kề hoặc ma trận trọng số trong file như sau: Dòng đầu tiên ghi lại số đỉnh của đồ thị; N dòng kế tiếp ghi lại ma trận kề của đồ thị. Hai phần tử khác nhau của ma trận kề được viết cách nhau một vài khoảng trống. Ví dụ ma trận kề gồm 6 đỉnh của Hình 2.1 được tổ chức trong file dothi.in như sau: dothi.in 5 0 1 0 0 0 0 1 0 1 1 1 0 1 1 0 1 0 0 0 1 0 1 0 1 0 0 1 1 0 0 2.2. Biểu diễn đồ thị bằng danh sách cạnh (cung ) Trong trường hợp đồ thị thưa (đồ thị có số cạnh m 6n), người ta thường biểu diễn đồ thị dưới dạng danh sách cạnh. Trong phép biểu diễn này, chúng ta sẽ lưu trữ danh sách tất cả các cạnh (cung) củaPTIT đồ thị vô hướng (có hướng). Mỗi cạnh (cung) e(x, y) được tương ứng với hai biến dau[e], cuoi[e]. Như vậy, để lưu trữ đồ thị, ta cần 2m đơn vị bộ nhớ. Nhược điểm lớn nhất của phương pháp này là để nhận biết những cạnh nào kề với cạnh nào chúng ta cần m phép so sánh trong khi duyệt qua tất cả m cạnh (cung) của đồ thị. Nếu là đồ thị có trọng số, ta cần thêm m đơn vị bộ nhớ để lưu trữ trọng số của các cạnh. 2.2.1. Biểu diễn đồ thị vô hướng bằng danh sách cạnh Đối với đồ thị vô hướng, mỗi cạnh là bộ không tính đến thứ tự các đỉnh. Ví dụ cạnh (u,v) và cạnh (v, u) được xem là một. Do vậy, trong khi biểu diễn đồ thị vô hướng bằng danh sách cạnh ta chỉ cần liệt kê các cạnh (u,v) mà không cần liệt kê cạnh (v,u). Để tránh nhầm lẫn, ta nên liệt kê các cạnh theo thứ tự tăng dần của đỉnh đầu mỗi cạnh. Trong 20
  21. trường hợp biểu diễn đa đồ thị vô hướng, ta bổ sung thêm một cột là số cạnh (socanh) nối giữa hai đỉnh của đồ thị. Hình 2.4 dưới đây mô tả chi tiết phương pháp biểu diễn đồ thị vô hướng bằng danh sách cạnh. Tính chất danh sách cạnh của đồ thị vô hướng: Đỉnh đầu nhỏ hơn đỉnh cuối mỗi cạnh. Số đỉnh có giá trị u thuộc cả vế phải và vế trái của danh sách cạnh là bậc của đỉnh u. Ví dụ giá trị u=1 xuất hiện 2 lần từ đó ta suy ra deg(1)=2, số 2 xuất hiện 4 lần vì vậy deg(2) = 4. 2 5 Đỉnh đầu Đỉnh cuối 1 2 1 3 2 3 1 6 2 4 2 5 3 4 4 5 3 4 4 6 5 6 Hình 2.4. Biểu diễn đồ thị vô hướng bằng danh sách cạnh. 2.2.2. Biểu diễn đồ thị có hướng bằng danh sách cạnh Trong trường hợp đồ thị có hướng, mỗi cạnh là bộ có tính đến thứ tự các đỉnh. Ví dụ cạnh (u,v) khác với cạnh (v, u). Do vậy, trong khi biểu diễn đồ thị vô hướng bằng danh sách cạnh ta đặc biệt chú ý đến hướng của các cạnh. Hình 2.5 dưới đây mô tả chi tiết phương pháp biểu diễn đồ thị có hướng bằng danh sách cạnh. 2 5 Đỉnh đầu Đỉnh Cuối 1 2 PTIT2 3 1 6 2 4 2 5 3 1 4 3 3 4 4 5 5 6 6 4 Hình 2.5. Biểu diễn đồ thị có hướng bằng danh sách cạnh. Tính chất danh sách cạnh của đồ thị vô hướng: Đỉnh đầu không nhất thiết phải nhỏ hơn đỉnh cuối mỗi cạnh. 21
  22. Số đỉnh có giá trị u thuộc cả vế phải các cạnh là deg+(u). Ví dụ giá trị u=1 xuất hiện 1 lần ở vế phải của tất cả các cạnh nên deg+(1) =1, giá trị u=2 xuất hiện 3 lần ở vế phải của tất cả các cạnh nên deg+(2) =3. Số đỉnh có giá trị u thuộc cả vế trái các cạnh là deg-(u). Ví dụ giá trị u=1 xuất hiện 1 lần ở vế trái của tất cả các cạnh nên deg-(1) =1, giá trị u=2 xuất hiện 1 lần ở vế trái của tất cả các cạnh nên deg-(2) =1. 2.2.3. Biểu diễn đồ thị trọng số bằng danh sách cạnh Trong trường hợp đồ thị có hướng (hoặc vô hướng) có trọng số, ta bổ sung thêm một cột là trọng số của mỗi cạnh. Hình 2.6 dưới đây mô tả chi tiết phương pháp biểu diễn đồ thị trọng số bằng danh sách cạnh. 2 3 5 Đỉnh đầu Đỉnh Cuối Trọng Số 4 1 2 5 5 2 3 8 6 5 1 8 6 2 4 6 2 5 3 2 3 3 1 2 4 3 7 3 7 4 4 5 5 5 6 4 6 4 3 Hình 2.6. Biểu diễn đồ thị có hướng bằng danh sách cạnh. Ưu điểm của danh sách cạnh: Trong trường hợp đồ thị thưa (m<6n), biểu diễn bằng danh sách cạnh tiết kiệm được không gian nhớ; Thuận lợi cho một số thuật toán chỉ quan tâm đến các cạnh của đồ thị. Nhược điểm của danh sáchPTIT cạnh: Khi cần duyệt các đỉnh kề với đỉnh u bắt buộc phải duyệt tất cả các cạnh của đồ thị. Điều này làm cho thuật toán có chi phí tính toán cao. 2.2.4. Qui ước khuôn dạng lưu trữ danh sách cạnh Để thuận tiện cho những nội dung kế tiếp, ta qui ước khuôn dạng dữ liệu biểu diễn đồ thị dưới dạng danh sách cạnh trong file như sau: Dòng đầu tiên ghi lại số N, M tương ứng với số đỉnh và số cạnh của đồ thị. Hai số được viết cánh nhau một vài khoảng trống; M dòng kế tiếp, mỗi dòng gi lại một cạnh của đồ thị, đỉnh đầu và đỉnh cuối mỗi cạnh được viết cách nhau một vài khoảng trống. 22
  23. Ví dụ với đồ thị trọng số cho bởi Hình 2.6 gồm 6 đỉnh và 9 cạnh được lưu trữ trong file dothi.in như sau: dothi.in 6 9 1 2 5 2 3 8 2 4 6 2 5 3 3 1 2 4 3 7 4 5 5 5 6 4 6 4 3 2.2.5. Cấu trúc dữ liệu biểu diễn danh sách cạnh Phương pháp tốt hơn cả để biểu diễn mỗi cạnh của đồ thị là sử dụng cấu trúc. Mỗi cấu trúc gồm có hai thành viên dau[e] và cuối cuoi[e]. Khi đó, danh sách cạnh của đồ thị dễ dàng được biểu diễn bằng mảng hoặc danh sách liên kết như dưới đây. Biểu diễn danh danh sách cạnh của đồ thị bằng mảng: typedef struct { //Định nghĩa một cạnh của đồ thị int dau; int cuoi; } Edge; Edge G[MAX]; //Danh sách các cạnh được biểu diễn trong mảng G. 2 5 Đỉnh đầu Đỉnh cuối 1 2 1 3 PTIT2 3 1 6 2 4 2 5 3 4 4 5 3 4 4 6 5 6 Hình 2.7. Biểu diễn đồ thị bằng danh sách cạnh 23
  24. Ví dụ với danh danh sách cạnh của đồ thị Hình 2.7, biểu diễn danh sách cạnh dựa vào mảng của đồ thị có dạng sau: Cạnh: G[1] G[2] G[3] G[4] G[5] G[6] G[7] G[8] G[9] G[i].dau 1 1 2 2 2 3 4 4 5 G[i].cuoi 2 3 3 4 5 4 5 6 6 Đối với đồ thị có hướng cũng được biểu diễn như trên nhưng ta cần chú ý đến hướng của mỗi cung. Đối với đồ thị trọng số ta chỉ cần bổ sung vào cấu trúc Edge một thành viên là trọng số của cạnh như sau: typedef struct { //Định nghĩa một cạnh có trọng số của đồ thị int dau; int cuoi; int trongso; } Edge; Edge G[MAX]; //Danh sách trọng số các cạnh biểu diễn trong mảng G. Biểu diễn danh danh sách cạnh của đồ thị bằng danh sách liên kết: typedef struct canh{ //Định nghĩa một cạnh của đồ thị int dau; int cuoi; struct node *next; } *Edge; Edge *G; //Các cạnh được của đồ thị biểu diễn bằng danh danh sách liên kết G. Ví dụ với danh danh sách cạnh của đồ thị Hình 2.7, biểu diễn danh sách cạnh dựa vào danh sách liên kết có dạng sau: 1 1 4 5 next next next Null 2 3 PTIT6 6 2.3. Biểu diễn đồ thị bằng danh sách kề Trong rất nhiều ứng dụng, cách biểu diễn đồ thị dưới dạng danh sách kề thường được sử dụng. Trong biểu diễn này, với mỗi đỉnh u của đồ thị chúng ta lưu trữ danh sách các đỉnh kề với nó mà ta ký hiệu là Ke(u), nghĩa là Ke(u) = { v V: (u, v) E}, Với cách biểu diễn này, mỗi đỉnh u của đồ thị, ta làm tương ứng với một danh sách tất cả các đỉnh kề với nó và được ký hiệu là List(u). Để biểu diễn List(u), ta có thể dùng các kiểu dữ liệu kiểu tập hợp, mảng hoặc danh sách liên kết. Hình 2.8 dưới đây đưa ra ví dụ chi tiết về biểu diễn đồ thị bằng danh sách kề. 24
  25. 2 5 Ke(1) = { 2, 3). Ke(2) = {1, 3, 4, 5}. 1 6 Ke(3) = {1, 2, 4}. Ke(4) = {2, 3, 5, 6}. Ke(5) = {2, 4, 6}. 3 4 Ke(6) = { 4, 5}. Hình 2.8. Biểu diễn đồ thị bằng danh sách kề. Ưu điểm của danh sách kề: Dễ dàng duyệt tất cả các đỉnh của một danh sách kề; Dễ dàng duyệt các cạnh của đồ thị trong mỗi danh sách kề; Tối ưu về phương pháp biểu diễn. Nhược điểm của danh sách kề: Khó khăn cho người đọc có kỹ năng lập trình yếu. 2.3.1. Biểu diễn danh sách kề dựa vào mảng Sử dụng một mảng để lưu trữ danh sách kề các đỉnh. Trong đó, mảng được chia thành n đoạn, đoạn thứ i trong mảng lưu trữ danh sách kề của đỉnh thứ i V. Ví dụ với đồ thị được cho trong Hình 2.8 ta tổ chức mảng A[] gồm 18 phần tử, trong đó mảng A[] được chia thành 6 đoạn, mỗi đoạn lưu trữ danh sách kề của đỉnh tương ứng như dưới đây. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 A[i]=? 2 3 1 3 4 5 1 2 4 2 3 5 6 2 4 6 4 5 Đoạn 1 Đoạn 2 PTITĐoạn 3 Đoạn 4 Đoạn 5 Đoạn 6 Để biết một đoạn thuộc mảng bắt đầu từ phần tử nào đến phần tử nào ta sử dụng một mảng khác dùng để lưu trữ vị trí các phần tử bắt đầu và kết thúc của đoạn. Ví dụ với danh sách kề gồm 6 đoạn như trên, ta cần xây dựng một mảng VT[6] = {0, 2, 6, 9, 13, 16, 18} để lưu trữ vị trí các đoạn trong mảng A[]. Dựa vào mảng VT[] ta có thể thấy: Ke(1) là A[1], A[2]; Ke(2) là A[3], A[4], A[5], A[6] 2.3.2. Biểu diễn danh sách kề bằng danh sách liên kết Với mỗi đỉnh u V, ta biểu diễn mỗi danh sách kề của đỉnh bằng một danh sách liên kết List(u). Ví dụ với đồ thị trong Hình 2.8 sẽ được biểu diễn bằng 6 danh sách liên kết List[1], List[2], , List[6] như dưới đây. 25
  26. List[1]: 2 3 Null List[2]: 1 3 4 5 Null List[3]: 1 2 4 Null List[4]: 2 3 5 6 Null List[5]: 2 4 6 Null Null List[6]: 4 5 2.3.3. Qui ước khuôn dạng lưu trữ danh sách kề: Dòng đầu tiên ghi lại số đỉnh của đồ thị; N dòng kế tiếp ghi lại danh sách kề của đỉnh tương ứng theo khuôn dạng: Phần tử đầu tiên là vị trí kết thúc của đoạn, tiếp đến là danh sách các đỉnh của danh sách kề. Các phần tử được ghi cách nhau một vài khoảng trống Ví dụ khuôn dạng lưu trữ danh sách kề của Hình 2.7 trong file dothi.in như sau: dothi.in 6 2 2 3 6 1 3 4 5 9 1 2 4 13 2 3 5 6 16 2 4 6 18 4 PTIT5 2.4. Những điểm cần ghi nhớ  Nắm vững và phân biệt rõ các loại đồ thị: đơn đồ thị, đa đồ thị, đồ thị vô hướng, đồ thị có hướng, đồ thị trọng số.  Nắm vững những khái niệm cơ bản về đồ thị: đường đi, chu trình, đồ thị liên thông.  Hiểu và nắm rõ bản chất của các phương pháp biểu diễn đồ thị trên máy tính. Phân tích ưu, nhược điểm của từng phương pháp biểu diễn.  Chuyển đổi các phương pháp biểu diễn qua lại lẫn nhau giúp ta hiểu được cách biểu diễn đồ thị trên máy tính. 26
  27. BÀI TẬP 1. Trong một buổi gặp mặt, mọi người đều bắt tay nhau. Hãy chỉ ra rằng số lượt người bắt tay nhau là một số chẵn. 2. Một đơn đồ thị với n đỉnh có nhiều nhất là bao nhiêu cạnh? 3. Hãy biểu diễn các đồ thị G1, G2, G3 dưới đây dưới dạng: ma trận kề, danh sách cạnh, danh sách kề. 2 5 2 5 1 4 7 1 4 7 a. Đồ thị vô hướng G1. b. Đồ thị có hướng G2. 3 6 3 6 B 8 E 5 3 7 4 A 2 PTIT D 9 G 1 6 5 9 C 4 F c. Đồ thị trọng số G3 4. Hãy tạo một file dữ liệu theo khuôn dạng như sau: a. Ma trận kề: - Dòng đầu tiên là số tự nhiên n là số các đỉnh của đồ thị. - N dòng kế tiếp là ma trận kề của đồ thị. b. Danh sách cạnh: 27
  28. - Dòng đầu tiên ghi lại số tự nhiên n và m là số các đỉnh và các cạnh của đồ thị. - M dòng kế tiếp ghi lại thứ tự đỉnh đầu, cuối của các cạnh. Hãy viết chương trình chuyển đổi một đồ thị cho dưới dạng ma trận kề thành một đồ thị cho dưới dạng danh sách cạnh và danh sách kề. Ngược lại, chuyển đổi một đồ thị cho dưới dạng danh sách cạnh thành đồ thị dưới dạng ma trận kề và danh sách cạnh. c. Danh sách kề: Dòng đầu tiên ghi lại số đỉnh của đồ thị; N dòng kế tiếp ghi lại danh sách kề của đỉnh tương ứng theo khuôn dạng: Phần tử đầu tiên là vị trí kết thúc của đoạn, tiếp đến là danh sách các đỉnh của danh sách kề. Các phần tử được ghi cách nhau một vài khoảng trống. M dòng kế tiếp ghi lại thứ tự đỉnh đầu, cuối của các cạnh. Hãy viết chương trình chuyển đổi một đồ thị cho dưới dạng ma trận kề thành một đồ thị cho dưới dạng danh sách cạnh và danh sách kề. Ngược lại, chuyển đổi một đồ thị cho dưới dạng danh sách cạnh thành đồ thị dưới dạng ma trận kề và danh sách cạnh. 5. Một bàn cờ 8 8 được đánh số theo cách sau: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50PTIT 51 52 53 54 55 56 57 58 59 60 61 62 63 64 Mỗi ô có thể coi là một đỉnh của đồ thị. Hai đỉnh được coi là kề nhau nếu một con vua đặt ở ô này có thể nhảy sang ô kia sau một bước đi. Ví dụ : ô 1 kề với ô 2, 9, 10, ô 11 kề với 2, 3, 4, 10, 12, 18, 19, 20. Hãy viết chương trình tạo ma trận kề của đồ thị, kết quả in ra file king.out. 6. Bàn cờ 8 8 được đánh số như bài trên. Mỗi ô có thể coi là một đỉnh của đồ thị . Hai đỉnh được gọi là kề nhau nếu một con mã đặt ở ô này có thể nhảy sang ô kia sau một nước đi. Ví dụ ô 1 kề với 11, 18, ô 11 kề với 1, 5, 17, 21, 26, 28. Hãy viết chương trình lập ma trận kề của đồ thị, kết quả ghi vào file matran.out. 28
  29. 7. Hãy biểu diễn đồ thị dưới đây dưới dạng ma trận kề, danh sách cạnh, danh sách kề. 2 6 9 10 1 3 5 11 13 4 8 7 12 8. Hãy biểu diễn đồ thị dưới đây dưới dạng ma trận kề, danh sách cạnh, danh sách kề. 6 10 3 9 1 12 2 13 5 4 8 11 7 9. Hãy biểu diễn đồ thị dưới đây dưới dạng ma trận trọng số, danh sách cạnh - trọng số. 2 9 7 2 8 2 12 2 2 2 1 2 9 9 9 1 8 3 8 6 8 9 8 13 6 6 6 7 1 1 6 7 4 7 5 7 10 7 11 10. Hãy biểu diễn đồ thị dưới đây dưới dạng ma trận kề, danh sách cạnh, danh sách kề. PTIT 6 10 3 9 1 12 2 13 5 4 8 11 7 29
  30. PTIT 30
  31. CHƯƠNG 3. TÌM KIẾM TRÊN ĐỒ THỊ Có nhiều thuật toán trên đồ thị được xây dựng để duyệt tất cả các đỉnh của đồ thị sao cho mỗi đỉnh được viếng thăm đúng một lần. Những thuật toán như vậy được gọi là thuật toán tìm kiếm trên đồ thị. Chúng ta cũng sẽ làm quen với hai thuật toán tìm kiếm cơ bản, đó là duyệt theo chiều sâu DFS (Depth First Search) và duyệt theo chiều rộng BFS (Breath First Search). Trên cơ sở của hai phép duyệt cơ bản, ta có thể áp dụng chúng để giải quyết một số bài toán quan trọng của lý thuyết đồ thị. Nội dung chính được đề cập trong chương này bao gồm:  Thuật toán tìm kiếm theo chiều sâu trên đồ thị.  Thuật toán tìm kiếm theo chiều rộng trên đồ thị.  Ứng dụng của thuật toán tìm kiếm theo chiều sâu.  Ứng dụng của thuật toán tìm kiếm theo chiều rộng. Bạn đọc có thể tìm hiểu sâu hơn về tính đúng đắn, độ phức tạp của các thuật toán trong các tài liệu [1, 2, 3]. 3.1. Thuật toán tìm kiếm theo chiều sâu (Depth First Search) Tư tưởng cơ bản của thuật toán tìm kiếm theo chiều sâu là bắt đầu tại một đỉnh v0 nào đó, chọn một đỉnh u bất kỳ kề với v0 và lấy nó làm đỉnh duyệt tiếp theo. Cách duyệt tiếp theo được thực hiện tương tự như đối với đỉnh v0 với đỉnh bắt đầu là u. Để kiểm tra việc duyệt mỗi đỉnh đúng một lần, chúng ta sử dụng một mảng chuaxet[] gồm n phần tử (tương ứng với n đỉnh), nếu đỉnh thứ u đã được duyệt, phần tử tương ứng trong mảng chuaxet[u] có giá trị FALSE. Ngược lại, nếu đỉnh chưa được duyệt, phần tử tương ứng trong mảng có giá trị TRUE. 3.1.1.Biểu diễn thuật toán DFS(u) Thuật toán DFS(u) có thểPTIT được mô tả bằng thủ tục đệ qui như sau: Thuật toán DFS (u): //u là đỉnh bắt đầu duyệt Begin ;//Duyệt đỉnh u chuaxet[u] := FALSE;//Xác nhận đỉnh u đã duyệt for each v ke(u) do //Lấy mỗi đỉnh v Ke(u). if (chuaxet[v] ) then //Nếu đỉnh v chưa duyệt DFS(v); //Duyệt theo chiều sâu bắt từ đỉnh v EndIf; EndFor; End. 31
  32. Thuật toán DFS(u) có thể khử đệ qui bằng cách sử dụng ngăn xếp như Hình 3.1 dưới đây: Thuật toán DFS(u): Begin Bước 1 (Khởi tạo): stack = ; //Khởi tạo stack là  Push(stack, u); //Đưa đỉnh u vào ngăn xếp ; //Duyệt đỉnh u chuaxet[u] = False; //Xác nhận đỉnh u đã duyệt Bước 2 (Lặp) : while ( stack  ) do s = Pop(stack); //Loại đỉnh ở đầu ngăn xếp for each t Ke(s) do //Lấy mỗi đỉnh t Ke(s) if ( chuaxet[t] ) then //Nếu t đúng là chưa duyệt ; // Duyệt đỉnh t chuaxet[t] = False; // Xác nhận đỉnh t đã duyệt Push(stack, s);//Đưa s vào stack Push(stack, t); //Đưa t vào stack break; //Chỉ lấy một đỉnh t EndIf; EndFor; EndWhile; Bước 3 (Trả lại kết quả): Return( ); End. Hình 3.1. Thuật toán DFS(u) dựa vào ngăn xếp. 3.1.2. Độ phức tạp thuật toán PTIT Độ phức tạp thuật toán DFS(u) phụ thuộc vào phương pháp biểu diễn đồ thị. Độ phức tạp thuật toán DFS(u) theo các dạng biểu diễn đồ thị như sau: Độ phức tạp thuật toán là O(n2) trong trường hợp đồ thị biểu diễn dưới dạng ma trận kề, với n là số đỉnh của đồ thị. Độ phức tạp thuật toán là O(n.m) trong trường hợp đồ thị biểu diễn dưới dạng danh sách cạnh, với n là số đỉnh của đồ thị, m là số cạnh của đồ thị. Độ phức tạp thuật toán là O(max(n, m)) trong trường hợp đồ thị biểu diễn dưới dạng danh sách kề, với n là số đỉnh của đồ thị, m là số cạnh của đồ thị. 32
  33. Bạn đọc tự chứng minh hoặc có thể tham khảo trong các tài liệu [1, 2, 3]. 3.1.3. Kiểm nghiệm thuật toán Ví dụ 1. Kiểm nghiệm thuật toán DFS(1) trên đồ thị gồm 13 đỉnh trong Hình 3.2 dưới đây? 2 6 8 7 1 4 5 3 10 11 9 13 12 Hình 3.2. Đồ thị vô hướng G. Đỉnh bắt đầu Các đỉnh đã duyệt: Các đỉnh chưa duyệt duyệt chuaxet[u]=False chuaxet[u]=True DFS(1) 1 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 DFS(2) 1, 2 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 DFS(4) 1, 2, 4 3, 5, 6, 7, 8, 9, 10, 11, 12, 13 DFS(3) 1,2,4, 3 5, 6, 7, 8, 9, 10, 11, 12, 13 DFS(6) 1,2,4,3, 6 5, 7, 8, 9, 10, 11, 12, 13 DFS(7) 1,2,4,3, 6,7PTIT 5, 8, 9, 10, 11, 12, 13 DFS(8) 1,2,4,3, 6,7,8 5, 9, 10, 11, 12, 13 DFS(10) 1,2,4,3, 6,7,8,10 5, 9, 11, 12, 13 DFS(5) 1,2,4,3, 6,7,8,10,5 9, 11, 12, 13 DFS(9) 1,2,4,3, 6,7,8,10,5,9 11, 12, 13 DFS(13) 1,2,4,3, 6,7,8,10,5,9,13 11, 12 DFS(11) 1,2,4,3, 6,7,8,10,5,9,13,11 12 DFS(11) 1,2,4,3, 6,7,8,10,5,9,13,11,12  Kết quả duyệt: 1, 2, 4, 3, 6, 7, 8, 10, 5, 9, 13, 11, 12 33
  34. Để bạn đọc làm quen với phương pháp kiểm nghiệm thuật toán dựa vào dữ liệu, chúng tôi sử dụng biểu diễn của đồ thị bằng ma trận kề như đã được trình bày trong Chương 2. Việc kiểm nghiệm thuật toán bằng các biểu diễn khác (danh sách cạnh, danh sách kề) xem như những bài tập để bạn đọc tự tìm ra lời giải. 0 1 1 1 0 0 0 0 0 0 0 0 0 Ví dụ 2. Cho đồ thị gồm 13 đỉnh được 1 0 1 1 0 1 0 0 0 0 0 0 0 biểu diễn dưới dạng ma trận kề như hình 1 1 0 1 1 0 0 0 0 0 0 0 0 bên ph ải. Hãy cho biết kết quả thực hiện 1 1 1 0 0 0 1 0 0 0 0 0 0 thuật toán trong Hình 3.1 bắt đầu tại đỉnh 0 0 1 0 0 1 1 1 0 0 0 1 0 u=1? Chỉ rõ trạng thái của ngăn xếp và 0 1 0 0 1 0 1 0 0 0 0 1 0 tập đỉnh được duyệt theo mỗi bước thực 0 0 0 1 1 1 0 1 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 1 0 hiện của thuật toán? 0 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 1 0 1 1 1 0 0 0 0 0 0 0 0 1 1 0 0 1 0 0 0 0 1 1 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 Lời giải. Trạng thái của ngăn xếp và tập đỉnh được duyệt theo thuật toán được thể hiện trong Bảng 3.1 dưới đây. Bảng 3.1. Kiểm nghiệm thuật toán DFS(1). STT Trạng thái stack Các đỉnh được duyệt 1 1 1 2 1, 2 1, 2 3 1, 2, 3 1, 2, 3 4 1, 2, 3, 4 1, 2, 3, 4 5 1, 2, 3, 4, 7 1, 2, 3, 4, 7 6 1, 2, 3, 4, 7, 5 1, 2, 3, 4, 7, 5 7 1, 2, 3, 4, 7, 5, 6 PTIT1, 2, 3, 4, 7, 5, 6 8 1, 2, 3, 4, 7, 5, 6, 12 1, 2, 3, 4, 7, 5, 6, 12 9 1, 2, 3, 4, 7, 5, 6, 12, 8 1, 2, 3, 4, 7, 5, 6, 12, 8 10 1, 2, 3, 4, 7, 5, 6, 12, 10 1, 2, 3, 4, 7, 5, 6, 12, 8, 10 11 1, 2, 3, 4, 7, 5, 6, 12, 10, 9 1, 2, 3, 4, 7, 5, 6, 12, 8, 10, 9 12 1, 2, 3, 4, 7, 5, 6, 12, 10, 9, 11 1, 2, 3, 4, 7, 5, 6, 12, 8, 10, 9, 11 13 1, 2, 3, 4, 7, 5, 6, 12, 10, 9, 11, 13 1, 2, 3, 4, 7, 5, 6, 12, 8, 10, 9, 11, 13 14  Kết quả duyệt DFS(1) = { 1, 2, 3, 4, 7, 5, 6, 12, 8, 10, 9, 11, 13}. 34
  35. Chú ý. Đối với đồ thị vô hướng, nếu DFS(u) = V ta có thể kết luận đồ thị liên thông. Đối với đồ thị có hướng, nếu DFS(u) = V ta có thể kết luận đồ thị liên thông yếu. 3.1.4. Cài đặt thuật toán Thuật toán được cài đặt theo khuôn dạng dữ liệu tổ chức trong file dothi.in được qui ước như được trình bày trong Mục 2.1.3 như sau: Dòng đầu tiên ghi lại số đỉnh của đồ thị; N dòng kế tiếp ghi lại ma trận kề của đồ thị. Hai phần tử khác nhau của ma trận kề được viết cách nhau một vài khoảng trống. Chương trình được thực hiện với các thủ tục như sau: Hàm Init() : đọc dữ liệu theo khuôn dạng từ file dothi.in và thiết lập mảng chuaxet[u] =True (u=1, 2, ,n). Hàm DFS_Dequi : Cài đặt thuật toán DFS(u) bằng đệ qui. Hàm DFS_Stack : Cài đặt thuật toán DFS(u) dựa vào stack. Ví dụ với file dothi.in dưới đây với u = 3 sẽ cho ta kết quả thực hiện chương trình như sau: dothi.in 10 0 1 1 1 0 0 0 0 0 0 1 0 0 1 1 0 0 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 PTIT1 1 0 0 1 0 0 1 0 1 0 0 0 1 0 0 0 0 1 1 0 0 1 0 0 0 0 0 0 0 0 1 0 0 1 1 0 0 0 0 1 0 0 0 1 1 0 0 0 1 0 0 1 1 0 1 0 0 0 0 0 0 1 1 1 0 DFS(3) = 3, 1, 2, 4, 5, 8, 9, 7, 6, 10. 35
  36. #include #include #include #define MAX 50 #define TRUE 1 #define FALSE 0 int A[MAX][MAX], n,chuaxet[MAX]; void Init(void){ int i,j;FILE *fp; fp=fopen("DOTHI.IN","r"); fscanf(fp,"%d",&n); printf("\n So dinh do thi:%d",n); for(i=1; i 0){ s=Stack[dau];dau ; for(t =1;t<=n; t++){ if(chuaxet[t] && A[s][t]){ printf("%3d",t); chuaxet[t] = FALSE; Stack[++dau]=s; Stack[++dau]=t;break; } } } } 36
  37. void main(void){ int u ;clrscr();Init(); cout >u; DFS_Stack(u); //DFS_Dequi(u); getch(); } 3.2. Thuật toán tìm kiếm theo chiều rộng (Breadth First Search) 3.2.1. Biểu diễn thuật toán Để ý rằng, với thuật toán tìm kiếm theo chiều sâu, đỉnh thăm càng muộn sẽ trở thành đỉnh sớm được duyệt xong. Đó là kết quả tất yếu vì các đỉnh thăm được nạp vào stack trong thủ tục đệ qui. Khác với thuật toán tìm kiếm theo chiều sâu, thuật toán tìm kiếm theo chiều rộng thay thế việc sử dụng stack bằng hàng đợi (queue). Trong thủ tục này, đỉnh được nạp vào hàng đợi đầu tiên là u, các đỉnh kề với u là ( v1, v2, . . ., vk) được nạp vào hàng đợi nếu như nó chưa được xét đến. Quá trình duyệt tiếp theo được bắt đầu từ các đỉnh còn có mặt trong hàng đợi. Để ghi nhận trạng thái duyệt các đỉnh của đồ thị, ta cũng vẫn sử dụng mảng chuaxet[] gồm n phần tử thiết lập giá trị ban đầu là TRUE. Nếu đỉnh u của đồ thị đã được duyệt, giá trị chuaxet[u] sẽ nhận giá trị FALSE. Thuật toán dừng khi hàng đợi rỗng. Hình 3.3. dưới đây mô tả chi tiết thuật toán BFS(u). Thuật toán BFS(u): Bước 1(Khởi tạo): Queue = ; Push(Queue,u); chuaxet[u] = False; Bước 2 (Lặp): while (Queue  ) do s = Pop(Queue);PTIT ; for each t Ke(s) do if ( chuaxet[t] ) then Push(Queue, t); chuaxet[t] = False; EndIf ; EndFor ; EndWhile ; Bước 3 (Trả lại kết quả) : Return( ) ; End. Hình 3.3. Thuật toán BFS(u). 37
  38. 3.2.2. Độ phức tạp thuật toán Độ phức tạp thuật toán BFS(u) phụ thuộc vào phương pháp biểu diễn đồ thị. Độ phức tạp thuật toán BFS(u) theo các dạng biểu diễn đồ thị như sau: Độ phức tạp thuật toán là O(n2) trong trường hợp đồ thị biểu diễn dưới dạng ma trận kề, với n là số đỉnh của đồ thị. Độ phức tạp thuật toán là O(n.m) trong trường hợp đồ thị biểu diễn dưới dạng danh sách cạnh, với n là số đỉnh của đồ thị, m là số cạnh của đồ thị. Độ phức tạp thuật toán là O(max(n, m)) trong trường hợp đồ thị biểu diễn dưới dạng danh sách kề, với n là số đỉnh của đồ thị, m là số cạnh của đồ thị. Bạn đọc tự chứng minh hoặc có thể tham khảo trong các tài liệu [1, 2, 3]. 3.2.3. Kiểm nghiệm thuật toán 0 1 1 1 0 0 0 0 0 0 0 0 0 Ví d ụ 3. Cho đồ thị gồm 13 đỉnh được 1 0 1 1 0 1 0 0 0 0 0 0 0 biểu diễn dưới dạng ma trận kề như hình 1 1 0 1 1 0 0 0 0 0 0 0 0 bên ph ải. Hãy cho biết kết quả thực hiện 1 1 1 0 0 0 1 0 0 0 0 0 0 thuật toán trong Hình 3.3 bắt đầu tại đỉnh 0 0 1 0 0 1 1 1 0 0 0 1 0 u=1? Chỉ rõ trạng thái của hàng đợi và tập 0 1 0 0 1 0 1 0 0 0 0 1 0 đỉnh được duyệt theo mỗi bước thực hiện 0 0 0 1 1 1 0 1 0 0 0 0 0 của thuật toán? 0 0 0 0 1 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 1 0 1 1 1 0 0 0 0 0 0 0 0 1 1 0 0 1 0 0 0 0 1 1 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 Lời giải. Trạng thái của hàng đợi và tập đỉnh được duyệt theo thuật toán được thể hiện trong Bảng 3.2 dưới đây. Bảng 3.2.PTIT Kiểm nghiệm thuật toán BFS(1). STT Trạng thái Queue Các đỉnh được duyệt 1 1  2 2, 3, 4 1 3 3, 4, 6 1, 2 4 4, 6, 5 1, 2, 3 5 6, 5, 7 1, 2, 3, 4 6 5, 7, 12 1, 2, 3, 4, 6 7 7, 12, 8 1, 2, 3, 4, 6, 5 8 12, 8 1, 2, 3, 4, 6, 5, 7 9 8, 10 1, 2, 3, 4, 6, 5, 7, 12 10 10 1, 2, 3, 4, 6, 5, 7, 12, 8 38
  39. 11 9, 11, 13 1, 2, 3, 4, 6, 5, 7, 12, 8,10 12 11, 13 1, 2, 3, 4, 6, 5, 7, 12, 8,10, 9 13 13 1, 2, 3, 4, 6, 5, 7, 12, 8,10, 9, 11 14  1, 2, 3, 4, 6, 5, 7, 12, 8,10, 9, 11, 13 Kết quả duyệt BFS(1) = { 1, 2, 3, 4, 6, 5, 7, 12, 8,10, 9, 11, 13}. Chú ý. Đối với đồ thị vô hướng, nếu BFS(u) = V ta có thể kết luận đồ thị liên thông. Đối với đồ thị có hướng, nếu BFS(u) = V ta có thể kết luận đồ thị liên thông yếu. 3.2.4. Cài đặt thuật toán Thuật toán được cài đặt theo khuôn dạng dữ liệu tổ chức trong file dothi.in được qui ước như được trình bày trong Mục 2.1.3 như sau: Dòng đầu tiên ghi lại số đỉnh của đồ thị; N dòng kế tiếp ghi lại ma trận kề của đồ thị. Hai phần tử khác nhau của ma trận kề được viết cách nhau một vài khoảng trống. Chương trình được thực hiện với các thủ tục như sau: Hàm Init() : đọc dữ liệu theo khuôn dạng từ file dothi.in và thiết lập mảng chuaxet[u] =True (u=1, 2, ,n). Hàm BFS_Dequi : Cài đặt thuật toán BFS(u) bằng hàng đợi. Ví dụ với file dothi.in dưới đây với u = 3 sẽ cho ta kết quả thực hiện chương trình như sau: dothi.in 10 0 1 1 1 0 0 0 0 0 0 1 0 0 1 PTIT1 0 0 0 0 0 1 0 0 1 0 1 0 0 0 0 1 1 1 0 1 1 0 0 1 0 0 1 0 1 0 0 0 1 0 0 0 0 1 1 0 0 1 0 0 0 0 0 0 0 0 1 0 0 1 1 0 0 0 0 1 0 0 0 1 1 0 0 0 1 0 0 1 1 0 1 0 0 0 0 0 0 1 1 1 0 BFS(3) = 3, 1, 4, 6, 2, 5, 9, 7, 8, 10. 39
  40. #include #include #define MAX 50 #define TRUE 1 #define FALSE 0 int A[MAX][MAX], n,chuaxet[MAX];FILE *fp; void Init(void){ int i,j; fp= fopen("dothi.in","r"); fscanf(fp,"%d",&n); printf("\n So dinh do thi:%d",n); for(i=1; i<=n; i++){ printf("\n");chuaxet[i]=TRUE; for(j=1; j<=n; j++){ fscanf(fp,"%d",&A[i][j]); printf("%3d",A[i][j]); } } } void BFS(int u){ int queue[MAX], low=1, high=1, v; queue[low]=u;chuaxet[u]=FALSE; printf("\n Ket qua:"); while(low<=high){ u = queue[low];low=low+1; printf("%3d", u); for(v=1; v<=n; v++){ if(A[u][v] && chuaxet[v]){ high = high+1; queue[high]=v; chuaxet[v]=FALSE; } PTIT } } } void main(void){ int u; Init(); printf("\n Dinh bat dau duyet:"); scanf("%d",&u); BFS(u); } 40
  41. 3.3. Ứng dụng của thuật toán DFS và BFS Có rất nhiều ứng dụng khác nhau của thuật toán DFS và BFS trên đồ thị. Trong khuôn khổ của giáo trình này, chúng tôi đề cập đến một vài ứng dụng cơ bản. Những ứng dụng cụ thể hơn bạn đọc có thể tìm thấy rất nhiều trong các tài liệu khác nhau hoặc Internet. Những ứng dụng cơ bản của thuật toán DFS và BFS được đề cập bao gồm: o Duyệt tất cả các đỉnh của đồ thị. o Duyệt tất cả các thành phần liên thông của đồ thị. o Tìm đường đi từ đỉnh s đến đỉnh t trên đồ thị. o Duyệt các đỉnh trụ trên đồ thị vô hướng. o Duyệt các đỉnh trụ trên đồ thị vô hướng. o Duyệt các cạnh cầu trên đồ thị vô hướng. o Định chiều đồ thị vô hướng. o Duyệt các đỉnh rẽ nhánh của cặp đỉnh s, t. o Xác định tính liên thông mạnh trên đồ thị có hướng. o Xác định tính liên thông yếu trên đồ thị có hướng. o Thuật toán tìm kiếm theo chiều rộng trên đồ thị. o Xây dựng cây khung của đồ thị vô hướng liên thông 3.3.1. Xác định thành phần liên thông của đồ thị a) Đặt bài toán Cho đồ thị đồ thị vô hướng G= , trong đó V là tập đỉnh, E là tập cạnh. Bài toán đặt ra là xác định các thành phần liên thông của G = ? b) Mô tả thuật toán PTIT Một đồ thị có thể liên thông hoặc không liên thông. Nếu đồ thị liên thông thì số thành phần liên thông của nó là 1. Điều này tương đương với phép duyệt theo thủ tục DFS(u) hoặc BFS(u) được gọi đến đúng một lần. Nói cách khác, DFS(u)=V và BFS(u)=V. Nếu đồ thị không liên thông (số thành phần liên thông lớn hơn 1) chúng ta có thể tách chúng thành những đồ thị con liên thông. Điều này cũng có nghĩa là trong phép duyệt đồ thị, số thành phần liên thông của nó bằng đúng số lần gọi tới thủ tục DFS() hoặc BFS(). Để xác định số các thành phần liên thông của đồ thị, chúng ta sử dụng thêm biến solt để nghi nhận các đỉnh cùng một thành phần liên thông. Khi đó, thuật toán xác định các thành phần liên thông của đồ thị được mô tả trong Hình 3.4. 41
  42. Thuật toán Duyet-TPLT: Bước 1 (Khởi tạo): solt = 0; //Khởi tạo số thành phần liên thông ban đầu là 0 Bước 2 (Lặp): for ( u =1; u n; u++) do //lặp trên tập đỉnh if (chuaxet[u] ) then solt = solt + 1; //Ghi nhận số thành phần liên thông ; BFS (u); //DFS(u); // endif; endfor; Bước 3 (Trả lại kết quả): Return(solt); end. Hình 3.4. Thuật toán duyệt các thành phần liên thông của đồ thị. c) Kiểm nghiệm thuật toán Ví dụ ta cần kiểm nghiệm thuật toán trên Hình 3.4 cho đồ thị được biểu diễn dưới dạng ma trận kề như dưới đây. 0 0 1 0 1 0 1 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 1 0 0 0 1 0 1 0 0 0 1 0 0 0 1 0 0 0 1 0 1 0 1 0 0 0 1 0 1 0 0 0 1 0 1 0 1 0 1 0 1 0 1 0 0 0 1 0 1 0 0 0 1 0 1 0 1 0 0 0 1 0 0 0 0 0 0 0 1PTIT 0 1 0 0 0 1 0 1 0 0 0 0 0 1 0 1 0 0 0 1 0 1 0 0 0 1 0 1 0 1 0 0 0 1 0 0 0 1 0 1 0 0 0 1 0 0 0 1 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 1 0 0 0 1 0 1 0 0 Thực hiện thuật toán DFS và BFS như được mô tả ở trên ta nhận được : Thành phần liên thông 1: BFS(1) = { 1, 3, 5, 7, 9, 11, 13}. Thành phần liên thông 2: BFS(2) = {2, 4, 6, 8, 10, 12}. 42
  43. d) Cài đặt thuật toán Chương trình duyệt các thành phần liên thông của đồ thị được cài đặt theo khuôn dạng dữ liệu biểu diễn dưới dạng ma trận kề trong Mục 2.3.1 với các thủ tục sau: Hàm Init() : Đọc dữ liệu theo khuôn dạng và khởi đầu mảng chuaxet[u] = True (1 i n). Hàm BFS (u), DFS(u) : Hai thuật toán duyệt theo chiều rộng và duyệt theo chiều sâu được sử dụng để xác định các thành phần liên thông. #include #include #define MAX 50 #define TRUE 1 #define FALSE 0 int A[MAX][MAX], n,chuaxet[MAX], solt=0; void Init(void){ int i,j;FILE *fp; fp=fopen("dothi.in","r"); fscanf(fp,"%d",&n); printf("\n So dinh do thi:%d",n); for(i=1; i<=n; i++){ printf("\n");chuaxet[i]=TRUE; for(j=1; j<=n; j++){ fscanf(fp,"%d",&A[i][j]); printf("%3d",A[i][j]); } } } void BFS(int u){ int queue[MAX],low=1,high=1, s,t; queue[low]=u;chuaxet[u]=FALSE; while(low<=high){PTIT s = queue[low];low=low+1; printf("%3d", s); for(t=1; t<=n; t++){ if(A[s][t] && chuaxet[t]){ high = high+1; queue[high]=t; chuaxet[t]=FALSE; } } } } 43
  44. void DFS(int u){ int v;printf("%3d",u); chuaxet[u]=FALSE; for(v=1; v (vô hướng hoặc có hướng), trong đó V là tập đỉnh, E là tập cạnh. Bài toán đặt ra là hãy tìm đường đi từ đỉnh s V đến đỉnh t V? b) Mô tả thuật toán Cho đồ thị G = , s, t là hai đỉnh thuộc V. Khi đó, dễ dàng nhận thấy, nếu t DFS(s) hoặc t BFS(s) thì ta có thể kết luận có đường đi từ s đến t trên đồ thị. Nếu t DFS(s) hoặc t BFS(s) thì ta có thể kết luận không có đường đi từ s đến t trên đồ thị. Vấn đề còn lại là ta ghi nhận thế nào đường đi từ s đến t? Để ghi nhận đường đi từPTIT s đến t dựa vào hai thuật toán DFS(s) hoặc BFS(s) ta sử dụng một mảng truoc[] gồm n phần tử (n=|V|). Khởi tạo ban đầu truoc[u]=0 với mọi u = 1, 2, , n. Trong quá trình thực hiện hai thuật toán DFS (s) và BFS(s), mỗi khi ta đưa đỉnh v Ke(s) vào ngăn xếp (trong trường hợp ta sử dụng thuật toán DFS) hoặc hàng đợi(trong trường hợp ta sử dụng thuật toán DFS) ta ghi nhận truoc[v] = s. Điều này có nghĩa, để đi được đến v ta phải qua đỉnh s. Khi hai thuật toán DFS và BFS duyệt đến đỉnh t thì truoc[t] sẽ nhận giá trị là một đỉnh nào đó thuộc V hay t DFS(s) hoặc t BFS(s). Trong trường hợp hai thủ tục DFS và BFS không duyệt được đến đỉnh t, khi đó truoc[t] =0 và ta kết luận không có đường đi từ s đến t. Hình 3.5 dưới đây mô tả thuật toán tìm đường đi từ đỉnh s đến đỉnh t trên đồ thị bằng thuât toán DFS. Hình 3.6 dưới đây mô tả thuật toán tìm đường đi từ đỉnh s đến đỉnh t trên đồ thị bằng thuât toán BFS. Hình 3.7 dưới đây mô tả thuật toán ghi nhận đường đi từ đỉnh s đến đỉnh t trên đồ thị. 44
  45. Thuật toán DFS(s): Begin Bước 1 (Khởi tạo): stack = ; //Khởi tạo stack là  Push(stack, s); //Đưa đỉnh s vào ngăn xếp chuaxet[s] = False; //Xác nhận đỉnh u đã duyệt Bước 2 (Lặp) : while ( stack  ) do u = Pop(stack); //Loại đỉnh ở đầu ngăn xếp for each v Ke(u) do //Lấy mỗi đỉnh u Ke(v) if ( chuaxet[v] ) then //Nếu v đúng là chưa duyệt chuaxet[v] = False; // Xác nhận đỉnh v đã duyệt Push(stack, u);//Đưa u vào stack Push(stack, v); //Đưa v vào stack truoc[v] = u; //Ghi nhận truoc[v] là u break; //Chỉ lấy một đỉnh t EndIf; EndFor; EndWhile; Bước 3 (Trả lại kết quả): Return( ); End. Hình 3.5. Thuật toán DFS tìm đường đi từ s đến t. Thuật toán BFS(s): Bước 1(Khởi tạo): Queue = ; Push(Queue,s); chuaxet[s] = False; Bước 2 (Lặp): while (Queue  ) do u = Pop(Queue);PTIT for each v Ke(u) do if ( chuaxet[v] ) then Push(Queue, v);chuaxet[v]=False;truoc[v]=u; EndIf ; EndFor ; EndWhile ; Bước 3 (Trả lại kết quả) : Return( ) ; End. Hình 3.6. Thuật toán BFS tìm đường đi từ s đến t. 45
  46. Thuật toán Ghi-Nhan-Duong-Di (s, t) { if ( truoc[t] == 0 ) { ; } else { ; //Đưa ra trước đỉnh t u = truoc[t]; //u là đỉnh trước khi đến được t while (u s ) { //Lặp nếu u chưa phải là s ; //Đưa ra đỉnh u u = truoc[u]; // Lần ngược lại đỉnh truoc[u]. } ; } } Hình 3.7. Thủ tục ghi nhận đường đi từ s đến t c) Kiểm nghiệm thuật toán Giả sử ta cần xác định đường đi từ đỉnh 1 đến đỉnh 13 trên đồ thị được biểu diễn dưới dạng ma trận kề. Khi đó, thứ tự các bước thực hiện theo thuật toán DFS được thể hiện trong Bảng 3.3, thứ tự các bước thực hiện theo thuật toán BFS được thể hiện trong Bảng 3.4. Bảng 3.3. Kiểm nghiệm thuật toán DFS(1). STT Trạng thái stack Truoc[s]=? 1 1 0 2 1, 2 truoc[2] =1 3 1, 2, 3 PTITtruoc[3] = 2 4 1, 2, 3, 4 truoc[4] =3 5 1, 2, 3, 4, 7 truoc[7] =4 6 1, 2, 3, 4, 7, 5 truoc[5] =7 7 1, 2, 3, 4, 7, 5, 6 truoc[6] =5 8 1, 2, 3, 4, 7, 5, 6, 12 truoc[12] =6 9 1, 2, 3, 4, 7, 5, 6, 12, 8 truoc[8] =12 10 1, 2, 3, 4, 7, 5, 6, 12, 10 truoc[10] =12 11 1, 2, 3, 4, 7, 5, 6, 12, 10, 9 truoc[9] =10 46
  47. 12 1, 2, 3, 4, 7, 5, 6, 12, 10, 9, 11 truoc[11] =9 13 1, 2, 3, 4, 7, 5, 6, 12, 10, 9, 11, 13 truoc[13] =11 14  Kết quả đường đi từ đỉnh 1 đến đỉnh 13: 13->11-9-10-12-6-5-7-4-3-2-1. Bảng 3.4. Kiểm nghiệm thuật toán BFS(1). STT Trạng thái Queue Truoc[s]=? 1 1 truoc[1]=0 2 2, 3, 4 truoc[2]=1; truoc[3]=1; truoc[4]=1; 3 3, 4, 6 truoc[6]= 2 4 4, 6, 5 truoc[5]=3 5 6, 5, 7 truoc[7]= 4 6 5, 7, 12 truoc[12]=6 7 7, 12, 8 truoc[8]=12 8 12, 8 9 8, 10 truoc[10]=12; 10 10 11 9, 11, 13 truoc[9]=10; truoc[11]=10; truoc[13]=10; 12 11, 13 13 13 14  Kết quả đường đi: 13-10-12-6-2-1. Chú ý. Đường đi từ đỉnh s đến đỉnh t theo thuật toán BFS đi qua ít nhất các cạnh của đồ thị (có độ dài nhỏ nhất). d) Cài đặt thuật toán #include PTIT #include #include #define MAX 50 #define TRUE 1 #define FALSE 0 int A[MAX][MAX], n,chuaxet[MAX], truoc[MAX], s, t; void Init(void){//Đọc dữ liệu và khởi đầu các biến int i,j;FILE *fp; fp=fopen("dothi.in","r"); fscanf(fp,"%d",&n); printf("\n So dinh do thi:%d",n); for(i=1; i<=n; i++){ 47
  48. printf("\n");chuaxet[i]=TRUE;truoc[i]=0; for(j=1; j<=n; j++){ fscanf(fp,"%d",&A[i][j]); printf("%3d",A[i][j]); } } } void DFS(int u){//Thuật toán DFS int v; printf("%3d",u);chuaxet[u]=FALSE; for(v=1; v<=n; v++){ if(A[u][v] && chuaxet[v]){ truoc[v]=u;DFS(v); } } } void BFS(int i){//Thuật toán BFS int queue[MAX], low=1, high=1, u, v; queue[low]=i;chuaxet[i]=FALSE; while(low<=high){ u = queue[low];low=low+1; for(v=1; v<=n; v++){ if(A[u][v] && chuaxet[v]){ high = high+1;queue[high]=v; truoc[v]=u; chuaxet[v]=FALSE; } } } } void Duongdi (void){ if (truoc[t]==0){PTIT printf("\n Khong ton tai duong di"); getch(); return; } printf("\n Duong di:"); int j = t;printf("%3d<=",j); while(truoc[j]!=s){ printf("%3d<=",truoc[j]);j=truoc[j]; } printf("%3d",s); getch(); } 48
  49. void main (void){ Init(); printf("\n Dinh dau:");scanf("%d",&s); printf("\n Dinh cuoi:");scanf("%d",&t); DFS(s); //BFS(s); Duongdi (); } 3.3.3. Tính liên thông mạnh trên đồ thị có hướng a) Đặt bài toán Đồ thị có hướng G= liên thông mạnh nếu giữa hai đỉnh bất kỳ của nó đều tồn tại đường đi. Cho trước đồ thị có hướng G = . Nhiệm vụ của ta là kiểm tra xem G có liên thông mạnh hay không? b) Mô tả thuật toán Đối với đồ thị vô hướng, nếu hai thủ tục DFS(u) = V hoặc BFS(u) = V thì ta kết luận đồ thị vô hướng liên thông. Đối với đồ thị có hướng, nếu DFS(u)=V hoặc BFS(u)=V thì ta mới chỉ có kết luận có đường đi từ u đến tất cả các đỉnh còn lại của đồ thị. Nhiệm vụ của ta là phải kiểm tra DFS(u)=V hoặc BFS(u)=V với mọi u V. Hình 3.8 dưới đây mô tả chi tiết thuật toán kiểm tra tính liên thông mạnh của đồ thị. Boolean Strong-Connective ( G = ) { ReInit(); //u V: chuaxet[u] = True; for each u V do {//Lấy mỗi đỉnh thuộc V if (DFS(u) V ) then //Nếu DFS(u) V hoặc BFS(u) V return(False); //Đồ thị không liên thông mạnh endif; ReInit();//KhởiPTIT tạo lại các mảng chuaxet[] endfor; return(True);//Đồ thị liên thông mạnh } Hình 3.8. Thuật toán kiểm tra tính liên thông mạnh. c) Kiểm nghiệm thuật toán Giả sử ta cần xác định đồ thị có hướng G = được biểu diễn dưới dạng ma trận kề dưới đây có liên thông mạnh hay không? Khi đó các bước thực hiện theo thuật toán Hình 3.8 được thực hiện theo Bảng 3.5 dưới đây. 49
  50. 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 Bảng 3.5. Kiểm nghiệm thuật toán kiểm tra tính liên thông mạnh. Đỉnh u V DFS(u)=?//BFS(u)=? DFS(u) =V? 1 V DFS(1) = 1, 6, 10, 2, 3, 9, 5, 7, 11, 8, 4, 12, 13 Yes 2 V DFS(2) = 2, 3, 9, 5, 7, 11, 8, 4, 1, 6, 10, 12, 13 Yes 3 V DFS(3) = 3, 9, 5, 7, 11, 2, 8, 4, 1, 6, 10, 12, 13 Yes 4 V DFS(4) = 4, 1, 6, 10, 2, 3, 9, 5, 7, 11, 8, 12, 13 Yes 5 V DFS(5) = 5, 7, 11, 2, 3, 9, 13, 8, 4, 1, 6, 10, 12 Yes 6 V DFS(6) = 6, 10, 2, 3, 9, 5, 7, 11, 8, 4, 1, 12, 13 Yes 7 V DFS(7) = 7, 11, 2, 3, 9, 5, 13, 8, 4, 1, 6, 10, 12 Yes 8 V DFS(8) = 8, 4, 1, 6, 10, 2, 3, 9, 5, 7, 11, 13, 12 Yes 9 V DFS(8) = 9, 5, 7, 11, 2, 3, 13, 8, 4, 1, 6, 10, 12 Yes 10 V DFS(10) = 10,PTIT 2, 3, 9, 5, 7, 11, 8, 4, 1, 6, 12, 13 Yes 11 V DFS(11) = 11, 2, 3, 9, 5, 7, 13, 8, 4, 1, 6, 10, 12 Yes 12 V DFS(12) = 12, 4, 1, 6, 10, 2, 3, 9, 5, 7, 11, 8, 13 Yes 13 V DFS(13) = 13, 9, 5, 7, 11, 2, 3, 8, 4, 1, 6, 10, 12 Yes Cột ngoài cùng của Bảng có DFS(u) = V với mọi u V nên ta kết luận G liên thông mạnh. Nếu tại một hàng nào đó có DFS(u) V thì ta kết luận đồ thị không liên thông mạnh và không cần phải kiểm tra tiếp các đỉnh còn lại. 50
  51. d) Cài đặt thuật toán Thuật toán được cài đặt theo khuôn dạng đồ thị được qui ước trong Mục 2.3.1 với các thủ tục sau: Thủ tục Read-Data() : Đọc ma trận kề biểu diễn đồ thị trong file dothi.in. Thủ tục ReInit() : Khởi tạo lại giá trị cho mảng chuaxet[]. Thủ tục DFS(u) : Thuật toán DFS bắt đầu tại đỉnh u. Thủ tục BFS(u) : Thuật toán BFS bắt đầu tại đỉnh u. Hàm Strong-Conective(): Kiểm tra tính liên thông mạnh của đồ thị. Chương trình kiểm tra tính liên thông mạnh của đồ thị được thể hiện như dưới đây. #include #include #include #define MAX 50 #define TRUE 1 #define FALSE 0 int A[MAX][MAX], n,chuaxet[MAX], solt=0; //Doc du lieu void Read_Data(void){ int i,j;FILE *fp; fp=fopen("dothi.IN","r"); fscanf(fp,"%d",&n); printf("\n So dinh do thi:%d",n); for(i=1; i<=n; i++){ printf("\n"); for(j=1; j<=n; j++){ fscanf(fp,"%d",&A[i][j]); printf("%3d",A[i][j]); } PTIT } } //Thuat toan BFS void BFS(int u){ int queue[MAX],low=1,high=1, s,t; queue[low]=u;chuaxet[u]=FALSE; while(low<=high){ s = queue[low];low=low+1; //printf("%3d", s); for(t=1; t<=n; t++){ if(A[s][t] && chuaxet[t]){ high = high+1; 51
  52. queue[high]=t; chuaxet[t]=FALSE; } } } } //Thuat toan DFS void DFS(int u){ int v;//printf("%3d",u); chuaxet[u]=FALSE; for(v=1; v 1? int Test_So_Lien_Thong(void) { for(int u=1; u<=n; u++) if(chuaxet[u]) return(1); return(0); } //Kiem tra tinh lien thong manh int Strong_Conective (void) { Read_Data(); ReInit(); for (int u=1; u<=n; u++){ chuaxet[u]=FALSE; if(u==1)PTIT DFS(2);//BFS(2); esle DFS(1); //BFS(1); if(Test_So_Lien_Thong()) return(0); ReInit(); } return(1); } void main(void){ if(Test_LT_Manh()) printf("\n Do thi lien thong manh"); else printf("\n Do thi khong lien thong manh"); } 52
  53. 3.3.4. Duyệt các đỉnh trụ a) Đặt bài toán Cho đồ thị vô hướng liên thông G = . Đỉnh u V được gọi trụ nếu loại bỏ đỉnh u cùng với các cạnh nối với u làm tăng thành phần liên thông của đồ thị. Bài toán đặt ra là xây dựng thuật toán duyệt các đỉnh trụ của đồ thị vô hướng G = cho trước? b) Mô tả thuật toán Không hạn chế tính tổng quát của bài toán ta có thể giả thiết đồ thị đã cho ban đầu là liên thông. Trong trường hợp đồ thị không liên thông, bài toán duyệt trụ thực chất giải quyết cho mỗi thành phần liên thông của đồ thị. Đối với đồ thị được biểu diễn dưới dạng ma trận kề, việc loại bỏ đỉnh u cùng với các cạnh nối với u tương ứng với việc loại bỏ hàng u và cột u tương ứng trong ma trận kề. Để thực hiện việc này trong các thủ tục DFS() hoặc BFS() ta chỉ cần thiết lập giá trị chuaxet[u] = False. Quá trình duyệt sẽ được thực hiện tại một đỉnh bất kỳ v u. Nếu DFS(v) = V\{u} hoặc BFS(v) = V\{u} thì đồ thị mới nhận được cũng chỉ có 1 thành phần liên thông và ta kết luận v không là trụ. Trường hợp DFS(v) V\{u} hoặc BFS(v) V\{u} thì v chính là trụ vì số thành phần liên thông của đồ thị đã tăng lên. Thuật toán duyệt các đỉnh trụ của đồ thị được mô tả chi tiết trong Hình 3.9. Thuật toán Duyet-Tru ( G = ) { ReInit(); //u V: chuaxet[u] = True; for each v V do {//Lấy mỗi đỉnh u tập đỉnh V chuaxet[v] = False; //Cấm DFS hoặc BFS duyệt vào đỉnh v if (DFS(u) V\{v} ) then //Duyệt DFS hoặc BFS tại đỉnh u v ; endif ; ReInit();//Khởi tạo lại các mảng chuaxet[] endfor; PTIT } Hình 3.9. Thuật toán duyệt các đỉnh trụ của đồ thị. c) Kiểm nghiệm thuật toán Giả sử ta cần xác định đồ thị vô hướng G = được biểu diễn dưới dạng ma trận kề dưới đây có những đỉnh trụ nào? Khi đó các bước thực hiện theo thuật toán Hình 3.9 được thực hiện theo Bảng 3.6 dưới đây. 53
  54. 0 1 1 1 0 0 0 0 0 0 0 0 0 1 0 1 1 0 0 0 0 0 0 0 0 0 1 1 0 1 1 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 1 1 1 0 0 0 0 0 0 0 0 1 0 1 0 1 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 0 1 0 1 0 1 0 0 0 0 0 0 0 0 1 1 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 1 1 1 0 0 0 0 0 0 0 0 0 1 0 1 1 0 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 0 1 1 1 0 Bảng 3.6. Kiểm nghiệm thuật toán duyệt các đỉnh trụ của đồ thị. Đỉnh v V DFS(u)=?//BFS(u)=? u v. DFS(u) V\v? 1 V DFS(2) = 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No 2 V DFS(1) = 1, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No 3 V DFS(1) = 1, 2, 4 Yes 4 V DFS(1) = 1, 2, 3, 5, 6, 7, 8, 9, 10, 11, 12, 13 No 5 V DFS(1) = 1, 2, 3, 4 Yes 6 V DFS(1) = 1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 12, 13 No 7 V DFS(1) = 1, 2, 3, 4, 5, 6, 9, 8, 10, 11, 12, 13 No 8 V DFS(1) = 1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 13 No 9 V DFS(1) = 1, 2, 3, 4, 5, 6, 7, 8 Yes 10 V DFS(1) = 1,PTIT 2, 3, 4, 5, 6, 7, 8, 9 Yes 11 V DFS(1) = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13 No 12 V DFS(1) = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13 No 13 V DFS(1) = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13 No Đỉnh ở hàng u có giá trị cột số 3 là Yes là những đỉnh trụ, các đỉnh có giá trị No không là đỉnh trụ. d) Cài đặt thuật toán Thuật toán được cài đặt theo khuôn dạng đồ thị được qui ước trong Mục 2.3.1 với các thủ tục sau: 54
  55. Thủ tục Read-Data() : Đọc ma trận kề biểu diễn đồ thị trong file dothi.in. Thủ tục ReInit() : Khởi tạo lại giá trị cho mảng chuaxet[]. Thủ tục DFS(u) : Thuật toán DFS bắt đầu tại đỉnh u. Thủ tục BFS(u) : Thuật toán BFS bắt đầu tại đỉnh u. Văn bản chương trình được thể hiện như dưới đây. #include #include #include #define MAX 50 #define TRUE 1 #define FALSE 0 int A[MAX][MAX], n,chuaxet[MAX], solt=0; //Doc du lieu void Read_Data(void){ int i,j;FILE *fp; fp=fopen("dothi.IN","r"); fscanf(fp,"%d",&n); for(i=1; i<=n; i++){ printf("\n"); for(j=1; j<=n; j++){ fscanf(fp,"%d",&A[i][j]); } } } //Thuat toan BFS void BFS(int u){ int queue[MAX],low=1,high=1, s,t; queue[low]=u;chuaxet[u]=FALSE; while(low<=high){PTIT s = queue[low];low=low+1; //printf("%3d", s); for(t=1; t<=n; t++){ if(A[s][t] && chuaxet[t]){ high = high+1; queue[high]=t; chuaxet[t]=FALSE; } } } } 55
  56. //Thuat toan DFS void DFS(int u){ int v;//printf("%3d",u); chuaxet[u]=FALSE; for(v=1; v 1? int Test_So_Lien_Thong(void) { for(int u=1; u . Cạnh e E được gọi là cầu nếu loại bỏ e làm tăng thành phần liên thông của đồ thị. Bài toán đặt ra là cho trước đồ thị vô hướng G = , hãy liệt kê tất cả các cạnh cầu của đồ thị. b) Mô tả thuật toán Không hạn chế tính tổng quát của bài toán ta cũng giả thiết đồ thị G= đã cho là liên thông. Trong trường hợp đồ thị không liên thông, bài toán duyệt cầu thực hiện trên mỗi thành phần liên thông của đồ thị. 56
  57. Trong trường hợp đồ thị được biểu diễn dưới dạng ma trận kề, ehi loại bỏ cạnh e=(u,v) E ra khỏi đồ thị ta thực hiện bằng cách cho các phần tử A[u][v]=0 và A[v][u]=0 (A[][] là ma trận kề biểu diễn đồ thị G). Đối với đồ thị được biểu diễn dưới dạng danh sách kề, danh sách kề của đỉnh u ta bớt đi đỉnh v và danh sách kề của đỉnh v ta bớt đi đỉnh u ( Ke(u) = Ke(u)\{v}, Ke(v) = Ke(v)\{u}). Thuật toán duyệt các cạnh cầu của đồ thị được ô tả chi tiết trong Hình 3.10. Thuật toán Duyet-Cau ( G = ) { ReInit(); //u V: chuaxet[u] = True; for each e E do {//Lấy mỗi cạnh thuộc đồ thị E = E\{e}; //Loại bỏ cạnh e ra khỏi đồ thị if (DFS(1) V ) then //Nếu tăng thành phần liên thông của đồ thị ; endif ; E = E{e} ; //Hoàn trả lại cạnh e ReInit();//Khởi tạo lại các mảng chuaxet[] endfor; } Hình 3.10. Thuật toán duyệt các cạnh cầu của đồ thị. c) Kiểm nghiệm thuật toán Giả sử ta cần xác định đồ thị vô hướng G = được biểu diễn dưới dạng ma trận kề dưới đây có những cạnh cầu? Khi đó các bước thực hiện theo thuật toán Hình 3.10 được thực hiện theo Bảng 3.7 dưới đây. 0 1 1 1 0 0 0 0 0 0 0 0 0 1 0 1 1 0 0 0 0 0 0 0 0 0 1 1 0 1 1 0 0 0 0 0 0 0 0 1 1 PTIT1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 1 1 1 0 0 0 0 0 0 0 0 1 0 1 0 1 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 0 1 0 1 0 1 0 0 0 0 0 0 0 0 1 1 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 1 1 1 0 0 0 0 0 0 0 0 0 1 0 1 1 0 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 0 1 1 1 0 57
  58. Bảng 3.7. Kiểm nghiệm thuật toán duyệt các cạnh cầu của đồ thị. Cạnh e E DFS(u)=?//BFS(u)=? DFS(u) V? (1,2) E DFS(1) = 1, 3, 2, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No (1,3) E DFS(1) = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No (1,4) E DFS(1) = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No (2,3) E DFS(1) = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No (2,4) E DFS(1) = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No (3,4) E DFS(1) = 1, 2, 3, 5, 6, 7, 8, 9, 10, 11, 12, 13, 4 No (3,5) E DFS(1) = 1, 2, 3, 4, 5 Yes (5,6) E DFS(1) = 1, 2, 3, 4, 5, 7, 6, 8, 9, 10, 11, 12, 13 No (5,7) E DFS(1) = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No (5,8) E DFS(1) = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No (5,9) E DFS(1) = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No (6,7) E DFS(1) = 1, 2, 3, 4, 5, 6, 9, 8, 7, 10, 11, 12, 13 No (6,9) E DFS(1) = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No (7,8) E DFS(1) =1, 2, 3, 4, 5, 6, 7, 9, 8, 10, 11, 12, 13 No (8,9) E DFS(1) =1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No (9,10) E DFS(1) =1, 2, 3, 4, 5, 6, 7, 8, 9 Yes (10,11) E DFS(1) =1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 11, 13 No (10,12) E DFS(1) =1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No (10,13) E DFS(1) =1,PTIT 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No (11,12) E DFS(1) =1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 12 No (11,13) E DFS(1) =1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No (12,13) E DFS(1) =1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 No Kết luận: cạnh (35), (9,10) là cầu d) Cài đặt thuật toán Thuật toán được cài đặt theo khuôn dạng đồ thị được qui ước trong Mục 2.3.1 với các thủ tục sau: 58
  59. Thủ tục Read-Data() : Đọc ma trận kề biểu diễn đồ thị trong file dothi.in. Thủ tục ReInit() : Khởi tạo lại giá trị cho mảng chuaxet[]. Thủ tục DFS(u) : Thuật toán DFS bắt đầu tại đỉnh u. Thủ tục BFS(u) : Thuật toán BFS bắt đầu tại đỉnh u. Chương trình kiểm tra tính liên thông mạnh của đồ thị được thể hiện như dưới đây. #include #include #include #define MAX 50 #define TRUE 1 #define FALSE 0 int A[MAX][MAX], n,chuaxet[MAX], solt=0; //Doc du lieu void Read_Data(void){ int i,j;FILE *fp; fp=fopen("dothi.IN","r"); fscanf(fp,"%d",&n); for(i=1; i<=n; i++){ printf("\n"); for(j=1; j<=n; j++){ fscanf(fp,"%d",&A[i][j]); } } } //Thuat toan BFS void BFS(int u){ int queue[MAX],low=1,high=1, s,t; queue[low]=u;chuaxet[u]=FALSE; while(low<=high){PTIT s = queue[low];low=low+1; //printf("%3d", s); for(t=1; t<=n; t++){ if(A[s][t] && chuaxet[t]){ high = high+1; queue[high]=t; chuaxet[t]=FALSE; } } } } 59
  60. //Thuat toan DFS void DFS(int u){ int v;//printf("%3d",u); chuaxet[u]=FALSE; for(v=1; v 1? int Test_So_Lien_Thong(void) { for(int u=1; u<=n; u++) if(chuaxet[u]) return(1); return(0); } //Duyệt cạnh cầu void main(void) { Read_Data(); ReInit(); for (int u=1; u<n; u++){ for(int v=u+1;v<=n; v++){ if(A[u][v]) { //Neu (u,v) la mot canh A[u][v]=0; A[v][u]=0;//Loai canh DFS(1);//BFS(1); if(Test_So_Lien_Thong()) printf("\n Canh %d%5d ",u, v); A[u][v]=1; A[v][u]=1; PTITReInit();//Khoi tao lai mang chuaxet } } } } 60
  61. 3.4. Một số bài toán quan trọng khác 2.4.1. Duyệt các thành phần liên thông mạnh của đồ thị Đối với đồ thị có hướng người ta quan tâm đến việc duyệt các thành phần liên thông mạnh của đồ thị. Mỗi thành phần liên thông mạnh của đồ thị là một đồ thị con của G mà giữa hai đỉnh bất kỳ của đồ thị con đều có đường đi. Bài toán đặt ra là hãy liệt kê tất cả các thành phần liên thông mạnh của đồ thị có hướng G= . Ví dụ với đồ thị trong Hình 3.11 dưới đây sẽ cho ta bốn thành phần liên thông mạnh. Thành phần liên thông mạnh 1: 7, 5, 6. Thành phần liên thông mạnh 2: 4, 3, 2. Thành phần liên thông mạnh 3: 11, 10, 9, 8. Thành phần liên thông mạnh 4: 1. Hình 3.11. Đồ thị có hướng G = 2.4.2. Bài toán định chiều đồ thịPTIT Một trong những ứng dụng quan trọng của đồ thị là biểu diễn đồ thị cho các hệ thống giao thông. Đối với hệ thống giao thông người ta quan tâm đến liệu hệ thống có thể định chiều được hay không. Định nghĩa. Phép định chiều đồ thị vô hướng liên thông là phép biến đổi đồ thị vô hướng liên thông thành đồ thị có hướng liên thông mạnh. Đồ thị vô hướng G = có thể dịch chuyển được thành đồ thị có hướng liên thông mạnh bằng cách định chiều mỗi cạnh vô hướng thành một cung có hướng được gọi là đồ thị định chiều được. Ví dụ đồ thị vô hướng trong Hình 3.12 dưới đây được gọi là định chiều được. 61
  62. 2 5 2 5 1 6 1 6 3 4 3 4 Hình 3.12. Phép định chiều đồ thị vô hướng liên thông. Định lý. Đồ thị vô hướng liên thông G = định chiều được khi và chỉ khi tất cả các cạnh e E của G đều không phải là cầu. Bạn đọc tự tìm hiểu cách chứng minh định lý trong các tài liệu [1, 2]. Bài toán. Cho đồ thị vô hướng liên thông G = . Hãy định chiều đồ thị G sao cho ta có thể nhận được đồ thị có hướng với ít nhất thành phần liên thông mạnh. 3.5. Một số điểm cần ghi nhớ Thuật toán duyệt theo chiều sâu bắt đầu tại đỉnh u V. Thuật toán duyệt theo rộng sâu bắt đầu tại đỉnh u V. Duyệt tất cả các đỉnh của đồ thị dựa vào DFS(u), BFS(u). Duyệt tất cả các thành phần liên thông của đồ thị dựa vào DFS(u), BFS(u). Tìm đường đi từ đỉn s đến t trên đồ thị dựa vào DFS(u), BFS(u). Kiểm tra tính liên thông mạnh của đồ thị dựa vào DFS(u), BFS(u). Duyệt các đỉnh trụ của đồ thị DFS(u), BFS(u). Duyệt các cạnh cầu của đồPTIT thị DFS(u), BFS(u). Một số ứng dụng quan trọng khác của DFS và BFS. 62
  63. BÀI TẬP 1. Cho đồ thị vô hướng được biểu diễn dưới 0 1 1 1 1 0 0 0 0 0 0 0 0 dạng ma trận kề như Hình bên phải. Hãy 1 0 1 1 0 0 0 0 0 0 0 0 0 thực hiện: 1 1 0 1 0 0 0 0 0 0 0 0 0 a) Trình bày thuật toán BFS(u)? 1 1 1 0 0 0 0 0 0 0 0 0 0 b) Kiểm nghiệm thuật toán BFS(u) bắt 1 0 0 0 0 1 1 1 1 0 0 0 0 đầu tại đỉnh u=1? Chỉ rõ kết quả trung 0 0 0 0 1 0 1 0 1 0 0 0 0 gian theo mỗi bước thực hiện của thuật 0 0 0 0 1 1 0 1 0 1 0 0 0 toán. 0 0 0 0 1 0 1 0 1 0 0 0 0 c) Kiểm nghiệm thuật toán BFS(u) bắt 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 1 1 đầu tại đỉnh u=7? Chỉ rõ kết quả trung 0 0 0 0 0 0 0 0 0 1 0 1 1 gian theo mỗi bước thực hiện của thuật 0 0 0 0 0 0 0 0 0 1 1 0 1 toán. 0 0 0 0 0 0 0 0 0 1 1 1 0 2. Cho đồ thị vô hướng được biểu diễn dưới 0 1 1 1 1 0 0 0 0 0 0 0 0 dạng ma trận kề như Hình bên phải. Hãy 1 0 1 1 0 0 0 0 0 0 0 0 0 thực hiện: 1 1 0 1 0 0 0 0 0 0 0 0 0 a) Trình bày thuật toán DFS(u)? 1 1 1 0 0 0 0 0 0 0 0 0 0 b) Kiểm nghiệm thuật toán DFS(u) bắt 1 0 0 0 0 1 1 1 1 0 0 0 0 đầu tại đỉnh u=1? Chỉ rõ kết quả trung 0 0 0 0 1 0 1 0 1 0 0 0 0 gian theo mỗi bước thực hiện của thuật 0 0 0 0 1 1 0 1 0 1 0 0 0 toán. 0 0 0 0 1 0 1 0 1 0 0 0 0 c) Kiểm nghiệm thuật toán DFS(u) bắt 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 1 1 đầu tại đỉnh u=7? Chỉ rõ kết quả trung 0 0 0 0 0 0 0 0 0 1 0 1 1 gian theo mỗi bước thực hiện của thuật 0 0 0 0 0 0 0 0 0 1 1 0 1 toán. 0 0 0 0 0 0 0 0 0 1 1 1 0 3. Cho đồ thị vô hướng được biểuPTIT diễn dưới 0 0 1 0 1 0 1 0 0 0 0 0 0 dạng ma trận kề như Hình bên phải. Hãy 0 0 0 1 0 1 0 0 0 0 0 0 0 thực hiện: 1 0 0 0 1 0 1 0 0 0 1 0 0 a) Trình bày thuật toán duyệt các thành 0 1 0 0 0 1 0 1 0 1 0 0 0 phần liên thông của đồ thị? 1 0 1 0 0 0 1 0 1 0 1 0 1 0 1 0 1 0 0 0 1 0 1 0 0 0 b) Kiểm nghiệm thuật toán trên đồ thị 1 0 1 0 1 0 0 0 1 0 0 0 0 đã cho? Chỉ rõ kết quả trung gian theo 0 0 0 1 0 1 0 0 0 1 0 1 0 mỗi bước thực hiện của thuật toán. 0 0 0 0 1 0 1 0 0 0 1 0 1 0 0 0 1 0 1 0 1 0 0 0 1 0 0 0 1 0 1 0 0 0 1 0 0 0 1 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 1 0 0 0 1 0 1 0 0 63
  64. 4. Cho đồ thị vô hướng được biểu diễn dưới 0 1 1 1 1 0 0 0 0 0 0 0 0 dạng ma trận kề như Hình bên phải. Hãy 1 0 1 1 0 0 0 0 0 0 0 0 0 thực hiện: 1 1 0 1 0 0 0 0 0 0 0 0 0 a) Dựa vào thuật toán BFS, xây dựng 1 1 1 0 0 0 0 0 0 0 0 0 0 thuật toán tìm đường đi từ đỉnh s đến 1 0 0 0 0 1 1 1 1 0 0 0 0 đỉnh t trên đồ thị? 0 0 0 0 1 0 1 0 1 0 0 0 0 b) Tìm đường đi từ đỉnh s=1 đến đỉnh t 0 0 0 0 1 1 0 1 0 1 0 0 0 =13 trên đồ thị đã cho? Chỉ rõ kết quả 0 0 0 0 1 0 1 0 1 0 0 0 0 trung gian theo mỗi bước thực hiện của 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 1 1 thuật toán. 0 0 0 0 0 0 0 0 0 1 0 1 1 c) Viết chương trình tìm đường đi từ s 0 0 0 0 0 0 0 0 0 1 1 0 1 đến t dựa vào biểu diễn đồ thị dưới dạng 0 0 0 0 0 0 0 0 0 1 1 1 0 ma trận kề. 5. Cho đồ thị vô hướng được biểu diễn dưới 0 1 1 1 1 0 0 0 0 0 0 0 0 dạng ma trận kề như Hình bên phải. Hãy 1 0 1 1 0 0 0 0 0 0 0 0 0 thực hiện: 1 1 0 1 0 0 0 0 0 0 0 0 0 a) Dựa vào thuật toán DFS, xây dựng 1 1 1 0 0 0 0 0 0 0 0 0 0 thuật toán tìm đường đi từ đỉnh s đến 1 0 0 0 0 1 1 1 1 0 0 0 0 đỉnh t trên đồ thị? 0 0 0 0 1 0 1 0 1 0 0 0 0 b) Tìm đường đi từ đỉnh s=1 đến đỉnh t 0 0 0 0 1 1 0 1 0 1 0 0 0 =13 trên đồ thị đã cho? Chỉ rõ kết quả 0 0 0 0 1 0 1 0 1 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 trung gian theo mỗi bước thực hiện của 0 0 0 0 0 0 1 0 0 0 1 1 1 thuật toán. 0 0 0 0 0 0 0 0 0 1 0 1 1 c) Viết chương trình tìm đường đi từ s 0 0 0 0 0 0 0 0 0 1 1 0 1 đến t dựa vào biểu diễn đồ thị dưới dạng 0 0 0 0 0 0 0 0 0 1 1 1 0 ma trận kề. 6. Cho đồ thị có hướng được biểu diễn dưới 0 0 0 0 0 1 0 0 0 0 0 0 0 dạng ma trận kề như Hình bên phải. Hãy 0 0 1 0 0 0 0 1 0 0 0 0 0 thực hiện: 0 0 0 0 0 0 0 0 1 0 0 0 1 PTIT1 0 0 0 0 1 0 0 0 0 0 0 0 a) Dựa vào thuật toán DFS, xây dựng 0 0 0 0 0 0 1 0 0 0 0 0 0 thuật toán kiểm tra tính liên thông mạnh 0 0 0 0 0 0 0 0 0 1 0 1 0 của đồ thị? 0 0 0 0 0 0 0 0 0 0 1 0 1 b) Kiểm nghiệm thuật toán trên đồ thị 0 0 0 1 0 0 0 0 0 0 0 1 0 đã cho? Chỉ rõ kết quả trung gian theo 0 0 0 0 1 0 1 0 0 0 0 0 0 mỗi bước thực hiện của thuật toán. 0 1 1 0 0 0 0 0 0 0 0 0 0 c) Viết chương trình kiểm tra tính liên 0 1 0 0 0 0 0 1 0 0 0 0 0 thông mạnh của đồ thị dựa vào biểu 0 0 0 1 0 0 0 0 0 1 0 0 0 diễn ma trận kề. 0 0 0 0 0 0 0 0 1 0 1 0 0 64
  65. 0 0 0 0 0 1 0 0 0 0 0 0 0 7. Cho đồ thị có hướng được biểu diễn dưới dạng ma trận kề như Hình bên phải. Hãy 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 thực hiện: a) Dựa vào thuật toán BFS, xây dựng 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 thuật toán kiểm tra tính liên thông mạnh 0 0 0 0 0 0 0 0 0 1 0 1 0 của đồ thị? 0 0 0 0 0 0 0 0 0 0 1 0 1 b) Kiểm nghiệm thuật toán trên đồ thị 0 0 0 1 0 0 0 0 0 0 0 1 0 đã cho? Chỉ rõ kết quả trung gian theo 0 0 0 0 1 0 1 0 0 0 0 0 0 mỗi bước thực hiện của thuật toán. 0 1 1 0 0 0 0 0 0 0 0 0 0 c) Viết chương trình kiểm tra tính liên 0 1 0 0 0 0 0 1 0 0 0 0 0 thông mạnh của đồ thị dựa vào biểu 0 0 0 1 0 0 0 0 0 1 0 0 0 diễn ma trận kề. 0 0 0 0 0 0 0 0 1 0 1 0 0 8. Cho đồ thị vô hướng được biểu diễn dưới 0 1 0 0 0 0 1 1 1 1 0 0 0 dạng ma trận kề như Hình bên phải. Hãy 1 0 1 0 0 0 1 0 1 0 0 0 0 thực hiện: 0 1 0 1 1 1 0 0 0 0 0 0 0 a) Dựa vào thuật toán BFS, xây dựng 0 0 1 0 1 1 0 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 thuật toán duyệt các đỉnh trụ của đồ thị? 0 0 1 1 1 0 0 0 0 0 0 0 0 b) Kiểm nghiệm thuật toán trên đồ thị 1 1 0 0 0 0 0 1 0 0 0 0 0 đã cho? Chỉ rõ kết quả trung gian theo 1 0 0 0 0 0 1 0 1 0 0 0 0 mỗi bước thực hiện của thuật toán. 1 1 0 0 0 0 0 1 0 0 0 0 0 c) Viết chương trình kiểm tra tính liên 1 0 0 0 0 0 0 0 0 0 1 1 1 thông mạnh của đồ thị dựa vào biểu 0 0 0 0 0 0 0 0 0 1 0 1 1 diễn ma trận kề. 0 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 0 1 1 1 0 9. Cho đồ thị vô hướng được biểu diễn dưới 0 1 0 0 0 0 1 1 1 1 0 0 0 dạng ma trận kề như Hình bên phải. Hãy 1 0 1 0 0 0 1 0 1 0 0 0 0 PTIT0 1 0 1 1 1 0 0 0 0 0 0 0 thực hiện: 0 0 1 0 1 1 0 0 0 0 0 0 0 a) Dựa vào thuật toán DFS, xây dựng 0 0 1 1 0 1 0 0 0 0 0 0 0 thuật toán duyệt các đỉnh trụ của đồ thị? 0 0 1 1 1 0 0 0 0 0 0 0 0 b) Kiểm nghiệm thuật toán trên đồ thị 1 1 0 0 0 0 0 1 0 0 0 0 0 đã cho? Chỉ rõ kết quả trung gian theo 1 0 0 0 0 0 1 0 1 0 0 0 0 mỗi bước thực hiện của thuật toán. 1 1 0 0 0 0 0 1 0 0 0 0 0 c) Viết chương trình kiểm tra tính liên 1 0 0 0 0 0 0 0 0 0 1 1 1 thông mạnh của đồ thị dựa vào biểu 0 0 0 0 0 0 0 0 0 1 0 1 1 diễn ma trận kề. 0 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 0 1 1 1 0 65
  66. 10. Cho đồ thị vô hướng được biểu diễn 0 1 0 0 0 0 1 1 1 1 0 0 0 1 0 1 0 0 0 1 0 1 0 0 0 0 dưới dạng ma trận kề như Hình bên phải. 0 1 0 1 1 1 0 0 0 0 0 0 0 Hãy thực hiện: a) Dựa vào thuật toán BFS, xây dựng 0 0 1 0 1 1 0 0 0 0 0 0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 thuật toán duyệt các cạnh cầu của đồ 0 0 1 1 1 0 0 0 0 0 0 0 0 thị? 1 1 0 0 0 0 0 1 0 0 0 0 0 b) Kiểm nghiệm thuật toán trên đồ thị 1 0 0 0 0 0 1 0 1 0 0 0 0 đã cho? Chỉ rõ kết quả trung gian theo 1 1 0 0 0 0 0 1 0 0 0 0 0 mỗi bước thực hiện của thuật toán. 1 0 0 0 0 0 0 0 0 0 1 1 1 c) Viết chương trình kiểm tra tính liên 0 0 0 0 0 0 0 0 0 1 0 1 1 thông mạnh của đồ thị dựa vào biểu 0 0 0 0 0 0 0 0 0 1 1 0 1 diễn ma trận kề. 0 0 0 0 0 0 0 0 0 1 1 1 0 11. Cho đồ thị vô hướng liên thông G = như dưới đây: Ke(1) = { 2, 3, 4}. Ke(5) = {3, 6, 7, 8, 12}. Ke(9) = {10, 11, 13}. Ke(2) = {1, 3, 4, 6}. Ke(6) = {2, 5, 7, 12}. Ke(10) = {9, 11, 12, 13}. Ke(3) = {1, 2, 4, 5}. Ke(7) = {4, 5, 6, 8}. Ke(11) = {9, 10, 13}. Ke(4) = {1, 2, 3, 7}. Ke(8) = {5, 7, 12}. Ke(12) = {5, 6, 8, 10}. Ke(13) = {9, 10, 11}. Hãy thực hiện: a) Tìm BFS(1) =? b) Tìm BFS(5) =? c) Tìm DFS(1) =? d) Tìm DFS(5) =? d) Tìm đường đi từ 1 đến 13 bằng thuật toán BFS? e) Tìm đường đi từ 1 đến 13 bằng thuật toán DFS? 12. Cho đồ thị vô hướng liên thông G = . Ta gọi đỉnh s V là đỉnh “thắt” của cặp đỉnh u, v V nếu mọi đường đi từ u đến v đều phải qua s. Dựa vào thuật toán duyệt theo chiều sâu (hoặc chiều rộng), hãyPTIT thực hiện: a) Xây dựng thuật toán tìm tất cả các đỉnh thắt s V của cặp đỉnh u, v V? b) Tìm tập đỉnh thắt s V của cặp đỉnh u=1, v=12 trên đồ thị đã cho, chỉ rõ kết quả theo mỗi bước thực hiện của thuật toán? c) Tìm tập đỉnh thắt s V của cặp đỉnh u =1, v =13 trên đồ thị được biểu diễn dưới dạng danh sách kề dưới đây, chỉ rõ kết quả theo mỗi bước thực hiện của thuật toán? Ke(1) = { 2, 3, 4}. Ke(5) = {3, 6, 7, 8, 12}. Ke(9) = {10, 11, 13}. Ke(2) = {1, 3, 4, 6}. Ke(6) = {2, 5, 7, 12}. Ke(10) = {9, 11, 12, 13}. Ke(3) = {1, 2, 4, 5}. Ke(7) = {4, 5, 6, 8}. Ke(11) = {9, 10, 13}. Ke(4) = {1, 2, 3, 7}. Ke(8) = {5, 7, 12}. Ke(12) = {5, 6, 8, 10}. Ke(13) = {9, 10, 11}. 66
  67. CHƯƠNG 4. ĐỒ THỊ EULER, ĐỒ THỊ HAMIL TON 4.1. Đồ thị Euler, đồ thị nửa Euler Định nghĩa. Chu trình đơn trong đồ thị G đi qua mỗi cạnh của đồ thị đúng một lần được gọi là chu trình Euler. Đường đi đơn trong G đi qua mỗi cạnh của nó đúng một lần được gọi là đường đi Euler. Đồ thị được gọi là đồ thị Euler nếu nó có chu trình Euler. Đồ thị có đường đi Euler được gọi là nửa Euler. Rõ ràng, mọi đồ thị Euler đều là nửa Euler nhưng điều ngược lại không đúng. Ví dụ 1. Xét các đồ thị G1, G2, G3 trong Hình 4.1. a b a b a b e e d c d c c d e G1 G2 G3 Hình 6.1. Đồ thị vô hướng G1, G2, G3. Đồ thị G1 là đồ thị Euler vì nó có chu trình Euler a, e, c, d, e, b, a. Đồ thị G3 không có chu trình Euler nhưng chứa đường đi Euler a, c, d, e, b, d, a, b vì thế G3 là nửa Euler. G2 không có chu trình Euler cũng như đường đi Euler. Ví dụ 2. Xét các đồ thị có hướng H1, H2, H3 trong Hình 4.2. a b a b a b c c d ePTIT d d c H1 H2 H3 Hình 4.2. Đồ thị có hướng H1, H2, H3. Đồ thị H2 là đồ thị Euler vì nó chứa chu trình Euler a, b, c, d, e, a vì vậy nó là đồ thị Euler. Đồ thị H3 không có chu trình Euler nhưng có đường đi Euler a, b, c, a, d, c nên nó là đồ thị nửa Euler. Đồ thị H1 không chứa chu trình Euler cũng như chu trình Euler. 4.2. Thuật toán tìm chu trình Euler Để tìm một chu trình Euler của đồ thị ta sử dụng kết quả của định lý sau. Định lý 1. Điều kiện cần và đủ để đồ thị G= là Euler. Đồ thị vô hướng liên thông G= là đồ thị Euler khi và chỉ khi mọi đỉnh của G đều có bậc chẵn. Đồ thị có 67
  68. hướng liên thông yếu G= là đồ thị Euler khi và chỉ khi tất cả các đỉnh của nó đều có bán đỉnh bậc ra bằng bán đỉnh bậc vào (điều này làm cho đồ thị là liên thông mạnh). 4.2.1. Chứng minh đồ thị là Euler Đối với đồ thị vô hướng, để chứng minh đồ thi có là Euler hay không ta chỉ cần thực hiện: Kiểm tra đồ thị có liên thông hay không? Điều này dễ dàng thực hiện bằng cách kiểm tra DFS(u) = V hoặc BFS(u) = V thì ta kết luận đồ thị là liên thông (u là đỉnh bất kỳ của đồ thị). Sử dụng tính chất của ma trận kề biểu đồ thị vô hướng để tính toán bậc của các đỉnh. - Vì BFS(1) = { 1, 2, 6, 3, 5, 7, 4, 11, 8, 0 1 0 0 0 1 0 0 0 0 0 0 0 10, 12, 9, 13} = V. Do vậy, G liên 1 0 1 0 1 1 0 0 0 0 0 0 0 0 1 0 1 1 0 0 0 0 0 1 0 0 thông. - Ta lại có : 0 0 1 0 0 0 1 1 0 0 1 0 0 0 1 1 0 0 1 1 0 0 0 0 0 0 deg(1) = deg(13) = 2. 1 1 0 0 1 0 1 0 0 0 0 0 0 deg (2) = deg(3) = 4 0 0 0 1 1 1 0 1 0 0 0 0 0 deg(4) = deg(5) = 4 0 0 0 1 0 0 1 0 1 1 0 0 0 deg(6) = deg(7) = 4 0 0 0 0 0 0 0 1 0 1 0 1 1 deg(8) = deg(9) = 4 0 0 0 0 0 0 0 1 1 0 1 1 0 deg(10) = deg(11) = deg(12) = 4 0 0 1 1 0 0 0 0 0 1 0 1 0 Chú ý: Tổng các phần tử của hàng u (cột 0 0 0 0 0 0 0 0 1 1 1 0 1 u) là bậc của đỉnh u. Ví dụ tổng các 0 0 0 0 0 0 0 0 1 0 0 1 0 phần tử của hàng 1 là 2 nên deg(1) = 2. Đối với đồ thị có hướng, để chứng minh đồ thi có là Euler hay không ta chỉ cần thực hiện: Kiểm tra đồ thị có liênPTIT thông yếu hay không? Điều này dễ dàng thực hiện bằng cách kiểm tra nếu tồn tại đỉnh u V để DFS(u) = V hoặc BFS(u) = V thì ta kết luận đồ thị là liên thông yếu. Sử dụng tính chất của ma trận kề biểu đồ thị có hướng để tính bán đỉnh bậc ra và bán đỉnh bậc vào của các đỉnh. Bán đỉnh bậc ra của đỉnh u là deg+(u) là số các số 1 của hàng u. Bán đỉnh bậc vào của đỉnh u là deg-(u) là số các số 1 của cột u. Ví dụ để chứng minh đồ thị có hướng được biểu diễn dưới dạng ma trận kề như dưới đây ta thực hiện như sau: 68
  69. + Vì BFS(1) = { 1, 2, 3, 5, 4, 11, 6, 7, 10, 12, 0 1 0 0 0 0 0 0 0 0 0 0 0 8, 9, 13} = V. Do vậy, G liên thông yếu. 0 0 1 0 1 0 0 0 0 0 0 0 0 + Ta lại có: 0 0 0 1 0 0 0 0 0 0 1 0 0 deg+(2)= deg-(2)=deg+(3)= eg-(3) =2 0 0 0 0 0 0 1 0 0 0 1 0 0 + - + - 0 0 1 0 0 1 0 0 0 0 0 0 0 deg (4)=deg (4)=deg (5)=deg (5) =2 1 1 0 0 0 0 0 0 0 0 0 0 0 + - + - deg (6)=deg (6)=deg (7)=deg (7) =2 0 0 0 0 1 1 0 0 0 0 0 0 0 + - + - deg (8)=deg (8)=deg (9)=deg (9) =2 0 0 0 1 0 0 1 0 0 0 0 0 0 deg+(10) = deg-(10) = 2 0 0 0 0 0 0 0 1 0 1 0 0 0 deg+(11) = deg-(11) =2 0 0 0 0 0 0 0 1 0 0 0 1 0 + - 0 0 0 0 0 0 0 0 0 1 0 1 0 deg (12) = deg (12) =2 0 0 0 0 0 0 0 0 1 0 0 0 1 + - - + deg (1)=deg (1)=deg (13)=deg (13) =1 0 0 0 0 0 0 0 0 1 0 0 0 0 G liên thông yếu và có bán đỉnh bậc ra bằng bán đỉnh bậc vào nên G là đồ thị Euler. 4.2.2. Biểu diễn thuật toán tìm chu trình Euler Để tìm một chu trình Euler trên đồ thị, ta thực hiện theo thuật toán trong Hình 4.3 dưới đây: Thuật toán Euler-Cycle(u): Bước 1 (Khởi tạo) : stack =  ; //Khởi tạo một stack bắt đầu là  CE =  ; //Khởi tạo mảng CE bắt đầu là  Push (stack, u) ; //Đưa đỉnh u vào ngăn xếp Bước 2 (Lặp ): while (stack  ) do { //Lặp cho đến khi stack rỗng s = get(stack); //Lấy đỉnh ở đầu ngăn xếp if ( Ke(s)  ) then { // Nếu danh sách Ke(s) chưa rỗng t = ; Push(stack, t)’ //Đưa t vào stack; PTITE = E \ (s,t); // Loại bỏ cạnh (s,t); } else { //Trường hợp Ke(s)= s = Pop(stack);// Đưa s ra khỏi ngăn xếp s CE; //Đưa s sang CE } } Bước 3 (Trả lại kết quả) : ; Hình 4.3. Thuật toán tìm chu trình Euler bắt đầu tại đỉnh u 69
  70. 4.2.3. Kiểm nghiệm thuật toán Ví dụ ta cần tìm một chu trình Euler bắt đầu tại đỉnh u=1 trên đồ thị G = được biểu diễn dưới dạng ma trận kề như dưới đây. Khi đó, các bước thực hiện của thuật toán được thực hiện như Bảng 4.1 (chú ý phần chứng minh ta đã thực hiện ở trên). Bước Trạng thái Stack Giá trị CE 1 1  2 1, 2  3 1, 2, 3  4 1, 2, 3, 4  5 1, 2, 3, 4,7  6 1, 2, 3, 4,7,5  7 1, 2, 3, 4,7,5,2  8 1, 2, 3, 4,7,5,2,6  9 1, 2, 3, 4,7,5,2,6,1  10 1, 2, 3, 4,7,5,2,6 1 11 1, 2, 3, 4,7,5,2,6,5 1 12 1, 2, 3, 4,7,5,2,6,5,3 1 13 1, 2, 3, 4,7,5,2,6,5,3,11 1 14 1, 2, 3, 4,7,5,2,6,5,3,11,4 1 15 1, 2, 3, 4,7,5,2,6,5,3,11,4,8 1 16 1, 2, 3, 4,7,5,2,6,5,3,11,4,8,7 1 17 1, 2, 3, 4,7,5,2,6,5,3,11,4,8,7,6 1 18 1, 2, 3, 4,7,5,2,6,5,3,11,4,8,7 1,6 19 1, 2, 3, 4,7,5,2,6,5,3,11,4,8 1,6,7 20 1, 2, 3, 4,7,5,2,6,5,3,11,4,8,9 1,6,7 21 1, 2, 3, 4,7,5,2,6,5,3,11,4,8,9,10 1,6,7 22 1, 2, 3, 4,7,5,2,6,5,3,11,4,8,9,10,8 1,6,7 23 1, 2, 3, 4,7,5,2,6,5,3,11,4,8,9,10 1,6,7,8 24 1, 2, 3, 4,7,5,2,6,5,3,11,4,8,9,10,11 1,6,7,8 25 1, 2, 3, 4,7,5,2,6,5,3,11,4,8,9,10,11,12PTIT 1,6,7,8 26 1, 2, 3, 4,7,5,2,6,5,3,11,4,8,9,10,11,12,9 1,6,7,8 27 1, 2, 3, 4,7,5,2,6,5,3,11,4,8,9,10,11,12,9,13 1,6,7,8 28 1, 2, 3, 4,7,5,2,6,5,3,11,4,8,9,10,11,12,9,13,12 1,6,7,8 29 1, 2, 3, 4,7,5,2,6,5,3,11,4,8,9,10,11,12,9,13,12,10 1,6,7,8 Đưa lần lượt các đỉnh trong Stack sang CE cho đến khi stack= 30 - CE = 1, 6, 7, 8, 10, 12, 13, 9, 12, 11, 10, 9, 8, 4, 11, 3, 5, 6, 2, 5, 7, 4, 3, 2, 1 Lật ngược lại các đỉnh trong CE ta được chu trình Euler 1- 2- 3- 4- 7-5-2-6-5-3-11-4-8-9-10-11-12-9-13-12-10-8-7-6-1 4.2.4. Cài đặt thuật toán 70
  71. Chương trình tìm một chu trình Euler của đồ thị bắt đầu tạo đỉnh u trên đồ thị vô hướng liên thông được cài đặt theo khuôn dạng đồ thị biểu diễn dưới dạng ma trận kề. Các thủ tục chính bao gồm: Thủ tục Init() : đọc dữ liệu theo khuôn dạng biểu diễn ma trận kề. Thủ tục Kiemtra(): Kiểm tra xem G có là Euler hay không. Thủ tục Euler-Cycle (u) : Xây dựng chu trình Euler bắt đầu tại đỉnh u. #include #include #include #include #define MAX 50 #define TRUE 1 #define FALSE 0 int A[MAX][MAX], n, u=1; void Init(void){ int i, j;FILE *fp; fp = fopen("CTEULER.IN", "r"); fscanf(fp,"%d", &n); printf("\n So dinh do thi:%d",n); printf("\n Ma tran ke:"); for(i=1; i 0) return(FALSE); return(TRUE); } void Euler-Cycle( int u){ 71
  72. int v, x, top, dCE; int stack[MAX], CE[MAX]; top=1; stack[top]=u;dCE=0; do { v = stack[top];x=1; while (x n) { dCE++; CE[dCE]=v; top ; } else { top++; stack[top]=x; A[v][x]=0; A[x][v]=0; } } while(top!=0); printf("\n Co chu trinh Euler:"); for(x=dCE; x>0; x ) printf("%3d", CE[x]); } void main(void){ clrscr(); Init(); if(Kiemtra()) Tim(); else printf("\n Khong co chu trinh Euler"); } 4.3. Thuật toán tìm đường đi Euler Một đồ thị không có chu trình Euler nhưng vẫn có thể có đường đi Euler. Để tìm một đường đi Euler trên đồ thị vô hướng ta sử dụng kết quả của định lý 2. Để tìm một đường đi Euler trên đồ thị có hướng ta sử dụng kết quả của định lý 3. Định lý 2. Đồ thị vô hưPTITớng liên thông G = là đồ thị nửa Euler khi và chỉ khi G có 0 hoặc 2 đỉnh bậc lẻ. Trong trường hợp G có hai đỉnh bậc lẻ, đường đi Euler xuất phát tại một đỉnh bậc lẻ và kết thúc tại đỉnh bậc lẻ còn lại. Trong trường hợp G có 0 đỉnh bậc lẻ G chính là đồ thị Euler. Định lý 3. Đồ thị có hướng liên thông yếu G = là đồ thị nửa Euler khi và chỉ khi tồn tại đúng hai đỉnh u, v V sao cho deg+(u) - deg-(u)= deg-(v) - deg+(v)=1, các đỉnh s u, s v còn lại có deg+(s) =deg-(s). Đường đi Euler sẽ xuất phát tại đỉnh u và kết thúc tại đỉnh v. 4.3.1. Chứng minh đồ thị là nửa Euler Để chứng tỏ đồ thị vô hướng G = là nửa Euler ta cần thực hiện: 72
  73. Chứng tỏ đồ thị đã cho liên thông. Điều này dễ ràng thực hiện được bằng cách sử dụng hai thủ tục DFS(u) hoặc BFS(u). Có 0 hoặc hai đỉnh bậc lẻ. Sử dụng tính chất của các phương pháp biểu diễn đồ thị để tìm ra bậc của mỗi đỉnh. Ví dụ. Chứng minh rằng, đồ thị vô hướng liên thông G = được biểu diễn dưới dạng ma trận kề dưới đây là đồ thị nửa Euler. Chứng minh. Theo tính chất của ma trận 0 1 0 0 1 1 0 0 0 0 0 0 0 1 0 1 0 1 1 0 0 0 0 0 0 0 kề, tổng các phần tử hàng u là bậc của đỉnh 0 1 0 1 1 0 0 0 0 0 1 0 0 u. Vì vậy ta có: 0 0 1 0 1 0 1 1 0 1 1 0 0 deg(1) = deg(13) = 3 1 1 1 1 0 1 1 0 0 0 0 0 0 deg (2) = deg(3) = deg(11) = 4 1 1 0 0 1 0 1 0 0 0 0 0 0 0 0 0 1 1 1 0 1 0 0 0 0 0 deg(12) = deg(6) = deg(7) = 4 0 0 0 1 0 0 1 0 1 1 0 0 0 deg(8) = deg(9) = 4 0 0 0 0 0 0 0 1 0 1 0 1 1 deg(5) = deg(4) = deg(10) = 6 0 0 0 1 0 0 0 1 1 0 1 1 1 0 0 1 1 0 0 0 0 0 1 0 1 0 G liên thông và có 2 đỉnh bậc lẻ u=1 và 0 0 0 0 0 0 0 0 1 1 1 0 1 u=13 nên G là nửa Euler. 0 0 0 0 0 0 0 0 1 1 0 1 0 Để chứng tỏ đồ thị có hướng G = là nửa Euler ta cần thực hiện: Chứng tỏ đồ thị đã cho liên thông yếu. Điều này dễ ràng thực hiện được bằng cách sử dụng hai thủ tục DFS(u) hoặc BFS(u). Có hai đỉnh u và v thỏa mãn deg+(u) - deg-(u)= deg-(v) - deg+(v)=1. Các đỉnh s u, s v còn lại có deg+(s) =deg-(s). Ví dụ. Chứng minh rằng, đồ thị có hướng liên thông yếu G = được biểu diễn dưới dạng ma trận kề dưới đây là đồ thị nửa Euler? Chứng minh. Theo tính chất của ma trận kề, 0 1 0 0 1 0 0 0 0 0 0 0 0 deg +(u) là tổng các phần tử hàngPTIT u, deg-(u) là 0 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 tổng các phần tử cột u. Vì vậy ta có: 0 0 0 0 0 0 1 0 0 1 1 0 0 + - + - deg (2) = deg (2) = deg (3) = deg (3) =2 0 0 1 1 0 1 0 0 0 0 0 0 0 deg+(6) = deg-(6) = deg+(7) = deg-(7) =2 1 1 0 0 0 0 0 0 0 0 0 0 0 + - + - 0 0 0 0 1 1 0 0 0 0 0 0 0 deg (8) = deg (8) = deg (9) = deg (9) =2 + - + - 0 0 0 1 0 0 1 0 0 0 0 0 0 deg (11) = deg (11) = deg (12) = deg (12) =2 0 0 0 0 0 0 0 1 0 1 0 0 0 + - + - deg (5) =deg (5)= deg (4) = deg (4) = 0 0 0 0 0 0 0 1 0 0 0 1 1 deg +(10) = deg-(10)=3 0 0 0 0 0 0 0 0 0 1 0 1 0 + - - + 0 0 0 0 0 0 0 0 1 0 0 0 1 deg (1) - deg (1) = deg (13) – deg (13) =1 0 0 0 0 0 0 0 0 1 0 0 0 0 G liên thông yếu và có 2 đỉnh u=1 và u=13 thỏa mãn điều kiện nên G là nửa Euler. 73
  74. 4.3.2. Thuật toán tìm đường đi Euler Thuật toán tìm đường đi Euler và chu trình Euler chỉ khác nhau duy nhất ở một điểm đó là đầu vào của thuật toán. Đối với thuật toán tìm chu trình Euler, đầu vào thuật toán là đỉnh u V bất kỳ. Đối với thuật toán tìm đường đi trình Euler, đầu vào thuật toán là đỉnh u V là đỉnh bậc lẻ đầu tiên trong trường hợp đồ thị vô hướng. Đối với đồ thị có hướng, đỉnh u V là đỉnh có deg+(u)-deg-(u)=1. Thuật toán tìm đường đi Euler trên đồ thị vô hướng hoặc có hướng được mô tả chi tiết trong Hình 4.4. Thuật toán Euler-Path (u ): - u là đỉnh bậc lẻ đầu tiên nếu G là đồ thị vô hướng - u là đỉnh có deg+(u) - deg-(u) =1. Bước 1 (Khởi tạo) : stack =  ; //Khởi tạo một stack bắt đầu là  dCE =  ; //Khởi tạo mảng dCE bắt đầu là  Push (stack, u) ; //Đưa đỉnh u vào ngăn xếp Bước 2 (Lặp ): while (stack  ) do { //Lặp cho đến khi stack rỗng s = get(stack); //Lấy đỉnh ở đầu ngăn xếp if ( Ke(s)  ) then { // Nếu danh sách Ke(s) chưa rỗng t = ; Push(stack, t)’ //Đưa t vào stack; E = E \ (s,t); // Loại bỏ cạnh (s,t); } else { //Trường hợp Ke(s)= s = Pop(stack);// Đưa s ra khỏi ngăn xếp s dCE; //Đưa s sang dCE } } Bước 3 (Trả lại kết quả) : ; Hình 4.4. Thuật toán tìm đường đi Euler trên đồ thị. 4.3.3. Kiểm nghiệm thuật toán Ví dụ ta cần tìm đường đi Euler trên đồ thị có hướng liên thông yếu được biểu diễn dưới dạng ma trận kề trong Mục 4.3.1. Khi đó, đỉnh u có deg+(u)-deg-(u)=1 là đỉnh 1. Kết quả thực hiện của thuật toán Hình 4.4 được thể hiện trong Bảng 4.2 dưới đây. 74