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

codegay

欧啦计划第七题10001st prime

#="""
julia解欧啦计划第七题10001st prime 
https://projecteuler.net/problem=7
2016年4月6日 05:59:30 codegay
"""=#
result=[]
for r in countfrom(3,2)
if isprime(r) 
append!(result,[r])
if length(result)==10000
println(r)
break
end
end
end
#104743
#[Finished in 4.4s]
复制代码

windows 版julia启动速度慢，3.2-4秒左右。拖慢了时间。
如果对julia熟悉的话，应该是可以写成一行。

codegay

回复 22# happy886rr


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

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

happy886rr

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

codegay

25th是什么鬼。

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

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

happy886rr

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

#  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)
复制代码

# 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)
复制代码

# 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

回复 18# happy886rr


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

happy886rr

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

codegay

回复 16# happy886rr

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

happy886rr

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

codegay

我还得跟你解释
我会gt tcc 。。一行命令把tcc下载下来
会存为xx.bat运行。。。

happy886rr

回复 13# codegay
tcc有个-run参数。我利用这个参数直接当脚本运行。也没发现tcc产生临时文件。tcc的官网也说它直接就支持不编译而运行，估计是读到内存里了。

codegay

乱码问题解决了。存为ANSI编码了。

codegay

回复 11# happy886rr


    我知道啊。
但是本质也是需要编译过的。
TCC就是编译器嘛。

happy886rr

回复 10# codegay
而且批处理代码和C代码融合的很好，相互可以调用和传递参数，基于的第三方C解释器tcc精简版才100多kb，运行非常迅捷。

codegay

回复 9# happy886rr


    我操太好用了！
只是我这里编译后运行，乱码

Problem7 -Project Euler
 * The 10001st prime number is 鈥?04743鈥?
复制代码

你那里正常？
win 7x64