A005836 -id:A005836 - cp's OEIS frontend

A032924 Numbers whose ternary expansion contains no 0.

1, 2, 4, 5, 7, 8, 13, 14, 16, 17, 22, 23, 25, 26, 40, 41, 43, 44, 49, 50, 52, 53, 67, 68, 70, 71, 76, 77, 79, 80, 121, 122, 124, 125, 130, 131, 133, 134, 148, 149, 151, 152, 157, 158, 160, 161, 202, 203, 205, 206, 211, 212, 214, 215, 229, 230, 232, 233, 238, 239

Offset: 1

Clark Kimberling

Complement of A081605. - Reinhard Zumkeller, Mar 23 2003
Subsequence of A154314. - Reinhard Zumkeller, Jan 07 2009
The first 28 terms are the range of A059852 (Morse codes for letters, when written in base 3) union {44, 50} (which correspond to Morse codes of Ü and Ä). Subsequent terms represent the Morse code of other symbols in the same coding. - M. F. Hasler, Jun 22 2020

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Robert Baillie and Thomas Schmelzer, Summing Kempner's Curious (Slowly-Convergent) Series, Mathematica Notebook kempnerSums.nb, Wolfram Library Archive, 2008.
David Garth and Adam Gouge, Affinely Self-Generating Sets and Morphisms, Journal of Integer Sequences, 10 (2007), Article 07.1.5., 1-13.
Clark Kimberling, Affinely recursive sets and orderings of languages, Discrete Math., 274 (2004), 147-160.
Index entries for 3-automatic sequences.

Cf. A005823, A005836, A065361, A007089, A081608, A132140, A132141.

Zeroless numbers in some other bases <= 10: A000042 (base 2), A023705 (base 4), A248910 (base 6), A255805 (base 8), A255808 (base 9), A052382 (base 10).

a032924 n = a032924_list !! (n-1)
a032924_list = iterate f 1 where
   f x = 1 + if r < 2 then x else 3 * f x'  where (x', r) = divMod x 3
-- Reinhard Zumkeller, Mar 07 2015, May 04 2012

f:= proc(n) local L,i,m;
   L:= convert(n,base,2);
   m:= nops(L);
   add((1+L[i])*3^(i-1),i=1..m-1);
end proc:
map(f, [$2..101]); # Robert Israel, Aug 04 2015

Select[Range@ 240, Last@ DigitCount[#, 3] == 0 &] (* Michael De Vlieger, Aug 05 2015 *)
Flatten[Table[FromDigits[#,3]&/@Tuples[{1,2},n],{n,5}]] (* Harvey P. Dale, May 28 2016 *)

apply( {A032924(n)=if(n<3,n,3*self()((n-1)\2)+2-n%2)}, [1..99]) \\ M. F. Hasler, Jun 22 2020

a(n) = fromdigits(apply(d->d+1,binary(n+1)[^1]), 3); \\ Kevin Ryde, Jun 23 2020

def a(n): return sum(3**i*(int(b)+1) for i, b in enumerate(bin(n+1)[:2:-1]))
print([a(n) for n in range(1, 61)]) # Michael S. Branicky, Aug 15 2022

def is_A032924(n):
    while n > 2:
       n,r = divmod(n,3)
       if r==0: return False
    return n > 0
print([n for n in range(250) if is_A032924(n)]) # M. F. Hasler, Feb 15 2023

def A032924(n): return int(bin(m:=n+1)[3:],3) + (3**(m.bit_length()-1)-1>>1) # Chai Wah Wu, Oct 13 2023

a(n) = A107680(n) + A107681(n). - Reinhard Zumkeller, May 20 2005
A081604(A107681(n)) <= A081604(A107680(n)) = A081604(a(n)) = A000523(n+1). - Reinhard Zumkeller, May 20 2005
A077267(a(n)) = 0. - Reinhard Zumkeller, Mar 02 2008
a(1)=1, a(n+1) = f(a(n)+1,a(n)+1) where f(x,y) = if x<3 and x<>0 then y, else if x mod 3 = 0 then f(y+1,y+1), else f(floor(x/3),y). - Reinhard Zumkeller, Mar 02 2008
a(2*n) = a(2*n-1)+1, n>0. - Zak Seidov, Jul 27 2009
A212193(a(n)) = 0. - Reinhard Zumkeller, May 04 2012
a(2*n+1) = 3*a(n)+1. - Robert Israel, Aug 05 2015
G.f.: x/(1-x)^2 + Sum_{m >= 1} 3^(m-1)*x^(2^(m+1)-1)/((1-x^(2^m))*(1-x)). - Robert Israel, Aug 04 2015
A065361(a(n)) = n. - Rémy Sigrist, Feb 06 2023
Sum_{n>=1} 1/a(n) = 3.4977362637842652509313189236131190039368413460747606236619907531632476445332666030262441154353753276457... (calculated using Baillie and Schmelzer's kempnerSums.nb, see Links). - Amiram Eldar, Apr 14 2025

A074940 Numbers having at least one 2 in their ternary representation.

Original entry on oeis.org

2, 5, 6, 7, 8, 11, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 29, 32, 33, 34, 35, 38, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 83, 86, 87, 88, 89, 92

Offset: 1

Benoit Cloitre and Reinhard Zumkeller, Oct 04 2002; revised Dec 03 2003

Also, numbers m such that 3 divides C(2m,m).
Also, numbers m such that the central trinomial coefficient A002426(m) == 0 (mod 3). - Emeric Deutsch and Bruce E. Sagan, Dec 04 2003
Also, numbers m such that A092255(m) == 0 (mod 3). - Benoit Cloitre, Mar 22 2004
Also, numbers m such that the coefficient of x^m equals 0 in Product_{k>=0} (1-x^(3^k)). - N. J. A. Sloane, Jun 01 2010

			12 is not in the sequence since it is 110_3, but 11 is in the sequence since it is 102_3. - _Michael B. Porter_, Jun 30 2016

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Emeric Deutsch and Bruce E. Sagan, Congruences for Catalan and Motzkin numbers and related sequences, arXiv:math/0407326 [math.CO], 2004.
Emeric Deutsch and Bruce E. Sagan, Congruences for Catalan and Motzkin numbers and related sequences, J. Num. Theory 117 (2006), 191-215.

Complement of A005836.
Cf. A006996, A007089, A081603, A081610, A081605, A081606.
A039966(a(n)) = 0.

a074940 n = a074940_list !! (n-1)
a074940_list = filter ((== 0) . a039966) [0..]
-- Reinhard Zumkeller, Jun 06 2012, Sep 29 2011

Select[Range@ 120, MemberQ[IntegerDigits[#, 3], 2] &] (* or *)
Select[Range@ 120, Divisible[Binomial[2 #, #], 3] &] (* Michael De Vlieger, Jun 29 2016 *)
Select[Range[100],DigitCount[#,3,2]>0&] (* Harvey P. Dale, Aug 25 2019 *)

is(n)=while(n,if(n%3==2,return(1));n\=3);0 \\ Charles R Greathouse IV, Aug 21 2011

from gmpy2 import digits
def A074940(n):
    def f(x):
        s = digits(x,3)
        for i in range(l:=len(s)):
            if s[i]>'1':
                break
        else:
            return n+int(s,2)
        return n+int(s[:i]+'1'*(l-i),2)
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    return m # Chai Wah Wu, Oct 29 2024

a(n) = n + O(n^0.631). - Charles R Greathouse IV, Aug 21 2011

More terms from Emeric Deutsch and Bruce E. Sagan, Dec 04 2003

A033042 Sums of distinct powers of 5.

Original entry on oeis.org

0, 1, 5, 6, 25, 26, 30, 31, 125, 126, 130, 131, 150, 151, 155, 156, 625, 626, 630, 631, 650, 651, 655, 656, 750, 751, 755, 756, 775, 776, 780, 781, 3125, 3126, 3130, 3131, 3150, 3151, 3155, 3156, 3250, 3251, 3255, 3256, 3275, 3276, 3280, 3281, 3750, 3751

Offset: 0

Clark Kimberling

Numbers without any base-5 digits larger than 1.
a(n) modulo 2 is the Prouhet-Thue-Morse sequence A010060. - Philippe Deléham, Oct 17 2011
Values of k where A008977(k) does not end with 0. - Henry Bottomley, Nov 09 2022

T. D. Noe, Table of n, a(n) for n = 0..1023
David Applegate, Omar E. Pol and N. J. A. Sloane, The Toothpick Sequence and Other Sequences from Cellular Automata, Congressus Numerantium, Vol. 206 (2010), 157-191. [There is a typo in Theorem 6: (13) should read u(n) = 4.3^(wt(n-1)-1) for n >= 2.]
K. Dilcher and L. Ericksen, Hyperbinary expansions and Stern polynomials, Elec. J. Combin, 22, 2015, #P2.24.
Hsien-Kuei Hwang, Svante Janson, and Tsung-Hsi Tsai, Identities and periodic oscillations of divide-and-conquer recurrences splitting at half, arXiv:2210.10968 [cs.DS], 2022, p. 45.
N. J. A. Sloane, Catalog of Toothpick and Cellular Automata Sequences in the OEIS

For generating functions Product_{k>=0} (1 + a*x^(b^k)) for the following values of (a,b) see: (1,2) A000012 and A000027, (1,3) A039966 and A005836, (1,4) A151666 and A000695, (1,5) A151667 and A033042, (2,2) A001316, (2,3) A151668, (2,4) A151669, (2,5) A151670, (3,2) A048883, (3,3) A117940, (3,4) A151665, (3,5) A151671, (4,2) A102376, (4,3) A151672, (4,4) A151673, (4,5) A151674.
Cf. A000695, A005836, A008977, A010060, A033043-A033052.
Row 5 of array A104257.

function a(n)
    m, r, b = n, 0, 1
    while m > 0
        m, q = divrem(m, 2)
        r += b * q
        b *= 5
    end
r end; [a(n) for n in 0:49] |> println # Peter Luschny, Jan 03 2021

a:= proc(n) local m, r, b; m, r, b:= n, 0, 1;
      while m>0 do r:= r+b*irem(m, 2, 'm'); b:= b*5 od; r
    end:
seq(a(n), n=0..100);  # Alois P. Heinz, Mar 16 2013

t = Table[FromDigits[RealDigits[n, 2], 5], {n, 1, 100}]
(* Clark Kimberling, Aug 02 2012 *)
FromDigits[#,5]&/@Tuples[{0,1},7] (* Harvey P. Dale, May 22 2018 *)

a(n) = subst(Pol(binary(n)), 'x, 5);
vector(50, i, a(i-1))  \\ Gheorghe Coserea, Sep 15 2015

a(n)=fromdigits(binary(n),5) \\ Charles R Greathouse IV, Jan 11 2017

def A033042(n): return int(bin(n)[2:],5) # Chai Wah Wu, Oct 30 2024

a(n) = Sum_{i=0..m} d(i)*5^i, where Sum_{i=0..m} d(i)*2^i is the base-2 representation of n.
Numbers j such that the coefficient of x^j is > 0 in Product_{k>=0} (1 + x^(5^k)). - Benoit Cloitre, Jul 29 2003
a(n) = A097251(n)/4.
a(2n) = 5*a(n), a(2n+1) = a(2n)+1.
a(n) = Sum_{k>=0} A030308(n,k)*5^k. - Philippe Deléham, Oct 17 2011
liminf a(n)/n^(log(5)/log(2)) = 1/4 and limsup a(n)/n^(log(5)/log(2)) = 1. - Gheorghe Coserea, Sep 15 2015
G.f.: (1/(1 - x))*Sum_{k>=0} 5^k*x^(2^k)/(1 + x^(2^k)). - Ilya Gutkovskiy, Jun 04 2017

Extended by Ray Chandler, Aug 03 2004

A039966 a(0) = 1; thereafter a(3n+2) = 0, a(3n) = a(3n+1) = a(n).

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 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, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0

Offset: 0

N. J. A. Sloane, Dec 11 1999

nonn

Number of partitions of n into distinct powers of 3.
Trajectory of 1 under the morphism: 1 -> 110, 0 -> 000. Thus 1 -> 110 ->110110000 -> 110110000110110000000000000 -> ... - Philippe Deléham, Jul 09 2005
Also, an example of a d-perfect sequence.
This is a composite of two earlier sequences contributed at different times by N. J. A. Sloane and by Reinhard Zumkeller, Mar 05 2005. Christian G. Bower extended them and found that they agreed for at least 512 terms. The proof that they were identical was found by Ralf Stephan, Jun 13 2005, based on the fact that they were both 3-regular sequences.

			The triples of elements (a(3k), a(3k+1), a(3k+2)) are (1,1,0) if a(k) = 1 and (0,0,0) if a(k) = 0.  So since a(2) = 0, a(6) = a(7) = a(8) = 0, and since a(3) = 1, a(9) = a(10) = 1 and a(11) = 0. - _Michael B. Porter_, Jul 11 2016

Reinhard Zumkeller, Table of n, a(n) for n = 0..10000
David Applegate, Omar E. Pol and N. J. A. Sloane, The Toothpick Sequence and Other Sequences from Cellular Automata, Congressus Numerantium, Vol. 206 (2010), 157-191. [There is a typo in Theorem 6: (13) should read u(n) = 4.3^(wt(n-1)-1) for n >= 2.]
D. Kohel, S. Ling and C. Xing, Explicit Sequence Expansions, in Sequences and their Applications, C. Ding, T. Helleseth, and H. Niederreiter, eds., Proceedings of SETA'98 (Singapore, 1998), 308-317, 1999.
N. J. A. Sloane, Catalog of Toothpick and Cellular Automata Sequences in the OEIS
Index entries for sequences that are fixed points of mappings
Index entries for characteristic functions

For generating functions Product_{k>=0} (1+a*x^(b^k)) for the following values of (a,b) see: (1,2) A000012 and A000027, (1,3) A039966 and A005836, (1,4) A151666 and A000695, (1,5) A151667 and A033042, (2,2) A001316, (2,3) A151668, (2,4) A151669, (2,5) A151670, (3,2) A048883, (3,3) A117940, (3,4) A151665, (3,5) A151671, (4,2) A102376, (4,3) A151672, (4,4) A151673, (4,5) A151674.
Cf. A062051, A000009, A000244, A004642.
Characteristic function of A005836 (and apart from offset of A003278).

a039966 n = fromEnum (n < 2 || m < 2 && a039966 n' == 1)
   where (n',m) = divMod n 3
-- Reinhard Zumkeller, Sep 29 2011

a := proc(n) option remember; if n <= 1 then RETURN(1) end if; if n = 2 then RETURN(0) end if; if n mod 3 = 2 then RETURN(0) end if; if n mod 3 = 0 then RETURN(a(1/3*n)) end if; if n mod 3 = 1 then RETURN(a(1/3*n - 1/3)) end if end proc; # Ralf Stephan, Jun 13 2005

(* first do *) Needs["DiscreteMath`Combinatorica`"] (* then *) s = Rest[ Sort[ Plus @@@ Table[UnrankSubset[n, Table[3^i, {i, 0, 4}]], {n, 32}]]]; Table[ If[ Position[s, n] == {}, 0, 1], {n, 105}] (* Robert G. Wilson v, Jun 14 2005 *)
CoefficientList[Series[Product[(1 + x^(3^k)), {k, 0, 5}], {x, 0, 111}], x] (* or *)
Nest[ Flatten[ # /. {0 -> {0, 0, 0}, 1 -> {1, 1, 0}}] &, {1}, 5] (* Robert G. Wilson v, Mar 29 2006 *)
Nest[ Join[#, #, 0 #] &, {1}, 5] (* Robert G. Wilson v, Jul 27 2014 *)

{a(n)=local(A,m); if(n<0, 0, m=1; A=1+O(x); while(m<=n, m*=3; A=(1+x)*subst(A,x,x^3)); polcoeff(A,n))} /* Michael Somos, Jul 15 2005 */

A039966(n)=vecmax(digits(n+!n,3))<2;
apply(A039966, [0..99]) \\ M. F. Hasler, Feb 15 2023

def A039966(n):
    while n > 2:
        n,r = divmod(n,3)
        if r==2: return 0
    return int(n!=2) # M. F. Hasler, Feb 15 2023

a(0) = 1, a(1) = 0, a(n) = b(n-2), where b is the sequence defined by b(0) = 1, b(3n+2) = 0, b(3n) = b(3n+1) = b(n). - Ralf Stephan
a(n) = A005043(n-1) mod 3. - Christian G. Bower, Jun 12 2005
a(n) = A002426(n) mod 3. - John M. Campbell, Aug 24 2011
a(n) = A000275(n) mod 3. - John M. Campbell, Jul 08 2016
Properties: 0 <= a(n) <= 1, a(A074940(n)) = 0, a(A005836(n)) = 1; A104406(n) = Sum(a(k), 1 <= k <= n). - Reinhard Zumkeller, Mar 05 2005
Euler transform of sequence b(n) where b(3^k) = 1, b(2*3^k) = -1 and zero otherwise. - Michael Somos, Jul 15 2005
G.f. A(x) satisfies A(x) = (1+x)*A(x^3). - Michael Somos, Jul 15 2005
G.f.: Product{k>=0} 1+x^(3^k). Exponents give A005836.

Entry revised Jun 30 2005

Offset corrected by John M. Campbell, Aug 24 2011

A020658 Lexicographically earliest increasing sequence of positive numbers that contains no arithmetic progression of length 7.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 15, 16, 17, 18, 19, 20, 22, 23, 24, 25, 26, 27, 29, 30, 31, 32, 33, 34, 36, 37, 38, 39, 40, 41, 50, 51, 52, 53, 54, 55, 57, 58, 59, 60, 61, 62, 64, 65, 66, 67, 68, 69, 71, 72, 73, 74, 75, 76, 78, 79, 80, 81, 82, 83, 85, 86, 87, 88, 89, 90, 99

Offset: 1

David W. Wilson

nonn

This is different from A047304: note the gap between 41 and 50. - M. F. Hasler, Oct 07 2014

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A047304.
Summary of increasing sequences avoiding arithmetic progressions of specified lengths (the second of each pair is obtained by adding 1 to the first):
3-term AP: A005836 (>=0), A003278 (>0);
4-term AP: A005839 (>=0), A005837 (>0);
5-term AP: A020654 (>=0), A020655 (>0);
6-term AP: A020656 (>=0), A005838 (>0);
7-term AP: A020657 (>=0), A020658 (>0);
8-term AP: A020659 (>=0), A020660 (>0);
9-term AP: A020661 (>=0), A020662 (>0);
10-term AP: A020663 (>=0), A020664 (>0).

Noap:= proc(N,m)
# N terms of earliest increasing seq with no m-term arithmetic progression
local A,forbid,n,c,ds,j;
A:= Vector(N):
A[1..m-1]:= <($1..m-1)>:
forbid:= {m}:
for n from m to N do
  c:= min({$A[n-1]+1..max(max(forbid)+1, A[n-1]+1)} minus forbid);
  A[n]:= c;
  ds:= convert(map(t -> c-t, A[m-2..n-1]),set);
  for j from m-2 to 2 by -1 do
    ds:= ds intersect convert(map(t -> (c-t)/j, A[m-j-1..n-j]),set);
    if ds = {} then break fi;
  od;
  forbid:= select(`>`,forbid,c) union map(`+`,ds,c);
od:
convert(A,list)
end proc:
Noap(100, 7); # Robert Israel, Jan 04 2016

Select[Range[0, 100], FreeQ[IntegerDigits[#, 7], 6]&] + 1 (* Jean-François Alcover, Aug 18 2023, after M. F. Hasler *)

a(n) = A020657(n)+1. - M. F. Hasler, Oct 07 2014

A020664 Lexicographically earliest increasing sequence of positive numbers that contains no arithmetic progression of length 10.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 12, 13, 14, 15, 16, 17, 18, 20, 21, 22, 23, 24, 25, 26, 27, 28, 31, 32, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 44, 45, 49, 50, 51, 52, 53, 54, 55, 58, 59, 60, 61, 62, 63, 64, 65, 66, 68, 69, 70, 71, 72, 73, 74, 75, 77, 78, 81, 82, 83, 84, 85, 88, 89, 96, 97

Offset: 1

David W. Wilson

nonn

Robert Israel, Table of n, a(n) for n = 1..10000

Summary of increasing sequences avoiding arithmetic progressions of specified lengths (the second of each pair is obtained by adding 1 to the first):
3-term AP: A005836 (>=0), A003278 (>0);
4-term AP: A005839 (>=0), A005837 (>0);
5-term AP: A020654 (>=0), A020655 (>0);
6-term AP: A020656 (>=0), A005838 (>0);
7-term AP: A020657 (>=0), A020658 (>0);
8-term AP: A020659 (>=0), A020660 (>0);
9-term AP: A020661 (>=0), A020662 (>0);
10-term AP: A020663 (>=0), A020664 (>0).

Noap:= proc(N,m)
# N terms of earliest increasing seq with no m-term arithmetic progression
local A,forbid,n,c,ds,j;
A:= Vector(N):
A[1..m-1]:= <($1..m-1)>:
forbid:= {m}:
for n from m to N do
  c:= min({$A[n-1]+1..max(max(forbid)+1, A[n-1]+1)} minus forbid);
  A[n]:= c;
  ds:= convert(map(t -> c-t, A[m-2..n-1]),set);
  for j from m-2 to 2 by -1 do
    ds:= ds intersect convert(map(t -> (c-t)/j, A[m-j-1..n-j]),set);
    if ds = {} then break fi;
  od;
  forbid:= select(`>`,forbid,c) union map(`+`,ds,c);
od:
convert(A,list)
end proc:
Noap(100,10); # Robert Israel, Jan 04 2016

A005838 Lexicographically earliest increasing sequence of positive numbers that contains no arithmetic progression of length 6.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 10, 12, 13, 14, 15, 17, 18, 19, 20, 22, 23, 24, 25, 26, 33, 34, 35, 36, 37, 39, 43, 44, 45, 46, 47, 49, 50, 51, 52, 59, 60, 62, 63, 64, 65, 66, 68, 69, 71, 73, 77, 85, 87, 88, 89, 90, 91, 93, 96, 97, 98, 99, 100, 103, 104, 107, 111, 114, 115, 117, 118, 120

Offset: 1

N. J. A. Sloane, Jeffrey Shallit

nonn

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Robert Israel, Table of n, a(n) for n = 1..10000
J. L. Gerver and L. T. Ramsey, Sets of integers with no long arithmetic progressions generated by the greedy algorithm, Math. Comp. 33 (1979), 1353-1359.

Summary of increasing sequences avoiding arithmetic progressions of specified lengths (the second of each pair is obtained by adding 1 to the first):
3-term AP: A005836 (>=0), A003278 (>0);
4-term AP: A005839 (>=0), A005837 (>0);
5-term AP: A020654 (>=0), A020655 (>0);
6-term AP: A020656 (>=0), A005838 (>0);
7-term AP: A020657 (>=0), A020658 (>0);
8-term AP: A020659 (>=0), A020660 (>0);
9-term AP: A020661 (>=0), A020662 (>0);
10-term AP: A020663 (>=0), A020664 (>0).

N:= 100: # to get a(1)..a(N)
A:= Vector(N):
A[1..5]:= <($1..5)>:
forbid:= {6}:
for n from 6 to N do
  c:= min({$A[n-1]+1::max(max(forbid)+1, A[n-1]+1)} minus forbid);
  A[n]:= c;
  ds:= convert(map(t -> c-t, A[4..n-1],set);
  if ds = {} then next fi;
  ds:= ds intersect convert(map(t -> (c-t)/4, A[1..n-4]),set);
  if ds = {} then next fi;
  ds:= ds intersect convert(map(t -> (c-t)/3, A[2..n-3]),set);
  if ds = {} then next fi;
  ds:= ds intersect convert(map(t -> (c-t)/2, A[3..n-2]),set);
  forbid:= select(`>`,forbid,c) union map(`+`,ds,c);
od:
convert(A,list); # Robert Israel, Jan 04 2016

A005838(n,show=1,i=1,o=6,u=0)={for(n=1,n,show&&print1(i,",");u+=1<M. F. Hasler, Jan 03 2016

a(n) = A020656(n) + 1.

Name and links/references edited by M. F. Hasler, Jan 03 2016

Further edited by N. J. A. Sloane, Jan 04 2016

A020655 Lexicographically earliest increasing sequence of positive numbers that contains no arithmetic progression of length 5.

Original entry on oeis.org

1, 2, 3, 4, 6, 7, 8, 9, 11, 12, 13, 14, 16, 17, 18, 19, 26, 27, 28, 29, 31, 32, 33, 34, 36, 37, 38, 39, 41, 42, 43, 44, 51, 52, 53, 54, 56, 57, 58, 59, 61, 62, 63, 64, 66, 67, 68, 69, 76, 77, 78, 79, 81, 82, 83, 84, 86, 87, 88, 89, 91, 92, 93, 94, 126, 127, 128, 129, 131, 132, 133

Offset: 1

David W. Wilson

nonn

Robert Israel, Table of n, a(n) for n = 1..10000

Summary of increasing sequences avoiding arithmetic progressions of specified lengths (the second of each pair is obtained by adding 1 to the first):
3-term AP: A005836 (>=0), A003278 (>0);
4-term AP: A005839 (>=0), A005837 (>0);
5-term AP: A020654 (>=0), A020655 (>0);
6-term AP: A020656 (>=0), A005838 (>0);
7-term AP: A020657 (>=0), A020658 (>0);
8-term AP: A020659 (>=0), A020660 (>0);
9-term AP: A020661 (>=0), A020662 (>0);
10-term AP: A020663 (>=0), A020664 (>0).

Noap:= proc(N,m)
# N terms of earliest increasing seq with no m-term arithmetic progression
local A,forbid,n,c,ds,j;
A:= Vector(N):
A[1..m-1]:= <($1..m-1)>:
forbid:= {m}:
for n from m to N do
  c:= min({$A[n-1]+1..max(max(forbid)+1, A[n-1]+1)} minus forbid);
  A[n]:= c;
  ds:= convert(map(t -> c-t, A[m-2..n-1]),set);
  for j from m-2 to 2 by -1 do
    ds:= ds intersect convert(map(t -> (c-t)/j, A[m-j-1..n-j]),set);
    if ds = {} then break fi;
  od;
  forbid:= select(`>`,forbid,c) union map(`+`,ds,c);
od:
convert(A,list)
end proc:
Noap(100, 5); # Robert Israel, Jan 04 2016

t = {1, 2, 3, 4}; Do[s = Table[Append[i, n], {i, Subsets[t, {4}]}]; If[! MemberQ[Table[Differences[i, 2], {i, s}], {0, 0, 0}], AppendTo[t, n]], {n, 5, 100}]; t (* T. D. Noe, Apr 17 2014 *)

A020656 Lexicographically earliest increasing sequence of nonnegative numbers that contains no arithmetic progression of length 6.

Original entry on oeis.org

0, 1, 2, 3, 4, 6, 7, 8, 9, 11, 12, 13, 14, 16, 17, 18, 19, 21, 22, 23, 24, 25, 32, 33, 34, 35, 36, 38, 42, 43, 44, 45, 46, 48, 49, 50, 51, 58, 59, 61, 62, 63, 64, 65, 67, 68, 70, 72, 76, 84, 86, 87, 88, 89, 90, 92, 95, 96, 97, 98, 99, 102, 103, 106, 110, 113, 114, 116, 117, 119, 121

Offset: 1

David W. Wilson

nonn

			5 is excluded since (0,1,2,3,4,5) would be a 6-term AP.
10 is excluded since (0,2,4,6,8,10) would be a 6-term AP.
Idem for 15 and 20.
25 is not excluded, but after 25 (1,6,11,16,21,26) would be a 6-AP, and similarly all of 26 through 32 are excluded.

David W. Wilson, Table of n, a(n) for n = 1..10000

Summary of increasing sequences avoiding arithmetic progressions of specified lengths (the second of each pair is obtained by adding 1 to the first):
3-term AP: A005836 (>=0), A003278 (>0);
4-term AP: A005839 (>=0), A005837 (>0);
5-term AP: A020654 (>=0), A020655 (>0);
6-term AP: A020656 (>=0), A005838 (>0);
7-term AP: A020657 (>=0), A020658 (>0);
8-term AP: A020659 (>=0), A020660 (>0);
9-term AP: A020661 (>=0), A020662 (>0);
10-term AP: A020663 (>=0), A020664 (>0).

A020656(n,show=1,i=0,o=6,u=0)={for(n=1,n,show&&print1(i,",");u+=1<M. F. Hasler, Jan 03 2016

a(n) = A005838(n) - 1.

Edited by M. F. Hasler, Jan 03 2016

Further edited (with new offset) by N. J. A. Sloane, Jan 04 2016

A020659 Lexicographically earliest increasing sequence of nonnegative numbers that contains no arithmetic progression of length 8.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 15, 16, 17, 18, 19, 20, 22, 23, 24, 25, 26, 27, 29, 30, 31, 32, 33, 34, 36, 37, 38, 39, 40, 41, 43, 44, 45, 46, 47, 48, 49, 58, 59, 60, 61, 62, 63, 64, 66, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 80, 83, 84, 86, 87, 88, 90, 91, 92, 94, 95, 96, 99

Offset: 1

David W. Wilson

nonn

David W. Wilson, Table of n, a(n) for n = 1..10000

Summary of increasing sequences avoiding arithmetic progressions of specified lengths (the second of each pair is obtained by adding 1 to the first):
3-term AP: A005836 (>=0), A003278 (>0);
4-term AP: A005839 (>=0), A005837 (>0);
5-term AP: A020654 (>=0), A020655 (>0);
6-term AP: A020656 (>=0), A005838 (>0);
7-term AP: A020657 (>=0), A020658 (>0);
8-term AP: A020659 (>=0), A020660 (>0);
9-term AP: A020661 (>=0), A020662 (>0);
10-term AP: A020663 (>=0), A020664 (>0).

cp's OEIS Frontend

A032924 Numbers whose ternary expansion contains no 0.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Python

Python

Python

Formula

A074940 Numbers having at least one 2 in their ternary representation.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Python

Formula

Extensions

A033042 Sums of distinct powers of 5.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Julia

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A039966 a(0) = 1; thereafter a(3n+2) = 0, a(3n) = a(3n+1) = a(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A020658 Lexicographically earliest increasing sequence of positive numbers that contains no arithmetic progression of length 7.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A020664 Lexicographically earliest increasing sequence of positive numbers that contains no arithmetic progression of length 10.

Views

Author