OFFSET
1,1
COMMENTS
EXAMPLE
MAPLE
with(numtheory):T:=array(1..111):k:=0:for n from 2 by 2 to 1000 do: p:=n^2+1:if type(p, prime)=true then k:=k+1:T[k]:=p:else fi:od:for i from 1 to k do:d0:=0:a:=sqrt(T[i]-1):b:=sqrt(T[i+1]-1):for j from a+1 to b-1 do:y:=factorset(j^2+1):n1:=nops(y):d:=y[n1]:if d>d0 then d0:=d:else fi:od: printf(`%d, `, d0):od:
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Dec 30 2013
STATUS
approved