# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a255898
Showing 1-1 of 1

%I A255898 #47 Nov 21 2018 11:39:22
%S A255898 23,5,3,7,2,3,43,47,3,3,7,11,17,11,3,29,3,11,3,109,11,43,71,19,71,11,
%T A255898 11,3,7,229,43,269,7,23,3,61,37,677,113,863,59,3,11,487,359,347,3,19,
%U A255898 53,173,3,127,229,7,3,3,13,3,241,41,79,79,3,83,23,31,71,31
%N A255898 Minimum prime p such that p^n is a concatenation of two primes.
%H A255898 Paolo P. Lava, <a href="/A255898/b255898.txt">Table of n, a(n) for n = 1..200</a>
%e A255898 23^1 = 23 = concat(2,3);
%e A255898 5^2 = 25 = concat(2,5);
%e A255898 3^3 = 27 = concat(2,7).
%p A255898 with(numtheory): P:= proc(q) local a,k,n,ok;
%p A255898 for a from 1 to q do for n from 1 to q do if isprime(n) then ok:=0;
%p A255898 for k from 1 to ilog10(n^a) do if isprime(trunc(n^a/10^k)) and isprime(n^a mod 10^k) then ok:=1; break; fi; od;
%p A255898 if ok=1 then lprint(a,n); break; fi; fi; od; od; end: P(10^9);
%t A255898 mp[n_]:=Module[{p=2},While[Count[PrimeQ[#]&/@Table[FromDigits/@ TakeDrop[ IntegerDigits[ p^n],i],{i,IntegerLength[p^n]}],{True,True}]== 0,p= NextPrime[ p]];p]; Array[mp,70] (* Requires Mathematica version 10 or later *) (* _Harvey P. Dale_, Oct 13 2016 *)
%o A255898 (PARI) a(n) = {forprime(p=2, , my(pn = p^n); for (k=1, #Str(pn), if (isprime(pn\10^k) && isprime(pn % 10^k), return (p));););} \\ _Michel Marcus_, Oct 22 2015
%Y A255898 Cf. A000040, A255579.
%K A255898 nonn,base,easy
%O A255898 1,1
%A A255898 _Paolo P. Lava_, Oct 21 2015

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE