[go: up one dir, main page]

login
A324870 revision #4

A324870
a(n) = A324863(n) - A252464(n).
5
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1
OFFSET
1
COMMENTS
Question: Are there any other terms than 0's and 1's ?
A324871 gives the numbers n where a(n) <> 0. The first such number which is not a square is 187 = 11*17.
FORMULA
a(n) = A324863(n) - A252464(n).
PROG
(PARI)
A061395(n) = if(1==n, 0, primepi(vecmax(factor(n)[, 1])));
A252464(n) = if(1==n, 0, (bigomega(n) + A061395(n) - 1));
A324861(n) = #binary(A324876(n)); \\ Needs also code from A324876.
A324863(n) = { my(m=0, w, c=0); fordiv(n, d, w=A324861(d); if(w>=m, if(w==m, c++, c=1; m=w))); (m); };
A324870(n) = (A324863(n) - A252464(n));
CROSSREFS
Sequence in context: A023975 A305821 A305892 * A011738 A011737 A011736
KEYWORD
nonn
AUTHOR
Antti Karttunen, Mar 21 2019
STATUS
editing