Metody programowania / Ćwiczenia 6

Zadania z programowania dynamicznego.

Oglądaj wskazówki i rozwiązania __SHOWALL__
Ukryj wskazówki i rozwiązania __HIDEALL__

Zadanie 1 (Odcinki)

Dana jest tablica A zawierająca posortowane rosnąco liczby rzeczywiste oznaczające współrzędne punktów na prostej. Oblicz maksymalną sumę długości rozłącznych odcinków (tzn. nie stykających się nawet końcami), które można otrzymać łącząc wybrane sąsiednie punkty. Nie ma obowiązku użycia wszystkich punktów i nie wolno zmieniać zawartości tablicy.

Wskazówka

Przeglądaj A od lewej od prawej, pamiętając dwie liczby: najlepszą sumę używającą ostatniego punktu i najlepszą sumę nieużywającą ostatniego punktu.

Rozwiązanie

function Odcinki(N:integer; var A:array[1..N] of real) : real;
// N>=0; A posortowana
var i : integer;
   zost, bezost, zost_pom, bezost_pom : real;
begin
  zost:=0;
  bezost:=0;
  // zost i bezost oznaczaja maksymalna sumę dlugosci z i bez punktu i-1
  for i:=2 to N do begin
    bezost_pom:=max(zost,bezost);
    zost_pom:=bezost+A[i]-A[i-1]; 
      // dokladamy odcinek A[i-1]..A[i]
      // możemy bo odcinki odpowiadające sumie bezost nie blokują i-1 
    bezost:=bezost_pom;
    zost:=zost_pom;
  end;
  // zost i bezost oznaczaja maksymalna sumę dlugosci z i bez punktu N
  // wystarczy wybrać tą większą
  Odcinki:=max(zost,bezost);
end // Odcinki

Koszt czasowy: liniowy względem N.

Koszt pamięciowy: stały

Poprawność rozwiązania: poprawność rozwiązania opiera się na poprawności następujących rekurencyjnych wzorów:

bezost[1]=0
zost[1]=0
bezost[i]=max(zost[i-1],bezost[i-1])
zost[i]=bezost[i-1]+A[i]-A[i-1]

gdzie bezost[i] (odp. zost[i]) oznaczaja największą sumę odcinków od początku tablicy A aż do punktu A[i], nieblokujących (odp. blokujących) punkt A[i]. Istotnie jeśli chcemy mieć największą sumę odcinków nieblokujących i to można wziąć albo blokujące i-1 albo nieblokujące. Jeśli chcemy największą sumę odcinków blokujących i, to do największej sumy odcinków nieblokujących i-1 należy dodać długość odcinka z i-1 do i. Wartości początkowe są oczywiście poprawne (budząca wątpliwość wartość zost[1] nie jest istotna, bo w następnej operacji zost[1] występuje tylko w wyrażeniu max(zost[1],bezost[1]), więc jego wartość 0 niczego nie popsuje).

Ponieważ, jak widać ze wzorów do obliczenia wartości dla i wystarczy pamiętać wartości dla i-1, rezygnujemy z pamiętania całych tablic zost i bezost.

Rozwiązanie 2

function Odcinki2(N:integer; var A:array[1..N] of real) : real;
// N>=0; A posortowana
var i : integer;
   mozezost, bezost, mozezost_pom, bezost_pom : real;
begin
  mozezost:=0;
  bezost:=0;
  // mozezost i bezost oznaczaja maksymalna sumę dlugosci odcinków 
  // mogących i niemogących blokować punkt i-1
  for i:=2 to N do begin
    bezost_pom:=mozezost;
    mozezost_pom:=max(bezost+A[i]-A[i-1],mozezost);
      // albo bieżemy odcinki nieblokujące i-1 oraz odcinek A[i-1]..A[i]
      // albo odcinki potencjalnie blokujące i-1
    bezost:=bezost_pom;
    mozezost:=mozezost_pom;
  end;
  Odcinki2:=mozezost;
end // Odcinki2

Koszt czasowy: liniowy względem $N$.

Koszt pamięciowy: stały

Poprawność rozwiązania: poprawność rozwiązania opiera się na poprawności następujących rekurencyjnych wzorów:

bezost[1]=0
mozezost[1]=0
bezost[i]=mozezost[i-1]
mozezost[i]=max(bezost[i-1]+A[i]-A[i-1],mozezost[i-1])

gdzie bezost[i] (odp. mozezost[i]) oznaczaja największą sumę odcinków od początku tablicy A aż do punktu A[i], które nie mogą blokować (odp. mogą blokować) punkt A[i]. W tym rozwiązaniu pozbywamy się wątpliwości związanych z zost[1].

Ćwiczenie

Jak zmodyfikować nasz program, aby zamiast liczenia maksymalnej sumy, wypisywał ciąg odcinków o maksymalnej sumie?

Zadanie 2 (Kwadrat)

Dana jest tablica A typu array[1..N,1..M] of integer. Oblicz rozmiar (długość boku) największego kwadratu w tej tablicy składającego się z samych zer.

Wskazówka 1

{{{3}}}

Rozwiązanie 1

{{{3}}}

Wskazówka 2

{{{3}}}

Rozwiązanie 2

{{{3}}}

Zadanie 3 (Prostokąt)

Dana jest tablica A typu array[1..N,1..M] of integer. Oblicz maksymalne pole prostokąta składającego się z samych zer w tej tablicy.

Wskazówka 1

{{{3}}}

Wskazówka 2

{{{3}}}

Rozwiązanie 1

{{{3}}}