OFFSET
0,3
COMMENTS
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..500
N. Heninger, E. M. Rains and N. J. A. Sloane, On the Integrality of n-th Roots of Generating Functions, arXiv:math/0509316 [math.NT], 2005-2006.
N. Heninger, E. M. Rains and N. J. A. Sloane, On the Integrality of n-th Roots of Generating Functions, J. Combinatorial Theory, Series A, 113 (2006), 1732-1745.
Ronald Orozco López, Generating Functions of Generalized Simplicial Polytopic Numbers and (s,t)-Derivatives of Partial Theta Function, arXiv:2408.08943 [math.CO], 2024. See p. 11.
Index entries for linear recurrences with constant coefficients, signature (3,6,-8).
FORMULA
a(n) = (2*4^n - (-2)^n - 1)/9;
a(n) = 3*a(n-1) + 6*a(n-2) - 8*a(n-3), a(0)=0, a(1)=1, a(2)=3.
G.f.: x/((1+2*x)*(1-x)*(1-4*x)).
E.g.f.: (2*exp(4*x) - exp(x) - exp(-2*x))/9.
a(n+1) - 4*a(n) = 1, -1, 3, -5, 11, ... = A001045(n+1) signed. - Paul Curtz, May 19 2008
a(n) = round(2^n/3) * round(2^(n+1)/3). - Gary Detlefs, Feb 10 2010
From Peter Bala, Mar 30 2015: (Start)
The shifted o.g.f. A(x) := 1/( (1 + 2*x)*(1 - x)*(1 - 4*x) ) = 1/(1 - 3*x - 6*x^2 + 8*x^3). Hence A(x) == 1/(1 - 3*x + 3*x^2 - x^3) (mod 9) == 1/(1 - x)^3 (mod 9). It follows by Theorem 1 of Heninger et al. that (A(x))^(1/3) = 1 + x + 4*x^2 + 10*x^3 + ... has integral coefficients.
MAPLE
for n from 1 to 25 do print(round(2^n/3)*round(2^(n+1)/3)) od; # Gary Detlefs, Feb 10 2010
MATHEMATICA
Table[(2*4^n -(-2)^n -1)/9, {n, 0, 30}] (* Vladimir Joseph Stephan Orlovsky, Feb 05 2011, modified by G. C. Greubel, Sep 21 2019 *)
LinearRecurrence[{3, 6, -8}, {0, 1, 3}, 25] (* Jean-François Alcover, Sep 21 2017 *)
PROG
(Sage) [gaussian_binomial(n, 2, -2) for n in range(1, 26)] # Zerinvary Lajos, May 28 2009
(Magma) [(2*4^n-(-2)^n-1)/9: n in [0..30]]; // Vincenzo Librandi, Jun 04 2011
(PARI) a(n)=(2*4^n-(-2)^n-1)/9 \\ Charles R Greathouse IV, Sep 24 2015
(GAP) List([0..30], n-> (2^(2*n+1) -(-2)^n -1)/9); # G. C. Greubel, Sep 21 2019
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Paul Barry, May 18 2003
STATUS
approved