Advertisement
Guest User

Untitled

a guest
Mar 26th, 2019
85
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.62 KB | None | 0 0
  1. AbsoluteTiming[
  2. f[n_, a_, b_, c_, k_] := Module[{x = IntegerDigits[n]},
  3. If[{a, k} != {1, 0},
  4. x[[a]] = x[[b]] = x[[c]] = k;
  5. FromDigits[x], 0]];
  6. g[n_, a_, b_, c_] := Length[Select[Flatten[
  7. Table[f[n, a, b, c, k], {k, 0, 9}]], PrimeQ]];
  8. h[n_] := Table[g[n, a, b, c], {a, 1, IntegerLength[n] - 3},
  9. {b, a + 1, IntegerLength[n] - 2}, {c, b + 1, IntegerLength[n] - 1}];
  10. nmin = SelectFirst[
  11. Table[Prime[k], {k, PrimePi[56003] + 1, PrimePi[10^6]}],
  12. Max[h[#]] == 8 &];
  13. lst = Accumulate@Flatten[Position[h[nmin], 8]];
  14.  
  15. ans = Min[
  16. Select[Table[f[nmin, lst[[1]], lst[[2]], lst[[3]], k], {k, 0, 9}],
  17. PrimeQ]]]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement