login
A341108
a(n) = A163176(n+1)*A003557(n+1).
3
1, 1, 4, 4, 48, 16, 576, 576, 11520, 768, 9216, 3072, 3870720, 552960, 442368, 442368, 26542080, 8847360, 2229534720, 445906944, 70071091200, 6370099200, 76441190400, 25480396800, 83473779916800, 1284211998720, 5136847994880, 81537269760, 4892236185600
OFFSET
0,3
FORMULA
a(n) = A195441(n)*A341107(n).
PROG
(Sage)
def A341108(n):
def L(n, p, r):
s, q = 0, p - r
while q <= n:
s += n // q
q *= p
return s
if n < 2: return 1
p = prod(p^(L(n, p, 1) - L(n+1, p, 0)) for p in primes(n+1))
q = prod(p for p in prime_divisors(n + 1))
return ((n + 1) // q) * p
print([A341108(n) for n in (0..28)])
(Julia)
A341108(n) = A163176(n+1)*A003557(n+1)
[A341108(n) for n in 0:30] |> println
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Feb 07 2021
STATUS
approved