OFFSET
0,3
COMMENTS
For n>1, a(n) is also the maximum value of the n-th difference of a permutation of 1..n. - Michel Marcus, Apr 15 2017
LINKS
Fung Lam, Table of n, a(n) for n = 0..3000
F. Disanto, A. Frosini, and S. Rinaldi, Square involutions, J. Int. Seq. 14 (2011) # 11.3.5.
FORMULA
a(n) = (n+1)*(2^(n-1)-binomial(n-1,n/2)) if n is even else ((n+1)/2)*(2^n-binomial(n,(n+1)/2)). - Vladeta Jovovic, Aug 23 2007
a(n) = (n+1)*(2^n-binomial(n,[n/2]))/2, where [x] is floor. - Graeme McRae, Jan 30 2012
G.f.: (1-sqrt((1-2*x)/(1+2*x)))/(2*(1-2*x)^2). - Vladeta Jovovic, Aug 24 2007
Asymptotics: a(n) ~ 2^(n-1)*(n+1-sqrt(2*n/Pi)). - Fung Lam, Mar 28 2014
D-finite with recurrence (n-1)*n*a(n) = 2*(n-1)*(n+1)*a(n-1) + 4*(n-2)*n*a(n-2) - 8*(n-1)*n*a(n-3). - Vaclav Kotesovec, Mar 28 2014
EXAMPLE
a(1)=1 because 0 1 has a first difference of 1;
a(2)=3 because 2 0 1 has a second difference of 3.
MAPLE
A130783:=n->(n+1)*(2^n-binomial(n, floor(n/2)))/2; seq(A130783(n), n=0..50); # Wesley Ivan Hurt, Nov 25 2013
MATHEMATICA
Table[(n + 1) (2^n - Binomial[n, Floor[n/2]])/2, {n, 0, 50}] (* Wesley Ivan Hurt, Nov 25 2013 *)
PROG
(PARI) a(n)=(n+1)*(2^n-binomial(n, n\2))/2 \\ Charles R Greathouse IV, Jan 30 2012
(Python)
from math import comb
def A130783(n): return (n+1)*((1<<n)-comb(n, n>>1))>>1 # Chai Wah Wu, Jun 04 2024
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
R. H. Hardin, Aug 19 2007
STATUS
approved