Bài 6: Thuật toán loang trên ma trận

Bài viết này là phần [part not set] trong 7 bài của Series Lý thuyết đồ thị căn bản

Thuật toán loang (Thuật toán vết dầu loang) là một trong những thuật toán được dùng khá nhiều trong tin học, điển hình là thuật toán loang trên ma trận này được ứng dụng để đếm số thành phần liên thông trên ma trận. Nó trong các trò chơi nổi tiếng như line 98, trò chơi dò mìn trên windows, tìm đường đi ngắn nhất trên ma trận…

Nội dung bài viết

1. Giới thiệu thuật toán vết dầu loang
2. Code loang đếm thành phần liên thông c++
3. Bài tập minh họa

1. Giới thiệu thuật toán vết dầu loang

Thuật toán loang có tên gọi như vậy vì nguyên lý hoạt động của nó giống vết dầu loang. Từ một vết dầu nhỏ sẽ được loang ra xung quanh. Thuật toán loang trên ma trận cũng vậy, bạn sẽ duyệt một ô trên ma trận và sau đó duyệt các điểm xung quanh nó và dần loang ra để giải quyết bài toán.

Như trong trò chơi dò mìn, để giải quyết được trò chơi, bạn phải chọn một điểm ban đầu, và xử lí các thông tin tại đó để có thể dò được các bãi mìn xung quanh, và bạn sẽ làm như vậy đến khi không còn loang được nữa.

Tư tưởng của thuật toán loang chỉ có vậy, trong bài viết này mình sẽ cố gắng giới thiệu bạn rõ hơn về việc ứng dụng nó qua một vài ví dụ.

VBGRASS spoj – Bãi cỏ ngon nhất

Bài trên bạn có thể dùng thuật toán BFS hoặc DFS để loang ra từng vùng có các kí tự giống nhau. cách tổ chức code khi loang ra bạn có thể tham khảo như bài này https://kienthuc24h.com/co-ban-ung-dung-bfs-de-giai-quyet-bai-tap-duong-di-cua-quan-ma-trong-thi/ bạn sẽ khai báo như sau

int cld[8]={1,-1,0,0};
int clc[8]={0,0,1,-1};

1 2	int cld[8]={1,-1,0,0}; int clc[8]={0,0,1,-1};

Và tiến hành loang và đếm thành phần liên thông theo hướng làm như sau:

2. Code loang đếm thành phần liên thông c++

res = 0;
for (int i=1; i<=r; i++)
    for (int j=1; j<=c; j++)
        if (chưa duyệt ô i,j) 
        {
            res++;
            bfs(i,j);
        }

res = 0;

for (int i=1; i<=r; i++)

for (int j=1; j<=c; j++)

if (chưa duyệt ô i,j)

{

res++;

bfs(i,j);

}

3. Bài tập minh họa

a. Đề bài toán thứ 2

Link đề bài http://www.spoj.com/THPTCBT/problems/MTKPMANG/

Khôi phục mảng

Chúng ta sẽ có một ma trận gồm các số 0,1. Và từ ma trận A chúng ta sẽ tạo được ma trận B với B[i,j] là số lượng ô có giá trị bằng 1 kề cạnh với A[i,j]

Cho trước mảng B, hãy khôi phục mảng A.

Nếu không thể khôi phục thì ghi ra dòng chữ “KHONG THE”

b. Example

Input:
4 4
1 2 1 1
2 2 3 1
3 2 2 2
0 3 1 1

Output:
1 0 0 1
1 1 0 1
0 1 1 0
1 0 1 0

Bài toán trên cũng là tư tưởng loang nhưng lúc này chúng ta không tiếp cận chúng từ 1 vị trí, mà dựa vào các đặc điểm tối thiểu của nó để kết luận kết quả bài toán.

Ví dụ trong bài toán trên, bạn có thể đi lên từ các trường hợp, số 0 ở các góc, khi đó bạn điền được 2 số 1, số 3 ở cạnh, trường hợp số 0 điền được 4 số 1,…. ý tưởng tổng quát là vậy.

c. Code lời giải

const   fi='';
        nmax=50;
        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,B:array[0..nmax+1,0..nmax+1] of data;
        n,m:data;
        kt,free:array[0..nmax+1,0..nmax+1] of boolean;
        sl:data;
 
procedure docfile;
var     i,j:data;
begin
        assign(f,fi); reset(f);
        readln(f,m,n);
        for i:=0 to m+1 do
                for j:=0 to n+1 do
                        b[i,j]:=2;
 
        for j:=0 to n+1 do
                begin
                        Free[0,j]:=true;
                        Free[m+1,j]:=true;
                        kt[0,j]:=true;
                        kt[m+1,j]:=true;
                        B[0,j]:=0;
                        B[m+1,j]:=0;
                end;
        for i:=0 to m+1 do
                begin
                        free[i,0]:=true;
                        free[i,n+1]:=true;
                        kt[i,0]:=true;
                        kt[i,n+1]:=true;
                        B[i,0]:=0;
                        B[i,n+1]:=0;
                end;
 
        for i:=1 to m do
                for j:=1 to n do
                        read(f,a[i,j]);
        close(f);
end;
 
procedure xuat;
var     i,j:data;
begin
        for i:=1 to m do
                begin
                        for j:=1 to n do
                                if b[i,j]=2 then
                                        write('. ')
                                else
                                write(b[i,j],' ');
                        writeln;
                end;
end;
 
procedure dem(i,j:data; var khong,mot,trong:data);
var     k,x,y:data;
begin
        khong:=0;
        mot:=0;
        for k:=1 to 4 do
                begin
                        x:=i+cld[k];
                        y:=j+clc[k];
                        if B[x,y]=0 then
                                inc(khong)
                        else
                                if b[x,y]=1 then
                                        inc(mot);
                end;
        trong:=4-khong-mot;
end;
 
procedure fill(i,j,c:data);
var     k,x,y:data;
begin 
        kt[i,j]:=true;
        for k:=1 to 4 do
                begin
                        x:=i+cld[k];
                        y:=j+clc[k];
                        if B[x,y]=2 then
                                begin
                                        B[x,y]:=c;
                                        inc(sl);
                                end;
                end;
end;
 
procedure update;
var     i,j,khong,mot,trong:data;
begin
        for i:=1 to m do
                for j:=1 to n do
                        if kt[i,j]=false then
                                begin
                                        dem(i,j,khong,mot,trong);
                                        if (a[i,j]=trong) and(khong=4-trong) then
                                                fill(i,j,1)
                                        else
                                                if trong+mot=a[i,j] then
                                                        fill(i,j,1)
                                                else
                                                        if (mot=khong+trong) and (mot=a[i,j]) then
                                                                fill(i,j,0)
                                                        else
                                                                if (mot=a[i,j]) then
                                                                        fill(i,j,0)
                                                                else
                                                                        if 4-khong=a[i,j] then
                                                                                fill(i,j,1); 
                                end;
end;
 
 
procedure duyet;
var     i,j:data;
        count:data;
begin
        sl:=0;
        repeat
                count:=sl;
                update;
        until count=sl;
        if sl=m*n then
        	xuat
        else
        	writeln('KHONG THE');
end;
 
 
begin
        docfile;
        duyet; 
end.

100

101

102

103

104

105

106

107

108

109

110

111

112

113

114

115

116

117

118

119

120

121

122

123

124

125

126

127

128

129

130

131

132

133

134

135

136

137

138

139

140

141

const fi='';

nmax=50;

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,B:array[0..nmax+1,0..nmax+1] of data;

n,m:data;

kt,free:array[0..nmax+1,0..nmax+1] of boolean;

sl:data;

procedure docfile;

var i,j:data;

begin

assign(f,fi); reset(f);

readln(f,m,n);

for i:=0 to m+1 do

for j:=0 to n+1 do

b[i,j]:=2;

for j:=0 to n+1 do

begin

Free[0,j]:=true;

Free[m+1,j]:=true;

kt[0,j]:=true;

kt[m+1,j]:=true;

B[0,j]:=0;

B[m+1,j]:=0;

end;

for i:=0 to m+1 do

begin

free[i,0]:=true;

free[i,n+1]:=true;

kt[i,0]:=true;

kt[i,n+1]:=true;

B[i,0]:=0;

B[i,n+1]:=0;

end;

for i:=1 to m do

for j:=1 to n do

read(f,a[i,j]);

close(f);

end;

procedure xuat;

var i,j:data;

begin

for i:=1 to m do

begin

for j:=1 to n do

if b[i,j]=2 then

write('. ')

else

write(b[i,j],' ');

writeln;

end;

procedure dem(i,j:data; var khong,mot,trong:data);

var k,x,y:data;

begin

khong:=0;

mot:=0;

for k:=1 to 4 do

begin

x:=i+cld[k];

y:=j+clc[k];

if B[x,y]=0 then

inc(khong)

else

if b[x,y]=1 then

inc(mot);

end;

trong:=4-khong-mot;

end;

procedure fill(i,j,c:data);

var k,x,y:data;

begin

kt[i,j]:=true;

for k:=1 to 4 do

begin

x:=i+cld[k];

y:=j+clc[k];

if B[x,y]=2 then

begin

B[x,y]:=c;

inc(sl);

end;

procedure update;

var i,j,khong,mot,trong:data;

begin

for i:=1 to m do

for j:=1 to n do

if kt[i,j]=false then

begin

dem(i,j,khong,mot,trong);

if (a[i,j]=trong) and(khong=4-trong) then

fill(i,j,1)

else

if trong+mot=a[i,j] then

fill(i,j,1)

else

if (mot=khong+trong) and (mot=a[i,j]) then

fill(i,j,0)

else

if (mot=a[i,j]) then

fill(i,j,0)

else

if 4-khong=a[i,j] then

fill(i,j,1);

end;

procedure duyet;

var i,j:data;

count:data;

begin

sl:=0;

repeat

count:=sl;

update;

until count=sl;

if sl=m*n then

xuat

else

writeln('KHONG THE');

end;

begin

docfile;

duyet;

end.

d. Một số bài tập khác liên quan

bài giải MTNTRAI spoj THPTCBT – 21697. Nông Trại

BCLKCOUN spoj PTIT – Đếm số ao

BCISLAND PTIT spoj – Nước biển

Bài 6: Thuật toán loang trên ma trận

1. Giới thiệu thuật toán vết dầu loang

2. Code loang đếm thành phần liên thông c++

3. Bài tập minh họa

a. Đề bài toán thứ 2

b. Example

c. Code lời giải

d. Một số bài tập khác liên quan

Bài viết liên quan

Trả lời Hủy