A269707 Decimal expansion of x = 3*Sum_{n in E} 1/10^n where E is the set of numbers whose base-4 representation consists of only 0's and 1's.
3, 3, 0, 0, 3, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 3, 0, 0, 3, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 3, 0, 0, 3, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 3, 0, 0, 3, 3, 0
Offset: 1
Examples
x = 3.3003300000000003300330000000000000000000000000000... 1/x = 0.303000003030000000000000000000003030000030300000...
References
- Daniel Duverney, Number Theory, World Scientific, 2010, 2.10 A striking number, pp. 19-20.
Programs
-
Maple
Digits:=200:nn:=5000:s:=0: for n from 0 to nn do: x:=convert(n,base,4):n0:=nops(x): it:=0:ii:=0: for k from 1 to n0 while(ii=0) do: if x[k]=0 or x[k]=1 then it:=it+1: else fi: od: if it=n0 then s:= s+evalf(1/10^n): else ii:=1:fi: od: print(3*s): print(1/(3*s)):
-
Mathematica
a[n_] := 3 * Boole[Max @ IntegerDigits[n-1, 4] <= 1]; Array[a, 100] (* Amiram Eldar, Aug 06 2021 *)
Extensions
Edited by Rick L. Shepherd, May 31 2016
Comments