最近在解 Project Euler 写了一段code
出现了一些error message
个人觉得这可能是Mathematica的bug
fexp[{n_Integer}] = 1;
fexp[exp_List] := fexp[exp] = Multinomial @@ exp;
f[i_] := fexp[Sort[FactorInteger[i]][[All, 2]]];
code蛮单纯的，就不解释了
执行
Sum[f[i], {i, 1000000}] // AbsoluteTiming
几秒钟后就计算完毕
但如果执行
Sum[f[i], {i, 1000001}] // AbsoluteTiming
开始出现错误讯息
FactorInteger::exact: Argument i in FactorInteger[i] is not an exact number.
Part::partd: "Part specification FactorInteger[i][[All,2]] is longer than
depth of object."
但事实上 单独执行f[1000001]是正常的
神奇的是 如果使用平行运算又没问题
ParallelSum[f[i], {i, 1000001}] // AbsoluteTiming
而且答案在Project Euler上是正确的
所以code应该是没错
使用Mathematica 10.2 跟 10.3 都是如此
想请教各位高手的是：
1) 这究竟是否是Bug，或其实是coding上有问题
2) 如果真是Bug，要怎么避开这个Bug
3) 当我丢平行运算的时候，各thread并无法共用内存，导致DP无法真正发挥效用
应该要如何设定才能够让各thread可以共用呢？