|
%I #16 Feb 27 2014 10:56:00
%S 23,37,53,73,113,137,173,193,197,211,223,229,233,241,271,283,293,311,
%T 331,337,347,353,359,367,379,383,389,397,433,523,541,547,571,593,613,
%U 617,673,677,719,733,743,761,773,977,1013,1033,1093,1097,1117,1123,1129
%N Primes which are the concatenation of two primes in exactly one way.
%C This is not a duplicate of A129800, which accepts "07" for example as the second prime.
%H Giovanni Resta, <a href="/A238056/b238056.txt">Table of n, a(n) for n = 1..5000</a>
%e 113 is in the sequence because 11 and 3 are both primes, but 1 and 13 are not both primes, so there is one way.
%t spl[n_] := Block[{d = IntegerDigits@n, c = 0, z}, z = Length@d; Do[If[PrimeQ@ FromDigits@ Take[d, k] && d[[k + 1]] > 0 && PrimeQ@ FromDigits@ Take[d, k - z], c++], {k, z - 1}]; c]; Select[ Prime@ Range@ 300, spl[#] == 1 &] (* _Giovanni Resta_, Feb 27 2014 *)
%o (Haskell)
%o a238056 n = a238056_list !! (n-1)
%o a238056_list = filter ((== 1) . length . f) a000040_list where
%o f x = filter (\(us, vs) ->
%o head vs /= '0' &&
%o a010051' (read us :: Integer) == 1 &&
%o a010051' (read vs :: Integer) == 1) $
%o map (flip splitAt $ show x) [1 .. length (show x) - 1]
%o -- _Reinhard Zumkeller_, Feb 27 2014
%Y Cf. A105184, A238057, A129800.
%Y Cf. A010051, A000040.
%K nonn,base
%O 1,1
%A _Colin Barker_, Feb 17 2014
|
