A154428 Primes of the form 50n^2 + 10n + 1.

61, 1301, 1861, 2521, 5101, 7321, 8581, 9941, 14621, 16381, 20201, 24421, 26681, 34061, 36721, 51521, 68821, 76441, 97241, 101701, 106261, 110921, 135721, 163021, 168781, 199081, 205441, 218461, 252761, 282001, 304981, 312841, 337021, 353641

Offset: 1

Author

Vincenzo Librandi, Jan 09 2009

Keywords

• nonn
• easy

Comments

Subsequence of A027862 associated with the values of A027861 that are multiples of 5. [R. J. Mathar, Jan 12 2009]

Links

• Vincenzo Librandi, Table of n, a(n) for n = 1..1000

Crossrefs

Cf. A017281, A048161.

Programs

GAP

Filtered(List([1..100],n->50*n^2+10*n+1),IsPrime); # Muniru A Asiru, Apr 25 2019

Magma

[a: n in [0..100] | IsPrime(a) where a is 50*n^2 + 10*n + 1]; // Vincenzo Librandi, Jul 23 2012

Maple

select(isprime,[50*n^2+10*n+1$n=1..100])[]; # Muniru A Asiru, Apr 25 2019

Mathematica

Select[Table[50n^2+10n+1,{n,0,200}],PrimeQ] (* Vincenzo Librandi, Jul 23 2012 *)

PARI

for (n=0, 100, if (isprime (k=50*n^2+10*n+1), print1 (k, ", "))); \\ Vincenzo Librandi, Jul 23 2012

Extensions

Replaced 13721 by 135721 - R. J. Mathar, Jan 12 2009