[rosalind] RSTR

Raku · Dec 7, 2012 · a19c36f · a19c36f
1 parent 85f27f2
commit a19c36f
Showing 1 changed file with 14 additions and 0 deletions.
diff --git a/rosalind/rstr-grondilu.pl b/rosalind/rstr-grondilu.pl
@@ -0,0 +1,14 @@
+use v6;
+
+my @data = $*IN.lines;
+my ($N, $gc-content) = @data.shift.split: " ";
+my $dna = @data.shift;
+
+sub prob(:$dna, :$gc-content) {
+    1/2**$dna.chars * [*] map { $_ eq 'G'|'C' ?? $gc-content !! (1 - $gc-content) }, $dna.comb
+}
+
+printf "%.3f",  1 - exp ($N * log(1-prob :$dna, :$gc-content));
+
+# vim: ft=perl6
+