OFFSET
1,1
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
PROG
(PARI) is(k)=my(m2, mm2, s); for(m=2, sqrtint(k\2), m2=2*m; mm2=2*m^2; forstep(n=if(m%2, 2, 1), min(k\m2-m, m-1), m%2+1, if(k%(mm2+m2*n)==0 && gcd(m, n)==1 && s++>1, return(1)))); 0 \\ Charles R Greathouse IV, Feb 03 2018
(PARI) list(lim)=my(v=List(), ct=Map(), m2, mm2, t); for(m=2, sqrtint(lim\2), m2=2*m; mm2=2*m^2; forstep(n=if(m%2, 2, 1), min(lim\m2-m, m-1), m%2+1, if(gcd(m, n)>1, next); t=mm2+m2*n; forstep(k=t, lim, t, mapput(ct, k, mapisdefined(ct, k))))); for(n=60, lim, if(mapisdefined(ct, n) && mapget(ct, n), listput(v, n))); ct=0; Vec(v) \\ Charles R Greathouse IV, Feb 03 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved