|
%I #20 Apr 01 2024 09:13:29
%S 23,37,53,73,113,131,137,173,179,197,311,313,317,373,379,419,431,479,
%T 613,617,619,673,719,797,971,1013,1019,1031,1097,1277,1373,1499,1571,
%U 1733,1811,1997,2113,2239,2293,2719,3079,3137,3313,3373,3491,3499,3593,3673,3677,3733
%N Consider a prime number p with m decimal digits. The sequence lists the primes p such that the prefix of length m-1 and the suffix of length m-1 are both prime numbers.
%C Let x(0)x(1)... x(q-1)x(q) denote the decimal expansion of a prime p. The sequence lists the primes p such that the prefix x(0)x(1)... x(q-1) and the suffix x(1)... x(q-1)x(q) are primes.
%C Superset of A051362; a(n) first differs from A051362 when n=12.
%H Paolo P. Lava, <a href="/A244282/b244282.txt">Table of n, a(n) for n = 1..1000</a>
%e The prime number 179 is in the sequence because 17 and 79 are primes.
%p with(numtheory):
%p for m from 1 to 200 do:
%p n:=ithprime(m):x:=convert(n, base, 10):n1:=nops(x):
%p s1:=sum('x[i]*10^(i-1) ', 'i'=1..n1-1):
%p s2:=(n-irem(n,10))/10:
%p if type(s1,prime)=true and type(s2,prime)=true
%p then
%p printf(`%d, `, n):
%p else
%p fi:
%p od:
%t Select[Prime[Range[1000]], (id = IntegerDigits[#]; PrimeQ[FromDigits[Take[id, {1, -2}]]] && PrimeQ[FromDigits[Take[id, {2, -1}]]]) &] (* _César Eliud Lozada_, Mar 31 2024 *)
%Y Cf. A000040, A051362, A244283.
%K nonn,base
%O 1,1
%A _Michel Lagneau_, Jun 25 2014
%E Inserted missing term a(49) and corrected a(50) by _Paolo P. Lava_, Dec 04 2017
|
