OFFSET
1,1
COMMENTS
The number of terms with 2, 4, 6, ... digits: 4, 92, 3223, 130607, 6350300, ..., . - Robert G. Wilson v, Dec 07 2008
REFERENCES
P. Giannopoulos, The brainteasers (unpublished)
LINKS
Robert G. Wilson v, Table of n, for n = 1..10000. [From Robert G. Wilson v, Dec 07 2008]
EXAMPLE
23 is a member because 2 and 3 are both primes.
MATHEMATICA
f[n_] := Block[{c = 0, lp = PrimePi[10^n] - PrimePi[10^(n - 1)], lq = PrimePi[10^n], lst = {}, pq, p = Prime@ Range[PrimePi[10^(n - 1)] + 1, PrimePi[10^n]], q = Prime@ Range[1, PrimePi[10^n]]}, Do[pq = p[[i]]*10^n + q[[j]]; If[PrimeQ@ pq, AppendTo[lst, pq]; c++ ], {i, lp}, {j, lq}]; lst]; Array[f, 2] // Flatten (* Robert G. Wilson v, Dec 07 2008 *)
pQ[n_]:=Module[{idn=IntegerDigits[n], len}, len=Length[idn]; EvenQ[len] && PrimeQ[FromDigits[Take[idn, len/2]]]&&PrimeQ[FromDigits[Take[idn, -len/2]]]]; Select[Prime[Range[500]], pQ] (* Harvey P. Dale, Nov 08 2011 *)
PROG
(PARI) t=1; forprime( p=2, 99, if( p>t, t*=10); forprime( q=3, t, isprime(p*t+q) & print1(p*t+q, ", "))) \\ M. F. Hasler
CROSSREFS
KEYWORD
base,easy,nonn
AUTHOR
P. Giannopoulos (pgiannop1(AT)yahoo.com), Mar 31 2003
EXTENSIONS
Corrected by Zak Seidov, Robert Israel, Farideh Firoozbakht and M. F. Hasler, Dec 06 2008 and Dec 07 2008
Edited by N. J. A. Sloane, Dec 07 2008
STATUS
approved