|
%I #47 Mar 27 2022 03:52:05
%S 0,15,54,117,204,315,450,609,792,999,1230,1485,1764,2067,2394,2745,
%T 3120,3519,3942,4389,4860,5355,5874,6417,6984,7575,8190,8829,9492,
%U 10179,10890,11625,12384,13167,13974,14805,15660,16539,17442,18369,19320,20295,21294,22317
%N Write 1,2,3,4,... counterclockwise in a hexagonal spiral around 0 starting left down, then a(n) is the sequence found by reading from 0 in the vertical upward direction.
%C Related to parity of Beatty sequences for exp(-(1/2)/n). Let f(k,n)=-sum(i=1,n,sum(j=1,i,(-1)^floor(j*exp(-(1/2)/n)))), then a(n)=Max{f(k,n) : 1<=k<=4*a(n)-2} and for 0<=i<=4*a(n)-3, f(i,n)=f(4*a(n)-2-i,n). - _Benoit Cloitre_, May 26 2004
%C Or, sum of multiples of 2 and 3 from 0 to 6n. - _Zak Seidov_, Aug 06 2016
%H Harry J. Smith, <a href="/A063436/b063436.txt">Table of n, a(n) for n = 0..1000</a>
%H Leo Tavares, <a href="/A063436/a063436.jpg">Illustration: Stellar Triangles</a>
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (3,-3,1).
%F a(n) = 3*n*(4*n+1) = 3*A007742(n).
%F a(n) = 24*n + a(n-1) - 9 (with a(0)=0). - _Vincenzo Librandi_, Aug 07 2010
%F From _Colin Barker_, Jul 07 2012: (Start)
%F a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
%F G.f.: 3*x*(5+3*x)/(1-x)^3. (End)
%F a(n) = A272399(n+1) - A003154(n+1). - _Leo Tavares_, Mar 24 2022
%F From _Amiram Eldar_, Mar 27 2022: (Start)
%F Sum_{n>=1} 1/a(n) = 4/3 - Pi/6 - log(2).
%F Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(3*sqrt(2)) + log(2)/3 + sqrt(2)*log(sqrt(2)+1)/3 - 4/3. (End)
%e The spiral begins:
%e .
%e 16--15--14
%e / \
%e 17 5---4 13
%e / / \ \
%e 18 6 0 3 12
%e / / / / /
%e 19 7 1---2 11 26
%e \ \ / /
%e 20 8---9--10 25
%e \ /
%e 21--22--23--24
%t a[n_] := 3*n*(4*n + 1); Array[a, 40, 0] (* _Amiram Eldar_, Mar 27 2022 *)
%o (PARI) { for (n=0, 1000, write("b063436.txt", n, " ", n*(12*n + 3)) ) } \\ _Harry J. Smith_, Aug 21 2009
%Y Cf. A062783, A000567.
%Y Cf. A272399, A003154.
%K nonn,easy
%O 0,2
%A _Floor van Lamoen_, Jul 21 2001
|
