A320017 a(1) = 1; for n > 1, a(n) = Product_{d|n} A019565(d)^[moebius(d) = +1].

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

Antti Karttunen, Table of n, a(n) for n = 1..8192

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

Adjacent sequences: A320014 A320015 A320016 * A320018 A320019 A320020

Keyword

nonn

Author

Antti Karttunen, Nov 08 2018

Status

approved