|
|
|
|
1, 2, 7, 16, 41, 98, 239, 576, 1393, 3362, 8119, 19600, 47321, 114242, 275807, 665856, 1607521, 3880898, 9369319, 22619536, 54608393, 131836322, 318281039, 768398400, 1855077841, 4478554082, 10812186007, 26102926096, 63018038201
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
The logarithmic derivative of this sequence is twice the g.f. of A113282, where a(2*n) = A113282(2*n), a(4*n+1) = A113282(4*n+1) - 3, a(4*n+3) = A113282(4*n+3) - 1.
Equals the self-convolution of integer sequence A113281. (End)
With an offset of 1, this sequence is the case P1 = 2, P2 = 0, Q = -1 of the 3-parameter family of 4th-order linear divisibility sequences found by Williams and Guy. - Peter Bala, Mar 19 2015
Floretion Algebra Multiplication Program, FAMP Code: -2ibaseiseq[B*C], B = - .5'i + .5'j - .5i' + .5j' - 'kk' - .5'ik' - .5'jk' - .5'ki' - .5'kj'; C = + .5'i + .5i' + .5'ii' + .5e
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (1+x^2)/((x-1)*(x+1)*(x^2+2*x-1).
a(n+2) - a(n+1) - a(n) = A100828(n+1).
a(n) = -(u^(n+1)-1)*(v^(n+1)-1)/2 with u = 1+sqrt(2), v = 1-sqrt(2). - Vladeta Jovovic, May 30 2007
a(n) = n * Sum_{k=1..n} Sum_{i=ceiling((n-k)/2)..n-k} binomial(i,n-k-i)*binomial(k+i-1,k-1)*(1-(-1)^k)/(2*k). - Vladimir Kruchinin, Apr 11 2011
exp( Sum_{n >= 1} 2*a(n-1)*x^n/n ) = 1 + 2*Sum_{n >= 1} Pell(n) *x^n. (End)
|
|
MATHEMATICA
|
a[n_] := n*Sum[ Sum[ Binomial[i, n-k-i]*Binomial[k+i-1, k-1], {i, Ceiling[(n-k)/2], n-k}]*(1-(-1)^k)/(2*k), {k, 1, n}]; Table[a[n], {n, 1, 29}] (* Jean-François Alcover, Feb 26 2013, after Vladimir Kruchinin *)
CoefficientList[Series[(1 + x^2) / ((x^2 - 1) (x^2 + 2 x - 1)), {x, 0, 30}], x] (* Vincenzo Librandi, Mar 20 2015 *)
LinearRecurrence[{2, 2, -2, -1}, {1, 2, 7, 16}, 30] (* Harvey P. Dale, Oct 10 2017 *)
|
|
PROG
|
(PARI) {a(n)=local(x=X+X*O(X^n)); polcoeff((1+x^2)/(1-x^2)/(1-2*x-x^2), n, X)} \\ Paul D. Hanna
(Maxima) a(n):=n*sum(sum(binomial(i, n-k-i)*binomial(k+i-1, k-1), i, ceiling((n-k)/2), n-k)*(1-(-1)^k)/(2*k), k, 1, n); /* Vladimir Kruchinin, Apr 11 2011 */
(Magma) [Floor((1+Sqrt(2))^(n+1)/2): n in [0..30]]; // Vincenzo Librandi, Mar 20 2015
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy,changed
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|