cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A121834 Primes p of the form 4*n^2 + 1 such that 4*p^2+1 is also prime.

This page as a plain text file.
%I A121834 #10 Nov 07 2016 17:40:22
%S A121834 5,37,677,1297,2917,8837,13457,50177,147457,156817,246017,341057,
%T A121834 414737,746497,1136357,1726597,1833317,2119937,2802277,2808977,
%U A121834 3013697,3549457,3865157,3896677,4104677,4384837,5354597,5410277,5779217,6031937,6635777,7001317
%N A121834 Primes p of the form 4*n^2 + 1 such that 4*p^2+1 is also prime.
%C A121834 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
%t A121834 fpQ[n_]:=PrimeQ[n]&&PrimeQ[4n^2+1]; Select[4Range[2000]^2+1,fpQ] (* _Harvey P. Dale_, Nov 07 2016 *)
%o A121834 (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*/
%Y A121834 Cf. A002496, A052291, A121326.
%K A121834 nonn
%O A121834 1,1
%A A121834 _Zak Seidov_, Aug 28 2006
%E A121834 More terms from _R. J. Mathar_, Sep 01 2006