|
%I #26 Sep 08 2022 08:46:05
%S 0,3,13,29,52,81,117,159,208,263,325,393,468,549,637,731,832,939,1053,
%T 1173,1300,1433,1573,1719,1872,2031,2197,2369,2548,2733,2925,3123,
%U 3328,3539,3757,3981,4212,4449,4693,4943,5200,5463,5733,6009
%N a(n) = floor(13*n^2/4).
%C This generalizes A032527, which uses 5, to 13 (the next prime 1 (mod 4)). The figures in A032527 use n/2 concentric dotted pentagons for even n, and (n-1)/2 concentric dotted pentagons plus an extra dot in the middle if n is odd. In the present case one can take n/2 concentric dotted 13-gons (the dot numbers of each side for these 13-gons are 2, 4, 6, ..., n) for even n>=2. There is no figure for n = 0. For odd n one has (n-1)/2 concentric dotted 13-gons (the dot numbers of each side for these 13-gons are 3, 5, 7, ..., n) and an extra dotted 3-gon in the middle. See the example section below for the counting.
%C a(n) = -N(-floor(n/2),n) with the N(a,b) = ((2*a+b)^2 - b^2*13)/4, the norm for integers a + b*omega(13), a, b rational integers, in the quadratic number field Q(sqrt(13)), where omega(13) = (1 + sqrt(13))/2.
%C a(n) = max({|N(a,n)|,a = -n..+n}) = |N(-floor(n/2),n)| = 3*n^2 + n*floor(n/2) - floor(n/2)^2 = floor(13*n^2/4) (the last eq. checks for even and odd n).
%C In the general case one has for primes 1 (mod 4), p(k) = A002144(k), k >= 1, a(p(k);n) = floor(p(k)*n^2/4) with o.g.f. G(p(k);x) = x*(A(k)*(1+x^2) + B(k)*x)/((1-x)^3*(1+x)), where A(k) = A005098(k) = (p(k)-1)/4 and B(k) = A119681(k) = (p(k)+1)/2. This follows from the alternative formula a(p(k),n) = p(k)*n^2/4 + ((-1)^n-1)/8, n >= 0 (which checks for even and odd n). Because the denominator of the o.g.f. is 1-2*x+2*x^3-x^4 the recurrence given by _Bruno Berselli_ below holds for all a(k;n) sequences with inputs for n = -1, 0, 1, 2 given by (p(k)-1)/4, 0, (p(k)-1)/4, p(k), respectively.
%C The dot counting in the concentric p(k)-gons is similar to the one described for p = 5 in A032527 and for p=13 here. For odd n one puts an additional dotted A(k)-gon into the center. - _Wolfdieter Lang_, Aug 08 2013
%H Ivan Panchenko, <a href="/A227541/b227541.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (2,0,-2,1).
%F a(n) = 13*n^2/4+((-1)^n-1)/8, n >= 0 (use even or odd n to prove it).
%F G.f.: x*(3+7*x+3*x^2)/((1-x)^3*(1+x)).
%F a(2*k) = k^2*13, k >= 0.
%F a(2*k+1) = 3 + k*(k+1)*13, k >= 0.
%F a(n) = a(-n) = 2*a(n-1) - 2*a(n-3) + a(n-4). - _Bruno Berselli_, Aug 08 2013
%F a(n) = Sum_{j=1..n} Sum_{i=1..n} ceiling((i+j-n+5)/2). - _Wesley Ivan Hurt_, Mar 12 2015
%e Counting dots in the concentric dotted 13-gons described above in a comment:
%e a(2*k), k >= 1: (2-1)*13 = 13, (1+(4-1))*13 = 52, (1+3+(6-1))*13 = 117, (1+3+5+7)*13 = 208, ... a(2*k+1), k >= 0: 3, 3+(3-1)*13 = 29, 3+(2+(5-1))*13 = 81, 3+2*(1+2+3)*13 = 159, ... (a dotted triangle is put into the middle of the k concentric 13-gons).
%p A227541:=n->floor(13*n^2/4); seq(A227541(n), n=0..50); # _Wesley Ivan Hurt_, Jun 09 2014
%t Table[Floor[13*n^2/4], {n, 0, 50}] (* _Wesley Ivan Hurt_, Jun 09 2014 *)
%o (Magma) [ Floor(13*n^2/4) : n in [0..50] ]; // _Wesley Ivan Hurt_, Jun 09 2014
%Y Cf. A032527 (case for prime 5).
%K nonn,easy
%O 0,2
%A _Wolfdieter Lang_, Aug 07 2013
|
