%I #24 Sep 08 2022 08:45:48
%S 4,5,9,11,15,17,19,23,25,29,33,35,39,41,45,47,51,57,67,69,71,73,77,81,
%T 83,85,107,113,131,133,145,149,155,157,165,167,169,173,175,179,187,
%U 193,195,203,215,217,225,229,233,239,245,247,269,291,293,295,329,331,337
%N Numbers k such that prime(k) - k - 1 is prime.
%H Harvey P. Dale, <a href="/A166562/b166562.txt">Table of n, a(n) for n = 1..1000</a>
%p isA166562 := proc(n) isprime(ithprime(n)-n-1) ; end: for n from 1 to 500 do if isA166562(n) then printf("%d,",n) ; fi; od: # _R. J. Mathar_, Oct 17 2009
%t Select[Range[400],PrimeQ[Prime[#]-#-1]&] (* _Harvey P. Dale_, Jan 19 2021 *)
%o (PARI) n=0;forprime(p=2,1e4,n++;if(isprime(p-n-1),print1(n", "))) \\ _Charles R Greathouse IV_, May 25 2011
%o (Magma) [ n: n in [1..340] | IsPrime(NthPrime(n)-n-1) ]; // _Bruno Berselli_, May 25 2011
%Y Cf. A065890. - _R. J. Mathar_, Oct 17 2009
%K nonn,less,easy
%O 1,1
%A _Juri-Stepan Gerasimov_, Oct 17 2009
%E 25 inserted, 39 inserted, 53 removed, etc. by _R. J. Mathar_, Oct 17 2009