Advertisement
Guest User

Untitled

a guest
Mar 29th, 2017
53
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.42 KB | None | 0 0
  1. In[1]:= EP146[maxn_] := Module[{a},
  2. HanteiQ[n_] :=
  3. Block[{m = n^2},
  4. PrimeQ[m + 1] && PrimeQ[m + 3] &&
  5. PrimeQ[m + 7] &&
  6. PrimeQ[m + 9] && PrimeQ[m + 13] &&
  7. NextPrime[m + 1, 5] == m + 27];
  8. SetSharedVariable[a];
  9. a = {};
  10. Parallelize[Scan[
  11. If[HanteiQ[#], AppendTo[a, #]] &, Range[10, maxn, 10]]];
  12. UnsetShared[a];
  13. Total[a]];
  14.  
  15. In[2]:= EP146[150000000] // Timing
  16.  
  17. Out[2]= {5.1562, 676333270}
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement