git.sesse.net Git - wloh/blob - www/index.pl

   1 #! /usr/bin/perl
   2 use strict;
   3 use warnings;
   4 no warnings qw(once);
   5 use CGI;
   6 use CGI::Carp qw( fatalsToBrowser );
   7 use DBI;
   8 use POSIX;
   9 use Devel::Peek;
  10 use HTML::Entities;
  11 use Encode;
  12 use utf8;
  13 use locale;
  14 require '../config.pm';
  15 require '../common.pm';
  16
  17 my $cgi = CGI->new;
  18
  19 my $dbh = DBI->connect($config::local_connstr, $config::local_username, $config::local_password)
  20         or die "connect: " . $DBI::errstr;
  21 $dbh->{AutoCommit} = 0;
  22 $dbh->{RaiseError} = 1;
  23
  24 my $trials = 25_000;
  25
  26 binmode STDOUT, ':utf8';
  27
  28 my %players = ();
  29 my %ratings = ();
  30 my %ratings_stddev = ();
  31 my @matches = ();
  32 my %parms = ();
  33 my $match_stddev;
  34
  35 sub sanitize {
  36         return HTML::Entities::encode_entities(shift);
  37 }
  38
  39 sub color {
  40         my $x = shift;
  41         return int(255.0 * ($x ** (1.80)));
  42 }
  43
  44 sub get_max_season {
  45         my $dbh = shift;
  46         my $ref = $dbh->selectrow_hashref('SELECT MAX(sesong) AS max_sesong FROM fotballserier');
  47         return $ref->{'max_sesong'};
  48 }
  49
  50 sub get_divisions {
  51         my ($dbh, $season) = @_;
  52
  53         my @divisions = ();
  54
  55         my $q = $dbh->prepare('SELECT DISTINCT(divisjon) FROM fotballserier WHERE sesong=? ORDER BY divisjon');
  56         $q->execute($season);
  57
  58         while (my $ref = $q->fetchrow_hashref) {
  59                 push @divisions, $ref->{'divisjon'};
  60         }
  61
  62         return @divisions;
  63 }
  64
  65 sub get_subdivisions {
  66         my ($dbh, $season, $division) = @_;
  67
  68         my @subdivisions = ();
  69
  70         my $q = $dbh->prepare('SELECT DISTINCT(avdeling) FROM fotballserier WHERE sesong=? AND divisjon=? ORDER BY avdeling');
  71         $q->execute($season, $division);
  72
  73         while (my $ref = $q->fetchrow_hashref) {
  74                 push @subdivisions, $ref->{'avdeling'};
  75         }
  76
  77         return @subdivisions;
  78 }
  79
  80 sub print_division_selector {
  81         my ($dbh, $divisions, $subdivisions, $division, $subdivision) = @_;
  82
  83         print <<"EOF";
  84     <form method="get" action="/">
  85 EOF
  86
  87         my $max_division = $divisions->[(scalar @$divisions) - 1];
  88
  89         print <<"EOF";
  90      <p>Divisjon:
  91         <select name="divisjon" onchange="form.submit();">
  92 EOF
  93
  94         for my $d (@$divisions) {
  95                 if ($d == $division) {
  96                         print "        <option value=\"$d\" selected=\"selected\">$d</option>\n";
  97                 } else {
  98                         print "        <option value=\"$d\">$d</option>\n";
  99                 }
 100         }
 101
 102         print <<"EOF";
 103         </select>
 104         Avdeling:
 105         <select name="avdeling" onchange="form.submit();">
 106 EOF
 107
 108         for my $sd (@$subdivisions) {
 109                 if ($sd == $subdivision) {
 110                         print "        <option value=\"$sd\" selected=\"selected\">$sd</option>\n";
 111                 } else {
 112                         print "        <option value=\"$sd\">$sd</option>\n";
 113                 }
 114         }
 115
 116         print <<"EOF";
 117         </select>
 118         <input type="submit" value="Vis" />
 119       </p>
 120     </form>
 121 EOF
 122 }
 123
 124 sub get_covariance_matrix {
 125         my ($dbh, @players) = @_;
 126
 127         my $player_sql = '{' . join(',', @players ) . '}';
 128         my $q = $dbh->prepare('SELECT * FROM covariance WHERE player1=ANY(?::smallint[]) AND player2=ANY(?::smallint[])', { pg_prepare_now => 0 });
 129         $q->execute($player_sql, $player_sql);
 130
 131         my $cov = {};
 132         while (my $ref = $q->fetchrow_hashref) {
 133                 $cov->{$ref->{'player1'}}{$ref->{'player2'}} = $ref->{'cov'};
 134         }
 135
 136         return $cov;
 137 }
 138
 139 sub write_parms_to_file {
 140         my ($used_ratings, $used_cov) = @_;
 141
 142         POSIX::setlocale(&POSIX::LC_ALL, 'C');
 143
 144         my $tmpnam = POSIX::tmpnam();
 145         open MCCALC, ">", $tmpnam
 146                 or die "$tmpnam: $!";
 147
 148         printf MCCALC "%f\n", $match_stddev;
 149         printf MCCALC "%d\n", scalar keys %players;
 150
 151         for my $id (keys %players) {
 152                 my $rating = $used_ratings->{$id} // 500.0;
 153                 printf MCCALC "%s %f\n", $id, $rating;
 154         }
 155
 156         # covariance matrix
 157         for my $id1 (keys %players) {
 158                 for my $id2 (keys %players) {
 159                         if ($id1 == $id2) {
 160                                 printf MCCALC "%f ", ($used_cov->{$id1}{$id2} // $parms{-3});
 161                         } else {
 162                                 printf MCCALC "%f ", ($used_cov->{$id1}{$id2} // 0.0);
 163                         }
 164                 }
 165                 printf MCCALC "\n";
 166         }
 167
 168         for my $match (@matches) {
 169                 printf MCCALC "%s %s %d %d\n", $match->[0], $match->[1], $match->[2], $match->[3];
 170         }
 171         close MCCALC;
 172
 173         POSIX::setlocale(&POSIX::LC_ALL, 'nb_NO.UTF-8');
 174
 175         return $tmpnam;
 176 }
 177
 178 sub make_table {
 179         my ($lowest_division, $used_ratings, $used_cov) = @_;
 180
 181         print <<"EOF";
 182
 183     <table>
 184       <tr>
 185         <th></th>
 186 EOF
 187
 188         my $tmpnam = write_parms_to_file($used_ratings, $used_cov);
 189         my %prob = ();
 190
 191         open MCCALC, "$config::base_dir/mcwordfeud $trials < $tmpnam |"
 192                 or die "mccalc: $!";
 193         while (<MCCALC>) {
 194                 chomp;
 195                 my @x = split /\s+/;
 196                 my $id = $x[0];
 197                 my $player = sprintf "%s (%.0f ± %.0f)", $players{$id}, ($ratings{$id} // 500.0), ($ratings_stddev{$id} // $parms{-3});
 198                 $prob{$player} = [ @x[1..$#x] ];
 199         }
 200         close MCCALC;
 201         #unlink $tmpnam;
 202
 203         my $num_games = scalar keys %prob;
 204         for my $i (1..$num_games) {
 205                 print "        <th>$i.</th>\n";
 206         }
 207         print "        <th>NEDRYKK</th>\n" unless ($lowest_division);
 208         print "      </tr>\n";
 209
 210         for my $player (sort { $a cmp $b } keys %prob) {
 211                 print "      <tr>\n";
 212                 print "        <th>$player</th>\n";
 213
 214                 for my $i (1..$num_games) {
 215                         my $pn = $prob{$player}->[$i - 1] / $trials;
 216
 217                         my $r = color(1.0 - $pn / 3);
 218                         my $g = color(1.0 - $pn / 3);
 219                         my $b = color(1.0);
 220
 221                         if ($i == 1) {
 222                                 ($g, $b) = ($b, $g);
 223                         } elsif ($i >= $num_games - 1 && !$lowest_division) {
 224                                 ($r, $b) = ($b, $r);
 225                         }
 226
 227                         printf "        <td style=\"background-color: rgb($r, $g, $b)\" class=\"num\">%.1f%%</td>\n", $pn * 100.0;
 228                 }
 229
 230                 unless ($lowest_division) {
 231                         my $pn = ($prob{$player}->[$num_games - 1] + $prob{$player}->[$num_games - 2]) / $trials;
 232
 233                         my $r = color(1.0);
 234                         my $g = color(1.0 - $pn / 3);
 235                         my $b = color(1.0 - $pn / 3);
 236                         printf "        <td style=\"background-color: rgb($r, $g, $b)\" class=\"num\">%.1f%%</td>\n", $pn * 100.0;
 237                 }
 238                 print "      </tr>\n";
 239         }
 240
 241         print << "EOF";
 242     </table>
 243 EOF
 244 }
 245
 246 # Get auxillary parameters
 247 my $q = $dbh->prepare('SELECT * FROM ratings WHERE id < 0');
 248 $q->execute;
 249
 250 while (my $ref = $q->fetchrow_hashref) {
 251         $parms{$ref->{'id'}} = $ref->{'rating'};
 252 }
 253 $match_stddev = $parms{-2} * sqrt(2.0);
 254
 255 my $season;
 256 my $division = $cgi->param('divisjon') // -1;
 257 my $subdivision = $cgi->param('avdeling') // -1;
 258 my $last_division = 0;
 259
 260 POSIX::setlocale(&POSIX::LC_ALL, 'nb_NO.UTF-8');
 261
 262 print $cgi->header(-type=>'text/html; charset=utf-8', -expires=>'now');
 263 printf <<"EOF", $match_stddev;
 264 <?xml version="1.0" encoding="UTF-8" ?>
 265 <!DOCTYPE
 266   html PUBLIC "-//W3C//DTD XHTML 1.1//EN"
 267   "http://www.w3.org/TR/xhtml11/DTD/xhtml11.dtd">
 268 <html xmlns="http://www.w3.org/1999/xhtml" xml:lang="no">
 269   <head>
 270     <title>WLoH-plasseringsannsynlighetsberegning</title>
 271     <link rel="stylesheet" href="/style" type="text/css" />
 272   </head>
 273   <body>
 274     <h1>WLoH-plasseringsannsynlighetsberegning</h1>
 275
 276     <p><em>Dette er et hobbyprosjekt fra tredjepart, og ikke en offisiell del av
 277       <a href="http://wordfeud.aasmul.net/">Wordfeud Leage of Honour</a>.</em></p>
 278
 279     <p>Beregningen tar ikke hensyn til ujevn spillestyrke, ting som er sagt i forumet e.l.;
 280       den antar at samtlige uspilte kamper trekkes fra en normalfordeling med standardavvik
 281       %.1f poeng. Sannsynlighetene kan summere til andre tall enn 100%% pga. avrunding.
 282       Tallene vil variere litt fra gang til gang fordi utregningen skjer ved randomisering.</p>
 283
 284     <p>Spillerne er sortert etter nick.</p>
 285 EOF
 286
 287 my $season = get_max_season($dbh);
 288 my @divisions = get_divisions($dbh, $season);
 289 $division = 1 if (!grep { $_ == $division } @divisions);
 290 my @subdivisions = get_subdivisions($dbh, $season, $division);
 291 $subdivision = 1 if (!grep { $_ == $subdivision } @subdivisions);
 292
 293 print_division_selector($dbh, \@divisions, \@subdivisions, $division, $subdivision);
 294
 295 # Get players and ratings
 296 my $sum_rating = 0.0;
 297
 298 $q = $dbh->prepare('SELECT fotballdeltagere.id,fotballdeltagere.navn,rating,rating_stddev FROM fotballdeltagere JOIN fotballserier ON fotballdeltagere.serie=fotballserier.nr LEFT JOIN ratings ON fotballdeltagere.id=ratings.id WHERE sesong=? AND divisjon=? AND avdeling=?');
 299 $q->execute($season, $division, $subdivision);
 300
 301 while (my $ref = $q->fetchrow_hashref) {
 302         my $id = $ref->{'id'};
 303         $players{$id} = sanitize(Encode::decode_utf8($ref->{'navn'}));
 304         $ratings{$id} = $ref->{'rating'};
 305         $ratings_stddev{$id} = $ref->{'rating_stddev'};
 306         $sum_rating += $ref->{'rating'};
 307 }
 308 $q->finish;
 309
 310 $q = $dbh->prepare('
 311 SELECT
 312   d1.id AS p1, d2.id AS p2, maalfor AS score1, maalmot AS score2
 313 FROM fotballresultater r
 314   JOIN fotballserier s ON r.serie=s.nr
 315   JOIN fotballdeltagere d1 ON r.lagrecno=d1.nr AND r.serie=d1.serie
 316   JOIN fotballdeltagere d2 ON r.motstander=d2.nr AND r.serie=d2.serie
 317 WHERE
 318   sesong=? AND divisjon=? AND avdeling=?
 319   AND lagrecno > motstander
 320 ');
 321 $q->execute($season, $division, $subdivision);
 322
 323 while (my $ref = $q->fetchrow_hashref) {
 324         push @matches, [ $ref->{'p1'}, $ref->{'p2'}, $ref->{'score1'}, $ref->{'score2'} ];
 325 }
 326 $q->finish;
 327
 328 my $cov = get_covariance_matrix($dbh, keys %players);
 329
 330 my $max_division = $divisions[$#divisions];
 331 my $lowest_division = ($division == $max_division);
 332 make_table($lowest_division, {}, {});
 333
 334 print <<"EOF";
 335     <p>Under er en variant som tar relativ spillestyrke med i beregningen;
 336       se <a href="/rating">ratingsiden</a>.</p>
 337 EOF
 338
 339 make_table($lowest_division, \%ratings, $cov);
 340
 341 my $avg_rating = $sum_rating / scalar keys %players;
 342 printf "    <p>Gjennomsnittlig rating i denne avdelingen er <strong>%.1f</strong>.</p>\n", $avg_rating;
 343
 344 wloh_common::output_last_sync($dbh);
 345
 346 print <<"EOF";
 347   </body>
 348 </html>
 349 EOF