[go: up one dir, main page]

login
A327154
a(n) = Product_{d|n, d>1} A008578(1+A286561(sigma(n),d)), where A286561(n,d) gives the highest exponent of d dividing n.
4
1, 1, 1, 1, 1, 12, 1, 1, 1, 2, 1, 6, 1, 5, 2, 1, 1, 2, 1, 2, 1, 3, 1, 48, 1, 2, 1, 80, 1, 45, 1, 1, 2, 2, 1, 1, 1, 3, 1, 8, 1, 44, 1, 6, 2, 5, 1, 6, 1, 1, 3, 2, 1, 20, 1, 20, 1, 2, 1, 80, 1, 11, 1, 1, 1, 63, 1, 2, 2, 7, 1, 2, 1, 2, 1, 6, 1, 20, 1, 2, 1, 2, 1, 264, 1, 3, 2, 6, 1, 48, 2, 10, 1, 7, 2, 108, 1, 1, 2, 1, 1, 125, 1, 2, 2
OFFSET
1,6
FORMULA
a(n) = Product_{d|n, d>1} A008578(1+A286561(sigma(n),d)), where sigma = A000203.
Other identities. For all n >= 1:
1+A001222(a(n)) = A073802(n).
PROG
(PARI) A327154(n) = { my(m=1, s=sigma(n), v); fordiv(n, d, if((d>1) && ((v = valuation(s, d))>0), m *= prime(v))); (m); };
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Sep 18 2019
STATUS
approved