From b058af01f9b352966240ad45788be00f3746ce11 Mon Sep 17 00:00:00 2001
From: Martin Mares <mj@ucw.cz>
Date: Sat, 15 Mar 2008 16:11:59 +0100
Subject: [PATCH] Fixed factorization.

---
 pyth.hs | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/pyth.hs b/pyth.hs
index 82e44fd..a59845b 100644
--- a/pyth.hs
+++ b/pyth.hs
@@ -7,7 +7,8 @@ sqr x = x*x
 -- Find (p,q) such that x=p^2*q and q is square-free
 sqf' :: Int -> Int -> (Int,Int)
 sqf' _ 1 = (1,1)
-sqf' d x | x `mod` (sqr d) == 0 = (d*p1,q1)
+sqf' d x | x < (sqr d)        = (1,x)
+         | x `mod` (sqr d) == 0 = (d*p1,q1)
          | x `mod` d == 0     = (p2,d*q2)
 	 | True               = sqf' (d+1) x
 	 where (p1,q1) = sqf' d (x `div` (sqr d))
-- 
2.39.2