OFFSET
1,1
PROG
(PARI) limit=100; v=vector(limit); for(n=0, limit, for(k=2, floor(limit^(1/2)), for(i=1, limit/k, for(j=1, limit/k, if(n==(k*i+1)*(k*j-1), v[n]=1, ))))); for(n=1, limit, if(v[n], print1(n", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Gionata Neri, Mar 11 2018
STATUS
approved