A055634 2-adic factorial function.

1, -1, 1, -3, 3, -15, 15, -105, 105, -945, 945, -10395, 10395, -135135, 135135, -2027025, 2027025, -34459425, 34459425, -654729075, 654729075, -13749310575, 13749310575, -316234143225, 316234143225, -7905853580625, 7905853580625, -213458046676875, 213458046676875

Offset

0,4

Comments

Also known as Morita's 2-adic gamma function. - Harry Richman, Jul 26 2023

References

Serge Lang, Cyclotomic Fields I and II, Springer-Verlag, 1990, p. 315.

Links

Kenny Lau, Table of n, a(n) for n = 0..806

Daniel Barsky, On Morita's p-adic Gamma function, Groupe de travail d'analyse ultramétrique, 5 (1977-1978), Talk no. 3, 6 p.

Wikipedia, P-adic gamma function.

Formula

a(2*n) = -a(2*n - 1) = (2*n - 1)!!

a(n) = (-1)^n*n!/A037223(n), A037223(n) = 2^floor(n/2)*floor(n/2)!. Exponential generating function: (1-x)*exp(x^2/2). - Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Nov 01 2002

Mathematica

a[ n_] := If[ n < 0, 0, n! (-1)^n / (n - Mod[n, 2])!!]; (* Michael Somos, Jun 30 2018 *)
4[ n_] := If[ n < 0, 0, n! SeriesCoefficient[ (1 - x) Exp[x^2/2], {x, 0, n}]]; (* Michael Somos, Jun 30 2018 *)

Prog

(PARI) {a(n) = if( n<1, 1, -if( n%2, n * a(n-1), a(n-1)))};
(PARI) a(n)=(-1)^n*(n=bitor(n-1, 1))!/(n\2)!>>(n\2) \\ Charles R Greathouse IV, Oct 01 2012
(Sage)
def Gauss_factorial(N, n): return mul(j for j in (1..N) if gcd(j, n) == 1)
def A055634(n): return (-1)^n*Gauss_factorial(n, 2)
[A055634(n) for n in (0..28)]  # Peter Luschny, Oct 01 2012
(Magma) /* Based on Gauss factorial n_2!: */ k:=2; [IsZero(n) select 1 else (-1)^n*&*[j: j in [1..n] | IsOne(GCD(j, k))]: n in [0..30]]; // Bruno Berselli, Dec 10 2013

Crossrefs

Cf. A000142, A006882, A001147, A133221.

Sequence in context: A290344 A217858 A185275 * A133221 A232097 A353584

Adjacent sequences: A055631 A055632 A055633 * A055635 A055636 A055637

Keyword

sign

Author

Michael Somos, Jun 06 2000

Status

approved