[go: up one dir, main page]

login
A339970
a(n) = A329697(A019565(2n)).
3
0, 1, 1, 2, 2, 3, 3, 4, 2, 3, 3, 4, 4, 5, 5, 6, 2, 3, 3, 4, 4, 5, 5, 6, 4, 5, 5, 6, 6, 7, 7, 8, 1, 2, 2, 3, 3, 4, 4, 5, 3, 4, 4, 5, 5, 6, 6, 7, 3, 4, 4, 5, 5, 6, 6, 7, 5, 6, 6, 7, 7, 8, 8, 9, 3, 4, 4, 5, 5, 6, 6, 7, 5, 6, 6, 7, 7, 8, 8, 9, 5, 6, 6, 7, 7, 8, 8, 9, 7, 8, 8, 9, 9, 10, 10, 11, 4, 5, 5, 6, 6, 7, 7, 8, 6, 7
OFFSET
0,4
FORMULA
a(n) = A329697(A019565(2*n)) = A329697(A003961(A019565(n))).
a(n) = A329697(A339971(n)) + A000120(n).
If 4n = 2^e1 + 2^e2 + ... + 2^ek [e1 ... ek distinct], then a(n) = A329697(A000040(e1)) + A329697(A000040(e2)) + ... + A329697(A000040(ek)).
PROG
(PARI)
A019565(n) = { my(m=1, p=1); while(n>0, p = nextprime(1+p); if(n%2, m *= p); n >>= 1); (m); };
A329697(n) = if(!bitand(n, n-1), 0, 1+A329697(n-(n/vecmax(factor(n)[, 1]))));
A339970(n) = A329697(A019565(2*n));
(PARI) A339970(n) = { my(s=0, p=2); while(n>0, p = nextprime(1+p); if(n%2, s += A329697(p)); n >>= 1); (s); };
CROSSREFS
Differs from A106486 for the first time at n=32, where a(32) = 1, while A106486(32) = 2.
Sequence in context: A260651 A116370 A261018 * A106486 A195743 A106494
KEYWORD
nonn
AUTHOR
Antti Karttunen, Dec 27 2020
STATUS
approved