A341108 a(n) = A163176(n+1)*A003557(n+1).

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

Author

Peter Luschny, Feb 07 2021

Keywords

• nonn

Crossrefs

Cf. A163176, A003557, A195441, A341107.

Programs

Julia

A341108(n) = A163176(n+1)*A003557(n+1)
[A341108(n) for n in 0:30] |> println

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)])

Formula

a(n) = A195441(n)*A341107(n).