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 A228456 #20 Dec 08 2018 11:50:41
%S A228456 1,3,41,2841,1010845,1790959059,15625174448153,684155685766047665,
%T A228456 153135670584610711281133,174306862433739623658457865851,
%U A228456 994128440039970508236558371536766105,28386094077591547319667447446929656332136825
%N A228456 Determinant of the (n+1) X (n+1) matrix with (i,j)-entry equal to A086618(i+j) for all i,j = 0,...,n.
%C A228456 Conjecture: (i) a(n) is always positive and odd, and not congruent to 7 modulo 8.
%C A228456 (ii) For any odd prime p, if p == 1 (mod 3) and p = x^2 + 3*y^2 with x == 1 (mod 3), then a(p-1) == (-1)^{(p-1)/2}*(2*x-p/(2*x)) (mod p^2); if p == 2 (mod 3) then a(p-1) == (-1)^{(p+1)/2}*3p/binomial((p+1)/2,(p+1)/6) (mod p^2).
%D A228456 Zhi-Wei Sun, Conjectures and results on x^2 mod p^2 with 4p = x^2+d*y^2, in: Number Theory and Related Area (eds., Y. Ouyang, C. Xing, F. Xu and P. Zhang), Higher Education Press & International Press, Beijing and Boston, 2013, pp. 147-195.
%H A228456 Zhi-Wei Sun, <a href="/A228456/b228456.txt">Table of n, a(n) for n = 0..25</a>
%H A228456 Zhi-Wei Sun, <a href="https://arxiv.org/abs/1103.4325">Conjectures and results on x^2 mod p^2 with 4p = x^2+d*y^2</a>, arXiv:1103.4325 [math.NT], 2011-2014.
%H A228456 Z.-W. Sun, <a href="http://dx.doi.org/10.1016/j.jnt.2012.05.014">On sums of Apéry polynomials and related congruences</a>, J. Number Theory 132(2012), 2673-2699.
%H A228456 Zhi-Wei Sun, <a href="https://doi.org/10.1016/j.jnt.2013.02.014">Connections between p = x^2 + 3*y^2 and Franel numbers</a>, J. Number Theory 133(2013), 2914-2928.
%H A228456 Zhi-Wei Sun, <a href="http://arxiv.org/abs/1308.2900">On some determinants with Legendre symbol entries</a>, preprint, arXiv:1308.2900 [math.NT], 2013-2018.
%e A228456 a(0) = 1 since A086618(0) = 1.
%t A228456 f[n_]:=Sum[Binomial[n,k]^2*Binomial[2k,k]/(k+1),{k,0,n}]
%t A228456 a[n_]:=Det[Table[f[i+j],{i,0,n},{j,0,n}]]
%t A228456 Table[a[n],{n,0,10}]
%o A228456 (PARI) f(n)=sum(k=0,n, binomial(n,k)^2*binomial(2*k,k)/(k+1))
%o A228456 a(n)=my(v=vector(2*n+1,k,f(k-1))); matdet(matrix(n+1,n+1,i,j,v[i+j-1])) \\ _Charles R Greathouse IV_, Jul 30 2016
%Y A228456 Cf. A086618.
%K A228456 nonn,nice
%O A228456 0,2
%A A228456 _Zhi-Wei Sun_, Aug 22 2013