串珠素数问题

无心人

import Primes readInt :: IO Int readInt = readLn getBlock b 0 = [10*a + b | a <- [0..10^(b-1)-1], b <- [1,3,7,9]] getBlock b 1 = [100*a + 10*b + c | a <- [0..10^(b-2)-1], b <- [1..9], c <-[1,3,7,9]] getBlock b p = [10^(p+1)*a + 10^p*b + 10*c + d | a <- [0..10^(b-1-p)-1], b <- [1..9], c <- [0..10^(p-1)-1], d <- [1,3,7,9]] multiBlock b n p = (p \`mod\` 10^(n \`mod\` b)) * 10^(b * (n \`div\` b))+ fromIntegral(sum (take (n \`div\` b) (iterate (10^b*) p))) isczPrime b n = [(ns, pp) | ns <- [b+1..n], ns \`mod\` b /= 0, let pps = map (multiBlock b ns) (getBlock b ((n \`mod\` b)-1)), pp <- pps, isPrime pp ] main = do putStrLn "Input Digits: " n <- readInt putStrLn "Input Block Size: " b <- readInt print (isczPrime b n)

gxqcn

这是 C++ 代码， 请先安装 HugeCalc 算法库，而后编译运行， 可以运行在任何 32 位或 64 位 windows 下。

无心人

21#程序改好了 不知道谁能写个别的语言的 相互印证 下面是编译后的程序 isczPrimen.rar (256.75 KB, 下载次数: 5)

2008-11-24 16:54 上传

点击文件名下载附件

gxqcn

原帖由 无心人 于 2008-11-24 16:51 发表 21#程序改好了 不知道谁能写个别的语言的 相互印证 下面是编译后的程序 595

下载了你的程序，可惜结果一闪而过，无法保留。

无心人

呵呵 没设定保存到文件 不过，很容易写保存程序 就增加一句 你重定向到文件吧 isczPrimen > czPrimes.txt 呵呵

gxqcn

我在 WinXP cmd 下，无论用 > 还是 >> 到文件，都不能成功。 因为无法输入参数了，得到一个0字节空文件。

mathe

不会的，不过输入数字参数需要注意。 比如程序名为test.exe,程序需要参数10,那么可以采用命令 test.exe 10 1>out 其中1>out表示将标准输出重定向到文件out. 如果省略了1而且10和>之间没有留空格，那么命令行会将数字10当成文件句柄从而没有输入参数了

gxqcn

学了一招。 但用在无心人提供的程序仍然不灵， 可能是他输入输出代码非标准问题。

无心人

可是我这里可以啊

无心人

import Primes readInt :: IO Int readInt = readLn getBlock b 0 = [10*a + b | a <- [0..10^(b-1)-1], b <- [1,3,7,9]] getBlock b 1 = [100*a + 10*b + c | a <- [0..10^(b-2)-1], b <- [1..9], c <-[1,3,7,9]] getBlock b p = [10^(p+1)*a + 10^p*b + 10*c + d | a <- [0..10^(b-1-p)-1], b <- [1..9], c <- [0..10^(p-1)-1], d <- [1,3,7,9]] multiBlock b n p = (p \`mod\` 10^(n \`mod\` b)) * 10^(b * (n \`div\` b))+ fromIntegral(sum (take (n \`div\` b) (iterate (10^b*) p))) isczPrime b n = [(ns, pp) | ns <- [b+1..n], ns \`mod\` b /= 0, let pps = map (multiBlock b ns) (getBlock b ((ns \`mod\` b)-1)), pp <- pps, isPrime pp ] main = do putStrLn "Input Digits: " n <- readInt putStrLn "Input Block Size: " b <- readInt let result = isczPrime b n let fileName = "P" ++ show(n) ++ "_" ++ show(b) ++ ".TXT" writeFile fileName (show result) ============================= 重新编译的，保存到文件 isczPrimen.rar (267.75 KB, 下载次数: 1)

2008-11-25 09:22 上传

点击文件名下载附件

本帖最后由 无心人 于 2008-11-25 09:22 编辑