OFFSET
1,4
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..1000
MATHEMATICA
GCD @@@ Partition[#, 2, 1] & [Array[Product[Quotient[#, k], {k, #}] &, 30]] (* Paolo Xausa, Feb 11 2025 *)
PROG
(Haskell)
a208448 n = a208448_list !! (n-1)
a208448_list = zipWith gcd a010786_list $ tail a010786_list
(PARI) f(n) = prod(k=1, n, n\k); \\ A010786
a(n) = gcd(f(n+1), f(n)); \\ Michel Marcus, Feb 03 2022
(Python)
from math import prod
from collections import Counter
from sympy import factorint
def A208448(n):
def f(n):
c, j = Counter(), 1
while j<=n:
k = n//j
a = -j+(j:=n//k+1)
c += {p:e*a for p, e in factorint(k).items()}
return c
x, y = f(n), f(n+1)
return prod(p**min(x[p], y[p]) for p in x) # Chai Wah Wu, May 06 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Feb 26 2012
STATUS
approved
