%I #30 Sep 08 2022 08:44:58
%S 1,2,3,4,5,7,8,9,10,11,13,14,15,16,17,19,22,23,25,26,27,28,29,31,32,
%T 33,34,35,37,38,39,41,43,44,45,46,47,49,50,51,52,53,55,56,57,58,59,61,
%U 62,64,65,67,68,69,70,71,73,74,75,76,77,79,81,82,83,85,86,87,88,89,91,92
%N Numbers n such that 2^n - 1 is squarefree.
%C Numbers n such that gcd(n, 2^n - 1) = 1 and n is not a multiple of A002326((q - 1)/2), where q is a Wieferich prime A001220. - _Thomas Ordowski_, Nov 21 2015
%C If n is in the sequence, then so are all divisors of n. - _Robert Israel_, Nov 23 2015
%H Max Alekseyev, <a href="/A049093/b049093.txt">Table of n, a(n) for n = 1..910</a>
%e a(7) = 8 because 2^8 - 1 = 255 = 3 * 5 * 17 is squarefree.
%p N:= 400: # to get all terms <= N
%p # This relies on the fact that the first N+1 members of A000225 have all been factored
%p # without any further Wieferich primes being found.
%p V:= Vector(N,1):
%p V[364 * [$1..N/364]]:= 0:
%p V[1755 * [$1..N/1755]]:= 0:
%p for n from 2 to N do
%p if V[n] = 0 then next fi;
%p if igcd(n, 2 &^n - 1 mod n) > 1 then
%p V[n * [$1..N/n]]:= 0
%p fi;
%p od:
%p select(t -> V[t] = 1, [$1..N]); # _Robert Israel_, Nov 23 2015
%t Select[Range@ 92, SquareFreeQ[2^# - 1] &] (* _Michael De Vlieger_, Nov 21 2015 *)
%o (PARI) isok(n) = issquarefree(2^n - 1); \\ _Michel Marcus_, Dec 19 2013
%o (Magma) [n: n in [1..100] | IsSquarefree(2^n-1)]; // _Vincenzo Librandi_, Nov 22 2015
%Y Complement of A049094.
%K nonn
%O 1,2
%A _Labos Elemer_
%E Terms a(73)-a(910) in b-file from _Max Alekseyev_, Nov 15 2014, Sep 28 2015