login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

Numbers n such that 10^n+1 is not squarefree.
9

%I #21 Nov 01 2024 09:35:24

%S 11,21,33,39,55,63,77,99,105,117,121,136,143,147,165,171,187,189,195,

%T 202,209,231,243,253,273,275,292,297,315,319,341,351,357,363,385,399,

%U 406,407,408,429,441,451,473,483,495,507,513,517,525,539,548,561,567

%N Numbers n such that 10^n+1 is not squarefree.

%C This sequence is the union of the collection of sequences formed from the nonzero terms of A086981 * A005408, the odd numbers. First occurrence of consecutive integers in sequence is 406,407,408.

%C From _Robert Israel_, Feb 13 2017: (Start)

%C Numbers n such that gcd(n, 10^n + 1) > 1 or n = k*m where k is odd and 2*m is the order of 10 modulo a member of A045616. [Corrected by _Jianing Song_, Nov 01 2024]

%C If n is in the sequence, then so is k*n for any odd k. (End)

%C Numbers of the form k*ord(10,p^2)/2, where k is an odd number and p is a prime such that ord(10,p) is even. Here ord(a,m) is the multiplicative order of a modulo m. Note that if p is not in A045616, then ord(10,p^2) = p*ord(10,p). - _Jianing Song_, Nov 01 2024

%H Robert Israel, <a href="/A086982/b086982.txt">Table of n, a(n) for n = 1..10000</a>

%p filter:= n -> (n mod 243 = 0 and (n/243)::odd) or igcd(n,(10 &^n +1 mod n)) > 1: # Note that this works if n < 28299156

%p select(filter, [$1..1000]); # _Robert Israel_, Feb 13 2017

%o (PARI) ord = [1, 486, 56598312]; \\ order of 10 modulo A045616

%o isA086982(n) = if(gcd(n, 10^n+1) > 1, return(1)); for(i=1, 3, if((ord[i] % 2 == 0) && (n % (ord[i]/2) == 0) && (n/(ord[i]/2) % 2 == 1), return(1))); return(0) \\ _Jianing Song_, Nov 01 2024, after _Robert Israel_'s comment; considering only the three currently-known terms of A045616

%Y Cf. A005408, A013929, A086981, A045616, A049096, A282344, A282345.

%K nonn

%O 1,1

%A _Ray Chandler_, Jul 27 2003