GnomSort

Aus ProgrammingWiki

(define listequal? (lambda (ls1 ls2) (cond [(and (null? ls1) (null? ls2)) #t] [(or (null? ls1) (null? ls2)) #f] [(and (list? (car ls1)) (list? (car ls1))) (and (listequal? (car ls1) (car ls2)) (listequal? (cdr ls1) (cdr ls2)))] [else (and (equal? (car ls1) (car ls2)) (listequal? (cdr ls1) (cdr ls2)))]))) (define rand 50) (define maxbreite (- 400 (* 2 rand))) (define maxhoehe (- 300 (* 2 rand))) (define minmax (lambda (ls) (letrec ([extremsuche (lambda (extrem ls op) (cond [(null? ls) extrem] [(op extrem (car ls)) (extremsuche (car ls) (cdr ls) op)] [else (extremsuche extrem (cdr ls) op)]))]) (cons (extremsuche 0 ls >) (extremsuche 0 ls <))))) (define balken (lambda (hoehe breite turtle) (forward turtle breite) (left turtle 90) (forward turtle hoehe) (left turtle 90) (forward turtle breite) (left turtle 90) (forward turtle hoehe) (left turtle 90) (forward turtle breite))) (define balkendiagramm (lambda (ls nr) (let* ([breite (/ maxbreite (length ls))] [extremwerte (minmax ls)] [min (car extremwerte)] [max (cdr extremwerte)] [canvas (eval (string->symbol (string-append "canvas" (number->string nr))))] [turtle (eval (string->symbol (string-append "turtle" (number->string nr))))]) (letrec ([zeichnen (lambda (min max ls) (if (null? ls) 'ok (begin (balken (floor (* (/ (car ls) (- max min)) maxhoehe)) (floor breite) turtle) (zeichnen min max (cdr ls)))))]) (clear canvas) (hide turtle) (move turtle rand (+ rand (* (/ max (- max min)) maxhoehe))) (right turtle 90) (zeichnen min max ls)))))

Inhaltsverzeichnis 1 Algorithmus 2 Effizienz von GnomSort 2.1 Average-Case-Experimente 2.2 Best- und Worst-Case-Experimente

Algorithmus

In einem Garten sind Bumentöpfe in einer Reihe aufgestellt. Gartenzwerge sollen sie nun der Größe nach aufsteigend sortieren. Sie benutzen dazu den folgenden Algorithmus:

• Beginne links und vergleiche die beiden ersten Töpfe.
• Stehen sie in der richtigen Reihenfolge gehe ein Schritt nach rechts.
• Vertausche die beiden Töpfe, wenn sie nicht in der korrekten Reihenfolge stehen, und gehe einen Schritt nach links.
  Ist der Schritt nach links nicht möglich, gehe ein Schritt nach rechts.
• Wiederhole das Sortieren bis zum letzten Blumentopf.

Die Größe der Blumentöpfe wird mit Zufallszahlen festgelegt:

(define zufallszahl (lambda (min max) (+ (randomint (- max min -1)) min))) (define zufallsliste (lambda (min max anzahl) (if (= anzahl 0) () (cons (zufallszahl min max) (zufallsliste min max (- anzahl 1))))))

Da keine Indizierung der Listenelemente möglich ist, wird der Algorithmus mit den beiden Listen ls1 und ls2 implementiert. Die Liste ls1 enthält die unsortierte Zahlenliste. Mit ls2 wird die sortierte Zahlenlise gebildet. Beim Schritt nach rechts wird das erste Listenelement von ls1 der Liste ls2 vorangestellt. Analog wird beim Schritt nach links das erste Listenelement von ls2 zurückgenommen und ls1 vorangestellt:

;GnomSort (define gnomsort (lambda (ls) (letrec ([sort (lambda (ls1 ls2) ...]) (reverse (sort ls ())))))

Ungeordnete Zufallszahlen:

Geordnete Zufallszahlen:

(define ls (zufallsliste 1 1000 100)) (balkendiagramm ls 1) (define ls-sort (gnomsort ls)) (balkendiagramm ls-sort 2)

ausführen

Hinweis: Um mehrere Zeichenflächen nutzen zu können, muss die Prozedur balkendiagramm mit einem zusätzlichen Parameter ausgestattet werden, der einer fortlaufenden Nummerierung entspricht:

(balkendiagramm
  <ls>       ; Liste mit Zahlenwerten   
  <nr>)      ; fortlaufende Nummerierung des Balkendiagramms

Effizienz von GnomSort

Average-Case-Experimente

Wir definieren unseren GnomSort-Algorithmus mit einem "eingebauten" (Zeittakt-)Zähler:

(define count 0) (define gnomsort-aufwand (lambda (ls) (letrec ([sort (lambda (ls1 ls2) (cond [(null? ls1) ls2] [(null? (cdr ls1)) (set! count (+ 1 count)) ;gib die einelementige Liste zurück (cons (car ls1) ls2)] [(> (car ls1) (cadr ls1)) (set! count (+ 1 count)) ;vertausche das n-te Paar und gehe nach links (bzw. nach rechts) (if (null? ls2) (sort (cons (car ls1) (cddr ls1)) (cons (cadr ls1) ls2)) (sort (cons (car ls2) (cons (cadr ls1) (cons (car ls1) (cddr ls1)))) (cdr ls2)))] [else (set! count (+ 1 count)) ;gehe nach rechts (sort (cdr ls1) (cons (car ls1) ls2))]))]) (reverse (sort ls ()))))) (define gnomsort-ausgabe (lambda (min max anzahl) (let ([ls (zufallsliste min max anzahl)]) (set! count 0) (display "Ungeordnete Zufallszahlen:") (newline) (display ls) (newline) (display "Geordnete Zufallszahlen:") (newline) (display (gnomsort-aufwand ls)) (newline) (display "Anzahl der Vertauschungen bzw. Bewegungen: ") (display count) (newline) 'ok)))

(gnomsort-ausgabe 1 1000 20)

ausführen

Im Weiteren wollen wir uns nur auf die Zählergebnisse unseres "Zeittaktmessers" konzentrieren.
Dazu sollen die Wiederholungen des Experimentes automatisiert werden:

(define gnomsort-mehrfach (lambda (min max anzahl whlg) (letrec ([summe (lambda (n) (if (= n 0) 0 (begin (set! count 0) (gnomsort-aufwand (zufallsliste min max anzahl)) (+ count (summe (- n 1))))))]) (exact->inexact (/ (summe whlg) whlg)))))

(gnomsort-mehrfach 1 1000 20 100)

ausführen

Schließlich sollen die Ergebnisse unserer empirischen Untersuchung wieder im Balkendiagramm visualisiert werden.

(define experiment (lambda (min max anzahl schrittweite whlg) (if (< anzahl 1) () (append (experiment min max (- anzahl schrittweite) schrittweite whlg) (list (gnomsort-mehrfach min max anzahl whlg))))))

Achtung! Mit kleinen Wiederholraten und etwas Geduld experimentieren!

(balkendiagramm (experiment 1 1000 200 5 10) 3)

ausführen

Die empirischen Untersuchungen zeigen im Mittel offensichtlich einen quadratischen Zeitaufwand von GnomSort.
Es gilt:

average case:	bzw.

Best- und Worst-Case-Experimente

Zur Untersuchung dieser Spezialfälle wollen wir die beiden Prozeduren gnomsort-mehrfach und experiment variabelstellig erweitern. Zusätzlich können nun beliebig viele Prozeduren übergeben, die sich auf Listen anwenden lassen. Damit kann beispielsweise die zufällige Zahlenliste "vorsortiert" oder vor dem Sortiervorgang ungekehrt werden:

(define gnomsort-spezial (lambda (min max anzahl whlg . proz-ls) (letrec ([anwenden (lambda (zuf-ls proz-ls) (if (null? proz-ls) zuf-ls (apply (car proz-ls) (list (anwenden zuf-ls (cdr proz-ls))))))] [summe (lambda (n) (if (= n 0) 0 (begin (set! count 0) (gnomsort-aufwand (anwenden (zufallsliste min max anzahl) proz-ls)) (+ count (summe (- n 1))))))]) (exact->inexact (/ (summe whlg) whlg))))) (define experiment-spezial (lambda (min max anzahl schrittweite whlg . proz-ls) (if (< anzahl 1) () (append (apply experiment-spezial min max (- anzahl schrittweite) schrittweite whlg proz-ls) (list (apply gnomsort-spezial min max anzahl whlg proz-ls))))))

GnomSort: best-case

GnomSort: worst-case

(balkendiagramm (experiment-spezial 1 1000 200 5 1 gnomsort) 4) (balkendiagramm (experiment-spezial 1 1000 200 5 5 reverse gnomsort) 5)

ausführen

Für die beiden Spezialfälle gilt:

best case:	bzw.
worst case:	bzw.

Die Effizienz von GnomSort verbessert sich also, wenn schon eine gewisse Ordnung ("Vorsortierung") vorliegt.
Im ungünstigsten Fall ist GnomSort etwas uneffizienter als MinSort.