login
A365426
a(n) = A336466(A163511(n)).
5
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 9, 1, 1, 1, 3, 1, 3, 3, 5, 1, 1, 1, 1, 1, 1, 1, 27, 1, 1, 1, 9, 1, 9, 9, 25, 1, 1, 1, 3, 1, 3, 3, 15, 1, 3, 3, 5, 3, 5, 5, 3, 1, 1, 1, 1, 1, 1, 1, 81, 1, 1, 1, 27, 1, 27, 27, 125, 1, 1, 1, 9, 1, 9, 9, 75, 1, 9, 9, 25, 9, 25, 25, 9, 1, 1, 1, 3, 1
OFFSET
0,16
LINKS
FORMULA
For all n >= 1, a(n) = a(2*n) = a(A000265(n)).
PROG
(PARI)
A000265(n) = (n>>valuation(n, 2));
A163511(n) = if(!n, 1, my(p=2, t=1); while(n>1, if(!(n%2), (t*=p), p=nextprime(1+p)); n >>= 1); (t*p));
A336466(n) = { my(f=factor(n)); prod(k=1, #f~, A000265(f[k, 1]-1)^f[k, 2]); };
CROSSREFS
Cf. also A365427.
Sequence in context: A204131 A317941 A318667 * A030598 A030395 A031243
KEYWORD
nonn
AUTHOR
Antti Karttunen, Sep 03 2023
STATUS
approved