login
A364571
a(n) = A297171(A163511(n)), where A297171 is the Möbius transform of the inverse permutation of A163511.
2
0, 1, 1, 3, 2, 2, 2, 7, 4, 4, 2, 4, 5, 3, 6, 15, 8, 8, 4, 8, 4, 4, 4, 8, 10, 10, 4, 5, 13, 11, 14, 31, 16, 16, 8, 16, 8, 8, 8, 16, 8, 8, 4, 8, 8, 8, 8, 16, 20, 20, 10, 20, 7, 9, 6, 9, 26, 26, 12, 21, 29, 27, 30, 63, 32, 32, 16, 32, 16, 16, 16, 32, 16, 16, 8, 16, 16, 16, 16, 32, 16, 16, 8, 16, 8, 8, 8, 16, 16, 16
OFFSET
0,4
FORMULA
a(n) = A297171(A163511(n)).
PROG
(PARI)
A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); t };
A054429(n) = ((3<<#binary(n\2))-n-1); \\ From A054429
A163511(n) = if(!n, 1, A005940(1+A054429(n)));
A243071(n) = if(n<=2, n-1, my(f=factor(n), p, p2=1, res=0); for(i=1, #f~, p = 1 << (primepi(f[i, 1]) - 1); res += (p*p2*(2^(f[i, 2]) - 1)); p2 <<= f[i, 2]); ((3<<#binary(res\2))-res-1)); \\ (Combining programs given in A156552 and A054429) - Antti Karttunen, Aug 05 2023
A297171(n) = sumdiv(n, d, moebius(n/d)*A243071(d));
CROSSREFS
Cf. also A364567.
Sequence in context: A130845 A134653 A090207 * A202538 A239736 A065437
KEYWORD
nonn,look
AUTHOR
Antti Karttunen, Aug 05 2023
STATUS
approved