From a02ee715c1b231e5dca96dbf633db3fdb25c0d17 Mon Sep 17 00:00:00 2001
From: Aaron Schaefer <aaron@elasticdog.com>
Date: Thu, 17 Jan 2008 12:55:02 -0500
Subject: [PATCH] Optimize solution to Project Euler problem 27

---
 extra/project-euler/027/027.factor | 9 +++++++--
 1 file changed, 7 insertions(+), 2 deletions(-)

diff --git a/extra/project-euler/027/027.factor b/extra/project-euler/027/027.factor
index abfe4d5704..3ce684549a 100644
--- a/extra/project-euler/027/027.factor
+++ b/extra/project-euler/027/027.factor
@@ -36,10 +36,15 @@ IN: project-euler.027
 ! SOLUTION
 ! --------
 
+! b must be prime since n = 0 must return a prime
+! a + b + 1 must be prime since n = 1 must return a prime
+! a < b
+
 <PRIVATE
 
 : source-027 ( -- seq )
-    1000 [ prime? ] subset dup [ neg ] map append dup cartesian-product ;
+    1000 [ prime? ] subset [ dup [ neg ] map append ] keep
+    cartesian-product [ first2 < ] subset ;
 
 : quadratic ( b a n -- m )
     dup sq -rot * + + ;
@@ -60,7 +65,7 @@ PRIVATE>
     source-027 max-consecutive drop product ;
 
 ! [ euler027 ] 100 ave-time
-! 1306 ms run / 58 ms GC ave time - 100 trials
+! 687 ms run / 23 ms GC ave time - 100 trials
 
 ! TODO: generalize max-consecutive/max-product (from #26) into a new word