OFFSET
0,3
PROG
(Python)
from sympy import primeomega
def A116452(n): return (1, 1, 2, 2, 1, 1, 2, 2, 1, 3, 1, 3, 2, 2, 1, 2, 2, 3, 3, 1, 3, 1, 2)[n] if n<23 else primeomega(n+1+(m:=(n+1).bit_length()+((n+1)//3).bit_length())+(n+1+m>=(1<<(m+1>>1) if m&1 else 3<<(m+1>>1)-1))) # Chai Wah Wu, Jan 25 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Feb 16 2006
STATUS
approved
