[go: up one dir, main page]

login
A324283
Multiplicative with a(p^e) = A276086(p^e).
2
1, 3, 6, 9, 18, 18, 10, 15, 30, 54, 90, 54, 50, 30, 108, 225, 450, 90, 250, 162, 60, 270, 2250, 90, 1250, 150, 3750, 90, 11250, 324, 14, 21, 540, 1350, 180, 270, 70, 750, 300, 270, 630, 180, 350, 810, 540, 6750, 3150, 1350, 1750, 3750, 2700, 450, 15750, 11250, 1620, 150, 1500, 33750, 78750, 972, 98, 42, 300, 441, 900
OFFSET
1,2
PROG
(PARI)
A276086(n) = { my(i=0, m=1, pr=1, nextpr); while((n>0), i=i+1; nextpr = prime(i)*pr; if((n%nextpr), m*=(prime(i)^((n%nextpr)/pr)); n-=(n%nextpr)); pr=nextpr); m; };
A324283(n) = { my(f=factor(n)); prod(i=1, #f~, A276086(f[i, 1]^f[i, 2])); };
CROSSREFS
KEYWORD
nonn,base,mult
AUTHOR
Antti Karttunen, Feb 20 2019
STATUS
approved