A136808 Numbers k such that k and k^2 use only the digits 0, 1 and 2.
0, 1, 10, 11, 100, 101, 110, 1000, 1001, 1010, 1011, 1100, 1101, 10000, 10001, 10010, 10011, 10100, 10110, 11000, 11001, 11010, 100000, 100001, 100010, 100011, 100100, 100101, 100110, 101000, 101001, 101100, 110000, 110001, 110010, 110100, 1000000, 1000001, 1000010, 1000011, 1000100
Offset: 1
Examples
101000100100001^2 = 10201020220210222010200200001.
Links
- Giovanni Resta, Table of n, a(n) for n = 1..10000 (first 1359 terms from Jonathan Wellons)
- Jonathan Wellons, Tables of Shared Digits [archived]
- Index to sequences related to squares.
Crossrefs
Programs
-
Maple
isA136808 := proc(n) local ndgs,n2dgs ; ndgs := convert(convert(n,base,10),set) ; n2dgs := convert(convert(n^2,base,10),set) ; if ( (ndgs union n2dgs) minus {0,1,2} ) = {} then true ; else false ; fi ; end: LtonRev := proc(L) local i ; add(op(i,L)*10^(i-1),i=1..nops(L)) ; end: A007089 := proc(n) convert(n,base,3) ; LtonRev(%) ; end: n := 1: for i from 0 do n3 := A007089(i) ; if isA136808(n3) then printf("%d %d ",n,n3) ; n := n+1 ; fi ; od: # R. J. Mathar, Jan 24 2008
-
Mathematica
Select[FromDigits/@Tuples[{0,1},7],Union[Take[DigitCount[#^2],{3,9}]]=={0}&] (* Harvey P. Dale, May 29 2013 *)
-
PARI
for(n=1,999,vecmax(digits((N=fromdigits(binary(n),10))^2))<3 && print1(N",")) \\ M. F. Hasler, Nov 15 2017
Comments