OFFSET
1,2
LINKS
Michel Lagneau, Table of n, a(n) for n = 1..5000
MAPLE
with(numtheory):T:=array(1..50000):T[1]:=2:kk:=1:nn:=10^5:
for i from 1 to nn do:
p:=4*i+1:
if type(p, prime)=true
then
kk:=kk+1:T[kk]:=p:
else
fi:
od:
for k from 1 to 5000 do:ii:=0:
y:=factorset(k^2+1):n2:=nops(y):t:=y[n2]:
for l from 1 to kk while(ii=0)do :
if t=T[l]
then
printf(`%d, `, l):
else
fi:
od:
od:
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Oct 13 2014
STATUS
approved