%I #26 Jun 14 2016 16:25:26
%S 1,5,7,13,23,53,125,215,373,1367,1373,1375,3551,4093,5471,5495,5503,
%T 30581,30589,32765,32767,56821,56831,89557,96119,96215,96223,97655,
%U 98135,98141,98143,98167,98293,98303,351743,352093,521599,521693,521717,521719,524119,524149,875893,875903,884725,884735
%N Numbers whose base-b expansions, for both b=3 and b=4, include no digits other than 1 and b-1.
%C 1, 7 and 32767 also share this property in base 2; their binary expansions consist only of a sequence of 1s.
%H Chai Wah Wu, <a href="/A262958/b262958.txt">Table of n, a(n) for n = 1..10000</a>
%e 53 is 1222 in base 3 and 311 in base 4; it only uses the digit 1 or the largest digit in the two bases and is therefore a term.
%e Similarly 215 is 21222 in base 3 and 3113 in base 4 so it is also a term.
%t Select[Range@ 1000000, Last@ DigitCount[#, 3] == 0 && Total@ Rest@ Drop[DigitCount[#, 4], {3}] == 0 &] (* _Michael De Vlieger_, Oct 05 2015 *)
%t Join[{1,5},Flatten[Table[Select[FromDigits[#,3]&/@Tuples[{1,2},n], Union[ IntegerDigits[ #,4]] =={1,3}&],{n,20}]]] (* _Harvey P. Dale_, Jun 14 2016 *)
%o (PARI) is(n)=!setsearch(Set(digits(n,3)),0) && #setintersect(Set(digits(n,4)),[0,2])==0 \\ _Charles R Greathouse IV_, Oct 12 2015
%o (Python)
%o from gmpy2 import digits
%o def f1(n):
%o s = digits(n,3)
%o m = len(s)
%o for i in range(m):
%o if s[i] == '0':
%o return(int(s[:i]+'1'*(m-i),3))
%o return n
%o def f2(n):
%o s = digits(n,4)
%o m = len(s)
%o for i in range(m):
%o if s[i] == '0':
%o return(int(s[:i]+'1'*(m-i),4))
%o if s[i] == '2':
%o return(int(s[:i]+'3'+'1'*(m-i-1),4))
%o return n
%o A262958_list = []
%o n = 1
%o for i in range(10**4):
%o m = f2(f1(n))
%o while m != n:
%o n, m = m, f2(f1(m))
%o A262958_list.append(m)
%o n += 1 # _Chai Wah Wu_, Oct 30 2015
%Y Cf. A207079, A258981, A261970.
%K nonn,base
%O 1,2
%A _Robin Powell_, Oct 05 2015
|