login
A381588
If n = Product (p_j^k_j) then a(n) = Product (lcm(p_j, k_j)), with a(1) = 1.
0
1, 2, 3, 2, 5, 6, 7, 6, 6, 10, 11, 6, 13, 14, 15, 4, 17, 12, 19, 10, 21, 22, 23, 18, 10, 26, 3, 14, 29, 30, 31, 10, 33, 34, 35, 12, 37, 38, 39, 30, 41, 42, 43, 22, 30, 46, 47, 12, 14, 20, 51, 26, 53, 6, 55, 42, 57, 58, 59, 30, 61, 62, 42, 6, 65, 66, 67, 34, 69, 70, 71
OFFSET
1,2
FORMULA
a(p) = p, for p prime.
EXAMPLE
a(18) = 12 because 18 = 2^1*3^2, lcm(2,1) = 2, lcm(3,2) = 6 and 2*6 = 12.
a(300) = 30 because 300 = 2^2*3^1*5^2, lcm(2,2) = 2, lcm(3,1) = 3, lcm(5,2) = 10 and 2*3*10 = 60.
MATHEMATICA
A381588[n_] := Times @@ LCM @@@ FactorInteger[n];
Array[A381588, 100]
PROG
(PARI) a(n) = my(f=factor(n)); prod(i=1, #f~, lcm(f[i, 1], f[i, 2])); \\ Michel Marcus, Mar 02 2025
CROSSREFS
Cf. A008473, A008477, A035306, A144338 (fixed points), A369008 (analogous for gcd).
Sequence in context: A308056 A336965 A381203 * A381201 A293303 A333569
KEYWORD
nonn,mult,easy,new
AUTHOR
Paolo Xausa, Feb 28 2025
STATUS
approved