|
|
A100763
|
|
a(1) = 2; for n>1, a(n) is the least prime with all prime digits which is == 1 (mod a(n-1)).
|
|
0
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
|
|
LINKS
|
|
|
EXAMPLE
|
a(4) = 337, a(5) = 32353 == 1 (mod 337).
|
|
MATHEMATICA
|
a[1] = 2; a[n_] := a[n] = Block[{k = 1, p = a[n - 1]}, While[ !PrimeQ[k*p + 1] || Union[ Join[{2, 3, 5, 7}, IntegerDigits[k*p + 1]]] != {2, 3, 5, 7}, k++ ]; k*p + 1]; Table[ a[n], {n, 6}] (* Robert G. Wilson v, Nov 27 2004 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
base,nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|