共有28种方法表一千为两素数之和

青青子衿

Select[IntegerPartitions[1000, 2], And @@ PrimeQ@# &] // Length
IntegerPartitions[1000, 2, Prime[Range[PrimePi[1000]]]] // Length

睡神

找到一个效率更高的代码
pr[n_]:=DeleteCases[If[PrimeQ[#] && PrimeQ[n-#], {#,n-#},0] & /@Range[3, IntegerPart[n/2], 2], 0]
m = 1000000;
Length[pr[m]]
运行时间：约1秒半钟
而楼上的运行时间需要约16分钟

hbghlyj

1. pr[n_]:=Select[Range[3,Floor[n/2],2],PrimeQ[#]&&PrimeQ[n-#]&];
2. m=1000000;
3. Length[pr[m]]

Copy the Code

Out[]= 5402

睡神

hbghlyj 发表于 2024-3-13 17:42
Compile

https://reference.wolfram.com/language/ref/ParallelMap.html

是不是哪里出错了？

hbghlyj

睡神 发表于 2024-3-13 13:17
是不是哪里出错了？

抱歉，我的意思是用ParallelMap替换/@，我的代码错误，已更正

睡神

hbghlyj 发表于 2024-3-13 21:58
抱歉，我的意思是用ParallelMap替换/@，我的代码错误，已更正

😂😂😂
一个钟了，还在运算，我只能重启电脑了...

hbghlyj

hbghlyj 发表于 2024-3-13 13:58
抱歉，我的意思是用ParallelMap替换/@，我的代码错误，已更正

文档说 ParallelMap 比 Map 更快。

1. pr[n_]:=DeleteCases[ParallelMap[If[PrimeQ[#]&&PrimeQ[n-#],{#,n-#},0]&,Range[3,IntegerPart[n/2],2]],0]
2. m=1000000;
3. Length[pr[m]]

Copy the Code

睡神

hbghlyj 发表于 2024-3-13 17:42
Out[]= 5402

为啥第一个我运行时是这样子的？电脑问题，还是版本问题？而第二个运算效率不算高，用时3分多钟

睡神

hbghlyj 发表于 2024-3-13 22:07
文档说 ParallelMap 比 Map 更快。

是不是我的版本问题？

hbghlyj

睡神 发表于 2024-3-13 14:18
是不是我的版本问题？

有没有可能pr没清除
Clear[pr]

睡神

hbghlyj 发表于 2024-3-13 22:43
有没有可能pr没清除
Clear[pr]

对的，没Clean[pr]。代码如下：
Clear[pr]
pr[n_] := DeleteCases[ParallelMap[If[PrimeQ[#] && PrimeQ[n - #], {#, n - #}, 0] &, Range[3, IntegerPart[n/2], 2]], 0]
m = 10000000;
Length[pr[m]]
约3.45s出结果，我的那个约5.45s出结果，效率提升了很多

hbghlyj

睡神 发表于 2024-3-13 15:26
对的，没Clean[pr]。

是不是拼写错误，好像是Clear[]而不是Clean[]