This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A174441 #9 May 08 2016 08:51:01 %S A174441 53,347,431,641,647,821,1709,1973,2081,2591,2657,2963,4073,4139,4643, %T A174441 4787,5039,5483,5657,6029,6791,6917,6959,7127,7673,8273,8693,8807, %U A174441 8849,9221,9311,9689,10139,10457,11423,12503,12743,13619,13913,14549 %N A174441 Primes p such that the concatenations p//1331 and 1331//p are both prime numbers (for naturals see A174355). %C A174441 See comments and references for A173836, A174213. %H A174441 Harvey P. Dale, <a href="/A174441/b174441.txt">Table of n, a(n) for n = 1..1000</a> %e A174441 531331 = prime(43928), 133153 = prime(12427) => p(1) = 53 = prime(16). %e A174441 3471331 = prime(248286), 1331347 = prime(102237) => p(2) = 347 = prime(69). %e A174441 139131331 = prime(7865788), 133113913 = prime(7544750) => p(39) = 13913 = prime(1645). %t A174441 Select[Prime[Range[2000]],AllTrue[{#*10^4+1331,1331*10^IntegerLength[ #]+#}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, May 08 2016 *) %o A174441 (PARI) isok(n) = isprime(n) && isprime(n*10^4 + 1331) && isprime(1331*10^(length(Str(n))) + n); \\ _Michel Marcus_, Aug 27 2013 %Y A174441 Cf. A168327, A168417, A173836, A174213, A174260, A174355. %K A174441 base,nonn %O A174441 1,1 %A A174441 Eva-Maria Zschorn (e-m.zschorn(AT)zaschendorf.km3.de), Mar 20 2010