cp's OEIS Frontend

A370486 Terms in A116260 that are not of the form 10^d-1.

%I A370486 #15 Feb 20 2024 08:24:29
%S A370486 36363636362,45454545453,54545454544,63636363635,72727272726,
%T A370486 81818181817,90909090908,428571428571428571427,571428571428571428570,
%U A370486 714285714285714285713,857142857142857142856,363636363636363636363636363636362,454545454545454545454545454545453
%N A370486 Terms in A116260 that are not of the form 10^d-1.
%C A370486 Since 10^d-1 for d>0 are terms of A116260, this sequence only includes terms of A116260 not of this form to highlight these nontrivial terms.
%H A370486 Chai Wah Wu, <a href="/A370486/b370486.txt">Table of n, a(n) for n = 1..44</a>
%o A370486 (Python)
%o A370486 from itertools import count, islice
%o A370486 from sympy import sqrt_mod
%o A370486 def A370486_gen(): # generator of terms
%o A370486     for j in count(0):
%o A370486         b = 10**j
%o A370486         a = b*10+1
%o A370486         for k in sorted(sqrt_mod(0,a,all_roots=True)):
%o A370486             if a*(b+4) <= k**2 < a*(a+3):
%o A370486                 yield k-2
%o A370486 A370486_list = list(islice(A370486_gen(),20))
%Y A370486 Cf. A116260.
%K A370486 nonn,base
%O A370486 1,1
%A A370486 _Chai Wah Wu_, Feb 19 2024