|
%I #20 Sep 08 2014 17:42:42
%S 2,3,4,5,7,14,21,22,26,27,32,43,51,58,62,65,82,131,148,207,229,249,
%T 257,320,334,386,423,440,481,747,823,1181,1314,1915,2025,2269
%N Numbers n such that numerator of Sum_{k=1..n} 1/(prime(k)-k) is prime.
%C Note that the definition here is subtly different from that of A092063.
%p count:= 0:
%p S:= 0: p:= 0;
%p for n from 1 to 2500 do
%p p:= nextprime(p);
%p S:= S + 1/(p - n);
%p if isprime(numer(S)) then
%p count:= count+1;
%p A[count]:= n;
%p fi
%p od:
%p seq(A[i],i=1..count); # _Robert Israel_, Sep 07 2014
%t f=0; Do[ p=Prime[n]; f=f+1/(p-n); g=Numerator[f]; If[ PrimeQ[g], Print[n]], {n,1,500} ]
%o (PARI) S=1;for(n=2,100,S=S+1/(prime(n)-n);if(isprime(numerator(S)),print1(n,","))) \\ _Edward Jiang_, Sep 08 2014
%Y Cf. A092066.
%K nonn
%O 1,1
%A Mohammed Bouayoun (mohammed.bouayoun(AT)sanef.com), Feb 20 2004; corrected Apr 24 2006
%E Sequence and Mathematica program corrected by _Alexander Adamchuk_, Jul 29 2007
%E a(30)-a(34) from _Vincenzo Librandi_, Nov 26 2012
%E a(35) and a(36) from _Robert Israel_, Sep 07 2014
|
