Curzon numbers: Difference between revisions
Content added Content deleted
Thundergnat (talk | contribs) m (syntax highlighting fixup automation) |
No edit summary |
||
Line 469: | Line 469: | ||
One thousandth: 46845 |
One thousandth: 46845 |
||
</pre> |
|||
=={{header|Haskell}}== |
|||
<syntaxhighlight lang=Haskell> |
|||
import Data.List.Split ( chunksOf ) |
|||
isGeneralizedCurzon :: Integer -> Integer -> Bool |
|||
isGeneralizedCurzon base n = mod ( base ^ n + 1 ) ( base * n + 1 ) == 0 |
|||
solution :: Integer -> [Integer] |
|||
solution base = take 50 $ filter (\i -> isGeneralizedCurzon base i ) [1..] |
|||
printChunk :: [Integer] -> String |
|||
printChunk chunk = foldl1 (++) $ map (\i -> (take ( 4 - (length $ show i) ) |
|||
$ repeat ' ' ) ++ show i ++ " ") chunk |
|||
prettyPrint :: [Integer] -> [String] |
|||
prettyPrint list = map printChunk $ chunksOf 10 list |
|||
oneThousandth :: Integer -> Integer |
|||
oneThousandth base = last $ take 950 $ filter (\i -> isGeneralizedCurzon base i ) |
|||
[(last $ solution base) + 1 ..] |
|||
printBlock :: Integer -> [String] |
|||
printBlock base = ["first 50 Curzon numbers using a base of " ++ show base ++ " :"] |
|||
++ (prettyPrint $ solution base) ++ ["one thousandth at base " ++ show base ++ |
|||
": " ++ (show $ oneThousandth base)] ++ [take 50 $ repeat '-'] |
|||
main :: IO ( ) |
|||
main = do |
|||
blocks <- return $ concat $ map (\i -> printBlock i ) [2 , 4 , 6 , 8 , 10] |
|||
mapM_ putStrLn blocks</syntaxhighlight> |
|||
{{out}} |
|||
<pre>first 50 Curzon numbers using a base of 2 : |
|||
1 2 5 6 9 14 18 21 26 29 |
|||
30 33 41 50 53 54 65 69 74 78 |
|||
81 86 89 90 98 105 113 114 125 134 |
|||
138 141 146 153 158 165 173 174 186 189 |
|||
194 198 209 210 221 230 233 245 249 254 |
|||
one thousandth at base 2: 8646 |
|||
-------------------------------------------------- |
|||
first 50 Curzon numbers using a base of 4 : |
|||
1 3 7 9 13 15 25 27 37 39 |
|||
43 45 49 57 67 69 73 79 87 93 |
|||
97 99 105 115 127 135 139 153 163 165 |
|||
169 175 177 183 189 193 199 205 207 213 |
|||
219 235 249 253 255 265 267 273 277 279 |
|||
one thousandth at base 4: 9375 |
|||
-------------------------------------------------- |
|||
first 50 Curzon numbers using a base of 6 : |
|||
1 6 30 58 70 73 90 101 105 121 |
|||
125 146 153 166 170 181 182 185 210 233 |
|||
241 242 266 282 290 322 373 381 385 390 |
|||
397 441 445 446 450 453 530 557 562 585 |
|||
593 601 602 605 606 621 646 653 670 685 |
|||
one thousandth at base 6: 20717 |
|||
-------------------------------------------------- |
|||
first 50 Curzon numbers using a base of 8 : |
|||
1 14 35 44 72 74 77 129 131 137 |
|||
144 149 150 185 200 219 236 266 284 285 |
|||
299 309 336 357 381 386 390 392 402 414 |
|||
420 441 455 459 470 479 500 519 527 536 |
|||
557 582 600 602 617 639 654 674 696 735 |
|||
one thousandth at base 8: 22176 |
|||
-------------------------------------------------- |
|||
first 50 Curzon numbers using a base of 10 : |
|||
1 9 10 25 106 145 190 193 238 253 |
|||
306 318 349 385 402 462 486 526 610 649 |
|||
658 678 733 762 810 990 994 1033 1077 1125 |
|||
1126 1141 1149 1230 1405 1422 1441 1485 1509 1510 |
|||
1513 1606 1614 1630 1665 1681 1690 1702 1785 1837 |
|||
one thousandth at base 10: 46845 |
|||
-------------------------------------------------- |
|||
</pre> |
</pre> |
||