Unelte utilizator

Unelte site


laboratoare:laborator-07

Diferențe

Aici sunt prezentate diferențele dintre versiunile selectate și versiunea curentă a paginii.

Link către această vizualizare comparativă

Ambele părți revizuirea anterioară Versiuni anterioare
Urmatoarea versiune
Versiuni anterioare
Ultima versiune Ambele părți următoarea reviziune
laboratoare:laborator-07 [2017/04/01 16:25]
mihai.iacov [3.1 Algoritmul lui Dijkstra]
laboratoare:laborator-07 [2017/04/07 16:55]
iulian.matesica [6. Exercițîi laborator]
Linia 1: Linia 1:
-====== Laborator 07: Parcurgerea grafurilor ======+====== Laborator 07: Drumuri de cost minim ======
  
 ===== 1.Obiective laborator ====== ===== 1.Obiective laborator ======
Linia 23: Linia 23:
  
 <note important>Pot exista mai multe drumuri de cost minim de la S la D.</note> <note important>Pot exista mai multe drumuri de cost minim de la S la D.</note>
 +
 +====2.3 Legătura muchii-arce: ====
 +<note tip>
 +Orice graf **neorientat** este echivalent cu un graf **orientat** dacă înlocuim fiecare **muchie** cu **două arce**(de acelaşi cost, câte un arc pentru fiecare sens).
 +
 +Algoritmii următori pot fi folosiţi(în limita observaţiilor finale) pe grafuri orientate la fel ca pe grafuri neorientate.
 +</note>
 +
 +
 ===== 3.Drumul de cost minim cu sursă unică ====== ===== 3.Drumul de cost minim cu sursă unică ======
 Următorii algoritmi caută drumurile de cost minim de la **un singur nod(sursă)** la **toate celelalte noduri**. Rezultatul acestor algoritmi este un **arbore cu drumuri de cost minim**, unde: Următorii algoritmi caută drumurile de cost minim de la **un singur nod(sursă)** la **toate celelalte noduri**. Rezultatul acestor algoritmi este un **arbore cu drumuri de cost minim**, unde:
Linia 44: Linia 53:
 <note tip>Construind astfel algoritmul, este **garantat** că, în momentul în care explorăm un nod(C), estimarea pentru costul drumului S-C este chiar **costul minim**.</note> <note tip>Construind astfel algoritmul, este **garantat** că, în momentul în care explorăm un nod(C), estimarea pentru costul drumului S-C este chiar **costul minim**.</note>
  
 +<note important>
 +Pentru grafuri orientate, ne referim la vecinii la care putem ajunge printr-un arc ce pleacă din C: (C,V).
 +</note>
 Paşii algoritmului lui Dijkstra: Paşii algoritmului lui Dijkstra:
 <code> <code>
Linia 67: Linia 79:
  
 ==== 3.2 Algoritmul Bellman-Ford ==== ==== 3.2 Algoritmul Bellman-Ford ====
 +Principii similare pentru algoritmul Bellman-Ford:
 +  * vom construi arborele, începând de la nodul S;
 +  * atribuim un posibil cost(o estimare a distanţei) pentru fiecare nod. (iniţial, S are costul 0, toate celelalte noduri au costul **infinit**);
 +  * la fiecare evaluare, dacă găsim o nouă estimare de cost mai bună decât cea precedentă, folosim, mai departe, noua estimare. Dacă dorim să ţinem evidenţa muchiilor folosite, actualizăm şi nodul părinte.
 +  * funcţia de estimare a costului este definită la fel ca la algoritmul lui Dijkstra(costul drumului de la S la nodul respectiv)
 +
 +<note tip>Rezultatul algoritmului va fi un arbore cu N noduri(unde N = |V|, numărul de noduri din graf). Prin urmare, **lungimea** oricărui drum de la S la alt nod va fi **maxim N-1**. (Presupunem că nu se repetă muchii)
 +</note>
 +
 +Diferenţa apare, în algoritmul Bellman-Ford, la alegerea nodurilor pentru care facem evaluarea:
 +  * Algoritmul nu are preferinţe pentru anumite noduri şi nu extrage, la fiecare pas, cel mai bun candidat.
 +  * În schimb, acest algoritm evaluează **toate muchiile** la un pas. Folosindu-se de principiul de mai sus, (N-1) astfel de paşi vor fi suficienţi.
 +
 +<note important>
 +În cazul grafurilor neorientate, evaluăm fiecare muchie în ambele sensuri.
 +</note>
 +
 +Paşii algoritmului Bellman-Ford:
 +<code>
 +1. Atribuim fiecărui nod o estimare iniţială a costului:
 + 0 pentru nodul sursă(S);
 + infinit pentru toate celelalte;
 +2. Executăm de N-1 ori:
 + 1. Pentru fiecare pereche (u, v) a.i. există muchie de la u la v
 + 1. calculăm noua estimare de cost = cost(drumul S-u) + cost(muchia (u,v));
 + 2. comparăm noua estimare cu vechiul cost(drumul S-v):
 + dacă noul cost e mai bun
 + 1. actualizăm cost(drumul S-v) = noul cost;
 + 2. actualizăm parinte(v) = u; (pentru păstrarea muchiei folosite)
 + altfel păstrăm vechiul cost
 +</code>
 +
 ===== 4.Drumul de cost minim între oricare 2 noduri ====== ===== 4.Drumul de cost minim între oricare 2 noduri ======
 Următorul algoritm caută cel mai scurt drum de la **orice nod sursă(S)** la **fiecare nod destinaţie(D)**. Următorul algoritm caută cel mai scurt drum de la **orice nod sursă(S)** la **fiecare nod destinaţie(D)**.
Linia 72: Linia 116:
  
 ==== 4.1 Algoritmul Floyd-Warshall ==== ==== 4.1 Algoritmul Floyd-Warshall ====
-===== 5.Exiciții laborator ======+Rezultatul algoritmului este o matrice N x N(unde N = |V|, numărul de noduri), iar valorea din matrice de la poziţia [i][j] va fi costul minim pentru drumul i-j. Fie această matrice numită **dist**. 
 + 
 +Pentru simplitate, algoritmul numerotează nodurile de la 1 la N şi foloseşte următoarea construcţie: 
 +  * defineşte o funcţie costMinim(i,j,k) = cel mai ieftin drum care: 
 +  *    -pleacă din i 
 +  *    -ajunge în j 
 +  *    -în afară de primul şi de ultimul nod, conţine numai noduri din {1,2,3,...,k}(primele k noduri). 
 + 
 +Algoritmul calculează costMinim(i,j,k) pentru toate perechile (i,j,k), folosind formula: **costMin(i,j,k+1) = min(costMin(i,j,k), costMin(i,k+1,k) + costMin(k+1,j,k))**; 
 + 
 +<note tip> 
 +Observaţii: 
 +  * costMin(i,j,0) = costMuchie(i,j) (dacă există muchie de la i la j) sau infinit(altfel); 
 +  * costMin(i,j,N) = drumul de cost minim de la i la j. 
 + 
 +</note> 
 + 
 +Paşii algoritmului Floyd-Warshall: 
 + 
 +<code> 
 +1. Declarăm matricile: 
 + dist, matrice N x N şi o iniţializăm dist[i][j] = infinit, pentru orice i şi j 
 + next, matrice N x N în care vom salva prima muchie din drumul i-j de cost minim 
 + //next este necesar doar în cazul în care ne interesează muchiile folosite 
 +//pasul k = 0 
 +2. Pentru fiecare nod v 
 + 1. dist[v][v] = 0; 
 +3. Pentru fiecare pereche (u,v) a.i. există muchie de la u la v 
 + 1. dist[u][v] = costMuchie(u,v); 
 + 2. next[u][v] = v; //pentru urmărirea muchiilor ce compun drumul 
 +//am terminat pasul k = 0 
 +4. Pentru fiecare pas k (de la 1 la N) 
 + 1. Pentru fiecare nod i (de la 1 la N) 
 + 1. Pentru fiecare nod j (de la 1 la N) 
 + 1. calculăm costul nou = dist[i][k] + dist[k][j]; 
 + 2. comparăm costul nou cu costul vechi = dist[i][j] 
 + dacă e mai bun costul nou: 
 + 1. dist[i][j] = costul nou; 
 + 2. next[i][j] = next[i][k]; //pentru urmărirea muchiilor 
 + altfel păstrăm costul vechi 
 +</code> 
 + 
 +Pentru obţinerea nodurilor ce formează drumul de cost minim de la u la v, putem folosi următoarea secvenţă: 
 +<code> 
 +funcţie afişareDrum(u, v) 
 +1. dacă u == v 
 + 1. print(v); STOP 
 +2. print(u); 
 +3. afişareDrum(next[u][v], v); 
 +</code> 
 + 
 +===== 5. Observaţii finale ====== 
 +<note tip> 
 +Cei 3 algoritmi funcţionează bine atunci când **toate muchiile** au **cost pozitiv**. 
 +</note> 
 +<note warning> 
 +Drumul de cost minim **NU** este bine definit atunci când există **cicluri cu cost negativ**. 
 +  * putem ajunge de la un nod la acelaşi nod, folosind acelaşi ciclu de oricâte ori => costMinim = -infinit); 
 +  * în aceste cazuri, nu putem pune problema găsirii drumurilor de cost minim. 
 + 
 +</note> 
 +<note warning> 
 +O **muchie cu cost negativ** este echivalentă cu **2 arce cu cost negativ** şi implică existenţa unui ciclu cu cost negativ. 
 +</note> 
 +<note tip> 
 +Algoritmii **Bellman-Ford** şi **Floyd-Warshall** pot detecta dacă un graf conţine **cicluri cu cost negativ**: 
 +  * Bellman-Ford: executăm încă o evaluare a tuturor muchiilor la final. Dacă găsim cel puţin o estimare nouă mai buna, graful conţine măcar un ciclu cu cost negativ. 
 +  * Floyd-Warshall: putem verifica, la fiecare pas, dacă avem o valoare negativă pe diagonala matricei dist. Dacă găsim cel puţin o valoare negativă, graful conţine măcar un ciclu cu cost negativ. 
 + 
 +</note> 
 +<note important> 
 +Pentru algoritmul lui **Dijkstra**, "garanţia" se pierde când există chiar şi **un singur arc cu cost negativ**. 
 +  * Algoritmul se bazează pe ideea că, odată explorat un nod, drumul de cost minim până la acel nod a fost găsit, ceea ce **NU** este mereu corect în prezenţa unui arc cu cost negativ. 
 +  * De aceea, algoritmul poate produce răspunsuri **greşite** în acest caz. 
 + 
 +</note> 
 +<note tip>În schimb, algoritmii **Bellman-Ford** şi **Floyd-Warshall** funcţionează pe grafurile cu **arce cu cost negativ**, atâta timp cât drumurile de cost minim sunt bine definite(**fără cicluri cu cost negativ**). 
 +</note> 
 + 
 +===== 6. Exerciții laborator ====== 
 + 
 +  - Daţi un exemplu de graf orientat cu un singur arc cu cost negativ pentru care algoritmul lui Dijkstra dă rezultate greşite. 
 +  - Cum putem folosi algoritmul lui Dijkstra sau algoritmul Bellman-Ford pentru a obţine aceleaşi rezultate ca algoritmul Floyd-Warshall(drumul de cost minim pentru toate perechile de noduri)? Ne limităm la grafurile pe care merg toţi cei trei algoritmi. 
 +  - Implementaţi unul din algoritmi pentru a calcula drumurile de cost minim de la un nod sursă la toate celelalte noduri într-un graf cu toate muchiile/arcele cu cost pozitiv. 
 +  - Verificaţi dacă un graf conţine cicluri negative. 
 + 
 +====6.1. Exerciții - schelet de laborator==== 
 +Pentru acest laborator puteți descărca scheletul de cod de [[http://elf.cs.pub.ro/sda-ab/wiki/_media/laboratoare/lab7_drumuri_minime-skel.zip|aici]]. Descărcați arhiva și dezarhivați-o.  
 + 
 +=== Linux=== 
 +Puteti folosi utilitarul ''%%wget%%'' pentru descarcare si utilitarul ''%%unzip%%'' pentru dezarhivare. 
 + 
 +  * ''%%wget http://elf.cs.pub.ro/sda-ab/wiki/_media/laboratoare/lab7_drumuri_minime-skel.zip%%'' 
 +  * ''%%unzip lab7_drumuri_minime-skel.zip%%'' 
 + 
 +Pentru compilare folositi comanda ''%%make%%''. Pentru rulare puteti folosi comanda ''%%make run%%'' sau ''%%./graph%%''
 + 
 + 
 +===Extra=== 
 + 
 +  - Se dă un graf pentru care se cunoaşte drumul de cost minim de la un nod sursă(S) la un nod destinaţie(D). Dacă adăugăm 100 la costul fiecărei muchii, se modifică drumul de cost minim? (dacă da, daţi un exemplu; dacă nu, explicaţi de ce). 
 +  - Aceeaşi întrebare dacă înmulţim fiecare cost cu 100. 
 +  - Cum găsiţi(mai rapid decât cu cei 3 algoritmi prezentaţi) drumul de cost minim de la S la D într-un graf în care toate muchiile au acelaşi cost(1)? Cum adaptaţi soluţia în cazul în care toate muchiile au costul 1 sau 2? 
 +  - Daţi un exemplu în care folosirea algoritmului lui Dijkstra(pentru a obţine drumul de cost minim pentru toate perechile de noduri) ar fi mai rapidă decât algoritmul Floyd-Warshall. 
  
-===== 6.Referințe ======+===== 7.Referințe ======
     - [[https://en.wikipedia.org/wiki/Dijkstra's_algorithm|Algoritmul lui Dijkstra]]     - [[https://en.wikipedia.org/wiki/Dijkstra's_algorithm|Algoritmul lui Dijkstra]]
     - [[https://en.wikipedia.org/wiki/Bellman%E2%80%93Ford_algorithm|Algoritmul lui Bellman Ford]]     - [[https://en.wikipedia.org/wiki/Bellman%E2%80%93Ford_algorithm|Algoritmul lui Bellman Ford]]
laboratoare/laborator-07.txt · Ultima modificare: 2018/02/25 22:45 de către mihai.iacov