OFFSET
1,2
COMMENTS
a(p-1) = 2*A000142(p-2) for odd primes.
a(n) = A060593((n - 1)/2) for odd n. - Gregory Gerard Wojnar, Jun 10 2021
LINKS
Robert Israel, Table of n, a(n) for n = 1..451
MAPLE
seq(denom((n+1)/(2*(n-1)!)), n=1..25); # Robert Israel, Oct 14 2018
MATHEMATICA
Table[Denominator[(n + 1) / (2 (n - 1)!)], {n, 25}] (* Vincenzo Librandi, Oct 15 2018 *)
PROG
(Magma) [Denominator(n + 1) / (2*Factorial(n - 1)): n in [1..30]]; // Vincenzo Librandi, Oct 15 2018
(Python)
from math import factorial
from sympy import isprime
def A090586(n): return (factorial(n)<<1)//n if n>1 and isprime(n+1) else (factorial(n)<<1)//(n*(n+1)) # Chai Wah Wu, Jun 20 2026
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Reinhard Zumkeller, Dec 03 2003
STATUS
approved
