use strict;
use warnings;
use List::Util qw(min max);

my $ranges148 = scan('148Nsorted.bed');
my $ranges162 = scan('162Nsorted.bed');
my $ranges174 = scan('174Nsorted.bed');
my $ranges175 = scan('175Nsorted.bed');


foreach my $chr (keys %$ranges148) {
    foreach my $range148 (@{$ranges148->{$chr}}) {
        foreach my $range162 (@{$ranges162->{$chr}}) {
            foreach my $range174 (@{$ranges174->{$chr}}) {
                foreach my $range175 (@{$ranges175->{$chr}}) {
                    my $minstart = min(
                        $range148->[0],
                        $range162->[0],
                        $range174->[0],
                        $range175->[0]);
                    my $maxstart = max(
                        $range148->[0],
                        $range162->[0],
                        $range174->[0],
                        $range175->[0]);
                    my $minend = min(
                        $range148->[1],
                        $range162->[1],
                        $range174->[1],
                        $range175->[1]);
                    my $maxend = max(
                        $range148->[1],
                        $range162->[1],
                        $range174->[1],
                        $range175->[1]);

                    if($maxstart <= $minend) {
                        print "overlapping $chr: $minstart, $maxend - ($range148->[0], $range148->[1]), ($range162->[0], $range162->[1]), ($range174->[0], $range174->[1]), ($range175->[0], $range175->[1])\n";
                    } else {
                        print "not overlapping: $chr: $minstart, $maxend - ($range148->[0], $range148->[1]), ($range162->[0], $range162->[1]), ($range174->[0], $range174->[1]), ($range175->[0], $range175->[1])\n";
                    }
                }
            }
        }
    }
}

exit(0);


sub scan {
    my ($file) = @_;

    open(my $fh, '<', $file) or die "$file: $!";

    my $ranges;
    foreach my $line (<$fh>) {
        chomp($line);
        my ($chr, $start, $end) = split(/\t/, $line);
        push(@{$ranges->{$chr}}, [ 0+$start, 0+$end ]);
    }
    return($ranges);
}