This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A066823 #9 Apr 24 2022 17:19:26 %S A066823 157,6011,55717,74636904300223,1524920857037139836753, %T A066823 533561095785570791522387803, %U A066823 960826779344998235261901369163166807780807 %N A066823 Primes in the main diagonal of the Stolarsky array (cf. A035506). %C A066823 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. %o A066823 (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)) %o A066823 for(x=1,5000,if(isprime(p=Stolarsky(x,x)),print1(p,","))) %Y A066823 Cf. A035506. %K A066823 easy,nonn %O A066823 1,1 %A A066823 _Randall L Rathbun_, Jan 25 2002 %E A066823 Edited by _Klaus Brockhaus_, May 26 2003 %E A066823 Edited by _T. D. Noe_, Oct 30 2008