|
|
A091137
|
|
The Hirzebruch numbers. a(n) = Product_{2 <= p <= n+1, p prime} p^floor(n / (p - 1)).
|
|
41
|
|
|
1, 2, 12, 24, 720, 1440, 60480, 120960, 3628800, 7257600, 479001600, 958003200, 2615348736000, 5230697472000, 31384184832000, 62768369664000, 32011868528640000, 64023737057280000, 51090942171709440000, 102181884343418880000, 33720021833328230400000, 67440043666656460800000
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
Largest number m such that number of times m divides k! is almost k/n for large k, i.e., largest m with A090624(m) = n.
This is always a relatively small multiple of n!, since the multiplicity with which a prime p divides n! is always <= n/(p-1); it is equal to floor(n/(p-1)) at least when n is a power of p. - Franklin T. Adams-Watters, May 31 2010
At least for most small n, a(n) = A002790(n) * n!; the first difference is n=15. It appears that A002790(n) * n! always divides a(n).
Conjecture: The denominators of the series reversion of the sequence with e.g.f. Polylog(2,x). - Benedict W. J. Irwin, Jan 05 2017
|
|
REFERENCES
|
P. Curtz, Integration numérique ..., Note 12, C.C.S.A., Arcueil, 1969; see pp. 36, 56.
F. Hirzebruch, Topological Methods in Algebraic Geometry, Springer, 3rd. ed., 1966; Lemma 1.7.3, p. 14. [From N. J. A. Sloane, Sep 06 2010]
|
|
LINKS
|
|
|
FORMULA
|
a(n) = Product_p {p prime} p^floor(n/(p-1)).
a(2n+1) = 2*a(2n).
a(n) = lcm_{p in P(n)} Product_{r in p}(r + 1), where P(n) are the partitions of n.
a(n) = lcm(A238963row(n)).
a(n) = A368116(1, n), seen as the lcm of the product of the 1-shifted partitions.
a(n) = A368093(1, n), seen as the cumulative product of the Clausen numbers A160014(1, n). (End)
|
|
EXAMPLE
|
Let n = 4. The partitions of 4 are [[4], [3, 1], [2, 2], [2, 1, 1], [1, 1, 1, 1]]. Thus a(4) = lcm([5, 4*2, 3*3, 3*2*2, 2*2*2*2]) = 720.
|
|
MAPLE
|
A091137 := proc(n) local a, i, p ; a := 1 ; for i from 1 do p := ithprime(i) ; if p > n+1 then break; fi; a := a*p^floor(n/(p-1)) ; od: a ; end: for n from 0 to 20 do printf("%a, ", A091137(n) ) ; od: # R. J. Mathar, Feb 23 2009
|
|
MATHEMATICA
|
A027760[n_] := Product[d, {d, Select[ Divisors[n] + 1, PrimeQ]}]; a[n_] := a[n] = A027760[n]*a[n-1]; a[0] = 1; Table[ a[n], {n, 0, 18}] (* Jean-François Alcover, Oct 04 2011 *)
|
|
PROG
|
(PARI) a(n) = local(r); r=1; forprime(p=2, n+1, r*=p^(n\(p-1))); r
(Python)
from math import prod
from sympy import primerange
def A091137(n): return prod(p**(n//(p-1)) for p in primerange(n+2))
(SageMath)
def a(n): return lcm(product(r + 1 for r in p) for p in Partitions(n))
# Or, more efficient:
from functools import cache
@cache
def a_rec(n):
if n == 0: return 1
p = mul(s for s in map(lambda i: i + 1, divisors(n)) if is_prime(s))
return p * a_rec(n - 1)
print([a_rec(n) for n in range(22)]) # Peter Luschny, Dec 12 2023
|
|
CROSSREFS
|
Starts similarly to A002207 especially for even n and all values of A002207 seen so far seem to divide a(n).
|
|
KEYWORD
|
nonn,changed
|
|
AUTHOR
|
|
|
EXTENSIONS
|
New name using a formula of the author by Peter Luschny, Dec 11 2023
|
|
STATUS
|
approved
|
|
|
|