|
%I #7 Apr 23 2013 03:05:44
%S 3,13,157,181,739,829,937,1009,1093,1483,1621,1879,2311,2503,2647,
%T 2719,3079,4969,4999,5209,5431,5569,6163,6961,8161,8329,9349,9631,
%U 10399,10459,10531,10657,11131,11953,13063,18523,20149,20731,21391,21589,26317,27481,28111,28351,29023
%N Primes p such that q = 2*p^2 - 1 and 2*p*q - 1 are also prime.
%C Subsequence of A106483, and more elementary version of A224614.
%t Select[Prime[Range[10000]], PrimeQ[2#^2 - 1] && PrimeQ[4#^3 - 2# - 1] &] (* _Alonso del Arte_, Apr 22 2013 *)
%o (PARI) forprime(p=1,3e4,isprime(r=2*p^2-1)&&isprime(2*p*r-1)&&print1(p","))
%K nonn
%O 1,1
%A _M. F. Hasler_, Apr 22 2013
|
