A067535 - OEIS

%I #22 Feb 16 2025 08:32:45

%S 1,2,3,5,5,6,7,10,10,10,11,13,13,14,15,17,17,19,19,21,21,22,23,26,26,

%T 26,29,29,29,30,31,33,33,34,35,37,37,38,39,41,41,42,43,46,46,46,47,51,

%U 51,51,51,53,53,55,55,57,57,58,59,61,61,62

%N Smallest squarefree number >= n.

%H Michael De Vlieger, <a href="/A067535/b067535.txt">Table of n, a(n) for n = 1..10000</a>

%H Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/Squarefree.html">Squarefree</a>.

%F a(n) = n + A081221(n). - _Amiram Eldar_, Oct 10 2023

%p A067535 := proc(n)

%p     for a from n do

%p         if issqrfree(a) then

%p             return a ;

%p         end if;

%p     end do:

%p end proc:

%p seq(A067535(n),n=1..100) ; # _R. J. Mathar_, May 31 2024

%t Table[k = n; While[! SquareFreeQ@ k, k++]; k, {n, 62}] (* _Michael De Vlieger_, Mar 18 2017 *)

%o (PARI) a(n) = while (! issquarefree(n), n++); n; \\ _Michel Marcus_, Mar 18 2017

%o (Python)

%o from itertools import count

%o from sympy import factorint

%o def A067535(n): return next(m for m in count(n) if max(factorint(m).values(),default=0)<=1) # _Chai Wah Wu_, Dec 04 2024

%Y Cf. A005117, A070321, A076260, A081221.

%K nonn,changed

%O 1,2

%A _Reinhard Zumkeller_, Jan 27 2002