Chu de 16 2015

6 điểm ĐƯỜNG THẲNG CHỨA 3 ĐIỂM Cho một dãy gồm N điểm khác nhau N 0 thì ghi ra SL bộ số để trong cặp ngoặc vuông, mỗi bộ số đó là chỉ số của các điểm nằm trên đường thẳng tương ứng.. 7 [r]

LUYỆN THI HỌC SINH GIỎI MÔN TIN HỌC Thời gian làm bài: 150 phút CHỦ ĐỀ 16 BÀI TOÁN HÌNH HỌC Tổng quan TT Tên bài Đường thẳng chứa điểm Điểm đa giác Diện tích phủ sóng Tệp chương trình Tệp liệu Tệp kết Thời Gian DT3D.PAS DT3D.INP DT3D.OUT 1s/test DTDG.PAS PSONG.PAS DTDG.INP PSONG.INP DTDG.OUT PSONG.OUT 1s/test 1s/test Bài (6 điểm) ĐƯỜNG THẲNG CHỨA ĐIỂM Cho dãy gồm N điểm khác (N <= 200) Yêu cầu: Tìm xem có bao nhiêu đường thẳng khác nhau, đường thẳng đó chứa ít điểm N điểm nói trên Dữ liệu: Đọc vào từ tệp DT3D.INP, gồm nhiều phương án, phương án trên dòng, trên dòng đó ghi N cặp số X, Y là tọa độ N điểm Kết quả: Ghi màn hình và ghi tệp DT3D.OUT, phương án trên dòng, trên dòng đó ghi số lượng (SL) các đường thẳng tìm (nếu không tìm đường thẳng nào thỏa mãn thì đặt SL = 0), tiếp đó, SL > thì ghi SL số để cặp ngoặc vuông, số đó là số các điểm nằm trên đường thẳng tương ứng Ví dụ: DT3D.INP DT3D.OUT 000102111222 [1 3] [1 6] [3 6] 00011112202526 [5 7] Bài (7 điểm) ĐIỂM TRONG ĐA GIÁC Cho đa giác không tự cắt có N đỉnh và điểm P =(px, py) Yêu cầu: Txác định vị trí điểm P đa giác đã cho Dữ liệu: Đọc từ tệp DTDG.INP, gồm nhiều phương án, phương án trên dòng, trên dòng đó ghi cặp số px, py, tiếp đó ghi N cặp số, cặp là tọa độ đỉnh đa giác Kết quả: Ghi màn hình và ghi tệp DTDG.OUT, phương án trên dòng, trên dòng đó ghi vị trí VTP điểm P đa giác với quy định sau: VTP = 1: P nằm đa giác, VTP = -1: P nằm ngoài đa giác, VTP = 0: P là đỉnh đa giác P nằm trên cạnh đa giác Ví dụ: DTDG.INP DTDG.OUT -2 -1 -1 -1 1 0 -1 -1 0 -1 -1 -1 1 0 -1 Bài (7 điểm) DIỆN TÍCH PHỦ SÓNG Một công ty kinh doanh dịch vụ điện thoại di động, lắp đặt N ăng ten thu phát sóng trên mặt đất Ăng ten thứ i lắt đặt vị trí (X[i], Y[i]), có công suất là R[i] và có diện tích phủ sóng là hình vuông với tâm là (X[i], Y[i]) và diện tích là 4R[i]2 Yêu cầu: Tính tổng diện tích phủ sóng N ăng ten đó Dữ liệu: Đọc từ tệp văn PSONG.INP, gồm nhiều phương án, phương án ghi trên dòng, trên dòng đó ghi N ba số X[i], Y[i], R[i], tương ứng với vị trí và công suất các ăng ten Các số là số thực Kết quả: Ghi màn hình và ghi tệp văn PSONG.OUT, phương án trên dòng, trên dòng đó ghi tổng diện tích phủ sóng tìm được, với chữ số thập phân Ví dụ: PSONG.INP 3.0 3.0 3.0 1.5 1.5 1.0 3.0 2.0 4.0 1.5 2.5 3.0 PSONG.OUT 36.00 67.00 TÓM TẮT LÝ THUYẾT Điểm, đường thẳng, đoạn thẳng, đa giác Điểm là đối tượng hình học sở, biểu diễn cặp số (x, y) Trong Pascal, có thể khai báo điểm là cặp số thông thường, có thể dùng ghi để khai báo các điểm: Type Point = record x,y: integer End; var P, Q, M: Point; Qua hai điểm khác luôn có và đường thẳng Phương trình đường thẳng qua hai điểm P(px,py) và Q(qx,qy) là: dy*(x-px) = dx*(y-py) đó dx = qx – px, dy = qy – py (Tỷ số dy/dx chính là tang góc tạo đường thẳng PQ và trục hoành) Trong Pascal, có thể khai báo đường thẳng cặp điểm thông thường, có thể dùng ghi để khai báo các đường thẳng: Type Line = record P1, P2: Point End; var L1, L2: Line; Đoạn thẳng biểu diễn cặp điểm (P,Q), là phần đường thẳng qua hai điểm đó, và gồm các điểm M nằm hai điểm đó Điều kiện cần và đủ để điểm M là điểm thuộc đoạn PQ là tồn số thực α với <= α <=1 Mx = (1- α) * Px + α * Qx và My = (1- α) * Py + α * Qy Trong Pascal, không có cách khai báo đoạn thẳng nào khác với cách khai báo đường thẳng nói trên Đa giác biểu diễn danh sách N điểm, đó hai điểm nối với đoạn thẳng (điểm cuối có điểm là điểm đầu) Trong Pascal, có thể khai báo đa giác theo nhiều cách khác nhau: a) Sử dụng hai mảng X, Y riêng biệt: var X, Y: array[1 NMax] of integer; {các mảng chứa NMax tọa độ} Cách này thường dùng cần xử lý các tọa độ riêng biệt b) Sử dụng mảng hai chiều: var P: array[1 Nmax, 2] of integer; {mảng chứa NMax điểm} Theo cách này điểm thứ I có tọa độ là Px = P[i,1] và Py = P[I,2] c) Sử dụng mảng các ghi: Type Diem = record x,y: integer End; var P: array[1 Nmax] of Diem; {mảng chứa NMax điểm} Theo cách này điểm thứ i có tọa độ là Px = P[i].x và Py = P[i].y Giao các đoạn thẳng Để xác định xem hai đoạn thẳng có giao (có điểm chung) hay không, người ta thường sử dụng hàm CCW sau đây để xác định hướng từ điểm (x1,y1) đến điểm (x2,y2) đến điểm (x3,y3): function CCW(x1,y1,x2,y2,x3,y3: integer): integer; Var dx1, dy1, dx2, dy2: integer; BEGIN dx1 := x2 – x1; dy1 := y2 – y1; dx2 := x3 – x1; dy1 := y3 – y1; if dx1*dy2 > dy1*dx2 then ccw := 1; {ngược chiều kim đồng hồ} if dx1*dy2 < dy1*dx2 then ccw := -1; {thuận chiều kim đồng hồ} if dx1*dy2 = dy1*dx2 then ccw := 0; {ba điểm thẳng hàng} END; Từ đó, giao hai đoạn thẳng (x1,y1)-(x2,y2) và (x3,y3)-(x4,y4) xác định hàm GiaoNhau sau: function GiaoNhau(x1,y1,x2,y2,x3,y3,x4,y4: integer): boolean; var c1, c2, c3, c4: integer; BEGIN c1 := CCW( x1,y1,x2,y2,x3,y3); {1 – - 3} c2 := CCW( x1,y1,x2,y2,x4,y4); {1 – - 4} c3 := CCW( x3,y3,x4,y4,x1,y1); {3 – - 1} c4 := CCW( x3,y3,x4,y4,x2,y2); {3 – - 2} GiaoNhau := (c1*c2 <= 0) and (c3*c4 <=0); END; Điểm nằm đa giác Để xác định vị trí điểm P đa giác GD ta có thể làm sau: 1) Nếu P là đỉnh P nằm trên cạnh đa giác thì đặt vtP = 2) Trái lại, tìm điểm M nằm ngoài đa giác cho đường thẳng MP không chứa đỉnh đa giác Tìm số điểm sc mà đoạn thẳng MP cắt các cạnh đa giác (sử dụng hàm CCW nói trên) Nếu sc lẻ vtP = (P nằm đa giác), trái lại đặt vtP = -1 (P nằm ngoài đa giác) GIỢI Ý CÁCH GIẢI Bài - Thủ tục TimSL để tìm số đường thẳng: Đặt SL = 0, n1 = N – 1, n2 = N - Khởi tạo mảng CD (0: chưa dùng, 1: đã dùng) Duyệt với i từ đến n2 {Điem thu 1} Đặt i1 = i + Duyệt với j = i1 đến n1 Nếu CD[j] = thì {Diem thu 2} Đặt j1 = j + Duyệt với k từ j1 đến n Nếu TH(i,j,k) thì {Diem thu 3} Tăng sl Đặt DT[sl,1] = i, DT[sl,2] = j Đặt CD[j] = 1, CD[k] = {Dánh dấu đã dùng} - Hàm TH(ii ,jj, kk: byte): boolean xác định điểm ii, jj, kk thẳng hàng Đặt dx = D[jj,1] - D[ii,1], dy = D[jj,2] - D[ii,2] if dx*(D[kk,2] - D[ii,2]) = dy*(D[kk,1] - D[ii,1]) thì TH = true trái lại TH = false Bài - Thủ tục TimVTP để tìm vị trí điểm P đa giác Đặt X[N+1] = X[1], Y[N+1] = Y[1], vtp = Nếu TrungDinh ThuocCanh) thì Thóat Trái lại TimM Đặt sdc = Duyệt với i từ đến N Nếu Cat(i) thì tăng sdc Nếu sdc mod = thì vtp = trái lại vtp = -1 - Hàm TrungDinh: xác định P trùng với đỉnh da giác Duyệt với i từ đến N Nếu (px = X[i]) và (py = Y[i]) thì Đặt TrungDinh = true Thoát TrungDinh = false - Hàm ThuocCanh: xác định P thuộc cạnh đa giác Duyệt với i từ đến N Nếu (px-X[i])*(Y[i+1] - Y[i]) = (py-Y[i])*(X[i+1] - X[i]) thì Nếu ((px - X[i])*(px-X[i+1])<=0) và ((py - Y[i])*(py-Y[i+1])<=0) thì Đặt ThuocCanh = true Thoát Đặt ThuocCanh = false - Thủ tục TimM: tìm điểm M nằm ngoài đa giác Đặt mx = -maxint, my = -maxint Duyệt với i từ đến N Nếu mx < X[i] thì mx = X[i] Nếu my < Y[i] thì my = Y[i] Tăng mx {Tránh trường hợp M trùng với đỉnh đa giác} Lặp lại cc = false Duyệt với i từ đến N Nếu (my-py)*(X[i] - px) = (mx-px)*(Y[i] - py) thì Tăng mx Đặt cc = true Thoát khỏi vòng lặp theo i; Cho đến cc = false - Hàm Cat(ii): xác định đoạn MP cắt cạnh ii đa giác Đặt x1 = X[ii], y1 = Y[ii], x2 = X[ii+1], y2 = Y[ii+1] Đặt cc1 = ccw(x1,y1,x2,y2,px,py) {1 P} Đặt cc2 = ccw(x1,y1,x2,y2,mx,my) {1 M} Đặt cc3 = ccw(px,py,mx,my,x1,y1) {P M 1} Đặt cc4 = ccw(px,py,mx,my,x2,y2) {P M 2} Trả Cat = (cc1*cc2 < 0) và (cc3*cc4 < 0) - Hàm ccw(x1,y1,x2,y2,x3,y3): xác định hướng từ 1-2-3 Đặt dx1 = x2 - x1, dx2 = x3 - x2, dy1 = y2 - y1, dy2 = y3 - y2 Nếu dx1*dy2 > dy1*dx2 thì ccw = {ngược chiều kim đồng hồ} Nếu dx1*dy2 < dy1*dx2 thì ccw := -1 {thuận chiều kim đồng hồ} Nếu dx1*dy2 = dy1*dx2 thì ccw:=0; {thẳng hàng} Bài - Thủ tục TinhDT: tính tổng diện tích phủ sóng SXTang(X, XX) SXTang(Y, YY) Đăt dt = 0, n21 = 2*N-1 Duyệt với i từ đến n21 Duyệt với j từ đến n21 Duyệt với k từ đến N Đặt x1 = XX[i], y1 = YY[j], x2 = XX[i+1], y2 = YY[j+1]; Đặt x3 = X[k]-r[k], y3 = Y[k]-r[k], x4 = X[k]+r[k], y4 = Y[k]+r[k] Nếu NT(x3, y3, x4, y4, x1, y1, x2, y2) thì dt = dt + abs((x1-x2)*(y1-y2)) thoát khỏi vòng lặp theo k - Thủ tục SXTang(mm, xt): tạo và xếp tăng dần mảng xt từ màng mm Duyệt với i từ đến N xt[2*i-1] = mm[i] - R[i], xt[2*i] = mm[i] + R[i]; n2 = 2*N, n21 = n2 -1 Duyệt với i từ đến n21 i1 = i+1; Duyệt với j từ i1 đến n2 Nếu (xt[i]>xt[j]) thì DoiCho(xt,i,j) - Thủ tục DoiCho(mm, ii,jj): đổi chỗ các phần tử ii, jj mảng mm t = mm[ii], mm[ii] = mm[jj], mm[jj] = t - Hàm NT(x1,y1,x2,y2,x3,y3,x4,y4): Xác định hcn 3-4 nằm hcn 1-2 Nếu (x1<=x3) và (x4<=x2) và (y1<=y3) và (y4<=y2) thì NT = true trái lại NT = false DỮ LIỆU ĐỂ KIỂM TRA BÀI Tệp DT3D.INP 0 0 1 0 1 0 1 0 0 1 0 1 0 1 1 1 1 1 1 2 2 2 2 2 2 0 0 2 2 2 2 5 5 5 3 3 3 Tệp DT3D.OUT [1 3] [1 6] [3 6] 3 3 4 4 4 4 4 5 4 6 4 7 7 5 7 9 1 3 [5 [1 [1 [1 [2 [1 [1 3 3 7] 8] 8] [2 9] 3] [1 10] [1 10] [1 10] [3 6] 10 11] [2 10 11] [2 10 11] 8] [2 12] [8 10 11] 8] [2 10 13] [2 10 13] [6 11 13] DỮ LIỆU ĐỂ KIỂM TRA BÀI Tệp DTDG.INP -2 -1 -1 -1 1 0 -1 0 -1 -1 -1 1 0 -1 -1 -2 -1 -1 1 0 -1 -5 -1 -1 -1 1 0 -1 -1 -2 -1 2 -1 -2 -2 -1 -1 3 -1 -1 -1 -2 -2 -1 1 0 -1 -5 -1 -1 1 0 -1 -1 -3 -5 -4 0 -1 -5 Tệp DTDG.OUT -1 -1 -1 -1 DỮ LIỆU ĐỂ KIỂM TRA BÀI Tệp PSONG.INP (Chú ý: liệu kéo dài đến dòng thục vào tiếp theo) 3.0 3.0 3.0 1.5 1.5 1.0 3.0 2.0 4.0 1.5 2.5 3.0 4.0 4.0 3.0 5.0 6.0 3.0 5.5 4.5 1.0 4.0 4.0 2.5 6.0 5.0 3.5 4.5 5.5 1.5 4.0 4.0 3.5 5.0 6.0 3.0 5.5 4.5 1.5 6.5 5.5 2.5 4.0 4.0 3.0 6.0 5.0 3.5 4.5 5.5 2.0 5.5 6.5 2.5 4.0 4.0 4.0 5.0 6.0 4.0 5.5 4.5 3.0 6.5 5.5 5.0 5.5 6.5 7.0 4.0 4.0 3.0 6.0 5.0 4.0 4.5 5.5 2.0 5.5 6.5 3.0 6.5 5.5 5.0 4.0 4.0 3.0 6.0 5.0 3.0 4.5 5.5 1.0 5.5 6.5 2.0 6.5 5.5 2.0 3.5 5.5 2.0 Tệp PSONG.OUT 36.00 67.00 52.00 54.00 65.25 62.75 196.00 103.00 54.75 MÃ NGUỒN BÀI (Tệp DT3D.PAS) {Phan khai bao chung} uses crt; Const tentep = 'DT3D'; {Ten tep} Nmax = 255; DTmax = 100; Var D: array[1 Nmax, 2] of integer; {Dau vao} sl: integer; {Dau ra} DT: array[1 Nmax, 2] of byte; {Dau ra} N: integer; {Bien phu} CD: array[1 nmax] of 1; {Mang phu} pa, sopa: byte; {QL phuong an} f: text; {Bien tep} {Dem so phuong an} procedure DemSoPA; BEGIN assign(f,tentep+'.INP'); reset(f); sopa:=0; while (not eof(f)) begin readln(f); inc(sopa); end; close(f); END; {Doc du lieu cua phuong an pa} procedure DocPA; var i:integer; BEGIN assign(f,tentep+'.INP'); reset(f); i:=1; while (i<pa) begin readln(f); inc(i); end; N := 0; while (not eoln(f)) begin inc(N); read(f,D[N,1],D[N,2]); end; close(f); END; {Xet diem, true: thang hang, false: khong thang hang} function TH(ii ,jj, kk: byte): boolean; var dx, dy: integer; BEGIN dx := D[jj,1] - D[ii,1]; dy := D[jj,2] - D[ii,2]; if dx*(D[kk,2] - D[ii,2]) = dy*(D[kk,1] - D[ii,1]) then TH := true else TH := false; END; {Tim so duong thang qua it nhat diem} Procedure TimSL; var n1, n2, i, i1, j, j1, k: integer; BEGIN SL := 0; n1 := N - 1; n2 := N - 2; fillchar(CD, sizeof(CD),0); for i := to n2 {Diem thu 1} begin i1 := i + 1; for j := i1 to n1 if CD[j] = then {Diem thu 2} begin j1 := j + 1; for k := j1 to n if TH(i,j,k) then {Diem thu 3} begin inc(sl); DT[sl,1] := i; DT[sl,2] := j; CD[j] := 1; CD[k] := 1; break; end; end; end; END; {Ghi ket qua man hinh va tep tin} procedure GhiKQ; var i, id, j: integer; BEGIN assign(f,tentep+'.OUT'); if pa=1 then rewrite(f) else append(f); write(SL); write(f,SL); if SL > then for i := to SL begin write(' [',DT[i,1],' ',DT[i,2]); write(f,' [',DT[i,1],' ',DT[i,2]); id := DT[i,2] + 1; for j := id to N if TH(DT[i,1], DT[i,2],j) then begin write(' ',j); write(f,' ',j); end; write(']'); write(f,']'); end; writeln; writeln(f); close(f); END; {Chuong trinh chinh} BEGIN clrscr; DemSoPA; for pa:=1 to sopa begin DocPA; TimSL; GhiKQ; end; writeln; write('Da xong Nhan ENTER de thoat'); Readln; END GIẢI THÍCH MÃ NGUỒN BÀI 1 Chương trình chính Sau gọi lệnh DemSoPA, vòng lặp For gọi các lệnh DocPA, TimSL, và GhiKQ Lệnh DemSoPA: đếm số phương án có tệp liệu Lệnh DocPA: đọc liệu phương án thứ pa Lệnh GhiKQ: ghi kết màn hình và tệp tin Lệnh TimSL: Tìm số lượng đường thẳng theo yêu cầu, thuật giải đã nêu phần Hướng dẫn cách giải MÃ NGUỒN BÀI (Tệp DTDG.PAS) {Phan khai bao chung} uses crt; const tentep = 'DTDG'; {Ten tep} var px, py: integer; {Dau vao} X,Y: Array[1 2000] of integer; {Dau vao} vtp: -1 1; {Dau ra} N, mx, my: integer; {Bien phu} sopa, pa: integer; f:text; var i: integer; BEGIN for i := to N if (px = X[i]) and (py = Y[i]) then begin TrungDinh := true; Exit; end; TrungDinh := false; END; {P thuoc canh cua da giac} function ThuocCanh: boolean; var i: integer; BEGIN for i := to N if (px-X[i])*(Y[i+1] - Y[i]) = (py-Y[i])*(X[i+1] X[i]) then {P, i, i+1 thang hang} if ((px - X[i])*(px-X[i+1])<=0) and ((py Y[i])*(py-Y[i+1])<=0) then {P thuoc canh i} begin ThuocCanh := true; Exit; end; ThuocCanh := false; END; {Tim diem M ngoai da giac va PM khong chua dinh cua da giac} procedure TimM; var i: integer; cc: boolean; {con chinh mx} BEGIN mx := -maxint; my := -maxint; for i := to N begin if mx < X[i] then mx := X[i]; if my < Y[i] then my := Y[i]; end; inc(mx); repeat cc := false; for i:=1 to N if (my-py)*(X[i] - px) = (mx-px)*(Y[i] - py) then begin {P, M, i thang hang} inc(mx); cc := true; break; end; until (not cc); END; {huong giua diem: 1: nguoc, -1: thuan, 0: thang hang} function ccw(x1,y1,x2,y2,x3,y3: integer): integer; var dx1, dx2, dy1, dy2: integer; BEGIN dx1 := x2 - x1; dx2 := x3 - x2; dy1 := y2 - y1; dy2 := y3 - y2; if dx1*dy2 > dy1*dx2 then ccw := 1; {nguoc} if dx1*dy2 < dy1*dx2 then ccw := -1; {thuan} if dx1*dy2 = dy1*dx2 then ccw:=0; {thang hang} END; {Doan thang MP cat canh thu ii} function Cat(ii: integer): boolean; var x1, y1, x2, y2, cc1, cc2, cc3, cc4: integer; BEGIN x1 := X[ii]; y1 := Y[ii]; x2 := X[ii+1]; y2 := Y[ii+1]; cc1 := ccw(x1,y1,x2,y2,px,py); {1 P} cc2 := ccw(x1,y1,x2,y2,mx,my); {1 M} cc3 := ccw(px,py,mx,my,x1,y1); {P M 1} cc4 := ccw(px,py,mx,my,x2,y2); {P M 2} Cat := (cc1*cc2 < 0) and (cc3*cc4 < 0); END; {vi tri cua diem P: 1: trong, -1: ngoai, 0: thuoc canh} procedure TimVTP; var i, sdc: integer; BEGIN X[N+1] := X[1]; Y[N+1] := Y[1]; vtp := 0; if (TrungDinh or ThuocCanh) then Exit else begin TimM; sdc := 0; for i := to N if Cat(i) then inc(sdc); if sdc mod = then vtp := else vtp := -1; end; END; {Ghi ket qua man hinh va tep tin} procedure GhiKQ; BEGIN Assign(f, tentep+'.OUT'); if (pa = 1) then Rewrite(f) else Append(f); writeln(vtp); writeln(f,vtp); Close(f); END; {Chuong trinh chinh} BEGIN clrscr; DemSoPA; for pa:=1 to soPA begin DocPA; TimVTP; GhiKQ; end; writeln; write('Da xong Nhan ENTER de thoat.'); readln; END GIẢI THÍCH MÃ NGUỒN BÀI Chương trình chính Sau gọi lệnh DemSoPA, vòng lặp For gọi các lệnh DocPA, TimVTP, và GhiKQ Lệnh DemSoPA: đếm số phương án có tệp liệu Lệnh DocPA: đọc liệu phương án thứ pa Lệnh GhiKQ: ghi kết màn hình và tệp tin Lệnh TimVTP: Tìm vị trí điểm P theo yêu cầu, thuật giải đã nêu phần Hướng dẫn cách giải MÃ NGUỒN BÀI (Tệp PSONG.PAS) {Phan khai bao chung} uses crt; const tentep = 'PSONG'; Nmax = 1000; {Ten tep} type MangST = array[1 Nmax] of real; var X, Y, R: MangST; dt: real; {Dau vao} {Dau ra} {hcn 1-2 chua hcn 3-4} function NT(x1,y1,x2,y2,x3,y3,x4,y4: real): boolean; BEGIN if (x1<=x3) and (x4<=x2) and (y1<=y3) and (y4<=y2) then NT := true else NT := false; END; {doi cho cac phan tu ii va jj cua mang mm} procedure DoiCho(var mm: MangST; ii,jj: integer); var t: real; BEGIN t := mm[ii]; mm[ii] := mm[jj]; mm[jj] := t; END; {tao va sap xep mang xx tang dan} procedure SXTang(mm: MangST; var xt: MangST); var i, i1, j: integer; BEGIN for i := to N begin xt[2*i-1] := mm[i] - R[i]; xt[2*i] := mm[i] + R[i]; end; n2 := 2*N; n21 := n2 -1; for i:= to n21 begin i1 := i+1; for j := i1 to n2 if (xt[i]>xt[j]) then DoiCho(xt,i,j); end; END; {Tinh dien tich phu song} Procedure TinhDT; var i, j, k: integer; x1, y1, x2, y2, x3, y3, x4, y4: real; BEGIN SXTang(X, XX); SXTang(Y, YY); dt :=0; n21 := 2*N-1; for i := to n21 for j := to n21 for k := to N begin x1 := XX[i]; y1 := YY[j]; x2 := XX[i+1]; y2 := YY[j+1]; x3 := X[k]-r[k]; y3 := Y[k]-r[k]; x4 := X[k]+r[k]; y4 := Y[k]+r[k]; if NT(x3, y3, x4, y4, x1, y1, x2, y2) then begin dt := dt + abs((x1-x2)*(y1-y2)); break; end; end; END; {Ghi ket qua man hinh va tep tin} procedure GhiKQ; BEGIN assign(f, tentep+'.OUT'); if pa = then rewrite(f) else append(f); writeln(dt:0:2); writeln(f,dt:0:2); close(f); END; {Chuong trinh chinh} BEGIN clrscr; DemSoPA; for pa:=1 to SoPA begin DocPA; TinhDT; GhiKQ; end; writeln; write('Da xong Nhan ENTER de thoat.'); readln; END GIẢI THÍCH MÃ NGUỒN BÀI Chương trình chính Sau gọi lệnh DemSoPA, vòng lặp For gọi các lệnh DocPA, TinhDT, và GhiKQ Lệnh DemSoPA: đếm số phương án có tệp liệu Lệnh DocPA: đọc liệu phương án thứ pa Lệnh GhiKQ: ghi kết màn hình và tệp tin Lệnh TinhDT: CTính tổng diện tích phủ sóng theo yêu cầu, thuật giải đã nêu phần Hướng dẫn cách giải