OFFSET
1,1
COMMENTS
Let p be the smallest prime that is coprime to n and let q be the second smallest prime that is coprime to n. Then a(n) = p*q.
Records in this sequence are set by n in A002110.
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
FORMULA
EXAMPLE
a(1) = 6 = 2*3, since p = 2, q = 3.
a(2) = 15 = 3*5, since p = 3, q = 5.
a(3) = 10 = 2*5, since p = 2, q = 5.
a(4) = 15 = 3*5, since p = 3, q = 5, a(2^i) = 15 for i > 0.
a(6) = 35 = 5*7, since p = 5, q = 7.
a(9) = 20 = 2*5, since p = 2, q = 5, a(3^i) = 10 for i > 0.
a(10) = 21 = 3*7, since p = 3, q = 7.
a(12) = 35 = 5*7, since p = 5, q = 7, a(k) = 35 for n in A033845 (i.e., n such that rad(n) = 6).
a(20) = 21 = 3*7, since p = 3, q = 7, a(k) = 21 for n in A033846 (i.e., n such that rad(n) = 10).
a(30) = 77 = 7*11, since p = 7, q = 11, etc.
MATHEMATICA
Table[c = 0; q = 2; Times @@ Reap[While[c < 2, While[Divisible[n, q], q = NextPrime[q]]; Sow[q]; q = NextPrime[q]; c++] ][[-1, 1]], {n, 120}]
PROG
(PARI) a(n) = my(k=2); while (isprime(k) || !issquarefree(k) || (gcd(k, n) != 1) , k++); k; \\ Michel Marcus, Apr 01 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michael De Vlieger, Mar 31 2025
STATUS
approved
