login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A146027 Numbers that can be written from base 2 to base 10 using only the digits 0 to 4. 6

%I #44 May 03 2019 08:40:21

%S 0,1,2,3,4,10,100,140004,140304,140312,1131032,1131033,1131034,1131040

%N Numbers that can be written from base 2 to base 10 using only the digits 0 to 4.

%C Conjectured to be complete.

%C Checked on Oct 31 2008 up to 2^16384 (or 1.2*10^4932). - _Daniel Mondot_, Jan 17 2016

%C No more entries < 10^78. - _Robert Israel_, Aug 31 2015

%C It is a plausible conjecture that there are no more terms, but this has not been proved. - _N. J. A. Sloane_, Nov 17 2017

%H Stuart A. Burrell, Han Yu, <a href="https://arxiv.org/abs/1905.00832">Digit expansions of numbers in different bases</a>, arXiv:1905.00832 [math.NT], 2019.

%p imax:= 20: # to consider numbers < 6^imax

%p L:= Matrix(5,imax):

%p Delta:= proc(L,b)

%p local i,j,m,Lloc;

%p if max(L) <= 4 then return 0 fi;

%p Lloc:= L;

%p m:= 0;

%p for j from 1 to imax while max(Lloc[j..imax]) > 4 do

%p m:= m + b^(j-1)*(b-Lloc[j]);

%p if j < imax then Lloc[j+1]:= Lloc[j+1]+1 fi

%p od;

%p m

%p end proc:

%p n:= 0: count:= 1: A[1]:= 0:

%p isdone:= false;

%p while max(L[..,imax]) < 5 and not isdone do

%p n:= n+1;

%p L[..,1]:= L[..,1]+<1,1,1,1,1>;

%p m:= max(seq(Delta(L[b-5,..],b),b=6..10));

%p while m > 0 and not isdone do

%p n:= n+m;

%p for b from 6 to 10 do

%p Lb:= convert(n,base,b);

%p if nops(Lb) > imax then isdone:= true; break fi;

%p L[b-5,1..nops(Lb)]:= Vector[row](Lb);

%p od:

%p m:= max(seq(Delta(L[b-5,..],b),b=6..10));

%p od;

%p if not isdone then

%p count:= count+1;

%p A[count]:= n;

%p fi

%p od:

%p seq(A[i],i=1..count); # _Robert Israel_, Aug 31 2015

%t f[n_] := Total[Total@ Drop[RotateRight[DigitCount[n, #]], 5] & /@ Range[6, 10]]; Select[Range[0, 1200000], f@ # == 0 &] (* Aug 29 2015, or *)

%t Select[Range[0, 1200000], Function[n, Times @@ Boole@ Map[Max@ IntegerDigits[n, #] <= 4 &, Range[2, 10]] > 0]] (* _Michael De Vlieger_, Aug 15 2016 *)

%o (PARI) isok(n) = if (n, for (b=6, 10, if (vecmax(digits(n,b))>4, return(0)))); 1; \\ _Michel Marcus_, Aug 30 2015

%Y Cf. A146025, A146026, A146028, A131646, A275600.

%K nonn,base

%O 1,3

%A _Daniel Mondot_, Oct 26 2008

%E Edited by _Charles R Greathouse IV_, Nov 01 2009

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 06:57 EDT 2024. Contains 371265 sequences. (Running on oeis4.)