OFFSET
1,2
COMMENTS
Squares n such that 10*n+1 is prime. - Robert Israel, Dec 09 2017
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
41, 4001, 6761 etc. are primes.
MAPLE
select(t -> isprime(10*t+1), map(`^`, [$1..300], 2)); # Robert Israel, Dec 09 2017
MATHEMATICA
Select[Range[150]^2, PrimeQ[10 # + 1] &] (* Michael De Vlieger, Dec 09 2017 *)
PROG
(MATLAB) A = []; count = 0; i = 1; while count < 60; s = i*i; if isprime(10*s + 1) A = [A s]; count = count + 1; end; i = i + 1; end; A
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Amarnath Murthy, Nov 23 2003
EXTENSIONS
Corrected and extended by David Wasserman, Feb 25 2004
STATUS
approved