login
Number of points of L1 norm 4 in cubic lattice Z^n.
6

%I #73 Mar 14 2024 17:20:54

%S 0,2,16,66,192,450,912,1666,2816,4482,6800,9922,14016,19266,25872,

%T 34050,44032,56066,70416,87362,107200,130242,156816,187266,221952,

%U 261250,305552,355266,410816,472642,541200,616962,700416,792066

%N Number of points of L1 norm 4 in cubic lattice Z^n.

%H Vincenzo Librandi, <a href="/A035598/b035598.txt">Table of n, a(n) for n = 0..1000</a>

%H J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (<a href="http://neilsloane.com/doc/Me220.pdf">pdf</a>).

%H Milan Janjic and Boris Petkovic, <a href="http://arxiv.org/abs/1301.4550">A Counting Function</a>, arXiv preprint arXiv:1301.4550 [math.CO], 2013. - _N. J. A. Sloane_, Feb 13 2013

%H Milan Janjic and Boris Petkovic, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL17/Janjic/janjic45.html">A Counting Function Generalizing Binomial Coefficients and Some Other Classes of Integers</a>, J. Int. Seq. 17 (2014), Article 14.3.5.

%H Joan Serra-Sagrista, <a href="http://dx.doi.org/10.1016/S0020-0190(00)00119-8">Enumeration of lattice points in l_1 norm</a>, Inf. Proc. Lett. 76 (1-2) (2000) 39-44.

%H <a href="/index/Rec#order_05">Index entries for linear recurrences with constant coefficients</a>, signature (5,-10,10,-5,1).

%F a(n) = 2*n^2*(n^2 + 2)/3. - _Frank Ellermann_, Mar 16 2002

%F G.f.: 2*x*(1+x)^3/(1-x)^5. - _Colin Barker_, Apr 15 2012

%F a(n) = 2*A014820(n-1). - _R. J. Mathar_, Dec 10 2013

%F a(n) = a(n-1) + A035597(n) + A035597(n-1). - _Bruce J. Nicholson_, Mar 11 2018

%F From _Shel Kaphan_, Feb 28 2023: (Start)

%F a(n) = 2*n*Hypergeometric2F1(1-n,1-k,2,2), where k=4.

%F a(n) = A001846(n) - A001845(n).

%F a(n) = A008412(n)*n/4. (End)

%F From _Amiram Eldar_, Mar 12 2023: (Start)

%F Sum_{n>=1} 1/a(n) = Pi^2/8 - 3*Pi*coth(sqrt(2)*Pi)/(8*sqrt(2)) + 3/16.

%F Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/16 + 3*Pi*cosech(sqrt(2)*Pi)/(8*sqrt(2)) - 3/16. (End)

%F E.g.f.: 2*exp(x)*x*(3 + 9*x + 6*x^2 + x^3)/3. - _Stefano Spezia_, Mar 14 2024

%p f := proc(d,m) local i; sum( 2^i*binomial(d,i)*binomial(m-1,i-1),i=1..min(d,m)); end; # n=dimension, m=norm

%t CoefficientList[Series[2*x*(1+x)^3/(1-x)^5,{x,0,40}],x] (* _Vincenzo Librandi_, Apr 22 2012 *)

%t LinearRecurrence[{5,-10,10,-5,1},{0,2,16,66,192},50] (* _Harvey P. Dale_, Dec 11 2019 *)

%o (PARI) a(n)=2*n^2*(n^2+2)/3 \\ _Charles R Greathouse IV_, Dec 07 2011

%o (Magma) [( 2*n^4 +4*n^2 )/3: n in [0..40]]; // _Vincenzo Librandi_, Apr 22 2012

%Y Cf. A035596, A035597, A035599, A035600, A035601, A035602, A035603, A035604, A035605, A035606.

%Y Cf. A001845, A001846, A008412, A014820.

%Y Column 4 of A035607, A266213, A343599.

%Y Row 4 of A113413, A119800, A122542.

%K nonn,easy

%O 0,2

%A _N. J. A. Sloane_