Nguồn đề bài: http://www.spoj.com/PTIT/problems/BCISLAND/
1. Đề bài BCISLAND spoj
Trái đất nóng lên kéo theo mực nước biển dâng. Hòn đảo nhỏ Gonnasinka thuê bạn để dự báo trước hiểm họa này. Cho trước 1 lưới tọa độ thể hiện cao độ của đảo, hãy giúp họ tính toán xem nước biển dâng cao bao nhiêu thì hòn đảo sẽ bị chia cắt.
Input
Input gồm nhiều bộ test, mỗi bộ test bao gồm:
- Dòng đầu ghi 2 số n, m là chiều dài và chiều rộng.
- Sau đó là n dòng, mỗi dòng gồm m số, mỗi số cho biết độ cao của ô đó, giá trị 0 chỉ mực nước biển. Những ô giá trị 0 dọc theo đường viền và những ô số 0 liền kề những ô này chỉ mặt biển. Những ô có giá trị 0 còn lại (được bao bọc bởi các số > 0) là đất liền bên trong đảo nhưng có độ cao ngang mặt nước biển. Hòn đảo lúc đầu chưa bị chia cắt. Số n và m không lớn hơn 100 và độ cao không lớn hơn 1000.
- Dòng cuối cùng của input chứa 2 số 0
Output
Với mỗi bộ test, in ra:
Case n: Island splits when ocean raises f feet. (Đảo bị chia khi nước biển dâng cao f feet)
Hoặc
Case n: Island never splits. (Đảo không bao giờ bị chia cắt)
Example
Input:
5 5
3 4 3 0 0
3 5 5 4 3
2 5 4 4 3
1 3 0 0 0
1 2 1 0 0
5 5
5 5 5 5 7
4 1 1 1 4
4 1 2 1 3
7 1 0 0 4
7 3 4 4 4
0 0
Output:
Case 1: Island never splits.
Case 2: Island splits when ocean rises 3 feet.
######################################
2. Gợi ý BCISLAND spoj
– Dùng thuật toán DFS để loang ra đếm số TLPT với mực nước là L nào đó. (L=[min(a[i,j])..max(a[i,j])])
– Nếu số thành phần liên thông lớn hơn 1 mà đó ko phải là nước biển thì đó là KQ bài toán.
Bài này mình ko biết phải viết hướng dẫn như thế nào nữa. 😀 các bạn coi code hoặc gợi ý đỡ vậy. 😀
3. code tham khảo BCISLAND spoj
const fi='';
nmax=100;
cld:array[1..4] of shortint=(-1,1,0,0);
clc:array[1..4] of shortint=(0,0,1,-1);
type data=longint;
var
f:text;
A:array[-1..nmax+2,-1..nmax+2] of data;
dd:array[-1..nmax+2,-1..nmax+2] of data;
n,m,L:data;
test,gtmin,gtmax:data;
procedure dfs(u,v:data);
var x,y,i:data;
begin
dd[u,v]:=1;
for i:=1 to 4 do
begin
x:=u+cld[i];
y:=v+clc[i];
if (dd[x,y]=0) and (a[x,y]<=L) then
DFS(x,y);
end;
end;
procedure dfs2(u,v:data);
var x,y,i:data;
begin
dd[u,v]:=2;
for i:=1 to 4 do
begin
x:=u+cld[i];
y:=v+clc[i];
if dd[x,y]=0 then
dfs2(x,y);
end;
end;
function check(x:data):boolean;
var i,j,stplt:data;
begin
L:=x;
for i:=0 to m+1 do
for j:=0 to n+1 do
dd[i,j]:=0;
dfs(0,0);
stplt:=0;
for i:=1 to m do
for j:=1 to n do
if dd[i,j]=0 then
begin
inc(stplt);
dfs2(i,j);
end;
if stplt>1 then
exit(true);
exit(false);
end;
procedure xuli;
var i,j,dau,cuoi,mid:data;
begin
// dien so 0 xung qanh mang A
for i:=0 to m+1 do
begin
a[i,0]:=0;
a[i,n+1]:=0;
end;
for j:=0 to n+1 do
begin
a[0,j]:=0;
a[m+1,j]:=0;
end;
// tao vien ko cho Loang ra ngoai bang
for i:=-1 to m+2 do
begin
dd[i,-1]:=1002;
dd[i,n+2]:=1002;
end;
for j:=-1 to n+2 do
begin
dd[-1,j]:=1002;
dd[m+2,j]:=1002;
end;
for i:=gtmin to gtmax do
if check(i) then
begin
writeln('Case ',test,': Island splits when ocean rises ',i,' feet.');
exit;
end;
writeln('Case ',test,': Island never splits.');
end;
function min(a,b:data):data;
begin
if a<b then exit(a); exit(b);
end;
function max(a,b:data):data;
begin
if a>b then exit(a); exit(b);
end;
procedure docfile;
var i,j:data;
begin
assign(f,fi); reset(f);
test:=0;
repeat
readln(f,m,n);
inc(test);
gtmax:=-10000;
gtmin:=10000;
if (n=0) and (m=0) then break;
for i:=1 to m do
begin
for j:=1 to n do
begin
read(f,a[i,j]);
gtmax:=max(gtmax,a[i,j]);
gtmin:=min(gtmin,a[i,j]);
end;
readln(f);
end;
xuli;
until false;
close(f);
end;
begin
docfile;
end.
