login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A216167 Composite numbers which yield a prime whenever a 5 is inserted anywhere in them, excluding at the end. 3

%I #19 Oct 02 2022 17:06:41

%S 9,21,57,63,69,77,87,93,153,231,381,407,413,417,501,531,581,651,669,

%T 741,749,783,791,987,1241,1551,1797,1971,2189,2981,3381,3419,3591,

%U 3951,4083,4503,4833,4949,4959,5049,5117,5201,5229,5243,5529,5547,5603,5691,5697

%N Composite numbers which yield a prime whenever a 5 is inserted anywhere in them, excluding at the end.

%H Michael S. Branicky, <a href="/A216167/b216167.txt">Table of n, a(n) for n = 1..1923</a> (terms 1..300 from Paolo P. Lava)

%e 4083 is not prime but 40853, 40583, 45083 and 54083 are all primes.

%p with(numtheory);

%p A216167:=proc(q,x)

%p local a,b,c,i,n,ok;

%p for n from 1 to q do

%p if not isprime(n) then

%p a:=n; b:=0; while a>0 do b:=b+1; a:=trunc(a/10); od; a:=n; ok:=1;

%p for i from 1 to b do c:=a+9*10^i*trunc(a/10^i)+10^i*x;

%p if not isprime(c) then ok:=0; break; fi;

%p od;

%p if ok=1 then print(n); fi;

%p fi;

%p od; end:

%p A216167(1000,5);

%t Select[Range[6000],CompositeQ[#]&&AllTrue[FromDigits/@Table[Insert[IntegerDigits[#],5,p],{p,IntegerLength[#]}],PrimeQ]&] (* _Harvey P. Dale_, Oct 02 2022 *)

%o (Magma) [n: n in [1..6000] | not IsPrime(n) and forall{m: t in [1..#Intseq(n)] | IsPrime(m) where m is (Floor(n/10^t)*10+5)*10^t+n mod 10^t}]; // _Bruno Berselli_, Sep 03 2012

%o (Python)

%o from sympy import isprime

%o def ok(n):

%o if n < 2 or n%10 not in {1, 3, 7, 9} or isprime(n): return False

%o s = str(n)

%o return all(isprime(int(s[:i] + '5' + s[i:])) for i in range(len(s)))

%o print(list(filter(ok, range(5698)))) # _Michael S. Branicky_, Sep 21 2021

%Y Cf. A068673, A068674, A068677, A068679, A069246, A215417, A215419-A215421, A216165, A216166, A216168, A216169.

%Y Subsequence of A053795.

%K nonn,look,base

%O 1,1

%A _Paolo P. Lava_, Sep 03 2012

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 11:59 EDT 2024. Contains 371254 sequences. (Running on oeis4.)