%I #28 Feb 12 2024 01:53:07
%S 2,3,5,7,11,31,71,13,73,17,37,97,79,101,701,311,131,941,151,751,761,
%T 971,181,191,991,113,313,733,743,353,953,373,383,983,107,907,727,337,
%U 937,347,157,757,167,967,787,797,709,919,929,739,149,359,769,179,389,199
%N Reverse digits of primes, append to sequence if result is a prime.
%C Conjecture: the Benford law limit is 2=Sum[N[Log[10, 1 + 1/d[[n]]]], {n, 1, Length[d]}]^2/(( #totalprimes/#totalPrimes)). At 50000 primes total it is 2.05931. - _Roger L. Bagula_ and _Gary W. Adamson_, Jul 02 2008
%C Presumably this does not satisfy Benford's law. - _N. J. A. Sloane_, Feb 09 2017
%H Indranil Ghosh, <a href="/A095180/b095180.txt">Table of n, a(n) for n = 1..25000</a> (terms 1..206 from Roger L. Bagula and Gary W. Adamson)
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/BenfordsLaw.html">Benford's Law</a>.
%H <a href="/index/Be#Benford">Index entries for sequences related to Benford's law</a>
%e The prime 107 in reverse is 701 which is prime.
%t b = Flatten[Table[If[PrimeQ[Sum[IntegerDigits[Prime[n]][[i]]*10^(i - 1), {i, 1, Length[IntegerDigits[Prime[n]]]}]], Sum[IntegerDigits[Prime[n]][[i]]*10^(i - 1), {i, 1, Length[IntegerDigits[Prime[n]]]}], {}], {n, 1,1000}]] (* _Roger L. Bagula_ and _Gary W. Adamson_, Jul 02 2008 *)
%t Select[FromDigits[Reverse[IntegerDigits[#]]]&/@Prime[Range[300]],PrimeQ] (* _Harvey P. Dale_, May 05 2015 *)
%o (PARI) r(n) = forprime(x=1,n,y=eval(rev(x));if(isprime(y),print1(y","))) \ Get the reverse of the input string rev(str) = { local(tmp,j,s); tmp = Vec(Str(str)); s=""; forstep(j=length(tmp),1,-1, s=concat(s,tmp[j])); return(s) }
%o (Haskell)
%o a095180 n = a095180_list !! (n-1)
%o a095180_list =filter ((== 1) . a010051) a004087_list
%o -- _Reinhard Zumkeller_, Oct 14 2011
%Y Cf. A007500.
%Y Cf. A004087, A010051.
%K base,easy,nonn
%O 1,1
%A _Cino Hilliard_, Jun 21 2004