OFFSET
1,2
COMMENTS
a(10) > 10^10. - Robert Price, Apr 02 2019
EXAMPLE
a(4)=33 because 331, 661, 991 and 1321 (1321=10*4*33+1) are all prime, but 1651 (1651=10*5*33+1) is not prime. - Robert Price, Apr 02 2019
MATHEMATICA
Table[k = 1;
While[! AllTrue[Table[10*r*k + 1, {r, 1, n}], PrimeQ] ||
PrimeQ[10*(n + 1)*k + 1], k++]; k, {n, 1, 9}] (* Robert Price, Apr 02 2019 *)
PROG
(PARI) isok(k, n) = {for (r=1, n, if (! isprime(10*r*k+1), return (0)); ); !isprime(10*(n+1)*k+1); }
a(n) = {my(k=1); while(! isok(k, n), k++); k; } \\ Michel Marcus, Apr 03 2019
CROSSREFS
KEYWORD
more,nonn
AUTHOR
Ray G. Opao, Jan 13 2005
EXTENSIONS
Corrected a(7) and added a(9) by Robert Price, Apr 02 2019
STATUS
approved