%I#20 2014年7月20日20:19:05

%S 1,2,0,0,6,2,14,4,10,14,16,4,12,6,18,22,18,14,14,14，24,24,34,20,40,18，

%电话：32,16,36,24,36,16,50,48,36,34,74,0,40,42,60,16,58,34,44,40,44,24,96，

%U 30,64,50,64,36,98,58,80,54,48,52124,36,72,64,74,60,66,52,80,60,92,52

%具有Gram矩阵的二次型N Theta级数[2，1，0，1；1，8，1，-3；0，1，10，2；1，-3，2，12]。

%这是四维Elkies_B晶格。

%H John Cannon，n的表，n=0..5000的a（n）</a>

%H.N.D.埃尔基，<a href=“http://www.math.harvard.edu/~elkies/modular.pdf“>有限域上的椭圆和模曲线及相关计算问题。

%e 1+2*q^2+6*q^8+2*q*10+14*q^12+4*q^14+10*q^16+14*q ^18+16*q^20+。。。

%o（PARI）B（x，y，z，w）=2*x^2+8*y^2+10*z^2+12*w^2+2*x*（y+w）+2*y*（z-3*w）+4*z*w；{θB（n，n，bx，by，bz，bw，ix，iy，iz，iw，nn）=n=2*n；bx=楼层（sqrt（n）*（1+1/sqrt（6）））；bz=楼层（n/7）；bw=楼层（mqrt（n/6））；by=楼层（m2（n/3））；

%o N=矢量（N/2+2）；对于（ix=-bx，bx，对于（iy=-by，by，对于（iz=-bz，bz，对于（iw=-bw，bw，nn=B（ix，iy，iz，iw）））；如果（nn<=n，n[1+nn/2]++）；））；N} a=θB（80）；

%o（PARI）{a（n）=如果（n<1，n==0，qfrep（[2，1，0，1；1，8，1，-3；0，1，10，2；1，-3，2，12]，n，1）[n]*2）}/*_Michael Somos_，2006年4月2日*/

%Y对偶格为A045867中的对偶格。

%K nonn公司

%0、2

%A _N.J.A.斯隆_

%E来自Kok Seng Chua（chuaks（AT）ihpc.nus.edu.sg）的更多条款，2000年6月22日