The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A053696 Numbers that can be represented as a string of three or more 1's in a base >= 2. 18

%I

%S 7,13,15,21,31,40,43,57,63,73,85,91,111,121,127,133,156,157,183,211,

%T 241,255,259,273,307,341,343,364,381,400,421,463,507,511,553,585,601,

%U 651,703,757,781,813,820,871,931,993,1023,1057,1093,1111,1123,1191

%N Numbers that can be represented as a string of three or more 1's in a base >= 2.

%C Numbers of the form (b^n-1)/(b-1) for n > 2 and b > 1. - _T. D. Noe_, Jun 07 2006

%C Numbers m that are nontrivial repunits for any base b >= 2. For k = 2 (I use k for the exponent since n is used as the index in a(n)) we get (b^k-1)/(b-1) = (b^2-1)/(b-1) = b+1, so every integer m >= 3 is a 2-digit repunit in base b = m-1. And for n = 1 (the 1-digit degenerate repunit) we get (b-1)/(b-1) = 1 for any base b >= 2. If we considered all k >= 1 we would get the sequence of all positive integers except 2 since it is the smallest uniform base used in positional representation (2 might be seen as the "repunit" in a nonpositional base representation such as the Roman numerals where 2 is expressed as II). - _Daniel Forgues_, Mar 01 2009

%C These repunits numbers belong to Brazilian numbers (A125134) (see Links: "Les nombres brĂ©siliens" - section IV, p. 32). - _Bernard Schott_, Dec 18 2012

%C The Brazilian primes (A085104) belong to this sequence. - _Bernard Schott_, Dec 18 2012

%H T. D. Noe and Charles R Greathouse IV, <a href="/A053696/b053696.txt">Table of n, a(n) for n = 1..10000</a> (first 1172 terms from Noe)

%H Bernard Schott, <a href="/A125134/a125134.pdf">Les nombres brĂ©siliens</a>, Quadrature, no. 76, avril-juin 2010, pages 30-38; included here with permission from the editors of Quadrature.

%F a(n) ~ n^2 since as n grows the density of repunits of degree 2 among all the repunits tends to 1. - _Daniel Forgues_, Dec 09 2008

%F A088323(a(n)) > 1. - _Reinhard Zumkeller_, Jan 22 2014

%e a(5) = 31 because 31 can be written as 111 base 5 (or indeed 11111 base 2).

%p N:= 10^4: # to get all terms <= N

%p V:= Vector(N):

%p for b from 2 while (b^3-1)/(b-1) <= N do

%p inds:= [seq((b^k-1)/(b-1), k=3..ilog[b](N*(b-1)+1))];

%p V[inds]:= 1;

%p od:

%p select(t -> V[t] = 1, [\$1..N]); # _Robert Israel_, Dec 10 2015

%t fQ[n_] := Block[{d = Rest@ Divisors[n - 1]}, Length@ d > 2 && Length@ Select[ IntegerDigits[n, d], Union@# == {1} &] > 1]; Select[ Range@ 1200, fQ]

%t lim=1000; Union[Reap[Do[n=3; While[a=(b^n-1)/(b-1); a<=lim, Sow[a]; n++], {b, 2, Floor[Sqrt[lim]]}]][[2, 1]]]

%t Take[Union[Flatten[With[{l=Table[PadLeft[{},n,1],{n,3,100}]}, Table[ FromDigits[#,n]&/@l,{n,2,100}]]]],80] (* _Harvey P. Dale_, Oct 06 2011 *)

%o (PARI) list(lim)=my(v=List(),e,t);for(b=2,sqrt(lim),e=3;while((t=(b^e-1)/(b-1))<=lim,listput(v,t);e++));vecsort(Vec(v),,8) \\ _Charles R Greathouse IV_, Oct 06 2011

%o (PARI) list(lim)=my(v=List(),e,t);for(b=2,lim^(1/3),e=4;while((t=(b^e-1)/(b-1))<=lim,listput(v,t);e++));vecsort(concat(Vec(v), vector((sqrtint (lim\1*4-3)-3)\2,i,i^2+3*i+3)),,8) \\ _Charles R Greathouse IV_, May 30 2013

%o a053696 n = a053696_list !! (n-1)

%o a053696_list = filter ((> 1) . a088323) [2..]

%o -- _Reinhard Zumkeller_, Jan 22 2014, Nov 26 2013

%Y Cf. A090503 (a subsequence), A119598 (numbers that are repunits in four or more bases), A125134, A085104.

%Y Cf. A108348.

%K nonn,base,easy

%O 1,1

%A _Henry Bottomley_, Mar 23 2000

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

Last modified April 5 06:13 EDT 2020. Contains 333238 sequences. (Running on oeis4.)