6 use CGI::Carp qw( fatalsToBrowser );
14 require '../config.pm';
15 require '../common.pm';
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;
26 binmode STDOUT, ':utf8';
30 my %ratings_stddev = ();
36 return HTML::Entities::encode_entities(shift);
41 return int(255.0 * ($x ** (1.80)));
46 my $ref = $dbh->selectrow_hashref('SELECT MAX(sesong) AS max_sesong FROM fotballserier');
47 return $ref->{'max_sesong'};
51 my ($dbh, $season) = @_;
55 my $q = $dbh->prepare('SELECT DISTINCT(divisjon) FROM fotballserier WHERE sesong=? ORDER BY divisjon');
58 while (my $ref = $q->fetchrow_hashref) {
59 push @divisions, $ref->{'divisjon'};
65 sub get_subdivisions {
66 my ($dbh, $season, $division) = @_;
68 my @subdivisions = ();
70 my $q = $dbh->prepare('SELECT DISTINCT(avdeling) FROM fotballserier WHERE sesong=? AND divisjon=? ORDER BY avdeling');
71 $q->execute($season, $division);
73 while (my $ref = $q->fetchrow_hashref) {
74 push @subdivisions, $ref->{'avdeling'};
80 sub print_division_selector {
81 my ($dbh, $divisions, $subdivisions, $division, $subdivision) = @_;
84 <form method="get" action="/">
87 my $max_division = $divisions->[(scalar @$divisions) - 1];
91 <select name="divisjon" onchange="form.submit();">
94 for my $d (@$divisions) {
95 if ($d == $division) {
96 print " <option value=\"$d\" selected=\"selected\">$d</option>\n";
98 print " <option value=\"$d\">$d</option>\n";
105 <select name="avdeling" onchange="form.submit();">
108 for my $sd (@$subdivisions) {
109 if ($sd == $subdivision) {
110 print " <option value=\"$sd\" selected=\"selected\">$sd</option>\n";
112 print " <option value=\"$sd\">$sd</option>\n";
118 <input type="submit" value="Vis" />
124 sub get_players_and_ratings {
125 my ($dbh, $season, $division, $subdivision) = @_;
127 my $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=?');
128 $q->execute($season, $division, $subdivision);
130 while (my $ref = $q->fetchrow_hashref) {
131 my $id = $ref->{'id'};
132 $players{$id} = sanitize(Encode::decode_utf8($ref->{'navn'}));
133 $ratings{$id} = $ref->{'rating'};
134 $ratings_stddev{$id} = $ref->{'rating_stddev'};
140 my ($dbh, $season, $division, $subdivision) = @_;
143 my $q = $dbh->prepare('
145 d1.id AS p1, d2.id AS p2, maalfor AS score1, maalmot AS score2
146 FROM fotballresultater r
147 JOIN fotballserier s ON r.serie=s.nr
148 JOIN fotballdeltagere d1 ON r.lagrecno=d1.nr AND r.serie=d1.serie
149 JOIN fotballdeltagere d2 ON r.motstander=d2.nr AND r.serie=d2.serie
151 sesong=? AND divisjon=? AND avdeling=?
152 AND lagrecno > motstander
154 $q->execute($season, $division, $subdivision);
156 while (my $ref = $q->fetchrow_hashref) {
157 push @matches, [ $ref->{'p1'}, $ref->{'p2'}, $ref->{'score1'}, $ref->{'score2'} ];
164 sub get_covariance_matrix {
165 my ($dbh, @players) = @_;
167 my $player_sql = '{' . join(',', @players ) . '}';
168 my $q = $dbh->prepare('SELECT * FROM covariance WHERE player1=ANY(?::smallint[]) AND player2=ANY(?::smallint[])', { pg_prepare_now => 0 });
169 $q->execute($player_sql, $player_sql);
172 while (my $ref = $q->fetchrow_hashref) {
173 $cov->{$ref->{'player1'}}{$ref->{'player2'}} = $ref->{'cov'};
179 sub write_parms_to_file {
180 my ($used_ratings, $used_cov) = @_;
182 POSIX::setlocale(&POSIX::LC_ALL, 'C');
184 my $tmpnam = POSIX::tmpnam();
185 open MCCALC, ">", $tmpnam
186 or die "$tmpnam: $!";
188 printf MCCALC "%f\n", $match_stddev;
189 printf MCCALC "%d\n", scalar keys %players;
191 for my $id (keys %players) {
192 my $rating = $used_ratings->{$id} // 500.0;
193 printf MCCALC "%s %f\n", $id, $rating;
197 for my $id1 (keys %players) {
198 for my $id2 (keys %players) {
200 printf MCCALC "%f ", ($used_cov->{$id1}{$id2} // $parms{-3});
202 printf MCCALC "%f ", ($used_cov->{$id1}{$id2} // 0.0);
208 for my $match (@matches) {
209 printf MCCALC "%s %s %d %d\n", $match->[0], $match->[1], $match->[2], $match->[3];
213 POSIX::setlocale(&POSIX::LC_ALL, 'nb_NO.UTF-8');
219 my ($lowest_division, $used_ratings, $used_cov) = @_;
228 my $tmpnam = write_parms_to_file($used_ratings, $used_cov);
231 open MCCALC, "$config::base_dir/mcwordfeud $trials < $tmpnam |"
237 my $player = sprintf "%s (%.0f ± %.0f)", $players{$id}, ($ratings{$id} // 500.0), ($ratings_stddev{$id} // $parms{-3});
238 $prob{$player} = [ @x[1..$#x] ];
243 my $num_games = scalar keys %prob;
244 for my $i (1..$num_games) {
245 print " <th>$i.</th>\n";
247 print " <th>NEDRYKK</th>\n" unless ($lowest_division);
250 for my $player (sort { $a cmp $b } keys %prob) {
252 print " <th>$player</th>\n";
254 for my $i (1..$num_games) {
255 my $pn = $prob{$player}->[$i - 1] / $trials;
257 my $r = color(1.0 - $pn / 3);
258 my $g = color(1.0 - $pn / 3);
263 } elsif ($i >= $num_games - 1 && !$lowest_division) {
267 printf " <td style=\"background-color: rgb($r, $g, $b)\" class=\"num\">%.1f%%</td>\n", $pn * 100.0;
270 unless ($lowest_division) {
271 my $pn = ($prob{$player}->[$num_games - 1] + $prob{$player}->[$num_games - 2]) / $trials;
274 my $g = color(1.0 - $pn / 3);
275 my $b = color(1.0 - $pn / 3);
276 printf " <td style=\"background-color: rgb($r, $g, $b)\" class=\"num\">%.1f%%</td>\n", $pn * 100.0;
286 sub find_avg_rating {
287 my ($ratings) = shift;
289 my $sum_rating = 0.0;
290 for my $r (values %$ratings) {
293 return $sum_rating / scalar keys %ratings;
296 # Get auxillary parameters
297 my $q = $dbh->prepare('SELECT * FROM ratings WHERE id < 0');
300 while (my $ref = $q->fetchrow_hashref) {
301 $parms{$ref->{'id'}} = $ref->{'rating'};
303 $match_stddev = $parms{-2} * sqrt(2.0);
306 my $division = $cgi->param('divisjon') // -1;
307 my $subdivision = $cgi->param('avdeling') // -1;
308 my $last_division = 0;
310 POSIX::setlocale(&POSIX::LC_ALL, 'nb_NO.UTF-8');
312 print $cgi->header(-type=>'text/html; charset=utf-8', -expires=>'now');
313 printf <<"EOF", $match_stddev;
314 <?xml version="1.0" encoding="UTF-8" ?>
316 html PUBLIC "-//W3C//DTD XHTML 1.1//EN"
317 "http://www.w3.org/TR/xhtml11/DTD/xhtml11.dtd">
318 <html xmlns="http://www.w3.org/1999/xhtml" xml:lang="no">
320 <title>WLoH-plasseringsannsynlighetsberegning</title>
321 <link rel="stylesheet" href="/style" type="text/css" />
324 <h1>WLoH-plasseringsannsynlighetsberegning</h1>
326 <p><em>Dette er et hobbyprosjekt fra tredjepart, og ikke en offisiell del av
327 <a href="http://wordfeud.aasmul.net/">Wordfeud Leage of Honour</a>.</em></p>
329 <p>Beregningen tar ikke hensyn til ujevn spillestyrke, ting som er sagt i forumet e.l.;
330 den antar at samtlige uspilte kamper trekkes fra en normalfordeling med standardavvik
331 %.1f poeng. Sannsynlighetene kan summere til andre tall enn 100%% pga. avrunding.
332 Tallene vil variere litt fra gang til gang fordi utregningen skjer ved randomisering.</p>
334 <p>Spillerne er sortert etter nick.</p>
337 my $season = get_max_season($dbh);
338 my @divisions = get_divisions($dbh, $season);
339 $division = 1 if (!grep { $_ == $division } @divisions);
340 my @subdivisions = get_subdivisions($dbh, $season, $division);
341 $subdivision = 1 if (!grep { $_ == $subdivision } @subdivisions);
343 print_division_selector($dbh, \@divisions, \@subdivisions, $division, $subdivision);
345 get_players_and_ratings($dbh, $season, $division, $subdivision);
347 my @matches = get_matches($dbh, $season, $division, $subdivision);
349 my $cov = get_covariance_matrix($dbh, keys %players);
351 my $max_division = $divisions[$#divisions];
352 my $lowest_division = ($division == $max_division);
353 make_table($lowest_division, {}, {});
356 <p>Under er en variant som tar relativ spillestyrke med i beregningen;
357 se <a href="/rating">ratingsiden</a>.</p>
360 make_table($lowest_division, \%ratings, $cov);
362 my $avg_rating = find_avg_rating(\%ratings);
363 printf " <p>Gjennomsnittlig rating i denne avdelingen er <strong>%.1f</strong>.</p>\n", $avg_rating;
365 wloh_common::output_last_sync($dbh);