use strict;
use warnings;

my $runs = 3;
my $range = 3;
my @X_info =
    ("1\t2\t3\t4", "1\t3\t6\t8", "1\t4\t9\t16", "1\t5\t12\t32");

# Outer loop: Repeat "$runs" times
for my $run (1 .. $runs) {
    my @tmp;
    my $pi_sum;
    my $L_sum;
    my $differences_sum;
    my $coverage_sum;

    push @tmp, $X_info[rand @X_info] for 1 .. $range;

    my @PAR1 = splice @tmp, 0, 3,;

    for my $k (0 .. @PAR1 - 1) {
        my @PAR1_info = $PAR1[$k];
        my @values = split '\t', $PAR1_info[0];

        $pi_sum          += $values[0];
        $L_sum           += $values[1];
        $differences_sum += $values[2];
        $coverage_sum    += $values[3];
    }

    my $PAR1_diversity =
        (($pi_sum / $L_sum) / ($differences_sum / $coverage_sum));

    printf "%2d: %6.3f\n", $run, $PAR1_diversity;
}