close
login
A086410
Smallest prime factor of 3-smooth numbers, with a(1)=1.
6
1, 2, 3, 2, 2, 2, 3, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2
OFFSET
1,2
LINKS
FORMULA
a(n) = A020639(A003586(n));
a(n) <= A086411(n) <= 3.
a(A033845(n)) = A086411(A033845(n))-1; a(A006899(n)) = A086411(A006899(n)). - Reinhard Zumkeller, Sep 25 2008
MATHEMATICA
Reap[For[n = 1, n <= 2*10^5, n++, If[EulerPhi[6*n] == 2*n, Sow[ FactorInteger[n][[1, 1]]]]]][[2, 1]] (* Jean-François Alcover, Sep 02 2016 *)
PROG
(Python)
from sympy import integer_log
from oeis_sequences.OEISsequences import bisection
def A086410(n):
if n == 1: return 1
def f(x): return n+x-sum((x//3**i).bit_length() for i in range(integer_log(x, 3)[0]+1))
return 2|(bisection(f, n, n)&1) # Chai Wah Wu, Mar 12 2026
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
Reinhard Zumkeller, Jul 18 2003
STATUS
approved