|
%I #12 Dec 19 2016 14:00:15
%S 2,13,23,29,37,41,79,107,127,149,211,239,263,313,383,389,397,439,467,
%T 509,547,631,673,743,757,773,787,827,829,907,997,1019,1061,1091,1231,
%U 1297,1367,1451,1459,1487,1543,1559,1601,1609,1613,1627,1637,1699,1721
%N Primes p such that the nearest integer to p * A is prime, where A is the Glaisher-Kinkelin constant.
%H Vincenzo Librandi, <a href="/A118796/b118796.txt">Table of n, a(n) for n = 1..1000</a>
%t lst={};Do[p=Prime[n];If[PrimeQ[Round[p*Glaisher]],AppendTo[lst,p]],{n,10^3}];lst
%t Select[Prime[Range[300]],PrimeQ[Round[Glaisher*#]]&] (* _Harvey P. Dale_, Dec 19 2016 *)
%o (PARI) A=exp(1/12-zeta'(-1)); forprime(p=2,1e4,if(isprime(round(p*A)), print1(p", "))) \\ _Charles R Greathouse IV_, Dec 12 2013
%Y Cf. A074962.
%K nonn,easy
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, Aug 25 2008
%E Edited by _Charles R Greathouse IV_, Nov 11 2009
|
