OFFSET
1,1
COMMENTS
Intersection of A001912 and A121326. Except for the first term all other terms are == 7 (mod 10). Also all the primes 4*p^2+1 are == 7 (mod 10). - Zak Seidov, Mar 05 2015
MATHEMATICA
fpQ[n_]:=PrimeQ[n]&&PrimeQ[4n^2+1]; Select[4Range[2000]^2+1, fpQ] (* Harvey P. Dale, Nov 07 2016 *)
PROG
(PARI) isA121834(n)={ if( (n-1) %4, return(0) ; ) ; if( issquare((n-1)/4), if( isprime(4*n^2+1), return(1), return(0) ), return(0) ; ) ; } { for(i=1, 1000000, p=prime(i) ; if( isA121834(p), print1(p, ", ") ; ) ; ) ; } /* R. J. Mathar, Sep 01 2006*/
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Aug 28 2006
EXTENSIONS
More terms from R. J. Mathar, Sep 01 2006
STATUS
approved