OFFSET
1,1
COMMENTS
Except for a(1) = 3, these numbers are not prime for n < 10000. See A089948.
LINKS
Robert Israel, Table of n, a(n) for n = 1..298
MAPLE
a[1]:= 3;
for n from 2 to 100 do
p:= ithprime(n+1);
a[n]:= 10^(1+ilog10(p))*a[n-1]+p
od:
seq(a[i], i=1..100); # Robert Israel, Aug 25 2016
MATHEMATICA
Table[FromDigits[Flatten[IntegerDigits[Prime[Range[n] + 1]]]], {n, 20}] (* Vincenzo Librandi, Nov 30 2015 *)
PROG
(PARI) concatprime(n) = { y=""; forprime(x=3, n, y=concat(Str(y), Str(x)); z=eval(y); print1(z", ") ) }
(Magma) [Seqint(Reverse(&cat[Reverse(Intseq(NthPrime(k+1))): k in [1..n]])): n in [1..17]]; // Vincenzo Librandi, Nov 30 2015
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Cino Hilliard, Jan 11 2004
EXTENSIONS
Offset changed by Robert Israel, Aug 25 2016
STATUS
approved