|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
FORMULA
|
a(1) = 1; for n > 1, a(n) = Product_{d|n} A019565(d)^[A008683(d) > 0].
|
|
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
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|