OFFSET
1,1
COMMENTS
After 3, 357111317192931414359 is the only prime in the sequence for n up to 10000.
Although 5 appears in two twin prime pairs (3, 5) and (5, 7), 5 is concatenated only once in the sequence. - Daniel Forgues, Aug 23 2016
a(n) == 0 mod 3 for n odd, a(n) == 2 mod 3 for n even. - Robert Israel, Sep 01 2016
LINKS
Robert Israel, Table of n, a(n) for n = 1..270
C. K. Caldwell, Twin Primes.
MAPLE
Primes:= select(isprime, {seq(i, i=1..100, 2)}):
T1:= Primes intersect map(`+`, Primes, 2):
Twins:= sort(convert(T1 union map(`-`, T1, 2), list)):
dcat:= (a, b) -> a*10^(1+ilog10(b))+b:
A[1]:= 3:
for n from 2 to nops(Twins) do A[n]:= dcat(A[n-1], Twins[n]) od:
seq(A[i], i=1..nops(Twins)); # Robert Israel, Sep 01 2016
MATHEMATICA
Table[FromDigits@ Flatten@ Map[IntegerDigits, Take[#, n]], {n, Length@ #}] &[Union@ Join[#, # + 2] &@ Select[Prime@ Range@ 17, NextPrime@ # - 2 == # &]] (* Michael De Vlieger, Sep 01 2016 *)
Module[{tps=Union[Flatten[Select[Partition[Prime[Range[50]], 2, 1], #[[2]]-#[[1]] == 2&]]]}, FromDigits[Flatten[IntegerDigits/@#]]&/@Table[Take[tps, n], {n, Length[tps]}]] (* Harvey P. Dale, Jun 16 2022 *)
PROG
(PARI) concattwprb(n) = { y=3; forprime(x=5, n, if(isprime(x+2) || isprime(x-2), y=eval(concat(Str(y), Str(x))); print1(y", ") ) ) }
CROSSREFS
KEYWORD
easy,nonn,base
AUTHOR
Cino Hilliard, Sep 08 2003
EXTENSIONS
Edited by N. J. A. Sloane, Jul 01 2008 at the suggestion of R. J. Mathar
STATUS
approved