OFFSET
1,3
COMMENTS
Generated with DrScheme.
A278038(18) = 10101, A136827(294) = 10110001101, A136831(1276) = 101100010001101 resp. A136836(1262) = 101090009991101 are the first terms from where on these four sequences differ from the present one. - M. F. Hasler, Nov 15 2017
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..10000 (first 1359 terms from Jonathan Wellons)
J. Wellons, Tables of Shared Digits [archived]
EXAMPLE
101000100100001^2 = 10201020220210222010200200001.
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 *)
PROG
(PARI) for(n=1, 999, vecmax(digits((N=fromdigits(binary(n), 10))^2))<3 && print1(N", ")) \\ M. F. Hasler, Nov 15 2017
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Jonathan Wellons (wellons(AT)gmail.com), Jan 22 2008
STATUS
approved