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

Author

Antti Karttunen, Nov 08 2018

Keywords

• nonn

Links

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

Crossrefs

Cf. A019565, A320018 (rgs-transform).

Cf. also A300831, A300832.

Programs

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));

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).