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

TOP

回复 16# happy886rr

windows版的问题。
linux版启动0.2秒
去学去写去用才有进步。安装python3代码存为xx.py 双击运行或右键用IDLE打开按F5运行

TOP

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

TOP

回复 18# happy886rr


    好的。但是我并不懂算法。
julia现在的primes() 函数还不是生成器,直接返回一个数组。所以会有大量内存消耗的问题。
去学去写去用才有进步。安装python3代码存为xx.py 双击运行或右键用IDLE打开按F5运行

TOP

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

三种不同算法,最后一个基于最好的数学理论.
#  Python解EulerPJ-7P:The 10001st prime number is?
# 常规解法用小于根号N的奇数去除
N=1;i=1;
while i<10001:
N+=2;v=1
for j in range(3,int(N**0.5)+1,2):
if N%j is 0:
v=0;
if v is 1:
i+=1
print(N)COPY
# Python解EulerPJ-7P:The 10001st prime number is?
# 构建素数数组p,只用素数数组中的数去除,按理说速度会提高,自己语法不太会,速度没提高很囧。
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 p[j]**2>N:
break
if v is 1:
i+=1;p.append(N)
print(N)COPY
# 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)COPY

TOP

25th是什么鬼。

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

数组长了,操作会变很慢的。
去学去写去用才有进步。安装python3代码存为xx.py 双击运行或右键用IDLE打开按F5运行

TOP

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

TOP

回复 22# happy886rr


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

网上说的字典和集合都是hash表存储的。缺点是默认不保持顺序。
去学去写去用才有进步。安装python3代码存为xx.py 双击运行或右键用IDLE打开按F5运行

TOP

返回列表