chiark / gitweb /
wip
[appendix-a6.git] / parse
diff --git a/parse b/parse
index a7c6bc0..13fe2f2 100755 (executable)
--- a/parse
+++ b/parse
@@ -1,50 +1,64 @@
 #!/usr/bin/perl -w
 use strict;
 use utf8;
+use autodie;
 
 use Data::Printer;
+use Graph::Directed;
 
-binmode STDIN, 'encoding(UTF-8)' or die;
-binmode STDOUT, 'encoding(UTF-8)' or die;
-binmode STDERR, 'encoding(UTF-8)' or die;
+binmode STDIN, 'encoding(UTF-8)';
+binmode STDOUT, 'encoding(UTF-8)';
+binmode STDERR, 'encoding(UTF-8)';
 
 our @choices;
 our %choices;
 our @invotes;
+our $defcho;
+our $quorum = 0;
 
 sub addchoice {
-    my $cho = shift @_;
-    $choices{$cho} = { @_, Index => (scalar @choices) };
-    push @choices, $cho;
+    my $choname = shift @_;
+    my $cho = $choices{$choname} = { @_, Index => (scalar @choices) };
+    push @choices, $choname;
+    return $cho;
 }
 
 while (<>) {
     s/\s+$//;
     next if m/^\s*\#/;
     next unless m/\S/;
+    next if m/^\s/;
     if (m/^([A-Z]+)\s*\=\s*(\S.*)$/) {
-       my ($cho, $desc) = ($1,$2);
-       my $smaj;
+       my ($choname, $desc) = ($1,$2);
+       my $cho = addchoice($choname, Desc => $desc);
        if ($desc =~ m/\[(\d+):(\d+)\]/) {
-           $smaj = [$1,$2];
-       } 
-       addchoice($cho, Desc => $desc, Smaj => $smaj);
+           $cho->{Smaj} = [$1,$2];
+       } elsif ($desc =~ m/\[default\]/) {
+           $defcho = $cho;
+       }
     } elsif (m/^V:\s+(\S+)\s+(\S.*)/) {
        push @invotes, [ $1, $2 ];
+    } elsif (m/^quorum = ([0-9.]+)$/) {
+       $quorum = $1+0.0;
     } else {
        die "invalid input";
     }
 }
 
-if (!$choices{FD}) {
-    addchoice('FD', Desc => "Further Discussion");
+$defcho ||= $choices{FD};
+if (!$defcho) {
+    $defcho = addchoice('FD', Desc => "Further Discussion");
 }
+my $defi = $defcho->{Index};
+die "FD has smaj?!" if $defcho->{Smaj};
 
-our $chofd = $choices{FD};
+our @vab;
+# $vab[$ia][$ib] = V(A,B)
+# Actually, it's a list of voters who prefer A to B (A.6(3)(1))
 
-die "FD has smaj?!" if $chofd->{Smaj};
+# Go through the voters and construct V(A,B)
 
-our @vab; # $vab[$ia][$ib] = V(A,B)
+print "\nPreference matrix\n";
 
 foreach my $iv (@invotes) {
     my ($votestr,$voter) = @$iv;
@@ -62,16 +76,191 @@ foreach my $iv (@invotes) {
            }
        }
        foreach my $ia (0..$#vs) {
-           foreach my $ib (0..$#vs) {
+           foreach my $ib ($ia+1..$#vs) {
                my $va = $vs[$ia];
                my $vb = $vs[$ib];
-               if ($va < $vb) { $vab[$ia][$ib]++; }
-               elsif ($vb < $va) { $vab[$ib][$ia]++; }
+               if ($va < $vb) { push @{ $vab[$ia][$ib] }, $voter }
+               elsif ($vb < $va) { push @{ $vab[$ib][$ia] }, $voter }
            }
        }
-    }
+    };
     die "voter $voter $@" if $@;
 }
 
-print p %choices;
-print p @vab;
+our @ch = map { $choices{$_} } @choices;
+
+# Print the counts V(A,B)
+foreach my $iy (-2..$#ch) {
+    foreach my $ix (-2..$#ch) {
+       if ($iy==-1) {
+           if ($ix==-1) {
+               printf "+";
+           } else {
+               printf "------";
+           }
+       } elsif ($ix==-1) {
+           printf "|";
+       } elsif ($ix==-2 && $iy==-2) {
+           printf "V(Y,X)";
+       } elsif ($iy==-2) {
+           printf "%5s ", $choices[$ix];
+       } elsif ($ix==-2) {
+           printf "%5s ", $choices[$iy];
+       } else {
+           my $v = \( $vab[$iy][$ix] );
+           $$v ||= [ ];
+           if (@$$v) {
+               printf "%5d ", (scalar @$$v);
+           } else {
+               printf "%5s ", "";
+           }
+       }
+    }
+    printf "\n";
+}
+
+sub drop ($$) {
+    my ($i,$why) = @_;
+    print "dropping $choices[$i]: $why\n";
+    $ch[$i]{Dropped} = $why;
+}
+
+print "\nQuorum A.6(2) (quorum is $quorum)\n";
+
+foreach my $i (0..$#choices) {
+    next if $ch[$i]{Dropped};
+    next if $i == $defi;
+    my $v = $vab[$i][$defi];
+    next if $v >= $quorum;
+    drop $i, "quorum ($v < $quorum)";
+}
+
+print "\nMajority ratio A.6(3)\n";
+
+foreach my $i (0..$#choices) {
+    next if $ch[$i]{Dropped};
+    next if $i == $defi;
+    my $majr = $ch[$i]{Smaj};
+    $majr ||= [1,1]; # A.6(3)(3)
+    my $vad = scalar @{ $vab[$i][$defi] };
+    my $vda = scalar @{ $vab[$defi][$i] };
+    next if $vad * $majr->[1] > $vda * $majr->[0];
+    drop $i, "majority ratio ($vad * $majr->[1] <= $vda * $majr->[0])";
+}
+
+print "\nDefeats A.6(4)\n";
+
+my $defeats = Graph::Directed->new;
+
+sub chremain () {
+    return grep { !$ch[$_]{Dropped} } (0..$#ch);
+}
+
+foreach my $ia (chremain()) {
+    $defeats->add_vertex($choices[$ia]);
+    foreach my $ib (chremain()) {
+       my $vab = scalar @{ $vab[$ia][$ib] };
+       my $vba = scalar @{ $vab[$ib][$ia] };
+       next unless $vab > $vba;
+       my $diff = $vab - $vba;
+       print "defeat: $choices[$ia] beats $choices[$ib]",
+           " ($vab > $vba: $diff)\n";
+       $defeats->add_edge($choices[$ia],$choices[$ib]);
+    }
+}
+
+sub chvab ($$) {
+    my ($ca,$cb) = @_;
+    my $v = $vab[ $choices{$ca}{Index} ][ $choices{$cb}{Index} ];
+    return scalar @$v;
+}
+
+sub weaker ($$) {
+    # A.6(7)(1)
+    my ($def1,$def2) = @_;
+    my ($ca,$cx) = @$def1;
+    my ($cb,$cy) = @$def1;
+    return 1 if chvab($ca, $cx) < chvab($cb, $cy);
+    return 1 if chvab($ca, $cx) == chvab($cb, $cy)
+       && chvab($cx, $ca) > chvab($cy, $cb);
+    return 0;
+}
+
+our $schwartz;
+
+for (;;) {
+    # loop from A6(5)
+
+    print "defeats graph: $defeats\n";
+
+    print "\nTransitive closure A.6(5)\n";
+
+    my $tdefeats = $defeats->transitive_closure();
+
+    print "closure graph: $tdefeats\n";
+
+    print "\nSchwartz set A.6(6)\n";
+    
+    $schwartz = $defeats->copy();
+
+    foreach my $ia (chremain()) {
+       foreach my $ib (chremain()) {
+           next if $tdefeats->has_edge($choices[$ia],$choices[$ib]);
+           next if !$tdefeats->has_edge($choices[$ib],$choices[$ia]);
+           print "not in Schwartz set: $choices[$ia] because $choices[$ib]\n";
+           $schwartz->delete_vertex($choices[$ia]);
+           last;
+       }
+    }
+
+    print "\nDropping weakest defeats A.6(7)\n";
+
+    our @weakest = ();
+
+    foreach my $edge ($schwartz->edges()) {
+#      print " considering @$edge\n";
+       if (!@weakest) {
+           # no weakest edges yet
+       } elsif (weaker($edge, $weakest[0])) {
+           # this edge is weaker than previous weakest, start new set
+           @weakest = ();
+       } elsif (weaker($weakest[0], $edge)) {
+           # weakest edge is weaker than this one, ignore this one
+           next;
+       } else {
+           # weakest edge is exactly as weak as this one, add this one
+       }
+       push @weakest, $edge;
+    }
+
+    last unless @weakest;
+
+    my $w = $weakest[0];
+    printf "weakest defeats %d > %d\n", 
+        chvab($w->[0], $w->[1]),
+       chvab($w->[1], $w->[0]);
+    foreach my $weakest (@weakest) {
+       my ($ca,$cb) = @$weakest;
+       print "weakest defeat $ca > $cb\n";
+       $defeats->delete_edge($ca,$cb);
+    }
+
+    print "defeats within the Schwartz set, round again\n";
+}
+
+print "no defeats within the Schwartz set\n";
+print "final schwartz set:\n\n";
+
+if ($schwartz->vertices() == 1) {
+    print "WINNER IS:\n";
+} else {
+    print "WINNER IS ONE OF (CASTING VOTE DECIDES):\n";
+}
+
+printf "    %-5s %s\n", $_, $choices{$_}{Desc}
+    foreach ($schwartz->vertices());
+
+print ".\n";
+
+#p %choices;
+#p @vab;