>>11 外側の正方形 Perl5 凸包の辺が正方形の辺に接するまたは対角二点で接する場合、>>135の露骨なﾊﾞｸﾞ一個修正
use List::Util qw{min max pairkeys pairvalues};
@s=qw{136 577 110 927 472 199 157 808 388 598 94 31 388 157 325 409 787 897 850 598};
@X = pairkeys @s; @Y = pairvalues @s;
sub sp {$_[0]*$_[2] + $_[1]*$_[3]}
sub rt {
　@e = (cos $th, -sin $th); @f = (sin $th, cos $th);
　my @x = map{sp @e, $X[$_], $Y[$_]} 0..$#X;
　my @y = map{sp @f, $X[$_], $Y[$_]} 0..$#Y;
　@x = (min(@x), max(@x)); @y = (min(@y), max(@y));
　$h = (max $x[1] - $x[0], $y[1] - $y[0]) / 2;
　$cx = ($x[1] + $x[0]) / 2; $cy = ($y[1] + $y[0]) / 2;
　@x = ($cx - $h, $cx + $h, $cx - $h, $cx + $h);
　@y = ($cy - $h, $cy - $h, $cy + $h, $cy + $h);
　($e[1], $f[0]) = (-$e[1], -$f[0]);
　@x = map{sp @e, $x[$_], $y[$_]} 0..3;
　@y = map{sp @f, $x[$_], $y[$_]} 0..3;
　(\@x, \@y, 2*$h) }
for $i (0..@X-2) { for $j ($i+1..$#X) {
　($dx, $dy) = ($X[$j] - $X[$i], $Y[$j] - $Y[$i]);
　($dx, $dy) = (-$dx, -$dy) if $dx < 0;
　$l = sqrt($dx*$dx + $dy*$dy);
　$th = $dx > abs($dy) ? -atan2($dy, $dx) : atan2($dx, $dy);
　($X, $Y, $w) = &rt;
　push @t, +{i,$i,j,$j,dx,$dx,dy,$dy,l,$l,th,$th,X,$X,Y,$Y,w,$w};
　$th += 3.14159265358979/4; ($X, $Y, $w) = &rt;
　push @t, +{i,$i,j,$j,dx,$dx,dy,$dy,l,$l,th,$th,X,$X,Y,$Y,w,$w};
} }
@t = sort{$$a{w}<=>$$b{w}} grep{0<=min@{$_->{X}}and max@{$_->{X}}<=999 and 0<=min@{$_->{Y}}and max@{$_->{Y}}<=999} @t;
do {@x = @{$t[$_]{X}}; @y = @{$t[$_]{Y}};
　　printf"%d: (%7.3f,%7.3f), (%7.3f,%7.3f), (%7.3f,%7.3f), (%7.3f,%7.3f): w=%7.3f, th=%7.3f°\n",
　　　$_+1,$x[0],$y[0],$x[1],$y[1],$x[2],$y[2],$x[3],$y[3],$t[$_]{w},$t[$_]{th}*180/3.14159265358979} for 0..4;