Hodnocení předmětu IB013 Logické programování Hana Rudová jaro 2013 Základní informace Přednáška: účast není povinná, nicméně ... Cvičení: účast povinná ■ individuální doplňující příklady za zmeškaná cvičení ■ nelze při vysoké neúčasti na cvičení ■ skupina 01, sudá středa, první cvičení 6.března, Hana Rudová ■ skupina 02, lichá středa, první cvičení 27.února, Adriana Strejčkova Předpoklad: znalost základů výrokové a predikátové logiky, např. z IB101 Web předmětu: interaktivní osnova v ISu ■ průsvitky dostupné postupně v průběhu semestru ■ sbírka příkladů včetně řešení (zveřejněna cca do 8.3.) ■ harmonogram výuky, předběžný obsah výuky během semestru ■ elektronicky dostupné materiály ■ informace o zápočtových projektech Hana Rudová, Logické programování 1,15. května 2013 3 Organizace předmětu ■ Průběžná písemná práce: až 30 bodů (základy programování v Prologu) ■ pro každého jediný termín: 26.března ■ alternativní termín pouze v případech závažných důvodů pro neúčast ■ vzor písemky na webu předmětu ■ Závěrečná písemná práce: až 1 50 bodů ■ vzor písemky na webu předmětu ■ opravný termín možnýjako ústní zkouška ■ Zápočtový projekt: celkem až 40 bodů ■ Hodnocení: součet bodů za projekt a za obě písemky ■ známka A za cca 1 75 bodů, známka F za cca 11 0 bodů ■ známka bude zapsána pouze těm, kteří dostanou zápočet za projekt ■ Ukončení předmětu zápočtem: zápočet udělen za zápočtový projekt Hana Rudová, Logickě programování 1,15. května 2013 2 Organizace předmětu Rámcový obsah předmětu Obsah přednášky ■ základy programování v jazyce Prolog ■ teorie logického programováni ■ logické programování s omezujícími podmínkami ■ DC gramatiky Obsah cvičení ■ zaměřeno na praktické aspekty, u počítačů ■ programování v Prologu ■ logické programování ■ logické programování s omezujícími podmínkami Hana Rudová, Logickě programování 1,15. května 2013 4 Organizace předmětu Literatura Bratko, I. Prolog Programming for Artificial Intelligence. Addison-Wesley, 2001. ■ prezenčně v knihovně Clocksin, W. F. - Mellish, Ch. S. Programming in Prolog. Springer, 1 994. Sterling, L. - Shapiro, E. Y. The art of Prolog : advanced programming techniques. MIT Press, 1 987. Nerode, A. - Shore, R. A. Logic for applications. Springer-Verlag, 1 993. ■ prezenčně v knihovně Dechter, R. Constraint Processing. Morgan Kaufmann Publishers, 2003. ■ prezenčně v knihovně - Elektronicky dostupné materiály (viz web předmětu) Hana Rudová, Logické programování 1,15. května 2013 5 Organizace předmětu Zápočtové projekty Zadání, dotazy, odevzdávání: Adriana Strejčkova Týmová práce na projektech, až 3 řešitelé ■ zápočtové projekty dostupné přes web předmětu Podrobné pokyny k zápočtovým projektům na webu předmětu ■ bodování, obsah předběžné zprávy a projektu ■ typ projektu: LP, CLP Předběžná zpráva ■ podrobné zadání ■ vjakém rozsahu chcete úlohu řešit ■ které vstupní informace bude program používat a co bude výstupem programu ■ scénáře použití programu (tj. ukázky dvojic konkrétních vstupů a výstupů) Průběžná písemná práce ■ Pro každého jediný termín 26. března ■ Alternativní termín pouze v závažných důvodech pro neúčast ■ Celkem až 30 bodů (1 50 závěrečná písemka, 40 projekt) ■ 3 příklady, 40 minut ■ Napsat zadaný predikát, porovnat chování programů ■ Obsah: první čtyři přednášky a první dvě cvičení ■ Oblasti, kterých se budou příklady zejména týkat ■ unifikace ■ seznamy ■ backtracking ■ optimalizace posledního volání ■ řez ■ aritmetika ■ Ukázka průběžné písemné práce na webu Hana Rudová, Logickě programování 1,15. května 2013 6 Organizace předmětu Časový harmonogram k projektům ■ Zveřejnění zadání (většiny) projektů: 26. února ■ Zahájení registrace řešitelů projektu: 13. března, 19:00 ■ Předběžná analýza řešeného problému: 12. dubna ■ Odevzdání projektů: 1 7. května ■ Předvádění projektů (po registraci): 23.května - 20.června Hana Rudová, Logickě programování I, 15. května 2013 7 Organizace předmětu Hana Rudová, Logickě programování I, 15. května 2013 8 Organizace předmětu Software: SICStus Prolog Doporučovaná implementace Prologu Dokumentace: http://www.fi.muni.cz/~hanka/sicstus/doc/html Komerční produkt ■ licence pro instalace na domácí počítače studentů Používané IDE pro SICStus Prolog SPIDER ■ dostupné od verze SICStus 4.1.3 ■ http://www.si cs.se/si cstus/spider ■ používá Eclipse SDK Podrobné informace dostupné přes web předmětu ■ stažení SICStus Prologu (sw + licenční klíče) ■ pokyny k instalaci (SICStus Prolog, Eclipse, Spider) Hana Rudová, Logické programování 1,15. května 2013 Organizace předmětu Úvod do Prologu SICStus IDE SPIDER 5 SICStus Debugging - My Prolog ProjecF/my.morMe.pro - Eclipse SDK File Edit SICStus Source Navigate Search Project Run Favori li ň * Bl in j « i * - O-*~ j j £j - Debug Si . Piclca Tcp-le -el Configuration Ĺi,.,;v:!-, Cer&K.ir'.or, T i J$ Prolog Target = call: suffK([3,_7551rc],_1810] = my.preolLlSlO) J Prolog Top-levíl Proces j IM-Variables £š~~\.__*b Breakpoints] nary (lists) , [f>qst suffix(?List, 'Suffix! itiue v/hen Lišt and Suffix art lištš and Suffix F-s a suffix of List. It termina =ple- íl i ;iiUšĚ'šr.iii-.i,:!C:--.f. nj!i*niie-Eiľ|:äEis|;!iišt:cn-: F'oico F'-cliít y „ uu k m - * - Hana Rudová, Logické programování 1,15. května 201 3 převzato z http://www.sics.se/sicstus/spider Organizace předmětu Prolog PROgramming in LOCic ■ část predikátové logiky prvního řádu Deklarativní programování ■ specifikační jazyk, jasná sémantika, nevhodné pro procedurální postupy ■ Co dělat namísto Jak dělat Základní mechanismy ■ unifikace, stromové datové struktury, automatický backtracking Hana Rudová, Logickě programování I, 15. května 2013 12 Úvod do Prologu Logické programování Historie ■ Rozvoj začíná po roce 1 970 ■ Robert Kowalski - teoretické základy ■ Alain Colmerauer, David Warren (Warren Abstract Machine) - implementace ■ SICStus Prolog vyvíjen od roku 1 985 ■ Logické programování s omezujícími podmínkami - od poloviny 80. let Aplikace ■ rozpoznávání řeči, telekomunikace, biotechnologie, logistika, plánování, data mining, business rules, ... ■ SICStus Prolog—the first 2 5 years, Mats Carlsson, Per Mildner. Theory and Practice of Logic Programming, 12 (1-2): 35-66, 2012. http://arxiv.org/abs/1011.5640. Hana Rudová, Logické programování I, 15. května 2013 13 Úvod do Prologu Program = fakta + pravidla ■ (Prologovský) program je seznam programových klauzulí ■ programové klauzule: fakt, pravidlo ■ Fakt: deklaruje vždy pravdivé věci ■ clovekC novak, 18, student ). ■ Pravidlo: deklaruje věci, jejichž pravdivost závisí na daných podmínkách ■ studujeC X ) :- clovekC X, _Vek, student ). ■ alternativní (obousměrný) význam pravidel pro každé X, pro každé X, X studuje, jestliže X je student, potom Xje student X studuje ■ pracujeC X ) :- clovekC X, _Vek, CoDela ), praceC CoDela ). ■ Predikát: seznam pravidel a faktů se stejným funktorem a aritou ■ značíme: clovek/3, student/l; analogie procedury v procedurálních jazycích, Hana Rudová, Logické programování I, 15. května 2013 14 Úvod do Prologu Komentáře k syntaxi ■ Klauzule ukončeny tečkou ■ Základní příklady argumentů ■ konstanty: (tomas , anna) ... začínají malým písmenem ■ proměnné ■X, Y ... začínají velkým písmenem ■ _, _A, _B ... začínají podtržítkem (nezajímá nás vracená hodnota) ■ Psaní komentářů clovekC novak, 18, student ). % komentář na konci řádku clovekC novotny, 30, učitel ). /* komentář '■'•'/ Hana Rudová, Logické programování 1,15. května 2013 15 Úvod do Prologu Dotaz ■ Dotaz: uživatel se ptá programu, zda jsou věci pravdivé ?- studujeC novak ). % yes splnitelný dotaz ?- studujeC novotny ). % no nesplnitelný dotaz ■ Odpověď na dotaz ■ positivní - dotaz je splnitelný a uspěl ■ negativní - dotaz je nesplnitelný a neuspěl ■ Proměnné jsou během výpočtu instanciovány (= nahrazeny objekty) ■ ?- clovekC novak, 18, Prace ). Prace = student ■ výsledkem dotazuje instanciace proměnných v dotazu ■ dosud nenainstanciovaná proměnná: volná proměnná ■ Prolog umí generovat více odpovědí, pokud existují ?- clovekC novak, Vek, Prace ). % všechna řešení přes ";" Hana Rudová, Logické programování I, 15. května 2013 16 Úvod do Prologu Klauzule = fakt, pravidlo, dotaz Rekurzivní pravidla Klauzule se skládá z hlavy a těla Tělo je seznam cílů oddělených čárkami, čárka = konjunkce Fakt: pouze hlava, prázdné tělo ■ rodic( pavla, robert ). Pravidlo: hlava i tělo ■ upracovany_c"lovekC X ) :- clovekC X, _Vek, Prace ), praceC Prace, tezka ). Dotaz: prázdná hlava, pouze tělo ■ ?- clovekC novak, Vek, Prace ). ?- rodicC pavla, Dite ), rodic( Dite, Vnuk ). predekC X, Z ) :- rodic( X, Z ). predekC X, Z ) predekC X, Z ) rodicC X, Y ), rodicC Y, Z ). rodicC X, Y ), predekC Y, Z ). % CD % (2) % C2") Hana Rudová, Logické programování I, 1 5. května 201 3 Úvod do Prologu Hana Rudová, Logické programování 1,15. května 201 3 Úvod do Prologu Příklad: rodokmen rodic( pavla, robert ). rodic( tomas, robert ). rodicC tomas, e "liská ). rodicC robert, anna ). rodicC robert, petr ). rodicC petr, ji rka ). Výpočet odpovědi na dotaz ?- predek(tomas,robert) rodi rodi rodi rodi rodi cC pavla, robert ) . cC tomas, robert ). c C tomas, eli ska ). cC robert, anna ). cC robert, petr ). rodicC petr, jirka ). predekC X, Z ) predekC X, Z ) predekC X, Z ) :- rodicC X, Z ). % CD rodicC X, Z ). rodicC X, Y ), predekC Y, Z ). predek(tomas,robert) dle (1) yes rodic(tomas, robert) % CD % C2') predekC X, Z ) :- rodicC X, Y ), predekC Y, Z ). % C2') Hana Rudová, Logické programování I, 1 5. května 201 3 Úvod do Prologu Hana Rudová, Logické programování 1,15. května 201 3 Uvod do Prologu Výpočet odpovědi na dotaz ?- predek(tomas, petr) Odpověď na dotaz s proměnnou predek(tomas, petr) dle (1) dle (2') no rodic( tomas, petr) rodic(tomas, Y) predek( Y, petr) rodicC tomas, robert ). rodicC tomas, eliska ). rodicC robert, petr ). Y=robert dle rodic(tomas, robert) predek( robert, petr) predekC X, Z ) predekC X, Z ) rodicC X, Z ). % (1) rodicC X, Y ), % (2') predekC Y, Z ). dle (1) rodic(robert, petr) yes Hana Rudová, Logické programování I, 1 5. května 2013 Úvod do Prologu rodicC pavla, robert ). rodicC tomas, robert ). rodicC tomas, eliska ). rodicC robert, anna ). rodicC robert, petr ). rodicC petr, jirka ). predekC X, Z ) :- rodicC X, Z ). predekC X, Z ) :- rodicC X, Y ), predekC Y, Z ). predekCpetr,Potomek) --> ??? predekCrobert,P) --> ??? Hana Rudová, Logické programování 1,15. května 201 3 1. P=anna, 2. Potomek=ji rka P=petr, 3. P=jirka Úvod do Prologu Syntaxe a význam Prologovských programů Syntaxe Prologovských programů ■ Typy objektů jsou rozpoznávány podle syntaxe ■ Atom ■ řetězce písmen, čísel, „_" začínající malým písmenem: pavel , pavel_novak, x25 ■ řetězce speciálních znaků: <-->, ====> ■ řetězce v apostrofech: 'Pavel', 'Pavel Novák' ■ Celá a reálná čísla: 0, -1056, 0.35 ■ Proměnná ■ řetězce písmen, čísel, „_" začínající velkým písmenem nebo „_" ■ anonymní proměnná: ma_diteCX) :- rodicC X, _ ). ■ hodnotu anonymní proměnné Prolog na dotaz nevrací: ?- rodicC X, _ ) ■ lexikální rozsah proměnné je pouze jedna klauzule: prvni CX,X,X). prvni CX,X,_). Hana Rudová, Logické programování 1,15. května 201 3 24 Syntaxe a význam Prologovských programů Termy ■ Term - datové objekty v Prologu: datum( 1, kveten, 2003 ) ■ funktor: datum ■ argumenty: 1, kveten, 2003 ■ arita - počet argumentů: 3 ■ Všechny strukturované objekty v Prologu jsou stromy ■ trojuhelnikC bod(4,2), bod(6,4), bod(7,l) ) ■ Hlavní funktor termu - funktor v kořenu stromu odpovídající termu ■ trojuhelnikje hlavní funktor v trojuhe"lnik( bod(4,2), bod(6,4), bod(7,l) ) Hana Rudová, Logické programování I, 1 5. května 201 3 25 Syntaxe a význam Prologovských programů Unifikace Termy S a T jsou unifikovatelné, jestliže 1. SaTjsou konstanty a tyto konstanty jsou identické; 2. S je proměnná a T cokoliv jiného - S je instanciována na T; Tje proměnná a S cokoliv jiného - Tje instanciována na S 3. S a T jsou termy ■ S a T mají stejný funktor a aritu a ■ všechny jejich odpovídající argumenty jsou unifikovatelné ■ výsledná substituce je určena unifikací argumentů Příklady: k = k ... yes, kl = k2 ... no, A = k(2,3) ... yes, k(s,a,1(1)) = A ... yes s(sss(2),B,ss(2)) = s(sss(2),4,ss(2),s(l))... no s(sssCA),4,ss(3)) = s(sss(2),4,ss(A))... no s(sss(A),4,ss(C)) = s(sss(t(B)),4,ss(A))... A=t(B),C=t(B)... yes Hana Rudová, Logické programování I, 1 5. května 201 3 27 Syntaxe a význam Prologovských programů Unifikace ■ Termy jsou unifikovatelné, jestliže ■ jsou identické nebo ■ proměnné v obou termech mohou být instanciovány tak, že termy jsou po substituci identické ■ datumC Dl, Ml, 2003 ) = datum( 1, M2, Y2) operátor = Dl = 1, Ml = M2, Y2 = 2003 ■ Hledáme nejobecnější unifikátor (most generál unifler (MGU) ■ jiné instanciace? ... Dl = 1, Ml = 5 , Y2 = 2003 ... není MCU ■ ?- datumC Dl, Ml, 2003 ) = datum( 1, M2, Y2), Dl = Ml. ■ Test výskytu (occurs check) ?- X=f(X). X = f(f(f(f(f(f(f(f(f(f(...)))))))))) Hana Rudová, Logické programování 1,15. května 201 3 26 Syntaxe a význam Prologovských programů Deklarativní a procedurální význam programů ■ p :- q, r. ■ Deklarativní: Co je výstupem programu? ■ p je pravdivé, jestliže q a r jsou pravdivé ■ Z q a r plyne p => význam mají logické relace ■ Procedurální: Jak vypočítáme výstup programu? ■ p vyřešíme tak, že nejprve vyřešíme q a pak r => kromě logických relací je významné i pořadí cílů ■ výstup ■ indikátor yes/no určující, zda byly cíle splněny ■ instanciace proměnných v případě splnění cílů Hana Rudová, Logické programování 1,15. května 201 3 28 Syntaxe a význam Prologovských programů Konjunce "," vs. disjunkce ";" cílů ■ Konjunce = nutné splnění všech cílů ■ p :- q, r. ■ Disjunkce = stačí splnění libovolného cíle ■ p :- q; r. p :- q. p :- r. ■ priorita středníku je vyšší (viz ekvivalentní zápisy): p :- q, r; s, t, u. P :- Cq, r) ; (s, t, u). P :- q, r. p : - s , t, u . Hana Rudová, Logické programování 1,15. května 2013 29 Syntaxe a význam Prologovských programu Pořadí klauzulí a cílů Ca) a(l). ?- a(l). a(X) :- b(X,Y), a(Y) . b(l,l). (b) a(X) :- b(X,Y), a(Y) . % změněné pořadí klauzulí v programu vzhledem k (a) a(l). b(l,l). % nenalezení odpovědi: nekonečný cyklus (c) a(X) :- b(X,Y), c(Y). ?- a(X). b(l,l). c(2). c(l). (d) a(X) :- c(Y), b(X,Y). % změněné pořadí cílů v těle klauzule vzhledem k (c) bCl.D. c(2). c(l) . % náročnější nalezení první odpovědi než u (c) V obou případech stejný deklarativní ale odlišný procedurální význam Hana Rudová, Logické programování 1,15. května 201 3 30 Syntaxe a význam Prologovských programů Pořadí klauzulí a cílů II. (1) a(X) :- c(Y), b(X,Y). (2) b(l,l). (3) c(2). (4) c Cl). ?- a(X). a(X) dle (1) c(Y), b(X,Y) dle (3j/Ý=2 dle (4)\Y=1 Vyzkoušejte si: (1) a(X) :- b(X,X), c(X). (3) a(X) :- b(Y,X), c(X), (4) b(2,2). (5) b(2,l). (6) c(l). b(X,2) no b(X,1) dle (2) X=1 yes Hana Rudová, Logické programování I, 1 5. května 201 3 Syntaxe a význam Prologovských programů Cvičení: průběh výpočtu a b b c. d. e e g- h. i. b,c,d. e,c,f,g. g,h. h. Jak vypadá průběh výpočtu pro dotaz ?- a. Hana Rudová, Logické programování 1,15. května 201 3 Syntaxe a význam Prologovských programu Operátory, aritmetika Typy operátoru ■ Typy operátorů ■ infixové operátory: xfx, xfy, yfx př. xfx = yfx - ■ prefixové operátory: fx, fy př. fx ?- fy- ■ postfixové operátory: xf, yf ■ x a y určují prioritu argumentu ■ x reprezentuje argument, jehož priorita musí být striktně menší než u operátoru ■ y reprezentuje argument, jehož priorita je menší nebo rovna operátoru ■ a-b-c odpovídá (a-b)-c a ne a-(b-c): „-" odpovídá yfx priorita: 500 Hana Rudová, Logické programování I, 15. května 2013 35 Operátory, aritmetika Operátory ■ Infixová notace: 2*a + b*c ■ Prefixová notace:+( *(2,a), *(b,c) ) priorita +: 500, priorita*: 400 ■ prefixovou notaci lze získat predikátem display/l :- dísplay((a:-s(0),b,c)). :-(a, ,(s(0), ,(b,c))) ■ Priorita operátorů: operátor s nejvyšší prioritou je hlavní funktor ■ Uživatelsky definované operátory: zna petr zna alese. zna( petr, alese). ■ Definice operátoru: :- op( 600, xfx, zna ). priorita: i ..i200 ■ :- opC 1100, xfy, ; ). nestrukturované objekty: 0 :- op( 1000, xfy, , ). p :- q, r; s, t. p :- (q, r) ; (s, t). ; má vyšší prioritu než , ■ :- op( 1200, xfx, :- ). :-má nejvyšší prioritu ■ Definice operátoru není spojena s datovými manipulacemi (kromě spec. případů) Hana Rudová, Logické programování 1,15. května 2013 34 Operátory, aritmetika Aritmetika ■ Předdefinované operátory + , -, *, /, '■'•"■'•'mocnina,// celočíselné dělení, mod zbytek po dělení ■ ?- X = 1 + 2. X=l + 2 = odpovídá unifikaci ■ ?- X i s 1 + 2. X= 3 „i s" je speciální předdefinovaný operátor, který vynutí evaluaci ■ porovnej: N = (1+1+1+1+1) N i s (1+1+1+1+1) ■ pravá strana musí být vyhodnotitelný výraz (bez proměnné) ■ výraz na pravé straně je nejdříve aritmeticky vyhodnocen a pak unifikován s levou stranou volání?- X i s Y + 1. způsobí chybu ■ Další speciální předdefinované operátory >, <, >=, =<, =:= aritmetická rovnost, =\= aritmetická nerovnost ■ porovnej: 1+2 =:= 2+1 1+2 = 2+1 ■ obě strany musí být vyhodnotitelný výraz: volání ?- 1 < A + 2. způsobí chybu Hana Rudová, Logické programování 1,15. května 2013 36 Operátory, aritmetika Různé typy rovností a porovnání X = Y X a Y jsou unifikovatelné X \= Y X a Y nejsou unifikovatelné, (také \+ X = Y) X == Y X a Y jsou identické porovnej: ?-A == B. ... no ?- A=B, A==B. ... B = A yes X \== Y X a Y nejsou identické porovnej: ?- A \== B. .. .yes ?- A=B, A \== B. ... A no X i s Y Y je aritmeticky vyhodnoceno a výsledek je přiřazen X X =:= Y X a Y jsou si aritmeticky rovny X =\= Y X a Y si aritmeticky nejsou rovny X < Y aritmetická hodnota X je menší než Y (=<, >, >=) X @< Y term X předchází term Y (@=<, @>, @>=) 1. porovnání termů: podle alfabetického n. aritmetického uspořádání 2. porovnání struktur: podle arity, pak hlavního funktoru a pak zleva podle argumentů ?- f( pavel, g(b)) @< f( pavel, h(a)). .. .yes Rez, negace Hana Rudová, Logické programování I, 1 5. května 201 3 Operátory, aritmetika Řez a upnutí Rez a ořezání f(X,0) f(X,2) f(X,4) X < 3, ! . 3 =< X, X < 6, ! 6 =< X. přidání operátoru řezu , , ! ' ' fCX.Y) s(X,Y) s(X,Y) s(X,Y). Y i s X + 1. Y i s X + 2. fCX.Y) s(X,Y) s(X,Y) s(X,Y), !. Y i s X + 1. Y i s X + 2. 6 ?- fCl,Y), Y>2. f(X,0) :- X < 3, !. %(1) f(X,2) :- X < 6, !. %(2) f(X,4). ?- fCl.Y). Smazání řezu v (1) a (2) změní chování programu Upnutí: po splnění podcílů před řezem se už další klauzule neuvažují ?- fCl,Z). Z = 2 ? ; Z = 3 ? ; no ?- fCl,Z). Z = 2 ? ; no Ořezání: po splnění podcílů před řezem se už neuvažuje další možné splnění těchto podcílů Smazání řezu změní chování programu Hana Rudová, Logické programování 1,15. května 2013 Hana Rudová, Logické programování 1,15. května 201 3 Chování operátoru řezu Předpokládejme, že klauzule H :- TI, T2.....Tm, !, ...Tn.je aktivována voláním cíle C, kterýje unifikovatelný s H. G=h(X,Y) V momentě, kdy je nalezen řez, existuje řešení cílů TI, ..., Tm X=l, Y=l Ořezání: při provádění řezu se už další možné splnění cílů TI, ..., Tm nehledá a všechny ostatní alternativy jsou odstraněny Y=2 Upnutí: dále už nevyvolávám další klauzule, jejichž hlava je také unifikovatelná s G X=2 ?- hCX.Y). h(l,Y) :- tl(Y), ! h(2,Y) :- a. tl(l) :- b. tl(2) :- c. h(X,Y) X=l / \ X=2 tl(Y) a (vynechej: upnutí) Y=l / \ Y=2 b c (vynechej: ořezání) / Hana Rudová, Logické programování I, 1 5. května 201 3 Řez, negace Rez: návrat na rodiče ?- a(X). Cl) aCX) :- hCX.Y). C2) aCX) :- d. C3) hCl.Y) :- tlCY), !, eCX). C4) hC2,Y) :- a. C5) tlCD :- b. C6) tlC2) :- c. C7) b :- c. C8) b :- d. C9) d. CIO) eCD . CH) eC2). a(x) h(X,Y) X/l/ tl(Y),!,e(X') upnutí Y/l / X b,!,e(X') ořezání c,!,e(X') d,!,e(X') I I no !,e(X') I ! e(X') X'/l/ \X72 □ □ d I □ 1 Po zpracování klauzule s řezem se vracím až na rodiče této klauzule, tj. a(X) Hana Rudová, Logické programování I, 15. května 2013 42 Řez, negace c(X) c(X) PCX). v(X). ?- c(X). X = 1 ? X = 2 ? X = 2 ? no P(D. P(2). ?- c(2). true ? ; %p(2) true ? ; %v(2) no %p(D %P(2) %v(2) Rez: příklad cl(X) :- p(X), ! cl(X) :- v(X). v(2). ?- cl(2). true ? ; %p(2) no ?- cl(X). X = 1 ? ; %p(l) no Rez: cvičení 1. Porovnejte chování uvedených programů pro zadané dotazy. a(X,X) :- b(X). a(X,Y) :- Y is X+l. b(X) :- X > 10. a(X,X) :- b(X),!. a(X,Y) :- Y is X+l. b(X) :- X > 10. a(X,X) :- b(X),c. a(X,Y) :- Y is X+l. b(X) :- X > 10. c : - ! . ?- a(X,Y). ?- a(l,Y). ?- a(ll,Y). 2. Napište predikát pro výpočet maxima max( X, Y, Max ) Hana Rudová, Logické programování I, 1 5. května 201 3 Řez, negace Hana Rudová, Logické programování I, 15. května 2013 Typy řezu Negace jako neúspěch Zlepšení efektivity programu: určíme, které alternativy nemá smysl zkoušet Poznámka: na vstupu pro X očekávám číslo Zelený řez: odstraní pouze neúspěšná odvození ■ f(X,l) :- X >= 0, !. f(X,-l) :- X < 0. bez řezu zkouším pro nezáporná čísla 2. klauzuli Modrý řez: odstraní redundantní řešení ■ f(X,l) :- X >= 0, !. f(0,l). f(X,-l) :- X < 0. bez řezu vrací f(0,l) 2x Červený řez: odstraní úspěšná řešení ■ f(X,l) :- X >= 0, !. f(_X,-l). bez řezu uspěje 2. klauzule pro nezáporná čísla Speciální cíl pro nepravdu (neúspěch) fail a pravdu true Xa Y nejsou unifikovatelné: different(X, Y) different( X, Y ) :- X = Y, !, fail. different( _X, _Y ). Xje muž: muz(X) muz( X ) :- zena( X ), !, fail. muz( _X ). Hana Rudová, Logické programování I, 15. května 2013 45 Rez, negace Negace jako neúspěch: operátor \+ different(X,Y) :-X=Y, !, fail. muz(X) :-zena(X), !, fail. different(_X,_Y). muz(_X). Unární operátor \+ P ■ jestliže P uspěje, potom \+ P neuspěje \+(P) :- P, !, fail. ■ v opačném případě \+ P uspěje \+(_) ■ different( X, Y ) :- \+ X=Y. muz( X ) :- \+ zena( X ). Pozor: takto definovaná negace \+P vyžaduje konečné odvození P Hana Rudová, Logické programováni 1,15. května 2013 47 Hana Rudová, Logické programování 1,15. května 201 3 Negace a proměnné \+(P) :- P, !, fail. % (I) \+(_). % (ID dobre( citroen ). dobre( bmw ). drahe( bmw ). rozumne( Auto ) :-\+ drahe( Auto ) . % (D % (2) % (3) % (4) ?- dobre( X ), rozumne( X ). dobre(X),rozumne(X) dle (1), X/citroen rozumne(citroen) dle (4) \+ drahe(citroen) \ dle (II) dle (I) drahe(citroen),!, fail □ yes no Hana Rudová, Logické programování I, 1 5. května 201 3 Negace a proměnny \+(P) :- P, !, fail. % (I) \+(_). % (ID dobre( citroen ). % (1) dobre( bmw ). % (2) drahe( bmw ). % (3) rozumne( Auto ) :- % (4) \+ drahe( Auto ). ?- rozumne( X ), dobre( X ). mne(X), dobre(X) dle (4) \+drahe(X),dobre(X) dle (I) drahe(X),!,fail,dobre(X) dle (3), X/bmw !, fail, dobre(bmw) fail,dobre(bmw) no Bezpečný cíl ?- \+ drahe( citroen ). yes ?- \+ drahe( X ). no ?- rozumne( citroen ). yes ?- rozumne( X ). no \+ P je bezpečný: proměnné P jsou v okamžiku volání P instanciovány ■ negaci používáme pouze pro bezpečný cíl P Hana Rudová, Logické programování 1,15. května 2013 49 Řez, negace Chování negace ■ ?- \+ drahe( citroen ). yes ?- \+ drahe( X ). no ■ ?- drahe( X ). PTÁME SE: existuje X takové, že drahe( X ) platí? ■ ?- \+ drahe( X ). \+ drahe( X ) :- drahe(X),!,fail. \+ drahe( X ). z definice \+ plyne: není dokazatelné, že existuje X takové, že drahe( X ) platí tj. pro všechna X platí \+ drahe( X ) ■ TEDY: pro cíle s negací neplatí existuje X takové, že \+ drahe( X ) => negace jako neúspěch není ekvivalentní negaci v matematické logice ■ Negace jako neúspěch používá předpoklad uzavřeného světa pravdivé je pouze to, co je dokazatelné Hana Rudová, Logické programování I, 15. května 2013 51 Řez, negace Hana Rudová, Logické programování I, 15. května 2013 50 Řez, negace Predikáty na řízení běhu programu I. ■ řez „!" ■ fai T cíl, který vždy neuspěje true: cíl, který vždy uspěje ■ \+P: negace jako neúspěch \+ P :- P, !, fail; true. ■ once(P): vrátí pouze jedno řešení cíle P once(P) :- P, !. ■ Vyjádření podmínky: P -> Q ; R ■ jestliže platí P tak Q (P -> Q ; R) : - P, ! , Q. ■ v opačném prípade R (P -> Q ; R) :- R. ■ príklad: min(X,Y,Z) : - X =< Y -> Z = X ; Z = Y. ■ P -> Q ■ odpovídá: (P -> Q; fail) ■ príklad: zaporne(X) :- number(X) -> X < 0. Hana Rudová, Logické programování I, 15. května 2013 52 Řez, negace Predikáty na řízení běhu programu II. call (P): zavolá cíl P a uspěje, pokud uspěje P nekonečná posloupnost backtrackovacích voleb: repeat repeat. repeat :- repeat. klasické použití: generuj akci X, proveď ji a otestuj, zda neskončit Hlava :- ... uloz_stav( StaryStav ), repeat, generuje X ), % deterministické: generuj, prováděj, testuj prováděj( X ), testuj( X ), i obnov_stav( StaryStav ) , Hana Rudová, Logické programování I, 15. května 2013 53 Řez, negao Reprezentace seznamu Seznam: [a, b, c], prázdný seznam [] Hlava (libovolný objekt), tělo (seznam): .(Hlava, Telo) ■ všechny strukturované objekty stromy - i seznamy ■ funktordva argumenty ■ .(a, .(b, .(c, []))) = [a, b, c] ■ notace: [ Hlava | Telo ] = [a | Tel o] Telo je v [a | Tel o] seznam, tedy píšeme [ a, b, c ] = [ a | [ b, c ] ] Lze psát i: [a, b | Tel o] ■ před " | "je libovolný počet prvků seznamu , za " | "je seznam zbývajících prvků ■ [a.b.c] = [a|[b,c]] = [a,b|[c]] = [a,b,c|[]] ■ pozor: [ [a,b] | [c] ] ŕ [ a,b | [c] ] Seznam jako neúplná datová struktura: [a,b,c|T] ■ Seznam = [a,b,c|T], T = [d,e|S], Seznam = [a,b,c,d,e|S] Hana Rudová, Logické programování I, 1 5. května 201 3 55 Seznamy Prvek seznamu 1 member( X, 5 ) 1 platí: member( b, [a,b,c] ). ■ neplatí: member( b, [[a,b]|[c]] ). 1 Xje prvek seznamu S, když ■ Xje hlava seznamu S nebo member( X, [ X | _ ] ). %(1) ■ Xje prvek těla seznamu S member( X, [ _ I Tel o ] ) :- memberC X, Telo ). %(2) ■ Příklady použití: ■ memberCl,[2,1,3]). ■ member(X,[1,2,3]). Hana Rudová, Logické programování 1,15. května 201 3 56 dle(1). □ yes dle (1), □ yes member(1,[2,1,3,1,4]) dle (2) member(1,[1,3,1,4]) | dle (2) member(1,[3,1,4]) dle (2) member(1,[1,4]) dle (2) member(1 ,[4]) | dle (2) member(1,[ ]) dle (2) no Spojení seznamů Cvičení: append/2 appencK LI, L2, L3 ) Platí: append( [a,b], [c,d], [a,b,c,d] ) Neplatí: append( [b,a], [c,d], [a,b,c,d] ), append( [a,[b]], [c,d], [a,b,c,d] ) Definice: ■ pokud je 1. argument prázdný seznam, pak 2. a 3. argument jsou stejné seznamy: append( [], S, S ). ■ pokud je 1. argument neprázdný seznam, pak má 3. argument stejnou hlavu jako 1. append( [X|S1], S2, [X|S3] ) :- append( SI, S2, S3). X S1 S2 S3 Hana Rudová, Logické programování I, 1 5. května 201 3 appendC [], S, S ). % CD appendC [X|S1], S2, [X|S3] ) :- appendC SI, S2, S3). % (2) :- append([l,2],[3,4],A). I (2) I A=[1|B] I :- append([2],[3,4],B) . I (2) I B=[2|C] => A=[1,2|C] I :- append([],[3,4],Q. I CD | C=[3,4] => A=[l,2,3,4], I yes Hana Rudová, Logické programování 1,15. května 201 3 58 Optimalizace posledního volání Last Call Optimization (LCO) Implementační technika snižující nároky na paměť Mnoho vnořených rekurzivních volání je náročné na paměť Použití LCO umožňuje vnořenou rekurzi s konstantními pamětovými nároky Typický příklad, kdyje možné použití LCO: ■ procedura musí mít pouze jedno rekurzivní volání: v posledním cíli poslední klauzule ■ cíle předcházející tomuto rekurzivnímu volání musí být deterministické ■ p( ...):- ... % žádné rekurzivní volám' v těle klauzule p( ...):- ... % žádné rekurzivní voláni v těle klauzule p(...) :- !, p( ... ). % řez zajišťuje determinismus Tento typ rekurze lze převést na iteraci LCO a akumulátor ■ Reformulace rekurzivní procedury, aby umožnila LCO ■ Výpočet délky seznamu length( Seznam, Délka ) lengthC [] , 0 ) . lengthC [ H | T ], Délka ) :- lengthC T, DelkaO ), Délka is 1 + DelkaO. ■ Upravená procedura, tak aby umožnila LCO: % lengthC Seznam, ZapocitanaDelka, CelkovaDelka ): % CelkovaDelka = ZapocitanaDelka + ,,počet prvků v Seznam'' lengthC Seznam, Délka ) :- lengthC Seznam, 0, Délka ). % pomocný predikát lengthC □ , Délka, Délka ). % celková délka = započítaná délka lengthC [ H | T ], A, Délka ) :- A0 is A + 1, lengthC T, A0, Délka ). ■ Přídavný argument se nazývá akumulátor Hana Rudová, Logické programování I, 15. května 2013 59 Seznamy Hana Rudová, Logické programování 1,15. května 2013 max_list s akumulátorem Výpočet největšího prvku v seznamu max_"li st (Seznam, Max) max_"list([X] , X) . max_"h'st([X|T] , Max) :-max_li s t(T,MaxT), ( MaxT >= X, !, Max = MaxT Max = X ). max_"list([H|T] ,Max) :- max_"list(T,H,Max) . max_listC[], Max, Max). max_list([H|T], CastecnyMax, Max) :-( H > CastecnyMax, !, max_list(T, H, Max ) max_1ist(T, CastecnyMax, Max) ). Hana Rudová, Logické programování 1,15. května 2013 61 Seznamy reverse/2: cvičení reverseC Seznam, OpacnySeznam ) :- % (1) reverseC Seznam, [], OpacnySeznam). reverseC [] , S, S ) . % (2) reverseC [ H I T ], A, Opacny ) :- % (3) reverseC T, [ H | A ], Opacny ). ? - reverse([l,2,3],0). reverse([l,2,3] ,0) -> CD reverse([l,2,3], [], 0) - (3) reverseC[2,3] , [1], 0) -~ (3) reverseC[3] , [2 ,1] , 0) -~ (3) reverseC[] , [3,2,1] , 0) ^ (2) yes 0=[3,2,1] Akumulátor jako seznam ■ Nalezení seznamu, ve kterém jsou prvky v opačném pořadí reverse( Seznam, OpacnySeznam ) ■ reverseC [], [] ) ■ reverseC [ H I T ], Opacny ) :-reverseC T, OpacnyT ), appendC OpacnyT, [ H ], Opacny ). ■ naivní reverse s kvadratickou složitosti ■ reverse pomocí akumulátoru s lineární složitostí *% reverseC Seznam, Akumulátor, Opacny ): % Opacny obdržíme přidáním prvků ze Seznam do Akumulátor v opačném poradi reversef Seznam, OpacnySeznam ) :- reversef Seznam, [], OpacnySeznam). reversef [], S, S ). reversef [ H | T ], A, Opacny ) :- reversef T, [ H | A ], Opacny ). % přidání H do akumulátoru ■ zpětná konstrukce seznamu (srovnej s předchozí dopřednou konstrukcí, např. append) Hana Rudová, Logické programování I, 15. května 2013 62 Seznamy Neefektivita při spojování seznamů ■ Sjednocení dvou seznamů ■ append( [], S, S ). append( [X|S1], S2, [X|S3] ) :- append( SI, S2, S3 ). ■ ?- append( [2,3], [1], S ). postupné volání cílů: append( [2,3], [1], S ) - append( [3], [1], S') - append( [], [1], S" ) ■ Vždyje nutné projít celý první seznam Hana Rudová, Logické programování I, 15. května 2013 63 Seznamy Hana Rudová, Logické programování 1,15. května 2013 Rozdílové seznamy Zapamatování konce a připojení na konec: rozdílové seznamy [a,b] = L1-L2 = [a,b|T]-T = [a,b,c|S]-[c|S] = [a,b,c]-[c] Reprezentace prázdného seznamu: L-L A1 Z1 A2 Z2 1 1 1 L1 L2 \ append( Al-Zl, Z1-Z2, A1-Z2 ). LI L2 L3 [2,3] [1] [2,3,1] [2,3|Z1]-Z1 [1|Z2]-Z2 [2,3,1|Z2]-Z2 L3 ?- append( [2,3|Z1]-Z1, [1|Z2]-Z2, S ). S = AI - Z2 = [2,3|Zl] - Z2 = [2,3| [1|Z2] ] - Z2 Zl = [1|Z2] S = [2,3,1|Z2]-Z2 Jednotková složitost, oblíbená technika ale není tak flexibilní Hana Rudová, Logické programování 1,15. května 2013 65 Vestavěné predikáty Akumulátor vs. rozdílové seznamy: reverse reverseC [], [] ) . reverseC [ H | T ], Opacny ) :- reverseC T, OpacnyT ), appendC OpacnyT, [ H ], Opacny ). kvadratická složitost reverseC Seznam, Opacny ) :- reverseOC Seznam, [], Opacny ). reverseOC [] , S, S ) . reverseOC [H | T ], A, Opacny ) :- reverseOC T, [ H | A ], Opacny ). akumulátor Oineárni) reverseC Seznam, Opacny ) :- reverseOC Seznam, Opacny-[]). reverseOC □ , S-S ). reverseOC [ H | T ], Opacny-OpacnyKonec ) :- reverseOC T, Opacny-[ H | OpacnyKonec] ). rozdílové seznamy Cl i neárni) Příklad: operace pro manipulaci s frontou ■ test na prázdnost, přidání na konec, odebrání ze začátku Hana Rudová, Logické programování 1,15. května 201 3 66 Vestavěné predikáty 1 Predikáty pro řízení běhu programu ■ fai 1 , true, ... 1 Různé typy rovností ■ unifikace, aritmetická rovnost, ... 1 Databázové operace ■ změna programu (programové databáze) zajeho běhu 1 Vstup a výstup 1 Všechna řešení programu 1 Testování typu termu ■ proměnná?, konstanta?, struktura?, ... 1 Konstrukce a dekompozice termu ■ argumenty?, funktor?, ... Hana Rudová, Logické programování 1,15. května 201 3 68 Vestavěné predikáty Databázové operace Příklad: databázové operace Databáze: specifikace množiny relací Prologovský program: programová databáze, kde jsou relace specifikovány explicitně (fakty) i implicitně (pravidly) Vestavěné predikáty pro změnu databáze během provádění programu: assert( Klauzule ) přidání Klauzule do programu asserta( Klauzule ) přidání na začátek assertz( Klauzule ) přidání na konec retract( Klauzule ) smazání klauzule unifikovatelné s Klauzule Pozor: nadměrné použití těchto operací snižuje srozumitelnost programu Hana Rudová, Logické programování 1,15. května 2013 Vestavěné predikáty Vstup a výstup uživatelsky terminal program může číst data ze vstupního proudu (input stream) program může zapisovat data do výstupního proudu (output stream) dva aktivní proudy ■ aktivní vstupní proud ■ aktivní výstupní proud uživatelský terminál - user ■ datový vstup z terminálu user chápán jako jeden ze vstupních proudů ■ datový výstup na terminál chápán jako jeden z výstupních proudů soubor 1 soubor 2 user program —ta — vstupni proudy vystupni proudy Hana Rudová, Logické programování I, 1 5. května 2013 Vestavěné predikáty Caching: odpovědi na dotazy jsou přidány do programové databáze ■ ?- solve( problem, Solution), asserta( solve( problem, Solution) ). ■ :- dynamic solve/2. % nezbytné při použití v SICStus Prologu Příklad: uloz_trojice( Seznámí, Seznam2 ) :-member( XI, Seznámí ), member( X2, Seznam2 ), spocitej_treti( XI, X2, X3 ), assertz( trojice( XI, X2, X3 ) ), f ai 1. uloz_trojice( _, _ ) :- !. Hana Rudová, Logické programování 1,15. května 201 3 Vestavěné predikáty Vstupní a výstupní proudy: vestavěné predikáty ■ změna (otevření) aktivního vstupního/výstupního proudu: see(S)/tell (S) ctěni( Soubor ) :- see( Soubor ), cteni_ze_souboru( Informace ), see( user ). ■ uzavření aktivního vstupního/výstupního proudu: seen/told ■ zjištění aktivního vstupního/výstupního proudu: seeing(S)/telling(S) ctěni( Soubor ) :- seeing( StarySoubor ), see( Soubor ), cteni_ze_souboru( Informace ), seen, see( StarySoubor ). Hana Rudová, Logické programování 1,15. května 201 3 Vestavěné predikáty Sekvenční přístup k textovým souborům Příklad čtení ze souboru čtení dalšího termu: read(Term) ■ při čtení jsou termy odděleny tečkou | ?- read(A), read( ahoj(B) ), read( [C,D] ). |: ahoj. ahoj( petre ). [ ahojC 'Petre!' ), jdeme ]. A = ahoj, B = petre, C = ahoj('Petre!'), D = jdeme ■ po dosažení konce souboru je vrácen atom end_of_file zápis dalšího termu: write(Term) ?- write( ahoj ). ?- write( 'Ahoj Petre!' ). nový řádek na výstup: nl N mezer na výstup: tab(N) čtení/zápis dalšího znaku: getO(Znak) , get(NeprazdnyZnak)/put(Znak) ■ po dosažení konce souboru je vrácena -1 Hana Rudová, Logické programování I, I 5. května 2013 Vestavěné predikáty Čtení programu ze souboru Interpretování kódu programu ■ ?- consult(program). ■ ?- consuTt('program.pl'). ■ ?- consult( [programl, 'program2.pl'] ). Kompilace kódu programu ■ ?- compile( [programl, 'program2.pl']). ■ ?- [program]. ■ ?- [user]. zadávání kódu ze vstupu ukončené CTRL+D ■ další varianty podobně jako u interpretování ■ typické zrychlení: 5 až 10 krát Hana Rudová, Logické programování 1,15. května 2013 Vestavěné predikáty process_file( Soubor ) :- seeing( StarySoubor ), see( Soubor ), repeat, read( Term ), process_term( Term ) , Term == end_of_file, % zjištěni aktivního proudu % otevřeni souboru Soubor % čtení termu Term % manipulace s termem % je konec souboru? seen, see( StarySoubor ) . % uzavření souboru % aktivace původního proudu repeat. repeat repeat. % opakování Hana Rudová, Logické programování 1,15. května 201 3 Vestavěné predikáty Všechna řešení Backtracking vrací pouze jedno řešení po druhém Všechna řešení dostupná najednou: bagof/3, setof/3, findall/3 bagof( X, P, S ): vrátí seznam S, všech objektů X takových, že P je splněno vek( petr, 7 ). vek( anna, 5 ). vek( tomas, 5 ). ?- bagofC Dite, vek( Dite, 5 ), Seznam ). Seznam = [ anna, tomas ] Volné proměnné v cíli P jsou všeobecně kvantifikovány ?- bagof( Dite, vek( Dite, Vek ), Seznam ). Vek = 7, Seznam = [ petr ]; Vek = 5, Seznam = [ anna, tomas ] Hana Rudová, Logické programování 1,15. května 201 3 Vestavěné predikáty Všechna řešení II. Existenční kvantifikátor „ Pokud neexistuje řešení bagof(X, P, S) neuspěje bagof: pokud nějaké řešení existuje několikrát, pak S obsahuje duplicity bagof, setof, findall: P je libovolný cíl vek( petr, 7 ). vek( anna, 5 ). vek( tomas, 5 ) . ?- bagofC Dite, ( vek( Dite, 5 ), Dite \= anna ), Seznam ). Seznam = [ tomas ] bagof, setof, findall: na objekty shromažďované vX nejsou žádná omezení: Xje term ?- bagofC Dite-Vek, vek( Dite, Vek ), Seznam ). Seznam = [petr-7,anna-5,tomas-5] Hana Rudová, Logické programování I, 15. května 2013 77 Vestavěné predikáty ■ Přidání existenčního kvantifikátoru ,," " => hodnota proměnné nemá význam ?- bagof( Dite, Vek"vek( Dite, Vek ), Seznam ). Seznam = [petr,anna,tomas] ■ Anonymní proměnné jsou všeobecně kvantifikovány, i když jejich hodnota není (jako vždy) vracena na výstup ?- bagof( Dite, vek( Dite, _Vek ), Seznam ). Seznam = [petr] ; Seznam = [anna,tomas] ■ Před operátorem ,," " může být i seznam ?- bagof( Vek ,[Jméno,Prijmeni]"vek( Jméno, Prijmeni, Vek ), Seznam ). Seznam = [7,5,5] Hana Rudová, Logické programování 1,15. května 2013 78 Vestavěné predikáty Všechna řešení III. setof( X, P, S ): rozdíly od bagof ■ S je uspořádaný podle @< ■ případné duplicity v S jsou eliminovány findalKX, P, S ): rozdíly od bagof ■ všechny proměnné jsou existenčně kvantifikovány ?- findallf Dite, vek( Dite, Vek ), Seznam ). svS jsou shromažďovány všechny možnosti i pro různá řešení => findall uspěje přesně jednou ■ výsledný seznam může být prázdný => pokud neexistuje řešení, uspěje a vrátí S = [] ■ ?- bagofC Dite, vek( Dite, Vek ), Seznam ). Vek = 7, Seznam = [ petr ]; Vek = 5, Seznam = [ anna, tomas ] ?- findallC Dite, vek( Dite, Vek ), Seznam ). Seznam = [petr,anna,tomas] Hana Rudová, Logické programování 1,15. května 2013 Vestavěné predikáty var(X) nonvar(X) Testování typu termu X je volná proměnná X není proměnná atom(X) X je atom (pavel, 'Pavel Novák', <-->) integer(X) X je integer float(X) X je float atomic(X) X je atom nebo číslo compound(X) X je struktura Hana Rudová, Logické programování 1,15. května 201 3 Vestavěné predikáty Určení počtu výskytu prvku v seznamu count( X, S, N ) :- count( X, S, 0, N ). count( _, [], N, N ). count( X, [X|S], NO, N) :- !, Nl is NO + 1, count( X, S, Nl, N). count( X, [_|S], NO, N) :- count( X, S, NO, N). :-? count( a, [a,b,a,a], N ) N=3 :-? count( a, [a,b,X,Y], N). N=3 count( _, [], N, N ). count( X, [Y|S], NO, N ) :- nonvar(Y), X = Y, !, Nl is NO + 1, count( X, S, Nl, N ). count( X, [_|S], NO, N ) :- count( X, S, NO, N ). Hana Rudová, Logické programování I, 1 5. května 201 3 Vestavěné predikáty Konstrukce a dekompozice atomu Atom (opakování) ■ řetězce písmen, čísel, „_" začínající malým písmenem: pavel, pavel_novak, x2, x4_34 ■ řetězce speciálních znaků:+, <->, ===> ■ řetězce v apostrofech: 'Pavel', 'Pavel Novák', 'prší', 'ano' ?- 'ano'=A. A = ano Řetězec znaků v uvozovkách ■ př. "ano", "Pavel" ?- A="Pavel". ?- A="ano". A = [80,97,118,101,108] A=[97,110,111] ■ př. použití: konstrukce a dekompozice atomu na znaky, vstup a výstup do souboru Konstrukce atomu ze znaků, rozložení atomu na znaky name( Atom, SeznamASCIIKodu ) Hana Rudová, Logické programování 1,15. května 201 3 nameC ano, [97,110,111] ) nameC ano, "ano" ) Vestavěné predikáty Konstrukce a dekompozice termu ■ Konstrukce a dekompozice termu Term =.. [ Funktor | SeznamArgumentu ] a(9,e) =.. [a,9,e] Cil =.. [ Funktor | SeznamArgumentu ], call( Cil ) atom =. . X => X = [atom] ■ Pokud chci znát pouze funktor nebo některé argumenty, pak je efektivnější: functor( Term, Funktor, Arita ) functorC a(9,e), a, 2 ) functor(atom,atom,0) functor(l,1,0) arg( N, Term, Argument ) argC 2, a(9,e), e) Rekurzivní rozklad termu ■ Term je proměnná (var/l), atom nebo číslo (atomic/1) => konec rozkladu ■ Term je seznam ([_|_]) => []... řešen výše jako atomi c procházení seznamu a rozklad každého prvku seznamu ■ Term je složený (=../2 , functor/3) => procházení seznamu argumentů a rozklad každého argumentu ■ Příklad: ground/1 uspěje, pokud v termu nejsou proměnné; jinak neuspěje ground(Term) :- atomic(Term), !. ground(Term) :- var(Term), !, fail. groundC[H|T]) :- !, ground(H), ground(T). ground(Term) :- Term =.. [ _Funktor | Argumenty ], groundC Argumenty ). ?- ground(s(2,[a(l,3),b,c],X)). ?- ground(s(2,[a(l,3),b,c])). no yes Hana Rudová, Logické programování I, 15. května 2013 83 Vestavěné predikáty Hana Rudová, Logické programování I, 15. května 2013 84 Vestavěné predikáty Príklad: dekompozice termu I. count_term( Integer, Term, N ) určí počet výskytů celého čísla v termu ■ ?- coLint_term( l, a(l, 2 , b(x, z(a, b, 1)) , Y) , N ). N=2 ■ count_term( X, T, N ) :- count_term( X, T, 0, N). count_term( X, T, NO, N ) :- integer(T), X = T, !, N i s NO + 1. count_term( _, T, N, N ) :- atomic(T), !. count_term( _, T, N, N ) :- var(T), !. count_term( X, T, NO, N ) :- T =.. [ _ | Argumenty ], count_arg( X, Argumenty, NO, N ). count_arg( _, [], N, N ). count_arg( X, [ H | T ], NO, N ) :- count_term( X, H, 0, NI), N2 i s NO + NI, count_arg( X, T, N2, N ). ■ ?- count_term( 1, [a,2,[b,c],[d,[e,f],Y]], N ). count_term( X, T, NO, N ) :- T = [_|_], !, count_arg( X, T, NO, N ). klauzuli přidáme před poslední klauzuli count_term/4 Hana Rudová, Logické programování I, 1 5. května 201 3 Vestavěné predikáty Technika a styl programování v Prologu Cvičení: dekompozice termu Napište predikát substitute( Podterm, Term, Podterml, Terml), který nahradí všechny výskyty Podterm v Term termem Podterml a výsledek vrátí v Terml Předpokládejte, že Term a Podterm jsou termy bez proměnných ?- substitute( sin(x), 2*sin(x)*f(sin(x)), t, F ). F=2*t*f(t) Hana Rudová, Logické programování 1,15. května 201 3 86 Vestavěné predikáty Technika a styl programování v Prologu Styl programování v Prologu ■ některá pravidla správného stylu ■ správný vs. špatný styl ■ komentáře Ladění Efektivita Hana Rudová, Logické programování 1,15. května 201 3 88 Technika a styl programování v Prologu Styl programování v Prologu I. ■ Cílem stylistických konvencí je ■ redukce nebezpečí programovacích chyb ■ psaní čitelných a srozumitelných programů, které se dobře ladí a modifikují ■ Některá pravidla správného stylu ■ krátké klauzule ■ krátké procedury; dlouhé procedury pouze s uniformní strukturou (tabulka) ■ klauzule se základními (hraničními) případy psát před rekurzivními klauzulemi ■ vhodnájmena procedur a proměnných ■ nepoužívat seznamy ([. . . ]) nebo závorky ({...}, (...)) pro termy pevné arity ■ vstupní argumenty psát před výstupními ■ struktura programu - jednotné konvence v rámci celého programu, např. ■ mezery, prázdné řádky, odsazení ■ klauzule stejné procedury na jednom místě; prázdné řádky mezi klauzulemi; každý cíl na zvláštním řádku Hana Rudová, Logické programování I, 1 5. května 201 3 89 Technika a styl programování v Prologu Špatný styl programování merge( SI, S2, S3 ) :- 51 =[],!, S3 = S2 ; % první seznam je prázdný 52 = [], !, S3 = SI; % druhý seznam je prázdný 51 = [X|TI], 52 = [Y|T2], ( X < Y, !, Z = X, % Z je hlava seznamu S3 merge( TI, S2, T3 ); Z = Y, merge( SI, T2, T3) ), 53 = [Z | T3 ]. Správný styl programování ■ konstrukce setříděného seznamu Seznam3 ze setříděných seznamů Seznámí, Seznam2: merge( Seznámí, Seznam2, Seznam3 ) ■ merge( [2,4,7], [1,3,4,8], [1,2,3,4,4,7,8] ) ■ merge( [], Seznam, Seznam ) :- ! . % prevence redundantních řešení merge( Seznam, [], Seznam ). merge( [X|Te"lol], [Y|Telo2], [X|Te"lo3] ) :-X < Y, !, merge( Telol, [Y|Telo2], Te"lo3 ). merge( Seznámí, [Y|Telo2], [Y|Te"lo3] ) :-merge( Seznámí, Te"lo2, Te"lo3 ). Hana Rudová, Logické programování 1,15. května 201 3 90 Technika a styl programování v Prologu Styl programování v Prologu II. ■ Středník „;" může způsobit nesrozumitelnost klauzule ■ nedávat středník na konec řádku, používat závorky ■ v některých případech: rozdělení klauzle se středníkem do více klauzulí ■ Opatrné používání operátoru řezu ■ preferovat použití zeleného řezu (nemění deklarativní sémantiku) ■ červený řez používat vjasně definovaných konstruktech negace: P, !, fail; true \+ P alternativy: Podmínka, !, Cill ; Cil2 Podmínka -> Cill ; Ci 12 ■ Opatrné používání negace „\+" ■ negace jako neúspěch: negace není ekvivalentní negaci v matematické logice ■ Pozor na assert a retract: snižuji transparentnost chování programu Hana Rudová, Logické programování I, 1 5. května 201 3 91 Technika a styl programování v Prologu Hana Rudová, Logické programování 1,15. května 201 3 92 Technika a styl programování v Prologu Dokumentace a komentáře Ladění co program dělá, jak ho používat (jaký cíl spustit a jaké jsou očekávané výsledky), příklad použití které predikáty jsou hlavní (top-level) jak jsou hlavní koncepty (objekty) reprezentovány doba výpočtu a paměťové nároky jaké jsou limitace programu zda jsou použity nějaké speciální rysy závislé na systému jaký je význam predikátů v programu, jaké jsou jejich argumenty, které jsou vstupní a které výstupní (pokud víme) ■ vstupní argumenty „+", výstupní, ■ ]menoPredikatu/Arita mergeC +Seznaml, +Seznam2, -Seznam3 ) merge/3 ■ algoritmické a implementační podrobnosti Hana Rudová, Logické programování I, 1 5. května 201 3 93 Technika a styl programování v Prologu Krabičkový (4-branový) model Vizualizace řídícího toku (backtrackingu) na úrovni predikátu ■ Call: volání cíle ■ Exit: úspěšné ukončení volání cíle ■ Fail: volání cíle neuspělo ■ Redo: jeden z následujících cílů neuspěl a systém backtrackuje, aby nalezl alternativy k předchozímu řešení Call 1 -----> + predekC X, I Z ) :- rodicC X, Z ). I 1 predekC X, Z ) :- rodicC X, Y ), <----------- ------+ predekC Y, Z ). Fail 1 + <- I Exit Redo Hana Rudová, Logické programování 1,15. května 2013 Technika a styl programování v Prologu Přepínače na trasování: trace/0, notrace/0 Trasování specifického predikátu: spy/1, nospy/1 ■ spy( merge/3 ) debug/0, nodebug/0: pro trasování pouze predikátů zadaných spy/1 Libovolná část programu může být spuštěna zadáním vhodného dotazu: trasování cíle ■ vstupní informace: jméno predikátu, hodnoty argumentů při volání ■ výstupní informace ■ při úspěchu hodnoty argumentů splňující cíl ■ při neúspěchu indikace chyby ■ nové vyvolání přes stejný cíl je volán při backtrackingu Hana Rudová, Logické programování 1,15. května 201 3 Technika a styl programování v Prologu Příklad: trasování a(X) :- nonvar(X). 1 ?- a(X). a(X) :- c(X). 1 1 Call a(_463) ? a(X) :- d(X). 2 2 Call nonvar(_463) ? c(l). 2 2 Fail nonvar(_463) ? d(2). 3 2 Call c(_463) ? 3 2 Exit c(l) ? ? 1 1 Exit a(l) ? Call I ------> + a(X) I a(X) <------+ a(X) Fai 1 I I Exi t nonvar(X).| ------> c(X). I d(X). + <------ I Redo 1 ? 1 1 Redo: a(l) ? 4 2 Call: d(_463) ? 4 2 Exit: d(2) ? 1 1 Exit: a(2) ? X = 2 ? no % trace I ?- Hana Rudová, Logické programování 1,15. května 201 3 96 Technika a styl programování v Prologu Efektivita ■ Čas výpočtu, paměťové nároky, a také časové nároky na vývoj programu ■ u Prologu můžeme častěji narazit na problémy s časem výpočtu a pamětí ■ Prologovské aplikace redukují čas na vývoj ■ vhodnost pro symbolické, nenumerické výpočty se strukturovanými objekty a relacemi mezi nimi ■ Pro zvýšení efektivity je nutno se zabývat procedurálními aspekty ■ zlepšení efektivity při prohledávání ■ odstranění zbytečného backtrackingu ■ zrušení provádění zbytečných alternativ co nejdříve ■ návrh vhodnějších datových struktur, které umožní efektivnější operace s objekty Hana Rudová, Logické programování I, 1 5. května 201 3 97 Technika a styl programování v Prologu Rezoluce v predikátové logice l.řádu Zlepšení efektivity: základní techniky ■ Optimalizace posledního volání (LCO) a akumulátory ■ Rozdílové seznamy při spojování seznamů ■ Caching: uložení vypočítaných výsledků do programové databáze ■ Indexace podle prvního argumentu ■ např. v SICStus Prologu ■ při volání predikátu s prvním nainstaniovaným argumentem se používá hašovací tabulka zpřístupňující pouze odpovídající klauzule ■ zamestnanecC Prijmeni, Krestni]meno, Odděleni, ...) ■ seznamyC [], ...) :- ... . seznamyC [H|T], ...) :- ... . ■ Determinismus: ■ rozhodnout, které klauzule mají uspět vícekrát, ověřit požadovaný determinismus Hana Rudová, Logické programování 1,15. května 201 3 98 Technika a styl programování v Prologu Rezoluce ■ rezoluční princip: z F v A, G v - formule je vždy pravdivá Hana Rudová, Logické programování I, 15. května 2013 100 Rezoluce v PL1 Formule ■ literál l ■ pozitivní literál = atomická formule p(h, ■ ■ ■ ,tn) ■ negativní literál = negace atomické formule ->p(t\, ■ ■ ■ ,tn) ■ ti, ■ ■ ■ , t„jsou termy ■ klauzule C = konečná množina literálů reprezentující jejich disjunkci ■ příklad: p (X) v q(a,f) v -ip(Y) notace: {p(X),q(a,f),-ip(Y)} ■ klauzule je pravdivá je pravdivý alespoň jeden z jejích literálů ■ prázdná klauzule se značí □ aje vždy nepravdivá (neexistuje v ní pravdivý literál) ■ formule F = množina klauzulí reprezentující jejich konjunkci ■ formule je v tzv. konjuktivní normální formě (konjunkce disjunkcí) ■ příklad: (p v q) a (-ip) a (p v -iq v r) notace: {{p,q}, {->p}, {p, -^q,r}} ■ formule je pravdivá <=> všechny klauzule jsou pravdivé ■ prázdná formule je vždy pravdivá (neexistuje klauzule, která by byla nepravdivá) ■ množinová notace: literál je prvek klauzule, klauzule je prvek formule, ... Hana Rudová, Logické programování I, 15. května 2013 101 Rezoluce v PL1 Rezoluční princip ve výrokové logice ■ Rezoluční princip = pravidlo, které umožňuje odvodit z klauzulí d u {1} a {-ii} u Cz klauzuli C\ u Cz Ciujil HluC2 QuC2 ■ Ci u C2 se nazývá rezolventou původních klauzulí ■ příklad: {p,r} {^r,s} (p v r) a (-if v í) {p,s} p v s obě klauzule (pvr)a (-ir v s) musí být pravdivé protože r nestačí k pravdivosti obou klauzulí, musí být pravdivé p (pokud je pravdivé - a zobrazením, které přiřadí konstantě c prvek T>, funkčnímu symbolu f In n-ární operaci v T> a predikátovému symbolu p In n-ární relaci. ■ příklad: F = {{f(a,b)=f(b,a)}, {f(f(a,a),b) = a}} interpretace 'h: D = Z,a:= l,b := -1,/ := " + " ■ Formule je splnitelná, existuje-li interpretace, pro kterou je pravdivá ■ formule je konjunkce klauzulí, tj. všechny klauzule musí být v dané interpretaci pravdivé ■ příklad (pokrač.): F je splnitelná (je pravdivá v íi) ■ Formule je nesplnitelná, neexistuje-li interpretace, pro kterou je pravdivá ■ tj. formule je ve všech iterpretacích nepravdivá ■ tj. neexistuje interpretace, ve které by byly všechny klauzule pravdivé ■ příklad: G = {{p(b)}, {p(a)}, {->p(«)}} je nesplnitelná ({p(a)\ a {-ip(a)} nemohou být zároveň pravdivé) Hana Rudová, Logické programování I, 15. května 2013 102 Rezoluce v PL1 Rezoluční důkaz ■ rezoluční důkaz klauzule C z formule F je konečná posloupnost Ci,..., Cn = C klauzulí taková, že Q je buď klauzule z F nebo rezolventa C,, Ct pro k,j < í. ■ příklad: rezoluční důkaz {p} z formule F = {{p,r}, {q, - -iF je nepravdivá ve všech interpretacích => F je vždy pravdivá začneme-li z klauzulí reprezentujících -iF, musíme postupným uplatňováním rezolučního principu dospět k prázdné klauzuli □ Příklad: F... -ifl v a G = -iF... a a -i a G = ->F... {{a}, {-"a}} Q = {a},C2 = {-ia} rezolventa C\ a C2 je □, tj. F je vždy pravdivá rezoluční důkaz □ z formule G se nazývá rezoluční vyvrácení formule G ■ a tedy C je nepravdivá ve všech interpretacích, tj. C je nesplnitelná Hana Rudová, Logické programování I, 15. května 2013 105 Rezoluce v PL1 Substituce co s proměnnými? vhodná substituce a unifikace ■ f(X,a,g(Y)) < l,f(h(c),a,Z) < 1, X = h(c), Z = g(Y) => f(h(c), a,g(Y)) < 1 substituce je libovolná funkce 9 zobrazující výrazy do výrazů tak, že platí ■ 9(E) = E pro libovolnou konstantu E ■ 0(f(Ei, ■ ■ ■ ,£„)) = f (9(Ei), ■ ■ ■ , e(En)) pro libovolný funkční symbol / ■ 9(p(Ei, ■ ■ ■ ,En)) = p(0(Ei), ■ ■ ■ , 9(En)) pro libovolný predik. symbol p substituce je tedy homomorfismus výrazů, který zachová vše kromě proměnných - ty lze nahradit čímkoliv substituce zapisujeme zpravidla ve tvaru seznamu [X\l^\, ■ ■ ■ ,Xn/%n] kde Xi jsou proměnné a §í substituované termy ■příklad: p(X)[X/f(a)] = p(f(a)) přejmenování proměnných: speciální náhrada proměnných proměnnými ■příklad: p(X)[X/Y] = p(Y) Hana Rudová, Logické programování 1,15. května 2013 Rezoluce v PL1 ■ strom rezolučního důkazu klauzule C z formule G je binární strom: ■ kořen je označen klauzulí C, ■ listy jsou označeny klauzulemi z G a ■ každý uzel, který není listem, ■ má bezprostředními potomky označené klauzulemi C\ a Ci ■ je označen rezolventou klauzulí C\ a Ci ■ příklad: G = {{p,r}, {q, -ir}, {-*q}, {^p}} C = D {p,r} {q,^r} {^q} {^p} strom rezolučního vyvrácení (rezoluční důkaz □ z G) příklad: {{p,r}, {q, ^r}, {^q}, {^p, t}, {^s}, {s, ^ŕ}} Hana Rudová, Logické programování 1,15. května 201 3 Rezoluce v PLI Unifikace Ztotožnění dvou literálů p, q pomocí vhodné substituce a takové, že pa = qa nazýváme unifikací a příslušnou substituci unifikátorem. Unifikátorem množiny S literálů nazýváme substituce 9 takovou, že množina S9 = {te\t e S} má jediný prvek. příklad: S = { datum( Dl, Ml, 2003 ), datum( 1, M2, Y2) ; unifikátor 9 = [Dl/l, Ml/2, M2/2, Y2/2003] se: ; datum( 1, 2, 2003 ) ; Unifikátor a množiny S nazýváme nejobecnějším unifikátorem (mgu - most generál unifier), jestliže pro libovolný unifikátor 9 existuje substituce A taková, že 9 = ctA. ■ příklad (pokrač.): nejobecnější unifikátor a = [Dl /l, Y2/2003, Ml /M2], A=[M2/2] Hana Rudová, Logické programování 1,15. května 201 3 Rezoluce v PL1 Rezoluční princip v PL1 základ: ....... , . . Ci u {1} HluC2 ■ rezoluční princip ve výrokové logice -—-—- Ci u c2 ■ substituce, unifikátor, nejobecnější unifikátor rezoluční princip v PL1 je pravidlo, které ■ připraví příležitost pro uplatnění vlastního rezolučního pravidla nalezením vhodného unifikátoru ■ provede rezoluci a získá rezolventu Ci u {A} {^B}uC2 C\pu u C2u ■ kde p je přejmenováním proměnných takové, že klauzule (C\ uA)p a {B} u C2 nemají společné proměnné ■ a je nejobecnější unifikátor klauzulí Ap a B Hana Rudová, Logické programování I, 15. května 2013 109 Rezoluce v PL1 Rezoluce v PL1 Obecný rezoluční princip v PL1 Ciu{4-Al {-.ďi,-■-,-.£„} uC2 C\pu u C2u ■ kde p je přejmenováním proměnných takové, že množiny klauzulí {Aip, ■ ■ ■ ,Amp, Cíp} a {Bi, ■ ■ ■ ,Bn, C2} nemají společné proměnné ■ ct je nejobecnější unifikátor množiny {A\p, ■ ■ ■ ,Amp,B\, ■ ■ ■ ,Bn} ■ příklad: Ai = a(X) vs. {-1B1, -iB2] = {^a(b),^a(Z)} v jednom kroku potřebuji vyrezolvovat zároveň B\ i B2 Rezoluce v PL1 ■ korektní: jestliže existuje rezoluční vyvrácení F, pak F je nesplnitelná ■ úplná: jestliže F je nesplnitelná, pak existuje rezoluční vyvrácení f Příklad: rezoluce v PL1 ■ příklad: d = {p(X,Y), q(Y)} C2 = {^q(a), s(X,W)} ■ přejmenování proměnných: p = [XIZ] Ci = {p(Z,Y), q(Y)} C2 = {^q(a), s(X,W)} ■ nejobecnější unifikátor: a = [Y/a] Ci = {p(Z,a), q(a)} C2 = {^q(a), s(X,W)} ■ rezoluční princip: C = {p(Z, a), s(X, W)} ■ vyzkoušejte si: d = {q(X), -.r(y), p(X,Y), p(f(Z)J(Z))} C2 = {n(Y), -^r(W), ->p(f(W),f(W)} Hana Rudová, Logické programování 1,15. května 2013 110 Rezoluce v Zefektivnění rezoluce ■ rezoluce je intuitivně efektivnější než axiomatické systémy ■ axiomatické systémy: který z axiomů a pravidel použít? ■ rezoluce: pouze jedno pravidlo ■ stále ale příliš mnoho možností, jak hledat důkaz v prohledávacím prostoru ■ problém SAT= {SIS je splnitelná } NP úplný, nicméně: menší prohledávací prostor vede k rychlejšímu nalezení řešení ■ strategie pro zefektivnění prohledávání => varianty rezoluční metody ■ vylepšení prohledávání ■ zastavit prohledávání cest, které nejsou slibné ■ specifikace pořadí, jak procházíme alternativními cestami Hana Rudová, Logické programování I, 1 5. května 2013 111 Rezoluce v PL1 Hana Rudová, Logické programování 1,15. května 2013 Varianty rezoluční metody Věta: Každé omezení rezoluce je korektní. ■ stále víme, že to, co jsme dokázali, platí T-rezoluce: klauzule účastnící se rezoluce nejsou tautologie ■ tautologie nepomůže ukázat, že formule je nesplnitelná úplná sémantická rezoluce: úplná zvolíme libovolnou interpretaci a pro rezoluci používáme jen takové klauzule, z nichž alespoň jedna je v této interpretaci nepravdivá ■ pokud jsou obě klauzule pravdivé, těžko odvodíme nesplnitelnost formule vstupní (input) rezoluce: neúplná alespoň jedna z klauzulí, použitá při rezoluci, je z výchozí vstupní množiny S " {{p,í?}, í-'P.q.}, {p,"■, ■ ■ ■ {Cn,Bn) taková, že C = Cn+i a ■ Co a každá Bi jsou prvky S nebo některé C,, j < i ■ každá Ci+i, i < n je rezolventa Q a Bi lineární vyvrácení S = lineární rezoluční důkaz □ z S co B o \/ C, Bj \/ C7 c Hana Rudová, Logické programování I, 1 5. května 201 3 Rezoluce a logické programování Rezoluce a logické programování Lineární rezoluce II. příklad: S = {A\, A2, A3, A4} Ai = {p,q} A2 = {p, -^q] Aj, = {^p,q} A4 = {^p, -^q] S: vstupní množina klauzulí (',: střední klauzule B{. boční klauzule C i B i {p,q}{p, -iq} \/ (pí l^p,q} \/ iq] {~r„} {jí} {-.Ti.....-.tu ■ Jí V -iTi V ■ ■ ■ V -ir„ Jí ".Ti V ■ ■ ■ V ->Tn Pravidlo: jeden pozitivní a alespoň jeden negativní literál ■ Prolog: Jí : - Ti, ■ ■ ■ , T„. Matematická logika: Jí <= Ti a ■ ■ ■ a Tn ■JÍ^T Hv-iľ Jí v -.Ti v ■ ■ ■ v -iTn Klauzule: {Jí, -.T,.....-iTn] Fakt: pouze jeden pozitivní literál ■ Prolog: Jí. Matematická logika: Jí Klauzule: {Jí} Cílová klauzule: žádný pozitivní literál ■ Prolog: : - Ti,... Tn. Matematická logika: -iTi v ■ ■ ■ v -iTn Klauzule: {-.Ti, ■ ■ ■ ->Tn} Hana Rudová, Logické programování I, 15. května 2013 117 Rezoluce a logické programování ■ Programová klauzule: právě jeden pozitivní literál (fakt nebo pravidlo) ■ Logický program: konečná množina programových klauzulí ■ Příklad: ■ logický program jako množina klauzulí: P = {Pi,P2,ft} Pi = {p], P2 = {p,^q}, P3={q] ■ logický program v prologovské notaci: V- v ■ -a- q. ■ cílová klauzule: G = {^q, -p} M Í^PJ M \/ □ {~pj (qJ f-ipl {p,->qt \/ {^q} Iq} □ Střední klauzule jsou cílové klauzule Hana Rudová, Logické programování I, 1 5. května 201 3 Rezoluce a logické programování Lineární vstupní rezoluce Vstupní rezoluce na P u {G} ■ (opakování:) alespoň jedna z klauzulí použitá při rezoluci je z výchozí vstupní množiny ■ začneme s cílovou klauzulí: Co = G ■ boční klauzule jsou vždy z P (tj. jsou to programové klauzule) (Opakování:) Lineární rezoluční důkaz C z P je posloupnost dvojic (C0,Ba), ... {Cn,Bn) taková, že C = Cn+i a ■ Co a každá J?í jsou prvky P nebo některé Cj,j < í ■ každá Cí+i, í < n je rezolventa C; a Bí Lineární vstupní (Linear ínput) rezoluce (Ll-rezoluce) C z P u {G} posloupnost dvojic (Co,Bo), ■ ■ ■ (Cn,Bn) taková, že C = Cn+\ a ■ Co = G a každá Bí jsou prvky P ■ každá Cf+i, i < n je rezolventa Q a Bí lineární rezoluce + vstupní rezoluce Hana Rudová, Logické programování 1,15. května 201 3 Rezoluce a logické programování Cíle a fakta při lineární rezoluci ■ Věta: Je-li S nesplnitelná množina Hornových klauzulí, pak S obsahuje alespoň jeden cíl a jeden fakt. ■ pokud nepoužiji cíl, mám pouze fakta (1 požit.literál) a pravidla (1 požit.literál a alespoň jeden negat. literál), při rezoluci mi stále zůstává alespoň jeden požit, literál ■ pokud nepoužiji fakt, mám pouze cíle (negat.literály) a pravidla (1 požit.literál a alespoň jeden negat. literál), v rezolventě mi stále zůstávají negativní literály ■ Věta: Existuje-li rezoluční důkaz prázdné množiny z množiny S Hornových klauzulí, pak tento rezoluční strom má v listech jedinou cílovou klauzuli. ■ pokud začnu důkaz pravidlem a faktem, pak dostanu zase pravidlo ■ pokud začnu důkaz dvěma pravidly, pak dostanu zase pravidlo ■ na dvou faktech rezolvovat nelze => dokud nepoužiji cíl pracuji stále s množinou faktů a pravidel ■ pokud použiji v důkazu cílovou klauzulí, fakta mi ubírají negat.literály, pravidla mi je přidávají, v rezolventě mám stále samé negativní literály, tj. nelze rezolvovat s dalším cílem Hana Rudová, Logické programování 1,15. května 2013 121 Rezoluce a logické programování Uspořádané klauzule (definite clauses) ■ Klauzule = množina literálů ■ Uspořádaná klauzule (definite clause) = posloupnost literálů ■ nelze volně měnit pořadí literálů ■ Rezoluční princip pro uspořádané klauzule: _{^Ao,..., ->An}_{B, -i/3q, ..., -^Bm}_ {-*A0,-lAf-i, -ifioP, ■ ■ ■, ^Bmp, -lAf+i,..., ^An}a ■ uspořádaná rezolventa: {-*Ao.....~vlí-i, -*Bop.....^Bmp, ~^4í+i.....~^An}a ■ p je přejmenování proměnných takové, že klauzule {Aq.....An\ a {B,Bq.....Bm}p nemají společné proměnné ■ Ll-rezoluce aplikovaná na (obecnou) formuli nezaručuje, že nalezeneme důkaz, i když formule platí! ■ Význam Ll-rezoluce pro Hornovy klauzule: ■ P = ÍPl.....Pn], G = {Gl.....Gm\ ■ Ll-rezolucí ukážeme nesplnitelnost P\ a ■ ■ ■ a Pn a (-iGi v ■ ■ ■ v -iGm) ■ pokud tedy předpokládáme, že program {P\.....Pn\ platí, tak musí být nepravdivá (-iGi v ■ ■ ■ v -iGm), tj. musí platit Gi a ■ ■ ■ a Gm Hana Rudová, Logické programování 1,15. května 2013 122 Rezoluce a logické programování Uspořádané klauzule II. ■ Uspořádané klauzule _{^A0,..., -^An}_{B, -i/3q, ..., -^Bm}_ {-*A0,-lAi-i, -ifioP, ■ ■ ■, -'Bmp, -*Ai+i,..., -^An}a Hornovy klauzule _: -A0,... ,An._B : -B0,... ,Bm._ : - (A0,...,Ai-i,B0p,...,Bmp,Ai+i,...,An)a. ■ Příklad: {^s(X), -.t(l), -^u(X)} {t(Z),^q(Z,X), ^r(3)} {-.í(A-),-.q(l,A),-.r(3),-.u(A-)} : -s(X), t(l),u(X). ř(Z) : -q(Z,X),r(3). : -s(X),q(l,A),r(3),u(X). p = [X/A] a=[Z/l] Hana Rudová, Logické programování 1,15. května 2013 124 Rezoluce a logické programování LD-rezoluce SLD-rezoluce LD-rezoluční vyvrácení množiny uspořádaných klauzulí p u {G} je posloupnost (Go, Co), ■ ■■, (Gn, Cn) taková, že ■ Gí,Cí jsou uspořádané klauzule ■ G = Go - Gn+1 = □ ■ Gi je uspořádaná cílová klauzule ■ Ci je přejmenování klauzule z p ■ Q neobsahuje proměnné, které jsou v Gj, j < í nebo v Q, fc < i ■ Gí+i,0 < i < n je uspořádaná rezolventa Gi a Cí LD-rezoluce: korektní a úplná Hana Rudová, Logické programování 1,15. května 2013 12 5 Rezoluce a logické programování ■ Lineární rezoluce se selekčním pravidlem = SLD-rezoluce (Selected Linear resolution for Definite clauses) ■ rezoluce ■ Selekční pravidlo ■ Lineární rezoluce ■ Definite (uspořádané) klauzule ■ vstupní rezoluce ■ Selekční pravidlo R je funkce, která každé neprázdné klauzuli C přiřazuje nějaký z jejích literálů R(C) e C ■ při rezoluci vybírám z klauzule literál určený selekčním pravidlem ■ Pokud se R neuvádí, pak se předpokládá výběr nejlevějšího literálu ■ nejlevější literál vybírá i Prolog Hana Rudová, Logické programování 1,15. května 2013 126 Rezoluce a logické programování Lineární rezoluce se selekčním pravidlem p= {{p}Ap,^q}Aq}}, G={^q,^p} {^q,^p} Iq} {^q,^p} (pI | y výběr | y výběr nejlevějšího {~p} M literálu výběr nejpravějšího í-"?? M literálu □ □ SLD-rezoluční vyvrácení p u {G} pomocí selekčního pravidla R je LD-rezoluční vyvrácení (Go, Co),..., (Gn, Cn) takové, že G = Go, Gn+\ = □ a R(Gí) je literál rezolvovaný v kroku i SLD-rezoluce - korektní, úplná Efektivita SLD-rezoluce je závislá na ■ selekčním pravidle R * způsobu výběru příslušné programové klauzule pro tvorbu rezolventy ■ v Prologu se vybírá vždy klauzule, která je v programu první Hana Rudová, Logické programování 1,15. května 2013 Rezoluce a logické programování Příklad: SLD-strom ř : -p,r. t: -s. P '■ -q,v. p : -u, w. q. s. u. (D (2) (3) (4) (5) (6) (7) fail fail Hana Rudová, Logické programování 1,15. května 201 3 Rezoluce a logické programování Strom výpočtu (SLD-strom) Příklad: SLD-strom a výsledná substituce SLD-strom je strom tvořený všemi možnými výpočetními posloupnostmi logického programu P vzhledem k cíli G kořenem stromu je cílová klauzule G v uzlech stromu jsou rezolventy (rodiče uzlu a programové klauzule) ■ číslo vybrané programové klauzule pro rezoluci je v příkladu uvedeno jako ohodnocení hrany listyjsou dvojího druhu: ■ označené prázdnou klauzulí - jedná se o úspěšné uzly (succes nodes) ■ označené neprázdnou klauzulí - jedná se o neúspěšné uzly (failure nodes) úplnost SLD-rezoluce zaručuje existenci cesty od kořene k úspěšnému uzlu pro každý možný výsledek příslušející cíli G (1) (2) (3) (4) (5) Hana Rudová, Logické programování I, 1 5. května 2013 Rezoluce a logické programování : -a(Z). a(X) : -b(X,Y),c(Y). a(X) : -c(X). b(2,3). b(l,2). c(2). Cvičení: p(B) : -q(A,E),r(E). p (A) : -q(A,A). q(a, a). q(a, b). r(b). Hana Rudová, Logické programování 1,15. května 201 3 :- a(Z). (1)/ \2) b(Z,Y), c(Y). :-c(Z). [Z/2.Y/3] (3)/ \(4) [Z/1 ,Y/2] \ (5) [Z/2] c(3). - c(2). (5) □ [Z/2] fail □ ÍZ/11 ve výsledné substituci jsou pouze proměnné z dotazu výsledné substituce jsou [Z/l] a [Z/2] nezajímá mě substituce [Y12] Rezoluce a logické programování Výsledná substituce {answer substitution) q(a). p(a,b). .-q(X), p(X,Y). X=a, Y=b q(X),p(X,Y). q(a). [X/a] :-piajy^pi^b). [Y/b] \/ □ [X/a, Y/b] Každý krok SLD-rezoluce vytváří novou unifikační substituci 0í => potenciální instanciace proměnné ve vstupní cílové klauzuli Výsledná substituce (answer substitution) 0 = 9q9i ■ ■ ■ 0n složení unifikací Význam SLD-rezolučního vyvrácení P u {G} Množina P programových klauzulí, cílová klauzule G Dokazujeme nesplnitelnost (DPA(VÍ)hGi(I) v ^G2(X) v ■ ■ ■ v -.G„(*)) kde G = {^Gi, -1G2, ■ ■ ■ , ^Gn} a X je vektor proměnných v G nesplnitelnost (1) je ekvivalentní tvrzení (2) a (3) (2) P h -G (3) Pr-(3Í)(Gi(í)A---AG„(í)) a jedná se tak o důkaz existence vhodných objektů, které na základě vlastností množiny P splňují konjunkci literálů v cílové klauzuli Důkaz nesplnitelnosti P u {G} znamená nalezení protipříkladu ten pomocí SLD-stromu konstruuje termy (odpověď) splňující konjunkci v (3) Hana Rudová, Logické programování 1,15. května 2013 Rezoluce a logické programování Hana Rudová, Logické programování 1,15. května 201 3 Rezoluce a logické programování Výpočetní strategie SLD-rezoluce v Prologu: úplnost Korektní výpočetní strategie prohledávání stromu výpočtu musí zaručit, že se každý (konečný) výsledek nalézt v konečném čase Korektní výpočetní strategie = prohledávání stromu do šířky ■ exponenciální paměťová náročnost ■ složité řídící struktury Použitelná výpočetní strategie = prohledávání stromu do hloubky ■ jednoduché řídící struktury (zásobník) ■ lineární paměťová náročnost ■ není ale úplná: nenalezne vyvrácení i když existuje ■ procházení nekonečné větve stromu výpočtu => na nekonečných stromech dojde k zacyklení ■ nedostaneme se tak najiné existující úspěšné uzly Hana Rudová, Logické programování 1,15. května 2013 Rezoluce a logické programování Test výskytu Kontrola, zda se proměnná vyskytuje v termu, kterým ji substituujeme ■ dotaz : -a(B,B). ■ logický program: a(X,f(X)). ■ by vedl k: [B/X], [X/f(X)] Unifikátorpro0(*i,...,*n) ag(f(X0,X0),f(X1,X1),...,f(Xn-1,Xn-1)) Xi=f(Xo,Xo), X2 = f(Xi,X\),..., Xn = f(Xn-i,Xn-i) X2=f(f(X0,Xo),f(X0,Xo)),... délka termu pro Xt exponenciálně narůstá => exponenciální složitost na ověření kontroly výskytu Test výskytu se při unifikaci v Prologu neprovádí Důsledek: ? - X = f(X) uspěje s X = /(/(/(/(/(/(/(/(/(/(...)))))))))) Hana Rudová, Logické programování 1,15. května 2013 Rezoluce a logické programování Prolog: prohledávání stromu do hloubky => neúplnost použité výpočetní strategie Implementace SLD-rezoluce v Prologu ■ není úplná logický program: q : -r. (1) r:-q. (2) q. (3) dotaz: : -q. q. (3) □ Hana Rudová, Logické programování 1,15. května 201 3 Rezoluce a logické programování SLD-rezoluce v Prologu: korektnost Implementace SLD-rezoluce v Prologu nepoužívá při unifikaci test výskytu => není korektní (1) t(X) : -p(X,X). P(XJ(X)). (2) t:-p(X,X). P(XJ(X)). :-t(X). X = /(/(/(/(...)))))))))) problém se projeví yes dokazovací systém nehledá unifikátor pro X a f(X) Řešení: problém typu (2) převést na problém typu (1) ? ■ každá proměnná v hlavě klauzule se objeví i v těle, aby se vynutilo hledání unifikátoru (přidáme X = X pro každou X, která se vyskytuje pouze v hlavě) ř : -p(X,X). P(XJ(X)):-X = X. ■ optimalizace v kompilátoru mohou způsobit opět odpověď „yes" Hana Rudová, Logické programování 1,15. května 201 3 Rezoluce a logické programování Řízení implementace: řez Příklad: řez ■ nemá ale žádnou deklarativní sémantiku ■ místo toho mění implementaci programu ■ p : -q,!, v. S ře^Še^fifeftfícÚy'chová jako kterýkoliv jiný literál ■ pokud uspěji => přeskočím řez a pokračuji jako by tam řez nebyl ■ pokud ale neuspěji (a tedy i při backtrackingu) a vracím se přes řez => vracím se až na rodiče : -p. a zkouším další větev => nezkouším tedy další možnosti, jak splnit p upnutí => a nezkouším ani další možnosti, jak splnit q v SLD-stromu ořezání Hana Rudová, Logické programování I, I 5. května 201 3 Rezoluce a logické programování Příklad: řez II a(X) :-b(X,Y),\,c(Y). (1) a(X):-c(X). (2) fc(2,3). (3) b(l,2). (4) c(2). (5) s(X) : -a(X). (6) s(X):-p(X). (7) p(B) : -q(A,B), r{B). (8) p(A) : -q(A,A). (9) q(a,a). (10) q(a,b). (11) r(b). (12) Hana Rudová, Logické programování I, 1 5. května 201 3 s(X). [X/2.Y/3] (3) (rez) J7) Rezoluce a logické programování ř : -p,r. t: -s. p : -q(X),\,v. p : -u, w. q(a). q(b). s. u. (D (2) (3) (4) (5) (6) (7) (8) [X/a] (i)/ \^) y X X' q(X),!,v,r. D (5) ^ !,v,r. (!) v, r. fail Hana Rudová, Logické programování 1,15. května 201 3 Rezoluce a logické programování Příklad: řez a(X) : -b(X,Y),c(Y),\. a(X) : -c(X). fc(2,3). b(l,2). c(2). s(X) : -a(X). s(X) : -p(X). p(B) : -q(A,B),r(B). p(A) : -q(A,A). q (a, a). q(a, b). r(b). (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) .- s(X). [X/2.Y/3] (3) Hana Rudová, Logické programování 1,15. května 201 3 Rezoluce a logické programování Operační sémantika Operační a deklarativní sémantika Operační sémantikou logického programu P rozumíme množinu 0{P) všech atomických formulí bez proměnných, které lze pro nějaký cíl g1 odvodit nějakým rezolučním důkazem ze vstupní množiny P u {Cj}. ^ímto výrazem jsou míněny všechny cíle, pro něž zmíněný rezoluční důkaz existuje. Deklarativní sémantika logického programu P 111 Opakování: interpretace Interpretace 1 jazyka £ je dána univerzem T> a zobrazením, které přiřadí konstantě c prvek T>, funkčnímu symbolu f/n n-ární operaci v D a predikátovému symbolu pln n-ární relaci. ■ příklad: F = {{f (a, b) = f (b,a)}, {f (f (a, a), b) = a}} interpretace 'h: D = Z, a := 1, b := -1,/ := " + " Interpretace se nazývá modelem formule, je-li v ní tato formule pravdivá ■ interpretace množiny N s obvyklými operacemi je modelem formule ( 0 + s(0) = s(0)) Hana Rudová, Logické programování 1,15. května 2013 Hana Rudová, Logické programování I, 15. května 2013 142 Sémantiky Herbrandovy interpretace 1 Omezení na obor skládající se ze symbolických výrazů tvořených z predikátových a funkčních symbolů daného jazyka ■ při zkoumání pravdivosti není nutné uvažovat modely nad všemi interpretacemi 1 Herbrandovo univerzum: množina všech termů bez proměnných, které mohou být tvořeny funkčními symboly a konstantami daného jazyka 1 Herbrandova interpretace: libovolná interpretace, která přiřazuje ■ proměnným prvky Herbrandova univerza ■ konstantám sebe samé ■ funkčním symbolům funkce, které symbolu f pro argumenty ti, ■ ■ ■ , tn přiřadí term f(ti, ■■■ ,tn) ■ predikátovým symbolům libovolnou funkci z Herbrand. univerza do pravdivostních hodnot 1 Herbrandův model množiny uzavřených formulí T: Herbrandova interpretace taková, že každá formule z P je v ní pravdivá. Hana Rudová, Logické programování 1,15. května 201 3 Specifikace Herbrandova modelu Příklad: Herbrandovy interpretace Herbrandovy interpretace mají předdefinovaný význam funktorů a konstant Pro specifikaci Herbrandovy interpretace tedy stačí zadat relace pro každý predikátový symbol Příklad: Herbrandova interpretace a Herbrandův model množiny formulí lichy(s(0)). % (1) lichyCs(sCX))) :- "lichy(X). % (2) 2i = 0 12 = {líchy(5(0))} % = {líchy (5(0)), líchy (5(5(5(0))))} 24 = {líchy(sn(0))\n e {1,3,5,7,...}} 25 = {líchy(sn(0))\n e N}} není model (1) není model (2) není model (2) Herbrandův model (1) i (2) Herbrandův model (1) i (2) Hana Rudová, Logické programování I, 15. května 2013 145 Sémantiky Deklarativní a operační sémantika Je-li S množina programových klauzulí a M libovolná množina Herbrandových modelů S, pak průnik těchto modelů je opět Herbrandův model množiny S. Důsledek: Existuje nejmenší Herbrandův model množiny S, který značíme M(S). Deklarativní sémantikou logického programu P rozumíme jeho minimální Herbrandův model M(P). Připomenutí: Operační sémantikou logického programu P rozumíme množinu O(P) všech atomických formulí bez proměnných, které lze pro nějaký cíl g1 odvodit nějakým rezolučním důkazem ze vstupní množiny P u {G.}. ^ímto výrazem jsou míněny všechny cíle, pro něž zmíněný rezoluční důkaz existuje. Pro libovolný logický program P platí M(P) = O(P) Hana Rudová, Logické programování 1,15. května 2013 rodic(a,b). rodi c(b,c). predek(X,Y) :- rodic(X,Y). predek(X,Z) :- rodic(X,Y), predek(Y,Z). 1\ = {rodičia, b),rodíc(b, c), predekia, h), predek(b, c), predekia, c)} 1z = {rodíc(a,b),rodíc(b,c), predekia, b), predekib, c), predekia, c), predekia, a)} 1\ i ?2 jsou Herbrandovy modely klauzulí Cvičení: Napište minimální Herbrandův model pro následující logický program. muz(petr). muz(pavel). zena(olga). zena(jitka). pary(X,Y) :- zena(X), muz(Y). Uveďte další model tohoto programu, který není minimální. Hana Rudová, Logické programování I, 15. května 2013 146 Séma Negace v logickém programování Negativní znalost logické programy vyjadřují pozitivní znalost negativní literály: pozice určena definicí Hornových klauzulí => nelze vyvodit negativní informaci z logického programu ■ každý predikát definuje úplnou relaci ■ negativní literál není logickým důsledkem programu relace vyjádřeny explicitně v nejmenším Herbrandově modelu ■ nad(X, Y) : -na(X, Y). na(c, b). nad(X, Y) : -na(X,Z),nad(Z, Y). na(b,a). ■ nej menší Herbrandův model: {na(b,a),na(c, b),nad(b, a),nad(c, b),nad(c,a)} ani program ani model nezahrnují negativní informaci ■ a není nad c, a není na c ■ i v realitě je negativní informace vyjádřena explicitně zřídka, např. jízdní řád Hana Rudová, Logické programování 1,15. května 2013 Negace v logickém programování Negace jako neúspěch (negation as failure) slabší verze CWA: definitivně neúspěšný (finitely failed) SLD-strom cíle : -A : -A má definitivně (konečně) neúspěšný SLD-strom (negation as failure, IMF normální cíl: cíl obsahující i negativní literály ■ : -nad(c,a), -inad(b,c). rozdíl mezi CWA a NF ■program nad(X, Y) : -nad(X,Y), cíl :--inad(b,c) ■ neexistuje odvození cíle podle NF, protože SLD-strom : -nad(b, c) je nekonečný ■ existuje odvození cíle podle CWA, protože neexistuje vyvrácení: -nad(b,c) CWA i NF jsou nekorektní: A není logickým důsledkem programu P řešení: definovat programy tak, aby jejich důsledkem byly i negativní literály zúplnění logického programu Hana Rudová, Logické programování 1,15. května 2013 Negace v logickém programování Předpoklad uzavřeného světa neexistence informace chápána jako opak: předpoklad uzavřeného světa (dosed world assumption, CWA) převzato z databází určitý vztah platí pouze když je vyvoditelný z programu. „odvozovací pravidlo" (A je (uzavřený) term): P if A ->A (CWA) pro SLD-rezoluci: P ý= nad(a, c), tedy lze podle CWA odvodit - 0): f(X\,.. .,Xm) ± giXi, ■ ■ ■, Yn) 7. pro každý f/m: f(Xu...,Xm) = f(Yu...,Ym) - Xi = Yi a ■ ■ ■ a Xm = Ym 8. pro každý term ř obsahující X jako vlastní podterm: ř ^ X X ± Y je zkrácený zápis -<(X = Y) Hana Rudová, Logické programování 1,15. května 2013 1 5 5 Negace v logickém programování na(X1,X2) : -3Y(X1 = c,X2 = b,f(Y)) v (X1 = b,X2 = a,g). ■ q/n predikátový symbol programu P na(c,b):-f(Y). na(b,a):-g. ■ Xi,..., Xn jsou „nové" proměnné, které se nevyskytují nikde v P ■ Nechť C je klauzule ve tvaru q(t\,..., tn) : -Li,..., Lm kde m > 0, ři,..., tn jsou termy a L\,... ,Lm jsou literály. Pak označme E(C) výraz 3Yi,Yu(X\ = t\,... ,Xn = tn,L\,... ,Lm) kde Yi,..., Yk jsou všechny proměnné v C. ■ Nechť de f (q/n) = {Ci,..., Cj}- Pak formuli l¥(q,P) získáme následujícím postupem: q(Xi,...,Xn) : -E(Ci) VE(C2) v ■ ■ ■ vE(Q) pro j > 0a q(X\,.. .,Xn) : - □ pro j = 0 [q/n není v programu P Hana Rudová, Logické programování 1,15. května 2013 1 54 Negace v logickém programování Korektnost a úplnost NF pravidla ■ Korektnost NF pravidla: Nechť P logický program a : -A cíl. Jestliže : -A má definitivně neúspěšný SLD-strom, pak V(-iA) je logickým důsledkem comp(P) (nebo-li comp(P) i= V(-iA)) ■ Úplnost NF pravidla: Nechť P je logický program. Jestliže comp(P) i= V(-iA), pak existuje definitivně neúspěšný SLD-strom : -A. ■ zůstává problém: není rozhodnutelné, zda daná atomická formule je logickým důsledkem daného logického programu. ■ teorém mluví pouze o existenci definitivně neúspěšného SLD-stromu ■ definitivně (konečně) neúspěšný SLD-strom sice existuje, ale nemusíme ho nalézt ■ např. v Prologu: může existovat konečné odvození, ale program přesto cyklí (Prolog nenajde definitivně neúspěšný strom) ■ Odvození pomocí NF pouze test, nelze konstruovat výslednou substituci ■ v (comp(P) i= V(-vl)) je A všeob. kvantifikováno, v V(-vi) nejsou volné proměnné Hana Rudová, Logické programování 1,15. května 2013 1 56 Negace v logickém programování Normální a stratifikované programy ■ normální program: obsahuje negativní literály v pravidlech ■ problém: existence zúplnění, která nemají žádný model ■ p : -->p. zúplnění: p «- -ip ■ rozdělení programu na vrstvy ■ vynucují použití negace relace pouze tehdy pokud je relace úplně definovaná ■ a. a. a: — ^b,a. a: —^b,a. b. b : --ia. stratifi kovaný není stratifi kovaný ■ normální program P je stratifikovaný: množina predikátových symbolů programu lze rozdělit do disjunktních množin So,...,Sm (Si = stratům) ■ p(...) : - ..., q(...), ... e P, p e Sk => q e S0 u ... u Sk ' p(...) : - (...), ... e P, p e Sk => q e S0 u ... u Sk-i Hana Rudová, Logické programování I, 15. května 2013 1 57 Negace v logickém programování Stratifi kované programy II ■ program je m-stratifikovaný m je nejmenší index takový, že So u ... u Sm je množina všech predikátových symbolů z P ■ Věta: Zúplnění každého stratifikovaného programu má Herbrandův model. ■ p :-~ úspěšné odvození Hana Rudová, Logické programování I, 15. května 2013 1 59 Negace v logickém programování SLDNF rezoluce: neúspěšné odvození ■ NF pravidlo: : _ C. má konečně neúspěšný SLD-strom ■ Pokud máme negativní podcíl - neúspěšné odvození Hana Rudová, Logické programování 1,15. května 2013 160 Negace v logickém programování SLDNF rezoluce: uvázlé odvození Cvičení: SLDNF odvození NF pravidlo: C. má konečně neúspěšný SLD-strom ■ Pokud máme negativní podcíl ->C v dotazu G, pak hledáme důkaz pro C ■ Pokud existuje vyvrácení C s neprázdnou substitucí (strom pro C je konečně úspěšný), pak je odvození G (i -,(Gi;Ci>,(G2;C2>,... kde v každém kroku m + l(m > 0), R vybírá pozitivní literál v Gm a dospívá k Gm+\ obvyklým způsobem. ■ konečné SLD+-odvození může být: 1. úspěšné: Gi = □ 2. neúspěšné 3. blokované: Gj je negativní (např. -,..., (Gi, e>,(: — Li,..., Lm_i, Lm+i,..., Ln) je (úspěšné) SLDNF odvození cíle Go ■ e označuje prázdnou cílovou substituci Hana Rudová, Logické programování 1,15. května 2013 165 Negace v logickém programování Korektnost a úplnost SLDNF odvození korektnost SLDNF-odvození: P normální program, : -G normální cíl a R je selekční pravidlo: je-li 9 cílová substituce SLDNF-odvození cíle : -G, pak G9 je logickým důsledkem comp(P) ■ implementace SLDNF v Prologu není korektní ■ Prolog neřeší uvázlé SLDNF-odvození (neprázdná substituce) ■ použití bezpečných cílů (negace neobsahuje proměnné) úplnost SLDNF-odvození: SLDNF-odvození není úplné ■ pokud existuje konečný neúspěšný strom : -A, pak -iA platí ale místo toho se odvozování: -A může zacyklit, tj. SLDNF rezoluce -iA neodvodí => -iA tedy sice platí, ale SLDNF rezoluce ho nedokáže odvodit Typy SLDNF odvození Konečné SLDNF-odvození může být: 1. úspěšné: Gi = □ 2. neúspěšné 3. uvázlé (flounder): Gj je negativní (-*A) a : -A je úspěšné s neprázdnou cílovou substitucí 4. blokované: Gi je negativní (-*A) a : -A nemá konečnou úroveň. Hana Rudová, Logické programování 1,15. května 201 3 1 66 Negace v logickém programování Logické programování s omezujícími podmínkami Constraint Logic Programming: CLP Hana Rudová, Logické programování I, 15. května 2013 167 Negace v logickém programování CP: elektronické materiály Dechter, R. Constraint Processing. Morgan Kaufmann Publishers, 2003. ■ http://www.ics.uci.edu/~dechter/books/materials.html průsvitky ke knize Barták R. Přednáška Omezující podmínky na MFF UK, Praha. ■ http://kti.ms.mff.cuni.cz/~bartak/podmi nky/i ndex.html SICStus Prolog User's Manual. Kapitola o CLP(FD). ■ http://www.fi.muni.cz/~hanka/sicstus/doc/html / Příklady v distribuci SICStus Prologu: cca 60 příkladů, zdrojový kód ■ 1 i b/sicstus-'-'-'/T i brary/cl pfd/examples/ Probírané oblasti ■ Obsah ■ úvod: od LP k CLP ■ základy programování ■ základní algoritmy pro řešení problémů s omezujícími podmínkami ■ Příbuzné přednášky na Fl ■ PA163 Programování s omezujícími podmínkami ■ viz interaktivní osnova IS ■ PA167 Rozvrhování ■ http://www.fi.muni.cz/~hanka/rozvrhováni ■ zahrnuty CP techniky pro řešení rozvrhovacích problémů Hana Rudová, Logické programování 1,15. května 2013 169 Logické programování s omezujícími podmínkami Omezení (constraint) Dána ■ množina (doménových) proměnných Y = {y\,yt] ■ konečná množina hodnot (doména) D = {D\,..., DjJ Omezení c na ľ je podmnožina D\ x ... x Dt ■ omezuje hodnoty, kterých mohou proměnné nabývat současně Příklad: ■ proměnné: A,B ■domény: {0,1} pro A {1,2} pro B ■omezení: A^B nebo (A,B) e {(0,1 ),(0,2),(1,2)} Omezení c definováno na y\,.. .yt je splněno, pokud pro di e D\, ... dt e Dt platí (d\,... dt) e c ■ příklad (pokračování): omezení splněno pro (0,1), (0,2), (1,2), není splněno pro (1,1) Hana Rudová, Logické programování I, 1 5. května 201 3 171 Logické programování s omezujícími podmínkami Hana Rudová, Logické programování 1,15. května 201 3 1 70 Logické programování s omezujícími podmínkám Problém splňování podmínek (CSP) ■ Dána ■ konečná množina proměnných X = {x\,xn} ■ konečná množina hodnot (doména) D = {D\,... ,Dn} ■ konečná množina omezení C = {c\,..., cm} ■ omezení je definováno na podmnožině X Problém splňování podmínek je trojice (X,D, C) (constraint satisfaction problem) ■ Příklad: ■ proměnné: A,B,C ■domény: {0,1} pro A {1,2} pro B {0,2} pro C ■ omezení: A^B, B^C Hana Rudová, Logické programování 1,15. května 201 3 1 72 Logické programování s omezujícími podmínkám Řešení CSP Příklad: jednoduchý školní rozvrh Částečné ohodnocení proměnných (d\,dk),k < n ■ některé proměnné mají přiřazenu hodnotu Úplné ohodnocení proměnných (d\,..., dn) ■ všechny proměnné mají přiřazenu hodnotu Řešení CSP ■ úplné ohodnocení proměnných, které splňuje všechna omezení ■ (di,..., dn) e Di x ... x Dn je řešení (X,D, C) ■ pro každé Cí e C na Xiv .. .Xík platí (d^,.. .dík) e Cí Hledáme: jedno nebo všechna řešení nebo optimální řešení (vzhledem k objektivní funkci) Hana Rudová, Logické programování 1,15. května 2013 Logické programování s omezujícími podmínkami CLPCFD) program Základní struktura CLP programu 1. definice proměnných a jejich domén 2. definice omezení 3. hledání řešení (1) a (2) deklarativní část ■ modelování problému ■ vyjádření problému splňování podmínek (3) řídící část ■ prohledávání stavového prostoru řešení ■ procedura pro hledání řešení (enumeraci) se nazývá labeling ■ umožní nalézt jedno, všechna nebo optimální řešení Hana Rudová, Logické programování 1,15. května 2013 Logické programování s omezujícími podmínkami ■ proměnné: Jan, Petr, ... ■ domény: {3,4, 5,6}, {3,4},... ■ omezení: al l_di sti nct( [Jan, Petr, . . . ]) ■ částečné ohodnocení: Jan=6, Anna=5, Marie=l ■ úplné ohodnocení: Jan=6, Petr=3, Anna=5, 0ta=2, Eva=4, Marie=6 ■ řešení CSP: Jan=6, Petr=3, Anna=5, 0ta=2, Eva=4, Marie=l ■ všechna řešení: ještě Jan=6, Petr=4, Anna=5, 0ta=2, Eva=3, Marie=l ■ optimalizace: ženy učí co nejdříve Anna+Eva+Marie #= Cena minimalizace hodnoty proměnné Cena optimální řešení: Jan=6, Petr=4, Anna=5, 0ta=2, Eva=3, Marie=l učitel min max Jan 3 6 Petr 3 4 Anna 2 5 Ota 2 4 Eva 3 4 Marie 1 6 Hana Rudová, Logické programování 1,15. května 201 3 Logické programování s omezujícími podmínkami Kód CLPCFD) programu % základní struktura CLP programu so"lve( Variables ) :- declare_variables( Variables ), post_constraints( Variables ), labelingC Variables ). domain([]an],3,6], ... all_di sti nct([]an,Petr, .]) % triviální labeling labeling( [] ). labeling( [Var|Rest] ) :- fd_min(Var,Min) , ( Var#=Min, labeling( Rest ) % výběr nejmenší hodnoty z domény Var#>Min , labeling( [Var|Rest] ) Hana Rudová, Logické programování 1,15. května 201 3 Logické programování s omezujícími podmínkami Príklad: algebrogram Od LP k CLP I. Přiřaďte cifry 0, ... 9 písmenům S, E, N, D, M, O, R, Y tak, aby platilo: ■ SEND + MORE = MONEY ■ různá písmena mají přiřazena různé cifry ■ S a M nejsou 0 domain([E,N,D,0,R,Y], 0, 9), domain([S,M],1,9) 1000*S + 100*E + 10*N + D + 1000*M + 100*0 + 10*R + E #= 10000*M + 1000*0 + 100*N + 10*E + Y all_distinct( [S,E,N,D,M,0,R,Y] ) TabelingC [S,E,N,D,M,0,R,Y] ) Hana Rudová, Logické programování I, 1 5. května 201 3 1 77 Logické programování s omezujícími podmínkami Od LP k CLP II. Pro řešení podmínek se používají konzistenční techniky ■ consistency techniques, propagace omezení (constraint propagation) * omezení: A in 0. .2 , B in 0. .2 , B #< A domény po propagaci omezení B #< A: A in 1..2, B in 0..1 Podmínky jsou deterministicky vyhodnoceny v okamžiku volání podmínky Prohledávání doplněno konzistenčními technikami ■A in 1..2, B in 0..1, B# Vyvolání dvou řešičů: unifikace + řešič omezení Hana Rudová, Logické programování 1,15. května 201 3 1 80 Logické programování s omezujícími podmínkám Operační sémantika CLP ■ CLP výpočet cíle G ■ Store množina aktivních omezení = prostor omezení (constraint store) ■ inicializace Store = 0 ■ seznamy cílů v G prováděny v obvyklém pořadí ■ pokud narazíme na cíl s omezením c: NewStore = Store u {c} ■ snažíme se splnit c vyvoláním jeho řešiče ■ při neúspěchu se vyvolá backtracking ■ při úspěchu se podmínky v NewStore zjednoduší propagací omezení ■ zbývající cíle jsou prováděny s upraveným NewStore ■ CLP výpočet cíle G je úspěšný, pokud se dostaneme z iniciálního stavu (G, 0} do stavu #= | #\= | #< | #=< | #> | #>= ■ A + B #=< 3, A #\= CC - 4) * C D - 5), A/2 #= 4 ■ POZOR: neplést #=< a #>= s operátory pro implikaci: #<= #=> ■ sum(Variables,RelOp,Suma) ■ domain([A,B,C,F],1,3), sumC[A,B,C],#= ,F) ■ Variables i Suma musí být doménové proměnné nebo celá čísla ■ scalar_product(Coeffs,Variables,RelOp,Seal arProduct) ■ domain([A,B,C,F],1,6), scalar_productC [1,2,3],[A,B,C],#= ,F) ■ Variables i Val ue musí být doménové proměnné nebo celá čísla, Coeff s jsou celá čísla ■ POZOR na pořadí argumentů, nejprve jsou celočíselné koeficienty, pak dom. proměnné ■ scalar_productCCoeffs, Variables, #= , Value, [consistencyCdomain)]) ■ silnější typ konzistence ■ POZOR: domény musí mít konečné hranice Hana Rudová, Logické programování 1,15. května 2013 188 CLP(f D) v SICStus Prologu Základní globální omezení Všechny proměnné různé all_distinct(List) ■ všechny proměnné různé cumulative(...) ■ disjunktivní a kumulativní rozvrhování cumulati ves(...) ■ kumulativní rozvrhování na více zdrojů all_distinct(Variables), all_different(Variables) Proměnné v seznamu Variables jsou různé all_distinct a all_different se liší úrovní propagace ■ all_distinct má úplnou propagaci ■ all_different má slabší (neúplnou) propagaci Příklad: učitelé musí učit v různé hodiny ■ all_di sti nct C [Jan, Petr, Anna, Ota, Eva, Marie]) Jan = 6, Ota = 2, Anna = 5, Marie = 1, Petr in 3..4, Eva in 3..4 ■ all_different([Jan , Petr,Anna,Ota, Eva,Marie]) Jan in 3..6, Petr in 3..4, Anna in 2..5, Ota in 2..4, Eva in 3..4, Marie in 1..6 učitel min max Jan 3 6 Petr 3 4 Anna 2 5 Ota 2 4 Eva 3 4 Marie 1 6 Hana Rudová, Logické programování I, 15. května 2013 189 CLP(fD) v SICStus Prologu Disjunktivní rozvrhování (unární zdroj) cumulative([task(Start, Duration, End, 1, Id) | Tasks]) Rozvržení úloh zadaných startovním a koncovým časem (Start, End), dobou trvání (nezáporné Duration) a identifikátorem (Id) tak, aby se nepřekrývaly ■ příklad s konstantami: cumulative([task(0,2,2,l ,1), task(3,l ,4,1,2), task(5,l ,6,1,3)]) 1 2 3 1 1 2 3 4 5 6 ■ příklad: vytvoření rozvrhu, za předpokladu, že doba trvání hodin není stejná JanE#= Jan+3, PetrE#= Petr+1, AnnaE#= Anna+2, ... cumulativeCtaskCJan,3,JanE,1,1),task(Petr,1,PetrE,1,2),task(Anna,2,AnnaE,1, task(0ta,2,OtaE,1,4),task(Eva,2,EvaE,1,5),task(Mari e,3,Mari eE,1,6)]) Hana Rudová, Logické programování I, 15. května 2013 191 CLP(fD) v SICStus Prologu Hana Rudová, Logické programování 1,15. května 2013 190 CLP(fD) v SICStus Prolog Kumulativní rozvrhování ■ cumulative([task(Start,Duration,End,Demand,TaskId) | Tasks], [li mi t (Limit)]) ■ Rozvržení úloh zadaných startovním a koncovým časem (Start, End), dobou trvání (nezáporné Duration), požadovanou kapacitou zdroje (Demand) a identifikátorem (Id) tak, aby se nepřekrývaly a aby celková kapacita zdroje nikdy nepřekročila Limit ■ Příklad s konstantami: cumulative([task(0,4,4,l,l),task(l,2,3,2,2),task(3,3,6,2,3),task(4,2,6,l,4)],[limit(3)]) 2 3 1 1 1 1 1 1 2 3 4 5 6 Hana Rudová, Logické programování 1,15. května 2013 1 92 CLP(fD) v SICStus Prolog Kumulativní rozvrhování s více zdroji Příklad: kumulativní rozvrhování ■ Rozvržení úloh tak, aby se nepřekrývaly a daná kapacita zdrojů nebyla překročena (limit zdroje chápán jako horní mez - bound (upper)) ■ cumulatives([tas k(Start,Duration,End,Demand,Machi neld)|Tasks], [machi ne(Id,Limit)|Machi nes],[bound(upper)]) ■ Úlohy zadány startovním a koncovým časem (Start, End), dobou trvání (nezáporné Duration), požadovanou kapacitou zdroje (Demand) a požadovaným typem zdroje (Machi neld) ■ Zdroje zadány identifikátorem (Id) a kapacitou (Limit) ■ Příklad: ?- domain([B,C],1,2), cumulatives([task(0,4,4,l,l),task(3,1,4,1,B), task(5,l,6,l,C)] [machine(l,l),machine(2,1)], [bound(upper)]). C in 1..2, B=2 Hana Rudová, Logické programování 1,15. května 2013 193 CLP(f D) v SICStus Prologu Řešení: kumulativní rozvrhování | ?- schedule(13, [16,6,13,7,5,18,4], [2,9,3,7,10,1,11], 69, Ss, End). Ss = [0,16,9,9,4,4,0], End = 22 ? schedule(l_imit, Ds, Rs, MaxCas, Ss, End) :-domainCSs, 0, MaxCas), End in 0..MaxCas, vytvor_u]ohy(Ss,Ds,Rs,1,Tasks), cumulative(Tasks, []i mi t(Li mit)]), afterCSs, Ds , End), % koncový čas appendCSs, [End], Vars), ]abe]ing([minimize(End)],Vars). vytvor_ulohy([] ,[],[] ,_Id, []) . vytvor_ulohy([S|Ss] , [D|Ds], [R|Rs], Id, [task(S,D,E,R,Id)|Tasks]):-Newld i s Id+1, E #= S+D, vytvor_u]ohy(Ss,Ds,Rs, Newld,Tasks). after([], [] , _) . after([S|Ss], [D|Ds], End) :- E #>= S+D, afterCSs, Ds, End). Hana Rudová, Logické programování I, 15. května 2013 195 CLP(fD) v SICStus Prologu ■ Vytvořte rozvrh pro následující úlohy, tak aby nebyla překročena kapacita 1 3 zdroje, a minimalizujte celkovou dobu trvání úloha doba trvání kapacita tl 16 2 t2 6 9 t3 13 3 t4 7 7 t5 5 10 t6 18 1 t7 4 11 Hana Rudová, Logické programování 1,15. května 2013 194 CLP(fD) v SICStus Prologu Vestavěné predikáty pro labeling ■ Instanciace proměnné Variable hodnotami v její doméně indomain( Variable ) hodnotyjsou instanciovány při backtrackingu ve vzrůstajícím pořadí ?- X in 4..5, indomain(X). X = 4 ? ; X = 5 ? labelingC [] ). labelingC [Var|Rest] ) :- % výběr nejlevější proměnné k instanciaci indomainC Var ), % výběr hodnot ve vzrůstajícím pořadí labelingC Rest ). ■ labe"ling( Options, Variables ) ?- A in 0..2, B in 0..2, B#< A, 1abeling([], [A,B]). Hana Rudová, Logické programování 1,15. května 2013 196 CLP(fD) v SICStus Prologu Uspořádání hodnot a proměnných Výběr hodnoty ■ Při prohledávání je rozhodující uspořádání hodnot a proměnných ■ Určují je heuristiky výběru hodnot a výběru proměnných labelingC [] ). labelingC Variables ) :- select_variable(Variables,Var,Rest), select_value(Var,Value), C Var #= Value, labelingC Rest ) Var #\= Value , % nemusí dojít k instanciaci Var labelingC Variables ) % proto pokračujeme se všemi proměnnými včetně Var ). ■ Statické uspořádání: určeno už před prohledáváním ■ Dynamické uspořádání: počítá se během prohledávání Hana Rudová, Logické programování 1,15. května 2013 197 CLP(f D) v SICStus Prologu Výběr proměnné ■ Obecný princip výběru proměnné: first-fail ■ výběr proměnné, pro kterou je nejobtížnější nalézt správnou hodnotu pozdější výběr hodnoty pro tuto proměnnou by snadněji vedl k failu ■ vybereme proměnnou s nejmenší doménou ■ ?- domai n( [A, B, C] , 1,3) , A#<3 , A#=B+C. nejlépe je začít s výběrem A ■ Parametry label i ng/2 ovlivňující výběr proměnné ■ leftmost: nejlevější (default) ■ ff: s (1) nejmenší velikostí domény fd_size(Var,Size) (2) (pokud s nejmenší velikostí domény více, tak) nejlevější z nich ■ f f c: s (1) nejmenší velikostí domény (2) největším množstvím omezením „čekajících" na proměnné fd_degree(Var,Size) (3) nejlevější z nich ■ mi n/max: s (1) nejmenší/největší hodnotou v doméně proměnné (2) nejlevnějšíz nich fd_min(Var,Min)/fd_max(Var,Max) Hana Rudová, Logické programování 1,15. května 2013 199 CLP(f D) v SICStus Prologu ■ Obecný princip výběru hodnoty: první úspěch (succeed first) ■ volíme pořadí tak, abychom výběr nemuseli opakovat ■ ?- domai n( [A, B, C] , 1,2) , A#=B+C. optimální výběr A=2,B=1 ,C=1 je bez backtrackingu ■ Parametry 1 abel i ng/2 ovlivňující výběr hodnoty př. labeling([down], Vars) ■ up: doména procházena ve vzrůstajícím pořadí (default) ■ down: doména procházena v klesajícím pořadí ■ Parametry 1 abel i ng/2 řídící, jak je výběr hodnoty realizován ■ step: volba mezi X #= M, X #\= M (default) ■ viz dřívější příklad u "Uspořádání hodnot a proměnných" ■ enum: vícenásobná volba mezi všemi hodnotami v doméně ■ podobně jako při indomain/1 Hana Rudová, Logické programování 1,15. května 2013 198 CLP(f D) v SICStus Prologu Hledání optimálního řešení (předpokládejme minimalizaci) ■ Parametry 1 abel i ng/2 pro optimalizaci: minimize(F)/maximize(F) ■ Cena #= A+B+C, labeling([minimize(Cena)], [A,B,C]) ■ Metoda větví a mezí (branch&bound) ■ algoritmus, který implementuje proceduru pro minimalizaci (duálně pro maximalizaci) ■ uvažujeme nejhorší možnou cenu řešení UB (např. cena už nalezeného řešení) ■ počítáme dolní odhad LB ceny částečného řešení Li? je tedy nejlepší možná cena pro rozšíření tohoto řešení ■ procházíme strom a vyžadujeme, aby prozkoumávaná větev měla cenu LB < UB pokud je LB > UB, tak víme, že v této větvi není lepší řešení a odřízneme ji ■ přidává se tedy inkrementálně omezení LB# 0 C4 : X2 + Xs - X% = 0 Hana Rudová, Logické programování 1,15. května 201 3 Algoritmy pro CSP Vrcholová a hranová konzistence Vrcholová konzistence (node consistency) NC ■ každá hodnota z aktuální domény V proměnné splňuje všechny unární podmínky s proměnnou Ví Hranová konzistence (arc consistency) AC pro binární CSP ■ hrana (Ví, V}) je hranově konzistentní, právě když pro každou hodnotu x z aktuální domény Dí existuje hodnota y tak, že ohodnocení [Ví = x,Vj = y] splňuje všechny binární podmínky nad Ví, V,. ■ hranová konzistence je směrová ■ konzistence hrany (VuVj) nezaručuje konzistenci hrany (Vj,Ví) A|3..7|ae= = H1-5| B A|3..4|<- - J| 1 ..5| B A|3..4|^=^|4..5| B A v ■ hranu (m, fe) vedoucí z proměnné Vm, která zmenšení domény způsobila, není třeba revidovat (změna sejí nedotkne) Hana Rudová, Logické programování 1,15. května 2013 206 Algoritmy pro CSP Algoritmus AC-3 procedure AC-3CC) Q := {(í.j) I Cíi j) e hrany(G), i 4= j} % seznam hran pro revizi while Q non empty do vyber a smaž (fe,m) z Q if revise((fc,m)) then % přidáni pouze hran, které Q := Q u {(í, k) e hrany(G), í 4 k, í4=in} % dosud nejsou ve fronte end while end AC-3 ®^®^® A řešení neexistuje ■ všechny domény jsou jednoprvkové => máme řešení ■ v obecném případě se alespoň zmenší prohledávaný prostor NE NE Hana Rudová, Logické programování 1,15. května 201 3 Algoritmy pro CSP k-konzistence Silná k-konzistence Mají NC a AC něco společného? ■ NC: konzistence jedné proměnné ■ AC: konzistence dvou proměnných "... můžeme pokračovat CSPje k-konzistentní právě tehdy, když můžeme libovolné konzistentní ohodnocení (k-1) různých proměnných rozšířit do libovolné k-té proměnné 4-konzistentní graf ■ Pro obecné CSP, tedy i pro nebinární podmínky Hana Rudová, Logické programování 1,15. května 2013 209 Algoritmy pro CSP 3-konzistentní graf (1.1) lze rozšířit na (1,1,1) (2.2) lze rozšířit na (2,2,2) 1,2 _ 1,2 _ 1 ?3 není 2-konzistentní (3) nelze rozšířit (1,3) ani (2,3) nejsou konzistentní dvojice (nerozšiřujeme je) ■ CSPje silně k-konzistentní právě tehdy, když je j-konzistentní pro každé j k-konzistence ■ Silná k-konzistence => j-konzistence V j < k ■ k-konzistence =fr silná k-konzistence ■ NC = silná 1-konzistence = 1-konzistence ■ AC = (silná) 2-konzistence Hana Rudová, Logické programování 1,15. května 201 3 Algoritmy pro CSP Konzistence pro nalezení řešení ■ Máme-li graf s n vrcholy, jak silnou konzistenci potřebujeme, abychom přímo našli řešení? ■ silná n-konzistence je nutná pro graf s n vrcholy ■ n-konzistence nestačí (viz předchozí příklad) ■ silná k-konzistence pro k min(A) = min(B)+min(C), max(A) = max(B)+max(C) min(B) = mi n (A)-max (C) , max(B) = max (A) -mi n (C) mi n(C) = mi n(A)-max(B) , max(C) = max(A)-mi n (B) ■ změna mi n (A) vyvolá pouze změnu mi n (B) a mi n (C) ■ změna max(A)vyvolá pouze změnu max(B) amax(C), ... ■ Příklad: A in 1..10, B in 1.. 10, A #= B + 2, A #> 5, A #\= 8 A #= B + 2 => min(A)=l+2, max(A)=10+2 => A in 3. . 10 => min(B)=l-2, max(B)=10-2 => B in 1..8 A #> 5 => min(A)=6 => A in 6.. 10 => min(B)=6-2 => B in 4.. 8 (nové vyvolání A #= B + 2) A #\= 8 => A i n (6..7) \/ (9.. 10) (meze stejné, k propagaci A #= B + 2 nedojde) ■ Vyzkoušejte si: A #= B - C, A #>= B + C Konzistence mezí ■ Bounds consistency BC: slabší než obecná hranová konzistence ■ podmínka má konzistentní meze (BC), právě když pro každou proměnnou Vj z této podmínky a každou hodnou x e D j existuje ohodnocení zbylých proměnných v podmínce tak, že je podmínka splněna a pro vybrané ohodnocení yi proměnné Ví platí min(Dí) < yi < max(Dj) ■ stačí propagace pouze při změně minimální nebo maximální hodnoty (při změně mezí) v doméně proměnné ■ Konzistence mezí pro nerovnice ■ A #> B=> min(A) = min(B)+l, max(B) = max(A)-l ■ příklad: A in 4. . 10, B in 6. . 18, A #> B mi n (A) = 6+1 => A in 7. . 10 max(B) = 10-1 => B in 6. . 9 ■ podobně: A #< B, A #>= B, A #=< B Hana Rudová, Logické programování 1,15. května 2013 214 Algoritmy pro CSP Globální podmínky ■ Propagace je lokální ■ pracuje se s jednotlivými podmínkami ■ interakce mezi podmínkami je pouze přes domény proměnných ■ Jak dosáhnout více, když je silnější propagace drahá? ■ Seskupíme několik podmínek do jedné tzv. globální podmínky ■ Propagaci přes globální podmínku řešíme speciálním algoritmem navrženým pro danou podmínku ■ Příklady: ■ all_distinct omezení: hodnoty všech proměnných různé ■ serialized omezení: rozvržení úloh zadaných startovním časem a dobou trvání tak, aby se nepřekrývaly Hana Rudová, Logické programování I, 15. května 2013 21 5 Algoritmy pro CSP Hana Rudová, Logické programování I, 15. května 2013 216 Algoritmy pro CSP Propagace pro all_ch"stinet U = {X2, X4, X5}, dom(U) = {2, 3, 4}: {2,3,4} nelze pro XI, X3, X6 XI in 5..6, X3 = 5, X6 in {1} \/ (5..6) Konzistence: V{Xi,.. .,Xk} c V : card{Di u ■ ■ ■ u Dk} > k stačí hledat Hallův interval J: velikost intervalu / je rovna počtu proměnných, jejichž doména je v J Inferenční pravidlo ■ U = {Xi,...,Xk}, dom(U) = {Di u ■ ■ ■ uDk] ' card(U) = card(dom(U)) ^Vve dom(U), VX e (V - U),X ± v ■ hodnoty v Hallově intervalu jsou pro ostatní proměnné nedostupné Složitost: 0(2") - hledání všech podmnožin množiny n proměnných (naivní) O(nlogn) - kontrola hraničních bodů Hallových intervalů (1 998) učitel min max Jan 3 6 Petr 3 4 Anna 2 5 Ota 2 4 Eva 3 4 Marie 1 6 Hana Rudová, Logické programování I, 1 5. května 201 3 Algoritmy pro CSP Prohledávání do hloubky Základní prohledávací algoritmus pro problémy splňování podmínek Prohledávání stavového prostoru do hloubky (depth first search) Dvě fáze prohledávání s navracením ■ dopředná fáze: proměnné jsou postupně vybírány, rozšiřuje se částečné řešení přiřazením konzistení hodnoty (pokud existuje) další proměnné ■ po vybrání hodnoty testujeme konzistenci ■ zpětná fáze: pokud neexistuje konzistentní hodnota pro aktuální proměnnou, algoritmus se vrací k předchozí přiřazené hodnotě Proměnné dělíme na ■ minulé - proměnné, které už byly vybrány (a mají přiřazenu hodnotu) ■ aktuální - proměnná, která je právě vybrána a je jí přiřazována hodnota ■ budoucí - proměnné, které budou vybrány v budoucnosti Hana Rudová, Logické programování 1,15. května 2013 Algoritmy pro CSP Prohledávání + konzistence Splňování podmínek prohledáváním prostoru řešení ■ podmínky jsou užívány pasivně jako test ■ přiřazuji hodnoty proměnných a zkouším co se stane ■ vestavěný prohledávací algoritmus Prologu: backtracking, triviální: generuj & testuj ■ úplná metoda (nalezneme řešení nebo dokážeme jeho neexistenci) ■ zbytečně pomalé (exponenciální): procházím i „evidentně" špatná ohodnocení Konzistenční (propagační) techniky ■ umožňují odstranění nekonzistentních hodnot z domény proměnných ■ neúplná metoda (v doméně zůstanou ještě nekonzistentní hodnoty) ■ relativně rychlé (polynomiální) Používá se kombinace obou metod ■ postupné přiřazování hodnot proměnným ■ po přiřazení hodnoty odstranění nekonzistentních hodnot konzistenčními technikami Hana Rudová, Logické programování 1,15. května 201 3 218 Algoritmy pro CSP Základní algoritmus prohledávání do hloubky 1 Pro jednoduchost proměnné očíslujeme a ohodnocujeme je v daném pořadí 1 Na začátku voláno jako 1 abel i ng (C, 1) procedure labeling(G,a) if a > |uzlyCC)| then return uzly(G) for V x e Da do if consistent(G,a) then % consistent(G,a) je nahrazeno FC(G,a), LA(G,<; R := Tabel ing(G,a +1) if R 4= fail then return R return fail end labeling Po přiřazení všech proměnných vrátíme jejich ohodnocení 1 Procedury consistent uvedeme pouze pro binární podmínky Hana Rudová, Logické programování 1,15. května 201 3 Algoritmy pro CSP Backtracking (BT) Príklad: backtracking Backtracking ověřuje v každém kroku konzistenci podmínek vedoucích z minulých proměnných do aktuální proměnné Backtracking tedy zajišťuje konzistenci podmínek ■ na všech minulých proměnných ■ na podmínkách mezi minulými proměnnými a aktuální proměnnou procedure BT(G,a) Q:={(Ví, V a) e hranyCO , i < a] % hrany vedoucí z minulých proměnných do aktuální Consistent := true while Q neni prázdná a Consistent do vyber a smaž libovolnou hranu (VJt,Vm) z Q Consistent := not reviseiV^, Vm) % pokud vyradíme prvek, bude doména prázdná return Consistent end BT Hana Rudová, Logické programování I, 15. května 2013 221 Algoritmy pro CSP Kontrola dopředu (FC - forward checking) FCje rozšíření backtrackingu FC navíc zajišťuje konzistenci mezi aktuální proměnnou a budoucími proměnnými, které jsou s ní spojeny dosud nesplněnými podmínkami procedure FC(G,a) Q:={(Ví, V a) e h rany (C) , í > a] % přidání hran z budoucích do aktuální proměnné Consistent := true while Q neni prázdná a Consistent do vyber a smaž libovolnou hranu (VJt,Vm) z Q if revi se((VJt, Vm)) then Consistent := (|Djt|>0) % vyprázdnění domény znamená nekonzistenci return Consistent end FC Hrany z minulých proměnných do aktuální proměnné není nutno testovat ■ Omezení: Vu V2, V3 in 1... 3, V1# = 3xV3 ■ Stavový prostor: 4 ■ červené čtverečky: chybný pokus o instanciaci, řešení neexistuje ■ nevyplněná kolečka: nalezeno řešení ■ černá kolečka: vnitřní uzel, máme pouze částečné přiřazení Hana Rudová, Logické programování 1,15. května 201 3 222 Příklad: kontrola dopředu ■ Omezení: Vlt V2, V3 in 1... 3, c:V1#=3x V3 ■ Stavový prostor: C => fail C => fail Hana Rudová, Logické programování I, 15. května 2013 223 Algoritmy pro CSP Hana Rudová, Logické programování 1,15. května 2013 224 Pohled dopředu (LA - looking aheacl) LAje rozšíření FC, navíc ověřuje konzistenci hran mezi budoucími proměnnými % začínáme s hranami do a procedure LA(G,a) Q := {(VuVa) e h rany (C), í > a] Consistent := true while Q neni prázdná a Consistent do vyber a smaž libovolnou hranu (Vk,Vm) z Q if revi seC(Vk, Vm)) then Q := Q u {(VuVk)\(VuVk) e hrany(G), í + k, í + m, í > a] Consistent := (\Dk\ > 0) return Consistent end LA ■ Hrany z minulých proměnných do aktuální proměnné opět netestujeme ■ Tato LA procedura je založena naAC-3, lze použít i jiné AC algoritmy LA udržuje hranovou konzistenci: protože ale LA(G,a) používá AC-3, musíme zajistit iniciální konzistenci pomocí AC-3 ještě před startem prohledávání Příklad: pohled dopředu (pomocí AC-3) Omezení: V1.y2.V3 in 1...4, cl:V1#>V2, c2 : V2# = 3 x V3 Stavový prostor (spouští se iniciální konzistence se před startem prohledávání) ^ d => V1 in 2..4 V2in 1..3 C2 => V2=3 V3= 1 d => V1= 4 Ví 41 V2 3, Vo 1j v3 o Hana Rudová, Logické programování I, 1 5. května 201 3 Algoritmy pro CSP Hana Rudová, Logické programování 1,15. května 201 3 Algoritmy pro CSP Přehled algoritmů Cvičení proměnné BT Backtracking (BT) kontroluje v kroku a podmínky c(Vi,Vfl),...,c(Vfl-i,Vfl) z minulých proměnných do aktuální proměnné Kontrola dopředu (FC) kontroluje v kroku a podmínky \^ c(V«+i,V«),...,c(V„,Vfl) z budoucích proměnných do aktuální proměnné Pohled dopředu (LA) kontroluje v kroku a podmínky ^ /I Vi(a < l < n), Vfe(a < fe 1 0. Implementace Prologu 2. Ukažte, jak je dosaženo hranové konzistence v následujícím příkladu: domain([X,Y,Z],l ,5]), X #< Y, Z#= Y+l . Literatura: ■ Matýska L., Toman D.: Implementační techniky Prologu, Informační systémy, (1 990), 21-59. http://www.i cs.muni.cz/people/matyska/vyuka/1 p/l p.html Hana Rudová, Logické programování I, 15. května 2013 229 Algoritmy pro CSP Opakování: základní pojmy ■ Konečná množina klauzulí Hlava :- Tělo tvoří program P. ■ Hlava je literál ■ Tělo je (eventuálně prázdná) konjunkce literálů Ti,... Ta, a > 0 ■ Literál je tvořen m-árním predikátovým symbolem (m/p) a m termy (argumenty) ■ Term je konstanta, proměnná nebo složený term. ■ Složený term s n termy na místě argumentů ■ Dotaz (cíl) je neprázdná množina literálů. Hana Rudová, Logické programování 1,15. května 2013 231 Implementace Prologu Interpretace Deklarativní sémantika: Hlava platí, platí-li jednotlivé literály těla. Procedurální (imperativní) sémantika: Entry: Hlava:: { call Ti call Ta } Volání procedury s názvem Hlava uspěje, pokud uspěje volání všech procedur (literálů) v těle. Procedurální sémantika = podklad pro implementaci Hana Rudová, Logické programování 1,15. května 2013 232 Implementace Prologu Abstraktní interpret Abstraktní interpret - pokračování Vstup: Logický program P a dotaz C. 1. Inicializuj množinu cílů S literály z dotazu C; S:=C 2. while ( S != empty ) do 3. Vyber AeS a dále vyber klauzuli A' : -Bi, . . . ,B„ (n > 0) z programu P takovou, že 3cr : Aer = A' cr; cr je nejobecnější unifikátor. Pokud neexistuje A' nebo cr, ukonči cyklus. 4. Nahraď A v S cíli Bi až B„. 5. Aplikuj a na C a S. 6. end while 7. Pokud S==empty, pak výpočet úspešne skončil a výstupem je C se všemi aplikovanými substitucemi. Pokud S!=empty, výpočet končí neúspěchem. Hana Rudová, Logické programování 1,15. května 2013 Implementace Prologu Nedeterminismus interpetu 1. Selekční pravidlo: výběr cíle A z množiny cílů S ■ neovlivňuje výrazně výsledek chování interpretu 2. Způsob prohledávání stromu výpočtu: výběr klauzule A' z programu P ■ je velmi důležitý, všechny klauzule totiž nevedou k úspěšnému řešení Vztah k úplnosti: 1. Selekční pravidlo neovlivňuje úplnost ■ možno zvolit libovolné v rámci SLD rezoluce 2. Prohledávání stromu výpočtu do šířky nebo do hloubky „Prozření" - automatický výběr správné klauzule ■ vlastnost abstraktního interpretu, kterou ale reálné interprety nemají Hana Rudová, Logické programování 1,15. května 2013 Implementace Prologu Kroky (3) až (5) představují redukci (logickou inferenci) cíle A. Počet redukcí za sekundu (LIPS) == indikátor výkonu implementace Věta Existuje-li instance C dotazu C, odvoditelná z programu P v konečném počtu kroků, pak bude tímto interpretem nalezena. Hana Rudová, Logické programování 1,15. května 201 3 Implementace Prologu Prohledávání do šířky 1. Vybereme všechny klauzule Ar, které je možno unifikovat s literálem A ■ nechť je těchto klauzulí q 2. Vytvoříme q kopií množiny S 3. V každé kopii redukujeme A jednou z klauzulí Ar. ■ aplikujeme příslušný nejobecnější unifikátor 4. V následujících krocích redukujeme všechny množiny Sí současně. 5. Výpočet ukončíme úspěchem, pokud se alespoň jedna z množin Sí stane prázdnou. ■ Ekvivalence s abstraktnímu interpretem ■ pokud jeden interpret neuspěje, pak neuspěje i druhý ■ pokud jeden interpret uspěje, pak uspěje i druhý Hana Rudová, Logické programování 1,15. května 201 3 Implementace Prologu Prohledávání do hloubky 1. Vybereme všechny klauzule A' i, které je možno unifikovat s literálem A. 2. Všechny tyto klauzule zapíšeme na zásobník. 3. Redukci provedeme s klauzulí na vrcholu zásobníku. 4. Pokud v nějakém kroku nenajdeme vhodnou klauzuli A', vrátíme se k předchozímu stavu (tedy anulujeme aplikace posledního unifikátoru cr) a vybereme ze zásobníku další klauzuli. 5. Pokud je zásobník prázdný, končí výpočet neúspěchem. 6. Pokud naopak zredukujeme všechny literály v S, výpočet končí úspěchem. ■ Není úplné, tj. nemusí najít všechna řešení ■ Nižší paměťová náročnost než prohledávání do šířky ■ Používá se v Prologu Hana Rudová, Logické programování I, 15. května 2013 237 Implementace Prologu Reprezentace objektů ■ Beztypový jazyk ■ Kontrola „typů" za běhu výpočtu ■ Informace o struktuře součástí objektu Typy objektů ■ Primitivní objekty: ■ konstanta ■ číslo ■ volná proměnná ■ odkaz (reference) ■ Složené (strukturované) objekty: ■ struktura ■ seznam Hana Rudová, Logické programování 1,15. května 2013 238 Implementace Prologu Reprezentace objektů II Příznaky (tags): Objekt Příznak volná proměnná FREE konstanta CONST celé číslo INT odkaz REF složený term FUNCT Obsah adresovatelného slova: hodnota a příznak. Primitivní objekty uloženy přímo ve slově Složené objekty ■ jsou instance termu ve zdrojovém textu, tzv. zdrojového termu ■ zdrojový term bez proměnných => každá instanciace ekvivalentní zdrojovému termu ■ zdrojový term s proměnnými => dvě instance se mohou lišit aktuálními hodnotami proměnných, jedinečnost zajišťuje kopírování struktur nebo sdílení struktur Hana Rudová, Logické programování 1,15. května 2013 Implementace Prologu Příklad: a(b(X),c(X,Y),d), Kopírování struktur FUNCT a/3 REF REF CONST d FUNCT c/2 REF FREE Y FUNCT b/1 FREE X Hana Rudová, Logické programování 1,15. května 201 3 Implementace Prologu Kopírování struktur II Term F s aritou A reprezentován A+l slovy: ■ funktor a arita v prvním slově ■ 2. slovo nese první argument (resp. odkaz najeho hodnotu) ■ ■ A+l slovo nese hodnotu A-tého argumentu Reprezentace vychází z orientovaných acyklických grafů: a/3 c/2 b/l X Vykopírována každá instance => kopírování struktur Termy ukládány na globální zásobník Hana Rudová, Logické programování 1,15. května 2013 Implementace Prologu Sdílení struktur II Příklad: a(b(X),c(X,Y),d) reprezentuje < a(b($l),c($l,$2),d) ; [FREE, FREE] > kde symbolem $i označujeme i-tou proměnnou. Implementace: < &kostra_termu; &rámec > (& vrací adresu objektu) Všechny instance sdílí společnou kostru_termu => sdílení struktur Sdílení struktur Vychází z myšlenky, že při reprezentaci je třeba řešit přítomnost proměnných Instance termu < kostra_termu; rámec > ■ kostra_termu je zdrojový term s očíslovanými proměnnými ■ rámec je vektor aktuálních hodnot těchto proměnných ■ i-tá položka nese hodnotu i-té proměnné v původním termu Hana Rudová, Logické programování 1,15. května 201 3 242 Implementace Prologu Srovnání: příklad Naivní srovnání: sdílení paměťově méně náročné Platí ale pouze pro rozsáhlé termy přítomné ve zdrojovém kódu Postupná tvorba termů: A = a(K,L,M), K = b(X), L = c(X,Y), M = d ■ Sdílení termů: kostra_a M:d kostra_b ^ kostra_c X —i Hana Rudová, Logické programování I, 15. května 2013 243 Implementace Prologu Hana Rudová, Logické programování 1,15. května 2013 244 Implementace Prologu Srovnání: příklad - pokračování Kopírování struktur: A = a(K,L,M), K = b(X), L = c(X,Y), M = d FUNCT a/3 REF - REF - CONST d FUNCT c/2 - REF FREE Y FUNCT b/1 FREE X tj. identické jako přímé vytvoření termu a(b(X) , c(X,Y) ,d) Hana Rudová, Logické programování I, 15. května 2013 245 Implementace Prologu Řízení výpočtu Dopředný výpočet ■ po úspěchu (úspěšná redukce) ■ jednotlivá volání procedur skončí úspěchem ■ klasické volání rekurzivních procedur Zpětný výpočet (backtracking) ■ po neúspěchu vyhodnocení literálu (neúspěšná redukce) ■ nepodaří se unifikace aktuálních a formálních parametrů hlavy ■ návrat do bodu, kde zůstala nevyzkoušená alternativa výpočtu ■ je nutná obnova původních hodnot jednotlivých proměnných ■ po nalezení místa s dosud nevyzkoušenou klauzulí pokračuje dále dopředný výpočet Srovnání II ■ Složitost algoritmů pro přístup k jednotlivým argumentům ■ sdílení struktur: nutná víceúrovňová nepřímá adresace ■ kopírování struktur: bez problémů ■ jednodušší algoritmy usnadňují i optimalizace ■ Lokalita přístupů do paměti ■ sdílení struktur: přístupy rozptýleny po paměti ■ kopírování struktur: lokalizované přístupy ■ při stránkování paměti - rozptýlení vyžaduje přístup k více stránkám ■ Z praktického hlediska neexistuje mezi těmito přístupy zásadní rozdíl Hana Rudová, Logické programování 1,15. května 2013 246 Implementace Prologu Aktivační záznam ■ Volání (=aktivace) procedury ■ Aktivace sdílí společný kód, liší se obsahem aktivačního záznamu ■ Aktivační záznam uložen na lokálním zásobníku ■ Dopředný výpočet ■ stav výpočtu v okamžiku volání procedury ■ aktuální parametry ■ lokální proměnné ■ pomocné proměnné ('a la registry) ■ Zpětný výpočet (backtracking) ■ hodnoty parametrů v okamžiku zavolání procedury ■ následující klauzule pro zpracování při neúspěchu Hana Rudová, Logické programování I, 15. května 2013 247 Implementace Prologu Hana Rudová, Logické programování I, 15. května 2013 248 Implementace Prologu Aktivační záznam a roll-back Neúspěšná klauzule mohla nainstanciovat nelokální proměnné ■ a(X) : - X = b(c,Y), Y = d. ?- W = b(Z,e), a(W) . (viz instanciace Z) Při návratu je třeba obnovit (roll-back) původní hodnoty proměnných Využijeme vlastností logických proměnných ■ instanciovat lze pouze volnou proměnnou ■ jakmile proměnná získá hodnotu, nelze ji změnit jinak než návratem výpočtu => původní hodnoty všech proměnných odpovídají volné proměnné Stopa (trail): zásobník s adresami instanciovaných proměnných ■ ukazatel na aktuální vrchol zásobníku uchováván v aktivačním záznamu ■ při neúspěchu jsou hodnoty proměnných na stopě v úseku mezi aktuálním a uloženým vrcholem zásobníku změněny na „volná" Globální zásobník: pro uložení složených termů ■ ukazatel na aktuální vrchol zásobníku uchováván v aktivačním záznamu ■ při neúspěchu vrchol zásobníku snížen podle uschované hodnoty v aktivačním záznamu Hana Rudová, Logické programování 1,15. května 2013 249 Implementace Prologu Řez Prostředek pro ovlivnění běhu výpočtu programátorem ■ a(X) :- b(X), !, c(X). a(3). b(l). b(2). c(l). c(2). Řez: neovlivňuje dopředný výpočet, má vliv pouze na zpětný výpočet Odstranění alternativních větví výpočtu => odstranění odpovídajících bodů volby ■ tj. odstranění bodů volby mezi současným vrcholem zásobníku a bodem volby procedury, která řez vyvolala (včetně bodu volby procedury s řezem) => změna ukazatele na „nejmladší" bod volby Vytváření deterministických procedur Optimalizace využití zásobníku Hana Rudová, Logické programování 1,15. května 2013 251 Implementace Prologu Okolí a bod volby Aktivační záznam úspěšně ukončené procedury nelze odstranit z lokálního zásobníku => rozdělení aktivačního záznamu: ■ okolí (environment) - informace nutné pro dopředný běh programu ■ bod volby (choice point) - informace nezbytné pro zotavení po neúspěchu ■ ukládány na lokální zásobník ■ samostatně provázány (odkaz na předchozí okolí resp. bod volby) Důsledky: ■ samostatná práce s každou částí aktivačního záznamu (optimalizace) ■ alokace pouze okolí pro deterministické procedury ■ možnost odstranění okolí po úspěšném vykonání (i nedeterministické) procedury (pokud okolí následuje po bodu volby dané procedury) ■ pokud je okolí na vrcholu zásobníku Hana Rudová, Logické programování 1,15. května 2013 250 Implementace Prologu Interpret Prologu Základní principy: ■ klauzule uloženy jako termy ■ programová databáze ■ pro uložení klauzulí ■ má charakter haldy ■ umožňuje modifikovatelnost prologovských programů za běhu (assert) ■ klauzule zřetězeny podle pořadí načtení ■ triviální zřetězení Vyhodnocení dotazu: volání procedur řízené unifikací Hana Rudová, Logické programování 1,15. května 2013 252 Implementace Prologu Interpret - Základní princip 1. Vyber redukovaný literál („první", tj. nejlevější literál cíle) 2. Lineárním průchodem od začátku databáze najdi klauzuli, jejíž hlava má stejný funktor a stejný počet argumentů jako redukovaný literál 3. V případě nalezení klauzule založ bod volby procedury 4. Založ dále okolí první klauzule (velikost odvozena od počtu lokálních proměnných v klauzuli) 5. Proveď unifikaci literálu a hlavy klauzule 6. Úspěch => přidej všechny literály klauzule k cíli („doleva", tj. na místo redukovaného literálu). Tělo prázdné => výpočet se s úspěchem vrací do klauzule, jejíž adresa je v aktuálním okolí. 7. Neúspěch unifikace => z bodu volby se obnoví stav a pokračuje se v hledání další vhodné klauzule v databázi. 8. Pokud není nalezena odpovídající klauzule, výpočet se vrací na předchozí bod volby (krátí se lokální i globální zásobník). 9. Výpočet končí neúspěchem: neexistuje již bod volby, k němuž by se výpočet mohl vrátit. 1 0. Výpočet končí úspěchem, jsou-li úspěšně redukovány všechny literály v cíli. Hana Rudová, Logické programování I, 1 5. května 2013 253 Implementace Prologu Optimalizace: Indexace ■ Zřetězení klauzulí podle pořadí načtení velmi neefektivní ■ Provázání klauzulí se stejným funktorem a aritou hlavy (tvoří jednu proceduru) ■ tj., indexace procedur ■ Hash tabulka pro vyhledání první klauzule ■ Možno rozhodnout (parciálně) determinismus procedury Interpret - vlastnosti ■ Lokální i globální zásobník ■ při dopředném výpočtu roste ■ při zpětném výpočtu se zmenšuje Lokální zásobník se může zmenšit při dopředném úspěšném výpočtu deterministické procedury. ■ Unifikace argumentů hlavy - obecný unifikační algoritmus Současně poznačí adresy instanciovaných proměnných na stopu. ■ „Interpret": interpret(Query, Vars) :- caVI(Query), success(Query, Vars). interpret(_,_) :- failure. ■ dotaz vsazen do kontextu této speciální nedeterministické procedury ■ tato procedura odpovídá za korektní reakci systému v případě úspěchu i neúspěchu Hana Rudová, Logické programování 1,15. května 2013 254 Implementace Prologu Indexace argumentů a(l) :- q Cl). a(a) :- b(X). a([A|T]) :- c(A,T). ■ Obecně nedeterministická ■ Při volání s alespoň částečně instanciovaným argumentem vždy deterministická (pouze jedna klauzule může uspět) ■ Indexace podle prvního argumentu Základní typy zřetězení: ■ podle pořadí klauzulí (aktuální argument je volná proměnná) ■ dle konstant (aktuální je argument konstanta) ■ formální argument je seznam (aktuální argument je seznam) ■ dle struktur (aktuální argument je struktura) Hana Rudová, Logické programování I, 1 5. května 2013 255 Implementace Prologu Hana Rudová, Logické programování I, 15. května 2013 256 Implementace Prologu Indexace argumentů II Tail Recursion Optimization, TRO Složitější indexační techniky ■ podle všech argumentů ■ podle nejvíce diskriminujícího argumentu ■ kombinace argumentů (indexové techniky z databází) ■ zejména pro přístup k faktům Iterace prováděna pomocí rekurze => lineární paměťová náročnost cyklů Optimalizace koncové rekurze (Tail Recursion Optimisation), TRO: Okolí se odstraní před rekurzivním voláním posledního literálu klauzule, pokud je klauzule resp. její volání deterministické. Řízení se nemusí vracet: ■ v případě úspěchu se rovnou pokračuje ■ v případě neúspěchu se vrací na předchozí bod volby („nad" aktuální klauzulí) ■ aktuální klauzule nemá dle předpokladu bod volby Rekurzivně volaná klauzule může být volána přímo z kontextu volající klauzule. Hana Rudová, Logické programování I, 15. května 2013 257 Implementace Prologu TRO - příklad Program: append([], L, L). append([A|X], L, [A|Y]) :- append(X, L, Y). Dotaz: ?- append([a,b,c], [x], L). append volán rekurzivně 4krát ■ bez TRO: 4 okolí, lineární paměťová náročnost ■ s TRO: 1 okolí, konstatní paměťová náročnost Hana Rudová, Logické programování I, 15. května 2013 259 Implementace Prologu Hana Rudová, Logické programování 1,15. května 2013 2 58 Implementace Prologu Optimalizace posledního volání TRO pouze speciální případ obecné optimalizace posledního volání (Last Call Optimization), LCO Okolí (před redukcí posledního literálu) odstraňováno vždy, když leží na vrcholu zásobníku. Nutné úpravy interpretu ■ disciplina směrování ukazatelů ■ vždy „mladší" ukazuje na „starší" („mladší" budou odstraněny dříve) ■ z lokálního do globálního zásobníku vyhneme se vzniku „visících odkazů" při předčasném odstranění okolí ■ „globalizace" lokálních proměnných: lokální proměnné posledního literálu ■ nutno přesunout na globální zásobník ■ pouze pro neinstanciované proměnné Hana Rudová, Logické programování 1,15. května 2013 260 Implementace Prologu Warrenův abstraktní počítač, WAM I. Navržen D.H.D. Warrenem v roce 1 983, modifikace do druhé poloviny 80. let Datové oblasti: ■ Oblast kódu (programová databáze) ■ separátní oblasti pro uživatelský kód (modifikovatelný) a vestavěné predikáty (nemění se) ■ obsahuje rovněž všechny statické objekty (texty atomů a funktorú apod.) ■ Lokální zásobník (Stack) ■ Stopa (Trail) ■ Globální zásobník n. ba\da(Heap) ■ Pomocný zásobník (Push Down List, PDL) ■ pracovní paměť abstraktního počítače ■ použitý v unifikaci, syntaktické analýze apod. Hana Rudová, Logické programování I, 15. května 2013 261 Implementace Prologu Registry WAM u ■ Stavové registry: P čitač adres (Program counter) CP adresa návratu (Continuation Pointer) E ukazatel na nejmladší okolí (Environment) B ukazatel na nejmladší bod volby (Backtrack point) TR vrchol stopy (TRail) H vrchol haldy (Heap) HB vrchol haldy v okamžiku založení posledního bodu volby (Heap on Backtrack point) S ukazatel, používaný při analýze složených termů (Structure pointer) CUT ukazatel na bod volby, na který se řezem zařízne zásobník ■ Argumentové registry: A1,A2 , . . . (při předávání parametrů n. pracovní registry) ■ Registry pro lokální proměnné: Y1.Y2, . . . ■ abstraktní znázornění lok. proměnných na zásobníku Hana Rudová, Logické programování 1,15. května 2013 263 Implementace Prologu Rozmístění datových oblastí ■ Příklad konfigurace Halda Stopa \ Zásobník , PDL > \ Oblast kodu ■ Halda i lokální zásobník musí růst stejným směrem ■ lze jednoduše porovnat stáří dvou proměnných srovnáním adres využívá se při zabránění vzniku visících odkazů Hana Rudová, Logické programování 1,15. května 2013 262 Implementace Prologu Typy instrukcí WAM u ■ put instrukce - příprava argumentů před voláním podcíle ■ žádná z těchto instrukcí nevolá obecný unifikační algoritmus ■ get instrukce - unifikace aktuálních a formálních parametrů ■ vykonávají činnost analogickou instrukcím unify ■ obecná unifikace pouze při get_va"lue ■ unify instrukce - zpracování složených termů ■ jednoargumentové instrukce, používají registr S jako druhý argument ■ počáteční hodnota S je odkaz na 1. argument ■ volání instrukce unify zvětší hodnotu S o jedničku ■ obecná unifikace pouze při unify_va"lue a unify_1oca"l_va"lue ■ Indexační instrukce - indexace klauzulí a manipulace s body volby ■ Instrukce řízení běhu - předávání řízení a explicitní manipulace s okolím Hana Rudová, Logické programování 1,15. května 2013 264 Implementace Prologu Instrukce put a get: příklad WAM - optimalizace Příklad: a(X,Y,Z) :- get_var A1.A5 get_var A2.A6 get_var A3.A7 put_const Al,f put_value A2.A5 put_value A3.A6 put_value A4.A7 execute b/4 b(f,X,Y,Z). 1. Indexace klauzulí 2. Generování optimální posloupnosti instrukcí WAMu 3. Odstranění redundancí při generování cílového kódu. ■ Příklad: a(X,Y,Z) :- b(f,X,Y,Z). naivní kód (vytvoří kompilátor pracující striktně zleva doprava) vs. optimalizovaný kód (počet registrů a tedy i počet instrukcí/přesunů v paměti snížen): get_var A1.A5 1 get_var A3, A4 get_var A2.A6 1 get_var A2.A3 get_var A3.A7 1 get_var A1.A2 put_const Al,f 1 put_const Al,f put_value A2.A5 1 execute b/4 put_value A3.A6 1 put_value A4.A7 1 execute b/4 1 Hana Rudová, Logické programování I, 1 5. května 201 3 Implementace Prologu Hana Rudová, Logické programování 1,15. května 201 3 Implementace Prologu Instrukce WAMu get instrukce put instrukce unify instrukce get_var Ai,Y put_var Ai ,Y unify_var Y get_value Ai ,Y put_value Ai ,Y unify_value Y get_const Ai,C put_unsafe_value Ai,Y unify_local_value Y get_nil Ai put_const Ai ,C unify_const C get_struct Ai,F/N put_ni1 Ai u n i fy_n i1 get_list Ai put_struct Ai,F/N unify_void N put_li st Ai instrukce řízení indexační instrukce allocate try_me_else Next try Next deallocate retry_me_else Next retry Next call Proc/N,A trust_me_else fail trust fail execute Proc/N proceed cut_last switch_on_ term Var,Const,List,Struct save_cut Y switch_on_ const Table load_cut Y switch_on_ struct Table Hana Rudová, Logické programování I, 15. května 2013 267 Implementace Prologu WAM - indexace ■ Provázání klauzulí: instrukce XX_me_el se: ■ první klauzule: try_me_else; založí bod volby ■ poslední klauzule: trust_me_else; zruší nejmladší bod volby ■ ostatní klauzule: retry_me_el se; znovu použije nejmladší bod volby po neúspěchu ■ Provázání podmnožiny klauzulí (podle argumentu): ■ try ■ retry ■ trust ■ „Rozskokové" instrukce (dle typu a hodnoty argumentu): ■ switch_on_term Var, Const, List, Struct výpočet následuje uvedeným návěstím podle typu prvního argumentu ■ switch_on_YY: hashovací tabulka pro konkrétní typ (konstanta, struktura) Hana Rudová, Logické programování 1,15. května 2013 268 Implementace Prologu Příklad indexace instrukcí Proceduře a(atom) :- bodyl. a(l) :- body2. a(2) :- body3. odpovídají instrukce a([X|Y]) a([X|Y]) aCs(N)) : a(fCN)) : - body4. - body5. body6. bodyľ. a: switch_on_term LI, L2, L3, L4 L5a: body2 L2: switch_on_const atom :Lla L6: retry. _me_ _el se L7 1 :L5a L6a: body3 2 :L6a L7: retry. _me_ _el se L8 l_3: try L7a L7a: body4 trust L8a L8: retry. _me_ _el se L9 L4: switch_on_struct s/l :L9a L8a: body5 f/1 :LlOa L9: retry. _me_ _el se L10 LI: try_me_else L5 L9a: body6 Lla: bodyl L10: trust. _me_ _el se fail L5: retry_me_else L6 LlOa body7 Hana Rudová, Logické programování 1,15. května 2013 269 WAM - řízení výpočtu execute Proč: ekvivalentní příkazu goto proceed: zpracování faktů allocate: alokuje okolí (pro některé klauzule netřeba, proto explicitně generováno) deallocate: uvolní okolí (je-li to možné, tedy leží-li na vrcholu zásobníku) cal 1 Proč, N: zavolá Proč, N udává počet lok. proměnných (odpovídá velikosti zásobníku) Možná optimalizace: vhodným uspořádáním proměnných lze dosáhnout postupného zkracování lokálního zásobníku a(A,B,C,D) :- b(D), c(A,C), d(B), e(A), f. generujeme instrukce allocate call b/1,4 call c/2,3 call d/1,2 call e/1,1 deallocate execute f/0 Implementace Prologu Hana Rudová, Logické programování 1,15. května 201 3 Implementace Prologu WAM - řez Implementace řezu (opakování): odstranění bodů volby mezi současným vrcholem zásobníku a bodem volby procedury, která řez vyvolala (včetně bodu volby procedury s řezem) Indexační instrukce znemožňují v době překladu rozhodnout, zda bude alokován bod volby ■ příklad: ?- a(X) . může být nedeterministické, ale ?- a(l) . může být deterministické cut_last: B := CUT save_cut Y: Y := CUT load_cut Y: B := Y Hodnota registru B je uchovávána v registru CUT instrukcemi cal 1 a execute. Je-li řez prvním predikátem klauzule, použije se rovnou cut_last. V opačném případě se použije jako první instrukce save_cut Y a v místě skutečného volání řezu se použije load_cut Y. Příklad: a(X,Z) :- b(X), !, c(Z). a(2,Z) :- !, c(Z). a(X,Z) :- d(X,Z). odpovídá save_cut Y2; get A2.Y1; call b/1,2; load_cut Y2; put Yl.Al; execute c/1 get_const AI,2; cut_last; put A2.A1; execute c/1 execute d/2 Hana Rudová, Logické programování I, 15. května 2013 271 Implementace Prologu