OFFSET
1,1
COMMENTS
The only unknown terms less than 10000, tested to 17500, are for n: 484, 1291, 2096, 2238, 3503, 3859, 6674, 7087, 7824, 8954.
LINKS
Vladimir Shevelev and Robert G. Wilson v, a(n) for n = 1..10000 with -1 for those entries where a(n) has not yet been found.
FORMULA
a(n)=k for the least k such that prime(n)*10^k+7*(10^k-1)/9 is prime, where prime(n) is the n-th prime.
MATHEMATICA
f[n_] := Block[{k = 1, p = Prime[n]}, While[ !PrimeQ[p*10^k + 7(10^k - 1)/9], k++]; k]; f[4] = 0; Array[f, 100]
PROG
(PARI) isok(k, dp) = ispseudoprime(fromdigits(concat(dp, vector(k, i, 7))));
a(n) = {if (prime(n) == 7, return(0)); my(k=1, p=prime(n)); while (!ispseudoprime(p*10^k+7*(10^k-1)/9), k++); k; } \\ Michel Marcus, Jan 20 2021
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Vladimir Shevelev and Robert G. Wilson v, Apr 24 2015
STATUS
approved