fork download
  1. ps = drop 2 $ takeWhile (<= 97) $ map head $ iterate (\x-> [ n | n<-x, mod n ( head x) /= 0 ] ) [ 2..]
  2.  
  3.  
  4.  
  5. main = do
  6. mapM_ print $ [ ( n,q,p,u,v, p == 24*n + 1, p==(4*u-1)*q-( div (4*u^2) v))| p<- map (^2) ps, u<-[1],v<-[-1], n<- [ div ( p - 1) 24 ], q<-[ div(p - 4) 3 ] ]
Success #stdin #stdout 0.01s 5448KB
stdin
Standard input is empty
stdout
(1,7,25,1,-1,True,True)
(2,15,49,1,-1,True,True)
(5,39,121,1,-1,True,True)
(7,55,169,1,-1,True,True)
(12,95,289,1,-1,True,True)
(15,119,361,1,-1,True,True)
(22,175,529,1,-1,True,True)
(35,279,841,1,-1,True,True)
(40,319,961,1,-1,True,True)
(57,455,1369,1,-1,True,True)
(70,559,1681,1,-1,True,True)
(77,615,1849,1,-1,True,True)
(92,735,2209,1,-1,True,True)
(117,935,2809,1,-1,True,True)
(145,1159,3481,1,-1,True,True)
(155,1239,3721,1,-1,True,True)
(187,1495,4489,1,-1,True,True)
(210,1679,5041,1,-1,True,True)
(222,1775,5329,1,-1,True,True)
(260,2079,6241,1,-1,True,True)
(287,2295,6889,1,-1,True,True)
(330,2639,7921,1,-1,True,True)
(392,3135,9409,1,-1,True,True)