login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A289623
a(n) = A055396(A048673(n)).
1
0, 1, 2, 3, 1, 1, 1, 1, 6, 5, 4, 9, 2, 7, 1, 13, 1, 1, 1, 1, 1, 1, 2, 1, 3, 1, 2, 1, 1, 16, 8, 1, 2, 10, 2, 30, 2, 3, 14, 3, 1, 23, 1, 17, 1, 1, 2, 4, 18, 1, 1, 4, 1, 1, 1, 35, 1, 15, 11, 1, 1, 1, 1, 3, 1, 1, 1, 1, 21, 1, 12, 1, 1, 1, 2, 1, 1, 1, 1, 1, 65, 3, 2, 1, 19, 20, 1, 1, 4, 56, 1, 32, 2, 1, 2, 1, 2, 1, 38, 6
OFFSET
1,3
COMMENTS
From the scatter plot it can be seen that the terms are grouped into two distinct populations by their magnitude, with significant gap between them.
LINKS
FORMULA
a(n) = A055396(A048673(n)).
PROG
(PARI)
A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ This function from Michel Marcus
A048673(n) = (A003961(n)+1)/2;
A055396(n) = if(n==1, 0, primepi(factor(n)[1, 1])); \\ This function from Charles R Greathouse IV, Apr 23 2015
CROSSREFS
Cf. A048673, A055396, A246263 (the positions of ones).
Sequence in context: A308638 A356320 A207572 * A086195 A086197 A139336
KEYWORD
nonn,look
AUTHOR
Antti Karttunen, Jul 16 2017
STATUS
approved