A300831 a(n) = Product_{d|n, d<n} A019565(d)^[moebius(n/d) = +1].

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 2, 2, 1, 1, 6, 1, 3, 2, 2, 1, 5, 1, 2, 1, 3, 1, 180, 1, 1, 2, 2, 2, 15, 1, 2, 2, 5, 1, 540, 1, 3, 6, 2, 1, 7, 1, 10, 2, 3, 1, 14, 2, 5, 2, 2, 1, 1575, 1, 2, 6, 1, 2, 756, 1, 3, 2, 900, 1, 35, 1, 2, 10, 3, 2, 1260, 1, 7, 1, 2, 1, 7875, 2, 2, 2, 5, 1, 44100, 2, 3, 2, 2, 2, 11, 1, 30, 6, 21, 1, 396, 1, 5, 1800

Offset

1,6

Links

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

Formula

a(n) = A293214(n) / (A300830(n)*A300832(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
A300831(n) = { my(m=1); fordiv(n, d, if((d < n)&&(1==moebius(n/d)), m *= A019565(d))); m; };

Crossrefs

Cf. A008683, A008966, A019565.

Cf. also A293214, A300830, A300832, A300833.

Sequence in context: A285572 A179926 A066882 * A068347 A284556 A025865

Adjacent sequences: A300828 A300829 A300830 * A300832 A300833 A300834

Keyword

nonn

Author

Antti Karttunen, Mar 16 2018

Status

approved