@x = values %ref_freq;

print"@x \n";

close(IN);


#my @x= map $_ ,$key;
#エントロピー計算
my $y=$count;
my $num;
my @pi = map $_ / $y, @x;
my $i;
print"@pi \n"; #テスト用プリント
#print"@x \n";
for($i=0;$i<=$#pi;$i++){
$num += -($pi[$i]*log($pi[$i]) / log(2));
}


#print"@pi\n\n";
print "-Σpi*log2(pi) = $num \n";