login
Numbers m such that the largest digit in the decimal expansion of 1/m is 6.
6

%I #18 Feb 17 2022 11:41:02

%S 6,15,16,24,39,60,64,88,96,150,156,160,165,219,240,246,273,275,375,

%T 378,384,390,399,462,600,606,615,624,625,640,792,822,858,880,888,956,

%U 960,975,984,1500,1515,1536,1554,1560,1584,1596,1600,1606,1626,1628,1638,1650,1665,1776,2145

%N Numbers m such that the largest digit in the decimal expansion of 1/m is 6.

%C If k is a term, 10*k is also a term.

%C First few primitive terms are 6, 15, 16, 24, 39, 64, 88, 96, 156, 165, ...

%C There is no prime up to 2.6*10^8 (see comments in A333237).

%C Subsequence: {6, 606, 60606, ...} = 6 * A094028.

%e 1/6 = 0.166666..., and 6 is the smallest number m such that the largest digit in the decimal expansion of 1/m is 6, so a(1) = 6.

%e As 1/39 = 0.025641025641..., 39 is a term.

%t f[n_] := Union[ Flatten[ RealDigits[ 1/n][[1]] ]]; Select[Range@1500000, Max@ f@# == 6 &]

%o (Python)

%o from itertools import count, islice

%o from sympy import n_order, multiplicity

%o def A351472_gen(startvalue=1): # generator of terms >= startvalue

%o for m in count(max(startvalue, 1)):

%o m2, m5 = multiplicity(2, m), multiplicity(5, m)

%o if max(str(10**(max(m2, m5)+n_order(10, m//2**m2//5**m5))//m)) == '6':

%o yield m

%o A351472_list = list(islice(A351472_gen(), 20)) # _Chai Wah Wu_, Feb 17 2022

%Y Similar with largest digit k: A333402 (k=1), A341383 (k=2), A350814 (k=3), A351470 (k=4), A351471 (k=5), this sequence (k=6), A351473 (k=7), A351474 (k=8), A333237 (k=9).

%Y Cf. A094028, A333236.

%K nonn,base

%O 1,1

%A _Bernard Schott_ and _Robert G. Wilson v_, Feb 16 2022