OFFSET
1,1
COMMENTS
EXAMPLE
MAPLE
with(numtheory):T:=array(1..200):k:=1:for p from 1 to 1000 do: if type(p, prime)=true
and irem(p, 4)=1 then T[k]:=p:k:=k+1:else fi:od:for q from 1 to k do:z:=T[q]:ind:=0:for n from 1 to 10000 while(ind=0) do: x:=n^2+1:y:=factorset(x):if z=y[1] then ind:=1:printf(`%d, `, n):else fi:od: od:
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Jan 18 2011
STATUS
approved