A049096 Numbers k such that 2^k + 1 is divisible by a square > 1.
3, 9, 10, 15, 21, 27, 30, 33, 39, 45, 50, 51, 55, 57, 63, 68, 69, 70, 75, 78, 81, 87, 90, 93, 99, 105, 110, 111, 117, 123, 129, 130, 135, 141, 147, 150, 153, 159, 165, 170, 171, 177, 182, 183, 189, 190, 195, 201, 204, 207, 210, 213, 219, 225, 230, 231, 234, 237, 243
Offset: 1
Keywords
Examples
9 is here because 2^9 + 1 = 513 is divisible by 9. 99 is here because 2^99 + 1 = 3^3*19*67*683*5347*20857*242099935645987 is divisible by 9, i.e. is not squarefree.
Links
- Robert Israel, Table of n, a(n) for n = 1..10000
- Robert Israel, When p^2 divides 2^n + 1
Crossrefs
Programs
-
Magma
[n: n in [3..220] | not IsSquarefree(2^n+1)]; // Vincenzo Librandi, Mar 08 2018
-
Maple
remove(n -> numtheory:-issqrfree(2^n+1), [$1..250]); # Robert Israel, Jul 07 2014
-
Mathematica
Select[Range[243], !SquareFreeQ[2^# + 1] &] (* Vladimir Joseph Stephan Orlovsky, Mar 18 2011*)
-
PARI
is(n)=!issquarefree(2^n+1) \\ Altug Alkan, Nov 20 2015
Formula
For any a(n+1) - a(n) <= 6 since numbers of form 3^a*(2k+1) a > 0, k >= 0, are in the sequence (2^(3*(2k+1) + 1 is divisible by 9). So are numbers of the form 20k + 10 since 2^(20k+10) + 1 is divisible by 25, 110k + 55 since 2^(110k+55) + 1 is divisible by 11^2, 78 + 156k since 2^(156k+78) + 1 is divisible by 13^2 ... - Benoit Cloitre, Aug 20 2002
Extensions
More terms from James Sellers, Dec 16 1999
More terms from Vladeta Jovovic, Apr 12 2002
Missing term 182 added by Rainer Rosenthal, Nov 01 2005
Comments