cp's OEIS Frontend

A108843 Palindromic primes in which all internal digits are 6.

%I A108843 #28 Jan 31 2023 08:29:21
%S A108843 16661,76667,7666667,1666666666661,16666666666666661,
%T A108843 1666666666666666661,1666666666666666666666666666666666661,
%U A108843 16666666666666666666666666666666666666666666666666661
%N A108843 Palindromic primes in which all internal digits are 6.
%C A108843 The external digits must either be ones or sevens. - _Harvey P. Dale_, Apr 06 2019
%C A108843 a(15) has 3385 digits. - _Michael S. Branicky_, Jan 27 2023
%H A108843 Michael S. Branicky, <a href="/A108843/b108843.txt">Table of n, a(n) for n = 1..14</a>
%t A108843 nn=80;With[{o=Table[FromDigits[Join[PadRight[{1},n,6],{1}]],{n,3,nn}], s= Table[ FromDigits[Join[PadRight[{7},n,6],{7}]],{n,3,nn}]}, Select[ Sort[ Join[o,s]],PrimeQ]] (* _Harvey P. Dale_, May 26 2014 *)
%o A108843 (PARI) n10np1(n,d) = { local(x,y,k); for(x=1,n, for(k=1,8, y=10^(x+1)*k+floor(10^x*d/9)*10+k; if(isprime(y),print1(y",")) ) ) }
%o A108843 (Python)
%o A108843 from sympy import isprime
%o A108843 from itertools import count, islice
%o A108843 def agen(): yield from (t for i in count(1) for f in "17" if isprime(t:=int(f + "6"*i + f)))
%o A108843 print(list(islice(agen(), 10))) # _Michael S. Branicky_, Jan 27 2023
%Y A108843 Similar sequences for digit d: A108845 (d=1), A108846 (d=2), A108841 (d=4), A108842 (d=5), A108844 (d=7), A108847 (d=8), A108848 (d=9).
%K A108843 easy,nonn,base
%O A108843 1,1
%A A108843 _Cino Hilliard_, Jul 11 2005
%E A108843 Name changed by _Arkadiusz Wesolowski_, Sep 07 2011