Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #12 Sep 08 2022 08:45:51
%S 2,107,137,167,257,269,311,317,557,593,761,773,809,911,1103,1151,1283,
%T 1289,1481,1487,1559,1709,1787,1931,2111,2141,2243,2339,2357,2657,
%U 2687,2777,2909,3137,3209,3251,3359,3371,3389,3449
%N Prime numbers p such that the concatenation p^3//1331 is a prime number.
%C See comments at A174213.
%C p^3//1331 is the concatenation of the cubes of two primes.
%C With the exception of a(1)=2, each term is necessarily of the form 6*k-1.
%H Harvey P. Dale, <a href="/A174409/b174409.txt">Table of n, a(n) for n = 1..1000</a>
%e The first prime is 2; 2^3 = 8, and 81331 = prime(7958), so a(1)=2.
%e The smallest prime p > 2 such that p^3//1331 yields a prime is p=107: 107^3 = 1225043, and 12250431331 = prime(552342812), so a(2)=107.
%t Select[Prime[Range[500]],PrimeQ[10000#^3+1331]&] (* _Harvey P. Dale_, May 30 2017 *)
%o (Magma) [p: p in PrimesUpTo(5000) | IsPrime(Seqint(Intseq(1331) cat Intseq(p^3)))]; // _Vincenzo Librandi_, Mar 05 2018
%Y Cf. A168327, A168417, A173836, A174213, A174260, A174229, A174355.
%K base,nonn
%O 1,1
%A Eva-Maria Zschorn (e-m.zschorn(AT)zaschendorf.km3.de), Mar 19 2010