OFFSET
1,1
LINKS
Michel Lagneau, Michael De Vlieger and Robert G. Wilson v, Table of n, a(n) for n = 1..1268
MAPLE
f:= proc(n)
local r, i, j, p;
r:= 7*(10^n-1)/9;
for p in sort([r, seq(seq(r + i*10^j, i=[$(-7)..(-1), 1, 2]), j=0..n-1)], `>`) do
if isprime(p) then return p fi
od;
error("no prime found")
end proc:
map(f, [$1..100]); # Robert Israel, Feb 24 2016
MATHEMATICA
f[n_] := Block[{k = 0, p = {}, r = 7 (10^n - 1)/9, s = Range@ 10 - 8}, While[k < n, AppendTo[p, Select[r + 10^k*s, PrimeQ]]; k++]; p = Max@ Flatten@ p]; Array[f, 20]
CROSSREFS
KEYWORD
nonn,base
AUTHOR
STATUS
approved