欧啦计划第七题10001st prime

happy886rr · 发表于 2016-4-7 23:19:34

回复 16# codegay
julia的计算能力更卓越，但不知为何你启动就要3秒，有那么费时吗。SSD 4k读写几十M应该是秒开的。

codegay · 发表于 2016-4-7 23:22:54

回复 16# happy886rr

windows版的问题。
linux版启动0.2秒

happy886rr · 发表于 2016-4-7 23:31:29

回复 17# codegay
另外你的解法是生成小于等于end的素数数组，但是你没有对end进行估值。
可以用素数定理逆推出end大体范围，从而减少不必要的素数生成。因为计算上亿级别的素数根本不是几秒的事。但如果用递增，就不用考虑这些。

codegay · 发表于 2016-4-7 23:38:16

回复 18# happy886rr

好的。但是我并不懂算法。
julia现在的primes() 函数还不是生成器，直接返回一个数组。所以会有大量内存消耗的问题。

happy886rr · 发表于 2016-4-10 23:36:05

本帖最后由 happy886rr 于 2016-4-11 09:34 编辑

三种不同算法，最后一个基于最好的数学理论.

复制代码

复制代码

# Python解EulerPJ-7P：The 10001st prime number is？
# 构建小于1000/3的素数数组p，只用小于333.3...的素数去除，理论上讲是最优解法，但是python数组效率没救了，居然用了0.7秒才出结果。脑汁已用尽！
N=1;i=1;p=[]
p.append(3)
while i<10001:
N+=2;v=1
for j in range(1,len(p)):
if N%p[j] is 0:
v=0
if v is 1:
i+=1
if N<333:
p.append(N)
print(N)

复制代码

codegay · 发表于 2016-4-11 02:52:34

25th是什么鬼。

当年我入门python可是用了半个月。。。

数组长了，操作会变很慢的。

happy886rr · 发表于 2016-4-11 09:36:20

回复 21# codegay
欧拉项目题多，搞错题号，代码已修改，又增加一个最优算法，见原楼。

codegay · 发表于 2016-4-11 16:17:26

回复 22# happy886rr

python很多对象都是可迭代的。不用len range index
for r in list\str\文件IO 都可以。

网上说的字典和集合都是hash表存储的。缺点是默认不保持顺序。

		自动登录	找回密码
密码			注册

[新手上路]批处理新手入门导读	[视频教程]批处理基础视频教程	[视频教程]VBS基础视频教程	[批处理精品]批处理版照片整理器
[批处理精品]纯批处理备份&还原驱动	[批处理精品]CMD命令50条不能说的秘密	[在线下载]第三方命令行工具	[在线帮助]VBScript / JScript 在线参考

[原创代码] 欧啦计划第七题10001st prime