OFFSET
1,1
COMMENTS
Corresponding row/column numbers of the array are 6, 12, 16, 57, 91, 117, 189, 597, 707, 1729, 3933, (through 5000). The values in row/column 597, 707, ... are probable primes.
PROG
(PARI) Stolarsky(r, c) = tau=(1+sqrt(5))/2; a=floor(r*(1+tau)-tau/2); b=round(a*tau); if(c==1, a, if(c==2, b, for(i=1, c-2, d=a+b; a=b; b=d; ); d))
for(x=1, 5000, if(isprime(p=Stolarsky(x, x)), print1(p, ", ")))
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Randall L Rathbun, Jan 25 2002
EXTENSIONS
Edited by Klaus Brockhaus, May 26 2003
Edited by T. D. Noe, Oct 30 2008
STATUS
approved