cp's OEIS Frontend

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

%I A341108 #7 Feb 18 2021 15:03:14
%S A341108 1,1,4,4,48,16,576,576,11520,768,9216,3072,3870720,552960,442368,
%T A341108 442368,26542080,8847360,2229534720,445906944,70071091200,6370099200,
%U A341108 76441190400,25480396800,83473779916800,1284211998720,5136847994880,81537269760,4892236185600
%N A341108 a(n) = A163176(n+1)*A003557(n+1).
%F A341108 a(n) = A195441(n)*A341107(n).
%o A341108 (Sage)
%o A341108 def A341108(n):
%o A341108     def L(n, p, r):
%o A341108         s, q = 0, p - r
%o A341108         while q <= n:
%o A341108             s += n // q
%o A341108             q *= p
%o A341108         return s
%o A341108     if n < 2: return 1
%o A341108     p = prod(p^(L(n, p, 1) - L(n+1, p, 0)) for p in primes(n+1))
%o A341108     q = prod(p for p in prime_divisors(n + 1))
%o A341108     return ((n + 1) // q) * p
%o A341108 print([A341108(n) for n in (0..28)])
%o A341108 (Julia)
%o A341108 A341108(n) = A163176(n+1)*A003557(n+1)
%o A341108 [A341108(n) for n in 0:30] |> println
%Y A341108 Cf. A163176, A003557, A195441, A341107.
%K A341108 nonn
%O A341108 0,3
%A A341108 _Peter Luschny_, Feb 07 2021