|
%I #41 Feb 14 2023 15:49:55
%S 0,1,7,12,26,35,57,70,100,117,155,176,222,247,301,330,392,425,495,532,
%T 610,651,737,782,876,925,1027,1080,1190,1247,1365,1426,1552,1617,1751,
%U 1820,1962,2035,2185,2262,2420
%N a(n) = (1/2) * ( 3*n^2 + n*(-1)^n ).
%C Previous name was: Pentagonal wave sequence of the first kind.
%C Odd-indexed terms = A033570, pentagonal numbers with odd index (1, 12, 35, 70, ...). Even-indexed terms = A049453, 2nd pentagonal numbers with even index (0, 7, 26, 57, 100, ...).
%C Companion sequence A104585 (Pentagonal wave sequence of the second kind), switches odd with even applications and vice versa. The pentagonal wave sequence triangle A104586 has A104584 in odd columns and A104585 in even columns.
%C Exponents of q in the identity Sum_{n >= 0} ( q^n*Product_{k = 1..n} (1 - q^(4*k-3)) ) = 1 + q - q^7 - q^12 + q^26 + q^35 - - + + .... Compare with Euler's pentagonal number theorem: Product_{n >= 1} (1 - q^n) = 1 - Sum_{n >= 1} ( q^n*Product_{k = 1..n-1} (1 - q^k) ) = 1 - q - q^2 + q^5 + q^7 - q^12 - q^15 + + - - .... - _Peter Bala_, Dec 03 2020
%H Vincenzo Librandi, <a href="/A104584/b104584.txt">Table of n, a(n) for n = 0..1000</a>
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Pentagonal_number_theorem">Pentagonal number theorem</a>.
%H <a href="/index/Rec#order_05">Index entries for linear recurrences with constant coefficients</a>, signature (1,2,-2,-1,1).
%F a(n) = a(n-1) + 2*a(n-2) - 2*a(n-3) - a(n-4) + a(n-5). - _Vincenzo Librandi_, Apr 04 2013
%F a(n) = (1/2) * (3*n^2 + n*(-1)^n ). - _Ralf Stephan_, May 20 2007
%F G.f. -x*(1+6*x+3*x^2+2*x^3) / ( (1+x)^2*(x-1)^3 ). - _R. J. Mathar_, Jan 10 2011
%F From _Amiram Eldar_, Feb 22 2022: (Start)
%F Sum_{n>=1} 1/a(n) = 6 - Pi/sqrt(3) - 4*log(2).
%F Sum_{n>=1} (-1)^(n+1)/a(n) = 2*Pi/sqrt(3) + 3*log(3) - 6. (End)
%e a(5) = 35 = A000326(5).
%e a(6) = 57 = A005449(6).
%t Table[(1/2) (3 n^2 + n (-1)^n), {n, 0, 100}] (* _Vincenzo Librandi_, Apr 04 2013 *)
%t LinearRecurrence[{1,2,-2,-1,1},{0,1,7,12,26},50] (* _Harvey P. Dale_, Feb 14 2023 *)
%o (Magma) I:=[0, 1, 7, 12, 26]; [n le 5 select I[n] else Self(n-1)+2*Self(n-2)-2*Self(n-3)-Self(n-4)+Self(n-5): n in [1..50]]; // _Vincenzo Librandi_, Apr 04 2013
%Y Cf. A000326, A005449, A049453, A033568, A104585, A104586.
%K nonn,easy
%O 0,3
%A _Gary W. Adamson_, Mar 17 2005
%E Better name, using formula from _Ralf Stephan_, _Joerg Arndt_, Sep 17 2013
|
