|
|
A084175
|
|
Jacobsthal oblong numbers.
|
|
15
|
|
|
0, 1, 3, 15, 55, 231, 903, 3655, 14535, 58311, 232903, 932295, 3727815, 14913991, 59650503, 238612935, 954429895, 3817763271, 15270965703, 61084037575, 244335800775, 977343902151, 3909374210503, 15637499638215, 62549992960455
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
Inverse binomial transform is A001019 doubled up.
|
|
LINKS
|
|
|
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
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.
Sum_{n >= 0} a(n+1)*x^n = exp( Sum_{n >= 1} J(3*n)/J(n)*x^n/n ), where J(n) = A001045(n) are the Jacobsthal numbers. Cf. A001656, A099930. (End)
|
|
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
|
|
|
PROG
|
(Sage) [gaussian_binomial(n, 2, -2) for n in range(1, 26)] # Zerinvary Lajos, May 28 2009
(GAP) List([0..30], n-> (2^(2*n+1) -(-2)^n -1)/9); # G. C. Greubel, Sep 21 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|