happy886rr

回复 4# codegay
julia在数值计算，绘图方面的优势很大，只是安装包体积太大。希望日后能出精简版。

happy886rr

回复 8# codegay

/*
TinyC script混编
move %0 "%~0.c"&cls&@tcc\tcc.exe -run "%~0.c"&rename "%~0.c" "%~nx0"
*/


#include<stdio.h>
main()
{
	int i,j,sum=0;
	for(i=1;i<=2000000000;i+=2){
		for(j=3;j*j<=i;j+=2){
			if(i%j==0) break;
		}
   		if(j*j>i) sum++;
		if(sum==10001) break;
	}
	printf("Problem7 -Project Euler\n");
	printf(" * The 10001st prime number is “%d”\n",i);
	getchar();
}
复制代码

happy886rr

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

happy886rr

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

happy886rr

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

happy886rr

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

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

happy886rr

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