OFFSET
1,1
LINKS
Jinyuan Wang, Table of n, a(n) for n = 1..10000
PROG
(PARI) A020482(k) = forprime(q=2, k, if(isprime(2*k-q), return(2*k-q)));
a(n) = {my(r=(ceil(sqrt(2*n+1)))^2-2*n+3); if(r%2==0, r=A020482(r/2), if(isprime(r-2), r-=2, r=A020482(r\2))); r; } \\ Jinyuan Wang, Jan 29 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Ray Chandler, Jan 10 2005
STATUS
approved