A322020 a(n) = Product_{d|n, gcd(d,n/d) is a prime power} gcd(d,n/d).

1, 1, 1, 2, 1, 1, 1, 4, 3, 1, 1, 4, 1, 1, 1, 16, 1, 9, 1, 4, 1, 1, 1, 16, 5, 1, 9, 4, 1, 1, 1, 64, 1, 1, 1, 36, 1, 1, 1, 16, 1, 1, 1, 4, 9, 1, 1, 256, 7, 25, 1, 4, 1, 81, 1, 16, 1, 1, 1, 16, 1, 1, 9, 512, 1, 1, 1, 4, 1, 1, 1, 144, 1, 1, 25, 4, 1, 1, 1, 256, 81, 1, 1, 16, 1, 1, 1, 16, 1, 81, 1, 4, 1, 1, 1, 4096, 1, 49, 9, 100, 1, 1, 1, 16, 1

Offset

1,4

Links

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

Formula

a(n) = Product_{d|n} gcd(d,n/d)^A069513(gcd(d,n/d)).

Prog

(PARI) A322020(n) = { my(m=1, p); fordiv(n, d, p = gcd(d, n/d); if(isprimepower(p), m *= p)); m; };

Crossrefs

Cf. A069513, A294876, A295666.

Sequence in context: A087653 A295666 A355003 * A294895 A285328 A321030

Adjacent sequences: A322017 A322018 A322019 * A322021 A322022 A322023

Keyword

nonn

Author

Antti Karttunen, Nov 28 2018

Status

approved