OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..422 (* First 118 terms from Paolo P. Lava *)
EXAMPLE
311683 is prime and also 3116873, 3116783, 3117683, 3171683 and 3711683.
MAPLE
with(numtheory);
A217065:=proc(q, x)
local a, b, c, i, n, ok;
for n from 5 to q do
a:=ithprime(n); b:=0; while a>0 do b:=b+1; a:=trunc(a/10); od; a:=ithprime(n); ok:=1;
for i from 1 to b-1 do
c:=a+9*10^i*trunc(a/10^i)+10^i*x; if not isprime(c) then ok:=0; break; fi; od;
if ok=1 then print(ithprime(n)); fi; od; end:
A217065(1000000, 7);
MATHEMATICA
Select[Prime[Range[5, 1600]], AllTrue[FromDigits/@Table[Insert[ IntegerDigits[ #], 7, i], {i, 2, IntegerLength[#]}], PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Oct 12 2016 *)
PROG
(PARI) is(n)=my(v=concat([""], digits(n))); for(i=2, #v-1, v[1]=Str(v[1], v[i]); v[i]=7; if(i>2, v[i-1]=""); if(!isprime(eval(concat(v))), return(0))); isprime(n) \\ Charles R Greathouse IV, Sep 26 2012
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Paolo P. Lava, Sep 26 2012
STATUS
approved