Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- In[1]:= EP146[maxn_] := Module[{a},
- HanteiQ[n_] :=
- Block[{m = n^2},
- PrimeQ[m + 1] && PrimeQ[m + 3] &&
- PrimeQ[m + 7] &&
- PrimeQ[m + 9] && PrimeQ[m + 13] &&
- NextPrime[m + 1, 5] == m + 27];
- SetSharedVariable[a];
- a = {};
- Parallelize[Scan[
- If[HanteiQ[#], AppendTo[a, #]] &, Range[10, maxn, 10]]];
- UnsetShared[a];
- Total[a]];
- In[2]:= EP146[150000000] // Timing
- Out[2]= {5.1562, 676333270}
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement