login
A320017
a(1) = 1; for n > 1, a(n) = Product_{d|n} A019565(d)^[moebius(d) = +1].
3
1, 2, 2, 2, 2, 30, 2, 2, 2, 42, 2, 30, 2, 210, 420, 2, 2, 30, 2, 42, 220, 330, 2, 30, 2, 462, 2, 210, 2, 132300, 2, 2, 52, 78, 156, 30, 2, 390, 780, 42, 2, 346500, 2, 330, 420, 2730, 2, 30, 2, 42, 1716, 462, 2, 30, 8580, 210, 4004, 6006, 2, 132300, 2, 30030, 220, 2, 68, 128700, 2, 78, 340, 343980, 2, 30, 2, 714, 420, 390, 2380, 2702700, 2, 42, 2
OFFSET
1,2
LINKS
FORMULA
a(1) = 1; for n > 1, a(n) = Product_{d|n} A019565(d)^[A008683(d) > 0].
For n >= 2, A048675(a(n)) = A318674(n).
PROG
(PARI)
A019565(n) = {my(j, v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ From A019565
A320017(n) = if(1==n, n, my(m=1); fordiv(n, d, if(1==moebius(d), m *= A019565(d))); (m));
CROSSREFS
Cf. A019565, A320018 (rgs-transform).
Cf. also A300831, A300832.
Sequence in context: A075182 A292023 A286851 * A289087 A154288 A225057
KEYWORD
nonn
AUTHOR
Antti Karttunen, Nov 08 2018
STATUS
approved