A123746 -id:A123746 - cp's OEIS frontend

A049606 Largest odd divisor of n!.

1, 1, 1, 3, 3, 15, 45, 315, 315, 2835, 14175, 155925, 467775, 6081075, 42567525, 638512875, 638512875, 10854718875, 97692469875, 1856156927625, 9280784638125, 194896477400625, 2143861251406875, 49308808782358125, 147926426347074375, 3698160658676859375

Offset: 0

N. J. A. Sloane, Feb 05 2000

Original name: Denominator of 2^n/n!.
For positive n, a(n) equals the numerator of the permanent of the n X n matrix whose (i,j)-entry is cos(i*Pi/3)*cos(j*Pi/3) (see example below). - John M. Campbell, May 28 2011
a(n) is also the number of binomial heaps with n nodes. - Zhujun Zhang, Jun 16 2019
a(n) is the number of 2-Sylow subgroups of the symmetric group S_n (see the Mathematics Stack Exchange link below). - Jianing Song, Nov 11 2022

			From _John M. Campbell_, May 28 2011: (Start)
The numerator of the permanent of the following 5 X 5 matrix is equal to a(5):
|  1/4  -1/4  -1/2  -1/4   1/4 |
| -1/4   1/4   1/2   1/4  -1/4 |
| -1/2   1/2    1    1/2  -1/2 |
| -1/4   1/4   1/2   1/4  -1/4 |
|  1/4  -1/4  -1/2  -1/4   1/4 | (End)

T. D. Noe, Table of n, a(n) for n = 0..100
Mathematics Stack Exchange, On the number of Sylow subgroups in Symmetric Group
Zhujun Zhang, A Note on Counting Binomial Heaps, ResearchGate, June 2019.

Numerators of 2^n/n! give A001316. Cf. A000680, A008977, A139541.
Factor of A160481. - Johannes W. Meijer, May 24 2009
Equals A003148 divided by A123746. - Johannes W. Meijer, Nov 23 2009
Different from A160624.
Cf. A011371.

[ Denominator(2^n/Factorial(n)): n in [0..25] ]; // Klaus Brockhaus, Mar 10 2011

f:= n-> n! * 2^(add(i,i=convert(n,base,2))-n); # Peter Luschny, May 02 2009
seq (denom (coeff (series(1/(tanh(t)-1), t, 30), t, n)), n=0..25); # Peter Luschny, Aug 04 2011
seq(numer(n!/2^n), n=0..100); # Robert Israel, Jul 23 2015

Denominator[Table[(2^n)/n!,{n,0,40}]] (* Vladimir Joseph Stephan Orlovsky, Apr 03 2011*)
Table[Last[Select[Divisors[n!],OddQ]],{n,0,30}] (* Harvey P. Dale, Jul 24 2016 *)
Table[n!/2^IntegerExponent[n!,2], {n,1,30}] (* Clark Kimberling, Oct 22 2016 *)

A049606(n)=local(f=n!);f/2^valuation(f,2); \\ Joerg Arndt, Apr 22 2011
(Python 3.10+)
from math import factorial
def A049606(n): return factorial(n)>>n-n.bit_count() # Chai Wah Wu, Jul 11 2022

a(n) = Product_{k=1..n} A000265(k).
a(n) = A000265(A000142(n)). - Reinhard Zumkeller, Apr 09 2004
a(n) = numerator(2*Sum_{i>=1} (-1)^i*(1-zeta(n+i+1)) * (Product_{j=1..n} i+j)). - Gerry Martens, Mar 10 2011
a(n) = denominator([t^n] 1/(tanh(t)-1)). - Peter Luschny, Aug 04 2011
a(n) = n!/2^A011371(n). - Robert Israel, Jul 23 2015
From Zhujun Zhang, Jun 16 2019: (Start)
a(n) = n!/A060818(n).
E.g.f.: Product_{k>=0} (1 + x^(2^k) / 2^(2^k - 1)).
(End)
log a(n) = n log n - (1 + log 2)n + Θ(log n). - Charles R Greathouse IV, Feb 12 2022

New name (from Amarnath Murthy) by Charles R Greathouse IV, Jul 23 2015

A003148 a(n+1) = a(n) + 2n(2n+1)a(n-1), with a(0) = a(1) = 1.

Original entry on oeis.org

1, 1, 7, 27, 321, 2265, 37575, 390915, 8281665, 114610545, 2946939975, 51083368875, 1542234996225, 32192256321225, 1114841223671175, 27254953356505875, 1064057291370698625, 29845288035840902625, 1296073464766972266375, 41049997128507054562875

Offset: 0

N. J. A. Sloane

Numerators of sequence of fractions with e.g.f. 1/((1-x)*(1+x)^(1/2)). The denominators are successive powers of 2.
a(n) is the coefficient of x^n in arctan(sqrt(2*x/(1-x)))/sqrt(2*x*(1-x)) multiplied by (2*n+1)!!.
This sequence is the linking pin between the a(n) formulas of the ED1, ED2, ED3 and ED4 array rows, see A167552, A167565, A167580 and A167591. - Johannes W. Meijer, Nov 23 2009

			arctan(sqrt(2*x/(1-x)))/sqrt(2*x*(1-x)) = 1 + 1/3*x + 7/15*x^2 + 9/35*x^3 + ...

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 0..100
P. S. Bruckman, An interesting sequence of numbers derived from various generating functions, Fib. Quart., 10 (1972), 169-181.
R. J. Mathar, Numerical Representation of the Incomplete Gamma Function of Complex Argument, arXiv:math/0306184 [math.NA], 2003-2004; cf. Eq. 22.

Cf. A002943, A046161, A049606, A077568, A084543, A091520, A123746.

Cf. A167552, A167565, A167580, A167591.

a003148 n = a003148_list !! n
a003148_list = 1 : 1 : zipWith (+) (tail a003148_list)
                          (zipWith (*) (tail a002943_list) a003148_list)
-- Reinhard Zumkeller, Nov 22 2011

[n le 2 select 1 else Self(n-1) + 2*(n-2)*(2*n-3)*Self(n-2): n in [1..30]]; // G. C. Greubel, Nov 04 2022

# double factorial of odd "l" df := proc(l) local n; n := iquo(l,2); RETURN( factorial(l)/2^n/factorial(n)); end: x := 1; for n from 1 to 15 do if n mod 2 = 0 then x := 2*n*x+df(2*n-1); else x := 2*n*x-df(2*n-1); fi; print(x); od; quit

a[n_] := a[n] = (-1)^n*(2n - 1)!! + 2n*a[n - 1]; a[0] = 1; Table[ a[n], {n, 0, 14}] (* Jean-François Alcover, Dec 01 2011, after R. J. Mathar *)
a[ n_] := If[ n < 0, 0, (2 n + 1)!! Hypergeometric2F1[ -n, 1/2, 3/2, 2]]; (* Michael Somos, Apr 20 2018 *)
a[ n_] := If[ n < 0, 0, n! SeriesCoefficient[ 1 / ((1 - 2 x) Sqrt[1 + 2 x]), {x, 0, n}]]; (* Michael Somos, Apr 20 2018 *)
RecurrenceTable[{a[0]==a[1]==1,a[n+1]==a[n]+2n(2n+1)a[n-1]},a,{n,20}] (* Harvey P. Dale, Jul 27 2019 *)

Vec(serlaplace(1/(sqrt(1+2*x + O(x^20))*(1-2*x)))) \\ Andrew Howroyd, Feb 05 2018

@CachedFunction
def a(n): return 1 if (n<2) else a(n-1) + 2*(n-1)*(2*n-1)*a(n-2) # a = A003148
[a(n) for n in range(31)] # G. C. Greubel, Nov 04 2022

a(n) = (-1)^n*(2n-1)!! + 2*n*a(n-1) with (2n-1)!! = 1*3*5*..*(2n-1) the double factorial. - R. J. Mathar, Jun 12 2003
a(n) = ((2*n+1)!!/4) * Integral_{-Pi..Pi} cos(x)^n * cos(x/2) dx. - R. J. Mathar, Jun 30 2003
a(n) = (2n+1)!! 2F1(-n, 1/2;3/2;2). - R. J. Mathar, Jun 30 2003
In terms of the (terminating) Gauss hypergeometric function/series, 2F1(., .; .; 2), a(n) is a special case of the family of integer sequences defined by a(m, n) = ((2*n+2*m+1)!!/(2*m+1)) * 2F1(-n, m+1/2; m+3/2; 2), for m >= 0, n >= 0. An integral form can be seen as a(m, n) = ((2*n+2*m+1)!!/4) * Integral_{-Pi..Pi} (sin(x/2))^(2*m) * (cos(x))^n * cos(x/2) dx. A recurrence property is 4*(n+1)*a(m, n) = (2*m-1)*a(m-1, n+1) + (-1)^n*(2*n+2*m+1)!!. Sequences that have these properties are a(0, n) = this sequence, a(1, n) = A077568, a(2, n) = A084543. - R. J. Mathar, Jun 30 2003
E.g.f.: 1/(sqrt(1+2*x)*(1-2*x)). - Vladeta Jovovic, Oct 12 2003
a(n) = (2^n)*n!*A123746(n)/A046161(n) = (2^n)*n!*Sum_{k=0..n} binomial(2*k,k)*(-1/4)^k. From the e.g.f. - Wolfdieter Lang, Oct 06 2008
a(n) = A049606(n)*A123746(n). - Johannes W. Meijer, Nov 23 2009
a(n) = A091520(n) * n! / 2^n. - Michael Somos, Mar 17 2011

a(16)-a(20) from Andrew Howroyd, Feb 05 2018

cp's OEIS Frontend

A049606 Largest odd divisor of n!.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

Extensions

A003148 a(n+1) = a(n) + 2n(2n+1)a(n-1), with a(0) = a(1) = 1.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

SageMath

Formula

Extensions

cp's OEIS Frontend

A049606 Largest odd divisor of n!.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

Extensions

A003148 a(n+1) = a(n) + 2n*(2n+1)*a(n-1), with a(0) = a(1) = 1.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

SageMath

Formula

Extensions

A003148 a(n+1) = a(n) + 2n(2n+1)a(n-1), with a(0) = a(1) = 1.