From 37c4f77fe7c41603b9410728ccd53b9021749717 Mon Sep 17 00:00:00 2001 From: Jesse Luehrs Date: Mon, 18 May 2009 00:04:45 -0500 Subject: problem 27 solution --- 027.pl | 34 ++++++++++++++++++++++++++++++++++ 1 file changed, 34 insertions(+) create mode 100755 027.pl diff --git a/027.pl b/027.pl new file mode 100755 index 0000000..b1f8808 --- /dev/null +++ b/027.pl @@ -0,0 +1,34 @@ +#!/usr/bin/perl +use strict; +use warnings; +use Math::Prime::XS qw/sieve_primes/; + +my @primes = sieve_primes(1e6); +my %primes = map { $_, 1 } @primes; + +sub is_prime { + my ($n) = @_; + if ($n <= 1e6) { + return exists $primes{$n}; + } + return Math::Prime::XS::is_prime($n); +} + +my $max = 0; +my $max_a = 0; +my $max_b = 0; +for my $a (-1000..1000) { + for my $b (-1000..1000) { + for my $n (0..1e6) { + if (!is_prime($n**2+$a*$n+$b)) { + if ($n - 1 > $max) { + $max = $n - 1; + $max_a = $a; + $max_b = $b; + } + last; + } + } + } +} +print "$max $max_a $max_b ", $max_a * $max_b, "\n"; -- cgit v1.2.3