OFFSET
1,1
COMMENTS
A002144 are the primes of the form 4q + 1.
EXAMPLE
MAPLE
with(numtheory):nn:=400:T:=array(1..nn):k:=1:for x from 1 to nn do: p:=4*x+1:if
type(p, prime)=true then T[k]:=p:k:=k+1:else fi:od:for n from 1 to k do: ind:=0:for
m from 1 to 500000 while(ind=0) do:y:=m^2+1:x:= factorset(y) : n1:=nops(x):n2
:=bigomega(y):if n1=2 and n2 = 3 and x[1]=T[n] and ind=0 then ind:=1:printf(`%d,
`, m):else fi:od:od:
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Jan 25 2011
STATUS
approved