Tối ưu hóa phân bổ và định giá đất đai theo thuật toán di truyền định hướng không gian

59 16 0
Tối ưu hóa phân bổ và định giá đất đai theo thuật toán di truyền định hướng không gian

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

Thông tin tài liệu

i ĐẠI HỌC THÁI NGUYÊN ĐẠI HỌC CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG TRẦN CAO BÁCH TỐI ƯU HÓA PHÂN BỔ VÀ ĐỊNH GIÁ ĐẤT ĐAI THEO THUẬT TOAN DI TRUYỀN ĐỊNH HƯỚNG KHÔNG GIAN LUẬN VĂN THẠC SĨ KHOA HỌC MÁY TÍNH THÁI NGUYÊN 2017 ii LỜI CẢM ƠN Đầu tiên, em xin gửi lời cảm ơn chân thành sâu sắc tới thầy Lê Hoàng Sơn, người đưa em đến với đề tài này, tạo điều kiện thuận lợi cho em nghiên cứu khoa học, nhiệt tình hướng dẫn để em hồn thành luận văn tốt nghiệp cách tốt Em xin gửi lời cảm ơn đến thầy cô khoa sau đại học trường đại học công nghệ thông tin truyền thông thái nguyên dạy bảo, giúp đỡ, tạo điều kiện cho em thời gian em học tập trường Xin gửi lời cảm ơn thầy cô, anh chị bạn Trung tâm Tính tốn Hiệu cao, Trường Đại học Khoa Học Tự Nhiên giúp đỡ em suốt trình học tập nghiên cứu trung tâm Em xin gửi lời cảm ơn chân thành tới công ty Tư vấn Geo Việt cung cấp số liệu, tài liệu liên quan nhiệt tình giúp đỡ em trình thực luận văn Xin gửi lời cảm ơn tới gia đình động viên, tạo điều kiện cho em thực luận văn, lời cảm ơn tới bạn bè giúp đỡ em học tập năm vừa qua Học viên Trần Cao Bách iii LỜI CAM ĐOAN Tơi xin cam đoan luận văn “Tối ưu hóa phân bổ định giá đất đai theo thuật toán di truyền định hướng không gian” thực hướng dẫn khoa học TS Lê Hoàng Sơn, số liệu kết nghiên cứu luận văn hoàn toàn trung thực chưa sử dụng để bảo vệ cơng trình khoa học nào, thơng tin, tài liệu trích dẫn luận văn rõ nguồn gốc phía cuối luận văn Mọi giúp đỡ cho việc hoàn thành luận văn cảm ơn Nếu có phát gian lận chép tài liệu, cơng trình nghiên cứu tác giả khác mà khơng ghi rõ tài liệu tham khảo, tơi hồn toàn chịu trách nhiệm Thái Nguyên, tháng năm 2017 Học viên Trần Cao Bách iv MỤC LỤC LỜI CẢM ƠN i LỜI CAM ĐOAN iii MỤC LỤC iv DANH MỤC CÁC HÌNH VẼ vi DANH MỤC CÁC THUẬT NGỮ vii MỞ ĐẦU .1 CHƯƠNG 1: BÀI TOÁN PHÂN BỔ VÀ ĐỊNH GIÁ ĐẤT 1.1 Tổng quan toán phân bổ định giá đất .3 1.1.1 Một số khái niệm 1.1.2 Các yếu tố ảnh hưởng trực tiếp đến giá đất 1.1.3 Đánh giá yếu tố ảnh hưởng tới giá đất 1.1.4 Nguyên tắc định giá đất 1.2 Các ứng dụng thực tiễn 12 1.3 Các nghiên cứu liên quan .13 1.4 Tổng kết chương 14 CHƯƠNG 2: THIẾT KẾ THUẬT TỐN DI TRUYỀN ĐỊNH HƯỚNG KHƠNG GIAN CHO BÀI TOÁN PHÂN BỔ VÀ ĐỊNH GIÁ ĐẤT 15 2.1 Tổng quan thuật tốn tối ưu tiến hóa 15 2.1.1 Cấu trúc liệu giải thuật tiến hóa 15 2.1.2 Công thức giải thuật di truyền 21 2.1.3 Các thành phần thuật toán .22 2.2 Thiết kế thuật tốn di truyền khơng gian cho toán định giá đất 23 2.2.1 Thuật toán di truyền .23 2.2.2 Thuật toán di truyền đinh hướng không gian (SGA) 26 2.3 Phân loại vùng không gian đồ theo giá đất hỗ trợ phân bổ đất đai thuật toán FCM 28 2.4 Tổng kết chương 30 v CHƯƠNG 3: XÂY DỰNG HỆ THỐNG THƠNG TIN CHO BÀI TỐN PHÂN BỔ VÀ ĐỊNH GIÁ ĐẤT 31 3.1 Giới thiệu hệ thống LIS 31 3.1.1 Khái niệm .31 3.1.2 Kiến thức chung 32 3.1.3 Chức vai trò hệ thông tin đất đai .34 3.2 Thiết kế hệ thống WebGIS cho toán phân bổ định giá đất 39 3.2.1 Giao diện chương trình 44 3.2.2 Thử nghiệm chương trình đánh giá kết 46 3.3 So sánh ưu nhược điểm thuật toán 48 3.4 Tổng kết chương 49 KẾT LUẬN .50 TÀI LIỆU THAM KHẢO 51 vi DANH MỤC CÁC HÌNH VẼ Hình 2.1: Sơ đồ thuật tốn di truyền [2] 20 Hình 2.2: Mô tả trọng số W 26 Hình 2.3: Mơ tả tìm hai hàng xóm 27 Hình 3.1: Các thành phần LIS[2] 35 Hình 3.2: Sơ đồ vận hành LIS 38 Hình 3.3: Các bước trình xây dựng ứng dụng .39 Hình 3.4: Sơ đồ hệ thống ứng dụng 40 Hình 3.5: Biểu đồ ứng dụng 41 Hình 3.6: Một phần sở liệu thuộc tính 42 Hình 3.7: Bảng sở liệu khơng gian số quận 43 Hình 3.8: Một phần sở liệu địa điểm khảo sát giá đất 43 Hình 3.9: Giao diện ban đầu hệ thống 44 Hình 3.10: Giao diện bật tắt Layer 44 Hình 3.11: Giao diện tắt Layers 45 Hình 3.12: Giao diện click vào điểm đồ .45 Hình 3.13: Giao diện nhập giá trị vào lọc 46 Hình 3.14: Bảng fitness .47 Hình 3.15: Giá hai hàng xóm 47 vii DANH MỤC CÁC THUẬT NGỮ Thuật ngữ Viết tắt Giải thích Spatial Genetic Algorithm SGA Thuật tốn di truyền khơng gian Genetic Algorithm GA Thuật tốn di truyền Real Estate BĐS Bất động sản Travelling Salesman Problem TSP Bài toán người du lịch Ant Colony Optimization ACO Thuật toán tối ưu đàn kiến Particle Swarm Optimization PSO Thuật toán tối ưu bầy đàn Chromosome NST Nhiễm sắc thể Geographic Information Systems GIS Hệ thống thông tin địa lý Land Information System LIS Hệ thống thông tin đất đai Data Exchange Center DEC Trung tâm trao đổi liệu Database Managemant System DBMS Hệ quản trị sở liệu Analyst Nhà phân tích System Administrator Người quản trị hệ thống Database Administrator Người quản trị sở liệu Data Server Máy chủ liệu Mean Square Error MSE Bình phương sai số trung bình Metadata Siêu liệu Web Server Máy chủ web HyperText Transfer Protocol HTTP Giao thức truyền tải siêu văn MỞ ĐẦU Bài toán phân bổ định giá đất đai vấn đề quan tâm nhiều ý nghĩa ứng dụng thực tiễn Bài tốn hiểu tái cấu lại đồ chuyên đề đất đai vùng địa lý cho trước theo số tiêu chí kinh tế - xã hội định, ví dụ: giá đất, tác động mơi trường Đây toán tối ưu đa mục tiêu với liệu gồm phần không gian liệu chun đề Chính việc tìm phương án tối ưu thể kết lên đồ chuyên đề thành cụm hay vùng khơng gian có mức độ khác theo chuyên đề có ý nghĩa việc hoạch định phát triển sách xã hội Tầm quan trọng ứng dụng tốn xem tài liệu [9] Trong thuật toán tối ưu đa mục tiêu, thuật toán di truyền [10] cơng cụ tốt cho phép mơ hình hóa tìm kiếm nghiệm không gian nghiệm vô lớn Thuật toán di truyền nhánh thuật toán tối ưu tiến hóa sử dụng chiến lược tìm nghiệm theo hành vi sinh học bầy đàn, quần thể, v.v Thuật toán ứng dụng rộng rãi nhiều toán khai phá liệu, nhận dạng mẫu công cụ tốt cho việc giải toán tối ưu với ràng buộc lớn Trong [9] tác giả đề xuất thuật tốn di truyền định hướng khơng gian cách tích hợp thơng tin khơng gian đối tượng đồ vào q trình tìm kiếm nghiệm Chính vậy, nghiệm tìm có độ xác chất lượng cao thuật tốn di truyển gốc thơng tin sử dụng bao gồm thông tin đồ chuyên đề Do vậy, đối tượng nghiên cứu chúng tơi luận văn Mục tiêu luận văn tìm hiểu thuật tốn di truyền định hướng không gian xây dựng ứng dụng phân bổ định giá đất đai hệ thống WebGIS Kết thu luận văn tài liệu chi tiết tiếp cận tối ưu tiến hóa, cụ thể thuật tốn di truyền định hướng khơng gian cho tốn phân bổ định giá đất đai, phần mềm ứng dụng thuật toán cho liệu thực tế đất đai Nội dung báo cáo gồm lời nói đầu chương chính: Chương 1: Bài tốn phân bổ định giá đất Chương đưa kiến thức tổng qua toán giá đất, ứng dụng thực tiễn, nghiên cứu liên quan mô tả liệu giá đất Hà Nội Chương 2: Thiết kế thuật tốn di truyền định hướng khơng gian cho toán phân bổ định giá đất Nội dung chương bao gồm nội dung sau: Tìm hiểu chung thuật tốn tối ưu tiến hóa, tìm hiểu cụ thể thuật toán di truyền thiết kế thuật tốn di truyền khơng gian cho tốn phân bổ định giá đất Chương 3: Xây dựng hệ thống thơng tin đất đai cho tốn phân bổ định giá đất Chương giới thiệu công nghệ WebGIS, đưa chức vai trò hệ thông tin đất đai, thiết kế hệ thống WebGIS cho toán định giá đất, xây dựng sở liệu, hiển thị giao diện chương trình kết đánh giá CHƯƠNG 1: BÀI TOÁN PHÂN BỔ VÀ ĐỊNH GIÁ ĐẤT 1.1 Tổng quan toán phân bổ định giá đất 1.1.1 Một số khái niệm Đất đai loại tài nguyên tự nhiên, loại tài nguyên sản xuất, bao gồm lục địa mặt nước bề mặt trái đất Theo quan điểm kinh tế học đất đai khơng bao gồm mặt đất cịn bao gồm tài nguyên lòng đất tất thứ sinh sơi mặt đất lịng đất không lao động người làm ra, tức bao gồm nước mặt đất nước ngầm, thổ nhưỡng, thực vật động vật Với nghĩa chung nhất, lớp bề mặt trái đất, bao gồm đồng ruộng, đồng cỏ, bãi chăn thả, rừng, bãi hoang, mặt nước, đầm lầy bãi đá v.v Với nghĩa hẹp đất đai biểu khối lượng tính chất quyền lợi tài sản mà người chiếm đất Nó bao gồm lợi ích đất mặt pháp lý quyền theo tập quán không thành văn [2] Đất đai tài sản đồng thời đất đai coi tài sản đặc biệt thân khơng lao động làm ra, mà lao động tác động vào đất đai để biến từ trạng thái hoang hố trở thành sử dụng vào đa mục đích giá đất giá trị quyền sử dụng đất tính đơn vị diện tích đất Định giá đất loại hoạt động chuyên mơn vừa mang tính kinh tế-kỹ thuật, tính pháp lý, vừa mang tính xã hội, tính nghệ thuật Hoạt động định giá hình thành, tồn phát triển gắn liền với hình thành, tồn phát triển thị trường Định giá tài sản hình thức phân tích kinh tế ứng dụng Nhiều khái niệm sử dụng lý thuyết định giá khái niệm kinh tế có số khác biệt Ngày nay, theo đà phát triển kinh tế thị trường, lý thuyết định giá ngày hoàn thiện trở thành ngành chuyên sâu kinh tế [2] 38 Thông tin kinh tế đất: phân hạng đất, giá trị đất cơng trình đất Các liệu quy hoạch đất, kế hoạch sử dụng đất Hình 3.2: Sơ đồ vận hành LIS c Chức năng, vai trò LIS Là phận GIS, LIS có chức hệ thống địa lí: nhập liệu, quản trị liệu,phân tích, mơ hình hóa liệu, hiển thị liệu Các chức đảm bảo khả thực vai trò LIS: - Tạo sở liệu địa lí đầy đủ thống (cho vùng hay lãnh thổ) để thể thơng tin có liên quan đên không gian - Là công cụ trực tiếp phục vụ cho việc hoạch định sách đất đai, cho quy định quy hoạch, kế hoạch sử dụng hợp lí hiệu đất đai cho mục tiêu quy hoạch phát triển kinh tế xã hội - Là cơng cụ để quản lí việc sử dụng đất đai đến đất - Là công cụ để quản lí thống liệu hồ sơ địa - Là cơng cụ hiệu cho việc cung cấp thông tin đất đai cho thị trường quyền sử dụng đất, thị trường bất động sản - Đáp ứng nhu cầu thông tin đất đai người nhu cầu chung phát triển xã hội nâng cao dân trí 39 3.2 Thiết kế hệ thống WebGIS cho toán phân bổ định giá đất a Các bước trình xây dựng ứng dụng: Hình 3.3: Các bước trình xây dựng ứng dụng - Bước 1: chạy thuật toán di truyền định hước không gian Bước cho phép đọc liệu đầu vào thuật toán - Bước 2: xây dựng phương án gồm giá đất ước tính Thực chuẩn hóa liệu dựa vào liệu để xây dựng hàm số phù hợp với tốn Qua q trình thuật tốn, xây dựng tham số phù hợp với fitness nhỏ - Bước 3: Sử dụng FCM để phân chia giá đất thành vùng từ thấp đến cao Thực phân chia để có vùng giá đất khác 40 - Bước 4: Ánh xạ hiển thị kết nhóm hệ thống thơng tin đất đai LIS Hiện thị kết lên đồ gồm điểm vùng giá đất b Thiết kế hệ thống Hệ thống cho phép người dùng ứng dụng xem chi tiết kết định giá đất Hình 3.4: Sơ đồ hệ thống ứng dụng c Chức ứng dụng Có chức sau: - Hiển thị đồ (có full screen, zoom in, zoom out) - Chú thích - Bộ lọc theo cụm - Bảng so sánh giá dự báo, MSE, giá chuẩn hóa GA với SGA - Tác nhân: người dùng - Đầu ra: thông tin phân bổ định giá đất - Mô tả chi tiết: + Người dùng click hay nhập giá vào điểm đồ hay lọc + Hệ thống xử lý nghiệp vụ hiển thị thông tin lên trình duyệt 41 - Biểu đồ trình tự: Hình 3.5: Biểu đồ ứng dụng d Xây dựng sở liệu - Dữ liệu thuộc tính Dữ liệu thu thập số quận Hà Nội (dữ liệu cung cấp Bộ Xây Dựng thông qua công ty Tư vấn Geo Việt) gồm nhiều trường thuộc tính liên quan trực tiếp đến giá đất Như mô tả chương có nhiều yếu tố ảnh hưởng trực tiếp tới giá đất, phần đặc tả liệu đề cập đến thuộc tính mà ta sử dụng để định giá đất ứng dụng mổ tả chương Từ liệu thu thập ta chia làm phần, phần dùng để training phần dùng để testing Cụ thể tập training có 627 điểm, testing có 269 điểm Trước tiên, từ liệu thu thập chuẩn hóa đoạn [0,1] cho hợp lý việc xử thuật tốn, sau ta xây dựng sở liệu thuộc tính tốn 42 Hình 3.6: Một phần sở liệu thuộc tính - Dữ liệu khơng gian Từ liệu shapefile thu thập được, ta xây dựng sở liệu khơng gian cho tốn.Dữ liệu khơng gian gồm có đồ lớp liệu điểm Bản đồ nền: Là đồ số quận huyện Hà Nội - Lớp liệu điểm: Mơ tả vị trí định giá đất, điểm gắn mã số 43 Sau bảng sở liệu khơng gian: Hình 3.7: Bảng sở liệu khơng gian số quận Hình 3.8: Một phần sở liệu địa điểm khảo sát giá đất 44 3.2.1 Giao diện chương trình a Giao diện ban đầu hệ thống Hình 3.9: Giao diện ban đầu hệ thống b Giao diện tắt lớp quận Hình 3.10: Giao diện bật tắt Layer 45 Mô tả: CheckBox dùng để bật tắt layer với lớp quận c Giao diện tắt Layer với liệu lớp Hình 3.11: Giao diện tắt Layers d Giao diện click vào điểm đồ Hình 3.12: Giao diện click vào điểm đồ 46 Mô tả: - ID: mã số điểm liệu đồ - Khi click vào điểm đồ hiển thị đầy đủ giá dự đoán GA giá dự đoán SGA e Giao diện nhập số liệu vào lọc Hình 3.13: Giao diện nhập giá trị vào lọc 3.2.2 Thử nghiệm chương trình đánh giá kết - Cơng thức tính bình phương sai số trung bình: MSE = với + Di giá trị thực + Fi giá trị tính theo hà mục tiêu 47 Hình 3.14: Bảng fitness Hình 3.15: Giá hai hàng xóm 48 Kết : Tập testing: GA SGA MSE 0.010 0.006 Thời gian chạy (s) 0.508 169.754 GA SGA MSE 0.009 0.002 Thời gian chạy (s) 0.943 1097.057 Tập training: 3.3 So sánh ưu nhược điểm thuật toán Từ bảng ta thấy tốc độ xử lý liệu thuật tốn SGA ln chậm thuật tốn GA thuật tốn SGA dựa tinh toán di truyền GA giá dự đốn GA SGA tiếp tục tìm hàng xóm gần tính tốn fitness dựa vào hàng xóm đưa giá dự đốn Nhưng độ xác thuật tốn SGA xác thuật tốn SGA có sai số bình phương MSE nhỏ GA Sau sử dụng thuật toán truyền định hướng khơng gian cho tốn phân bổ định giá đất ta nhận thấy + Thời gian để tính toán đưa giá dự báo nhiều thời gian + Thuật toán đưa giá dự đốn tương đối (là cơng cụ hỗ trợ) khơng thể đưa giá xác số liệu hạn chế để định giá 49 khơng phụ thuộc vào yếu tố cố định mà cịn nhiều yếu tố khơng cố định thay đổi hàng năm ví dụ diện tích mặt tiếp xúc, diện tích v.v 3.4 Tổng kết chương Chương mô tả tổng thể hệ thống WebGIS hệ thống xử lí ứng dụng định giá đất dựa hệ thống thơng tin địa lí Một số hình ảnh minh họa rõ ràng hoạt động chương trình đưa kết chạy thử nghiệm 50 KẾT LUẬN Với nhiều ý nghĩa thực tế, vấn đề định giá đất quan tâm nhiều nhà kinh tế học, vấn đề định giá đất phức tạp, phụ thuộc vào nhiều yếu tố: kinh tế, xã hội, vị trí, điều kiện phát triển,v.v Mặt khác yếu tố biến động theo thời gian nên cần phân tích kỹ lưỡng định giá Xuất phát từ thực tiễn đó, với quan tâm, hứng thú tới vấn đề giá đất, giá đất, luận văn nghiên cứu trình bày số vấn đề sau: - Tìm hiểu tổng quan giá đất yếu tố ảnh hưởng trực tiếp tới giá đất phương pháp định giá đất Từ đặt toán phân bổ định giá đất với phương pháp định giá đất hợp lý với biến động yếu tố ảnh hưởng trực tiếp tới giá đất - Tìm hiểu thuật tốn di truyền định hướng khơng gian áp dụng thuật tốn vào toán phân bổ định giá đất Thiết kế thuật toán cách hợp lý đề phù hợp với toán - Xây dựng hệ thống định giá đất WebGIS người dùng dễ sử dụng trực quan Dựa vào kết đạt được, tương lai, đề tài phát triển theo hướng sau: - Tiếp tục cải tiến để phù hợp với biến động yếu tố ảnh hưởng đến giá đất, nâng cao độ xác cho tốn - Mở rộng vùng định giá đất luận văn nghiên cứu số quận Hà Nội - Tích thêm thuộc tính khơng gian để tính tốn cho giá trị MSE nhỏ 51 TÀI LIỆU THAM KHẢO I Tài liệu tiếng Việt [1] Lê Hoàng Sơn, Nguyễn Thọ Thơng (2015) Lập trình ứng dụng webgis Hà Nội: NXB giáo dục Việt Nam [2] Hồ Thị Lam Trà (ngày 15 tháng 12 năm 2005) Giáo trình định giá đất [3] Trần Kim Hương, Nguyễn Thị Ngọc Chi Giải thuật di truyền (GAs) ứng dụng [4] Nguyễn Thị Thúy Hoài Giải thuật di truyền ứng dụng.Tuyển tập báo cáo “ hội nghị sinh viên nghiên cứu khoa học” lần thứ Đại học Đà Nẵng - 2008 II Tài liệu tiếng Anh [5] Bezdek, J.C., R Ehrlich, et al (1984), "FCM: the fuzzy c-means clustering algorithm", Computers and Geosciences, 10, pp.191-203 [6] Dawkins, R (2006) The Selfish Gene: 30th Anniversary Edition Oxford University Press: 21(4),pp.211-217 [7] Bagley, J D (1967) The Behavior of Adaptive Systems Which Employ Genetic and Correlative Algorithms [8].Charbonneau, P (1988) An introduction to genetic algorithms for numerical optimization [9].Demetriou, D (2013) The development of an integrated planning and decision support system (IPDSS) for land consolidation Springer Science & Business Media [10] Goldberg, D E (1988) Genetic algorithms and machine learning Machine learning 3(2), 95-99 [11] Hartigan J A and Wong.M A Algorithm AS 136: A K-Means Clustering Algorithm Vol 28, No (1979), pp 100-108, ournal of the Royal Statistical Society 52 [12] Jaco F Schutte EGM 6365 The Particle Swarm Optimization Algorithm - Structural Optimization Fall 2005 [13] Bäck, T., Hoffmeister, F and Schwefel (H P (1991).A survey of evolution strategies) Proceedings of the 4th International Conference on Genetic Algorithms Morgan Kaufmann, San Mateo, CA [14] Booker, L ((1987) Improving search in genetic algorithms) In Genetic Algorithms and Simulating Annealing ed L Davis, pp 61–73: Morgan Kaufmann, San Mateo, CA [15] Davis, LD (1991) Handbook of Genetic Algorithms Van Nostrand Reinhold, New York [16] Eshelman, L (1991) The CHC adaptive search algorithm In Foundations of Genetic Algorithms ed G Rawlins, pp 256–283: Morgan Kaufmann, San Mateo, CA [17] Freisleben, P Mers (1996) New genetic local search operator for the travelling University of Siegen [18] Goldberg, D (1987) Simple genetic algorithms and the minimal, deceptive problem In Genetic Algorithms and Simulated Annealing ed L Davis Pitman, London [19]S.V Sivanandam, S.N Deepa (2008) Introduction to Genetic Algorithms, Spinger [20] Wright, A (1991) Genetic algorithms for real parameter optimization In Foundations of Genetic Algorithms ed G Rawlins Morgan Kaufmann, San Mateo, CA ... đến giá 2.3 Phân loại vùng không gian đồ theo giá đất hỗ trợ phân bổ đất đai thuật toán FCM Dựa vào kết SGA ta có liệu khơng gian điểm giá đất từ ta xây dựng toán phân bổ đất đai sử dụng thuật toán. .. chạy theo số hệ, tiếp đến đánh giá cho chạy theo kết ngược lại.” 2.2 Thiết kế thuật toán di truyền khơng gian cho tốn phân bổ định giá đất 2.2.1 Thuật toán di truyền Bài toán định giá đất dựa vào... chung thuật tốn tối ưu tiến hóa, tìm hiểu cụ thể thuật toán di truyền thiết kế thuật tốn di truyền khơng gian cho tốn phân bổ định giá đất Chương 3: Xây dựng hệ thống thơng tin đất đai cho tốn phân

Ngày đăng: 24/02/2021, 19:31

Từ khóa liên quan

Tài liệu cùng người dùng

  • Đang cập nhật ...

Tài liệu liên quan