cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 51-60 of 238 results. Next

A020660 Lexicographically earliest increasing sequence of positive numbers that contains no arithmetic progression of length 8.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 14, 16, 17, 18, 19, 20, 21, 23, 24, 25, 26, 27, 28, 30, 31, 32, 33, 34, 35, 37, 38, 39, 40, 41, 42, 44, 45, 46, 47, 48, 49, 50, 59, 60, 61, 62, 63, 64, 65, 67, 69, 70, 71, 72, 74, 75, 76, 77, 78, 79, 81, 84, 85, 87, 88, 89, 91, 92, 93, 95, 96, 97
Offset: 1

Views

Author

David W. Wilson

Keywords

Links

Crossrefs

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).

Programs

  • Maple
    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, 8); # Robert Israel, Jan 04 2016

A020661 Lexicographically earliest increasing sequence of nonnegative numbers that contains no arithmetic progression of length 9.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 14, 15, 16, 18, 19, 20, 21, 22, 23, 26, 27, 28, 29, 30, 31, 32, 33, 36, 37, 38, 39, 40, 42, 43, 44, 45, 46, 47, 48, 49, 52, 54, 55, 56, 57, 58, 59, 63, 64, 65, 66, 67, 68, 69, 70, 77, 78, 79, 80, 81, 82, 83, 84, 86, 87, 90, 91, 93, 94, 95, 96, 97
Offset: 1

Views

Author

David W. Wilson

Keywords

Links

Crossrefs

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).

A020662 Lexicographically earliest increasing sequence of positive numbers that contains no arithmetic progression of length 9.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 16, 17, 19, 20, 21, 22, 23, 24, 27, 28, 29, 30, 31, 32, 33, 34, 37, 38, 39, 40, 41, 43, 44, 45, 46, 47, 48, 49, 50, 53, 55, 56, 57, 58, 59, 60, 64, 65, 66, 67, 68, 69, 70, 71, 78, 79, 80, 81, 82, 83, 84, 85, 87, 88, 91, 92, 94, 95, 96, 97
Offset: 1

Views

Author

David W. Wilson

Keywords

Links

Crossrefs

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).

Programs

  • Maple
    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,9); # Robert Israel, Jan 04 2016

A020663 Lexicographically earliest increasing sequence of nonnegative numbers that contains no arithmetic progression of length 10.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 16, 17, 19, 20, 21, 22, 23, 24, 25, 26, 27, 30, 31, 32, 33, 34, 35, 36, 37, 38, 40, 41, 42, 43, 44, 48, 49, 50, 51, 52, 53, 54, 57, 58, 59, 60, 61, 62, 63, 64, 65, 67, 68, 69, 70, 71, 72, 73, 74, 76, 77, 80, 81, 82, 83, 84, 87, 88, 95, 96
Offset: 1

Views

Author

David W. Wilson

Keywords

Links

Crossrefs

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).

A185256 Stanley Sequence S(0,3).

Original entry on oeis.org

0, 3, 4, 7, 9, 12, 13, 16, 27, 30, 31, 34, 36, 39, 40, 43, 81, 84, 85, 88, 90, 93, 94, 97, 108, 111, 112, 115, 117, 120, 121, 124, 243, 246, 247, 250, 252, 255, 256, 259, 270, 273, 274, 277, 279, 282, 283, 286, 324, 327, 328, 331, 333, 336, 337, 340, 351, 354, 355, 358, 360, 363
Offset: 1

Views

Author

N. J. A. Sloane, Mar 19 2011

Keywords

Comments

Given a finite increasing sequence V = [v_1, ..., v_k] containing no 3-term arithmetic progression, the Stanley Sequence S(V) is obtained by repeatedly appending the smallest term that is greater than the previous term and such that the new sequence also contains no 3-term arithmetic progression.

Examples

			After [0, 3, 4, 7, 9] the next term cannot be 10 or we would have the 3-term A.P. 4,7,10; it cannot be 11 because of 7,9,11; but 12 is OK.

References

  • R. K. Guy, Unsolved Problems in Number Theory, E10.

Links

Crossrefs

For other examples of Stanley Sequences see A005487, A005836, A187843, A188052, A188053, A188054, A188055, A188056, A188057.
See also A004793, A033160, A033163.

Programs

  • Maple
    # Stanley Sequences, Discrete Math. vol. 311 (2011), see p. 560
ss:=proc(s1,M) local n,chvec,swi,p,s2,i,j,t1,mmm; t1:=nops(s1); mmm:=1000;
s2:=Array(1..t1+M,s1); chvec:=Array(0..mmm);
for i from 1 to t1 do chvec[s2[i]]:=1; od;
# Get n-th term:
for n from t1+1 to t1+M do # do 1
# Try i as next term:
for i from s2[n-1]+1 to mmm do # do 2
swi:=-1;
# Test against j-th term:
for j from 1 to n-2 do # do 3
p:=s2[n-j];
if 2*p-i < 0 then break; fi;
if chvec[2*p-i] = 1 then swi:=1; break; fi;
od; # od 3
if swi=-1 then s2[n]:=i; chvec[i]:=1; break; fi;
od; # od 2
if swi=1 then ERROR("Error, no solution at n = ",n); fi;
od; # od 1;
[seq(s2[i],i=1..t1+M)];
end;
ss([0,3],80);
  • Mathematica
    ss[s1_, M_] := Module[{n, chvec, swi, p, s2, i, j, t1, mmm}, t1 = Length[s1]; mmm = 1000; s2 = Table[s1, {t1 + M}] // Flatten; chvec = Array[0&, mmm];
For[i = 1, i <= t1, i++, chvec[[s2[[i]] ]] = 1];
(* get n-th term *)
For[n = t1+1, n <= t1 + M, n++,
(* try i as next term *)
For[i = s2[[n-1]] + 1, i <= mmm, i++, swi = -1;
(* test against j-th term *)
For[j = 1, j <= n-2, j++, p = s2[[n - j]]; If[2*p - i < 0, Break[] ];
If[chvec[[2*p - i]] == 1, swi = 1; Break[] ] ];
If[swi == -1, s2[[n]] = i; chvec[[i]] = 1; Break[] ] ];
If[swi == 1, Print["Error, no solution at n = ", n] ] ];
Table[s2[[i]], {i, 1, t1 + M}] ];
ss[{0, 3}, 80] (* Jean-François Alcover, Sep 10 2013, translated from Maple *)
  • PARI
    A185256(n,show=1,L=3,v=[0,3], D=v->v[2..-1]-v[1..-2])={while(#v1||next(2), 2); break)); if(type(show)=="t_VEC", v, v[n])} \\ 2nd (optional) arg: zero = silent, nonzero = verbose, vector (e.g. [] or [1]) = get the whole list [a(1..n)] as return value, else just a(n). - M. F. Hasler, Jan 18 2016

A117940 a(0)=1, thereafter a(3n) = a(3n+1)/3 = a(n), a(3n+2)=0.

Original entry on oeis.org

1, 3, 0, 3, 9, 0, 0, 0, 0, 3, 9, 0, 9, 27, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 9, 0, 9, 27, 0, 0, 0, 0, 9, 27, 0, 27, 81, 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, 9, 0, 9, 27, 0, 0, 0, 0, 9, 27, 0, 27, 81, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Paul Barry, Apr 05 2006

Keywords

Comments

a(n) = a(3n)/a(0) = a(3n+1)/a(1). a(n) mod 2 = A039966(n). Row sums of A117939.
Observation: if this is written as a triangle (see example) then at least the first five row sums coincide with A002001. - Omar E. Pol, Nov 28 2011

Examples

			Contribution from Omar E. Pol, Nov 26 2011 (Start):
When written as a triangle this begins:
1,
3,0,
3,9,0,0,0,0,
3,9,0,9,27,0,0,0,0,0,0,0,0,0,0,0,0,0,
3,9,0,9,27,0,0,0,0,9,27,0,27,81,0,0,0,0,0,0,0,0,0,0,0,...
(End)

Links

Crossrefs

For generating functions Prod_{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.

Formula

G.f.: Product{k>=0, 1+3x^(3^k)}; a(n)=sum{k=0..n, sum{j=0..n, L(C(n,j)/3)*L(C(n-j,k)/3)}} where L(j/p) is the Legendre symbol of j and p.

A005839 Lexicographically earliest increasing nonnegative sequence that contains no 4-term arithmetic progression.

Original entry on oeis.org

0, 1, 2, 4, 5, 7, 8, 9, 14, 15, 16, 18, 25, 26, 28, 29, 30, 33, 36, 48, 49, 50, 52, 53, 55, 56, 57, 62, 64, 65, 66, 79, 86, 87, 88, 90, 93, 98, 101, 104, 105, 108, 109, 110, 121, 125, 135, 144, 148, 150, 151, 159, 162, 166, 168, 169, 170, 173, 175, 176, 182
Offset: 1

Views

Author

N. J. A. Sloane

Keywords

Comments

a(n) = A005837(n) - 1. - Alois P. Heinz, Jan 31 2014

References

  • R. K. Guy, Unsolved Problems in Number Theory, E10.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Links

Crossrefs

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).

Programs

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

Extensions

More terms from Jeffrey Shallit, Aug 15 1995.
Edited (with new offset, etc.) by N. J. A. Sloane, Jan 04 2016

A033052 a(1) = 1, a(2n) = 16a(n), a(2n+1) = a(2n)+1.

Original entry on oeis.org

0, 1, 16, 17, 256, 257, 272, 273, 4096, 4097, 4112, 4113, 4352, 4353, 4368, 4369, 65536, 65537, 65552, 65553, 65792, 65793, 65808, 65809, 69632, 69633, 69648, 69649, 69888, 69889, 69904, 69905, 1048576, 1048577, 1048592, 1048593, 1048832
Offset: 0

Views

Author

Clark Kimberling

Keywords

Comments

Numbers whose set of base 16 digits is {0,1}.
a(n) = Xpower(n,4). - Antti Karttunen, Apr 26 1999
Sums of distinct powers of 16.
For every nonnegative n, A000695(n) is a unique sum of the form a(k) + 4a(l). Thus every nonnegative n is a unique sum of the form a(p) + 2a(q) + 4a(r) + 8a(s). This gives a one-to-one map of the set N_0 of all nonnegative integers to (N_0)^4. Furthermore, if, for a fixed positive integer m, to consider all sums of distinct powers of 4^m, then one can obtain a one-to-one map of the set N_0 to (N_0)^(2^m). - Vladimir Shevelev, Nov 14 2008

Links

Crossrefs

Cf. A000695, A005836, A033042-A033051.
Column 4 of A048723. Row 15 of array A104257.

Programs

Formula

a(n) = Sum_{i=0..m} d(i)*16^i, where Sum_{i=0..m} d(i)*2^i is the base-2 representation of n.
a(n) = A097262(n)/15.
a(2n) = 16*a(n), a(2n+1) = a(2n)+1.
a(n) = Sum_{k>=0} A030308(n,k)*16^k. - Philippe Deléham, Oct 19 2011
G.f.: (1/(1 - x))*Sum_{k>=0} 16^k*x^(2^k)/(1 + x^(2^k)). - Ilya Gutkovskiy, Jun 04 2017
a(n) = A000695(A000695(n)). - Alan Michael Gómez Calderón, Mar 23 2025

Extensions

Extended by Ray Chandler, Aug 03 2004
Simpler definition from Ralf Stephan, Jun 18 2005

A062050 n-th chunk consists of the numbers 1, ..., 2^n.

Original entry on oeis.org

1, 1, 2, 1, 2, 3, 4, 1, 2, 3, 4, 5, 6, 7, 8, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17
Offset: 1

Views

Author

Marc LeBrun, Jun 30 2001

Keywords

Comments

a(k) is the distance between k and the largest power of 2 not exceeding k, where k = n + 1. [Consider the sequence of even numbers <= k; after sending the first term to the last position delete all odd-indexed terms; the final term that remains after iterating the process is the a(k)-th even number.] - Lekraj Beedassy, May 26 2005
Triangle read by rows in which row n lists the first 2^(n-1) positive integers, n >= 1; see the example. - Omar E. Pol, Sep 10 2013

Examples

			From _Omar E. Pol_, Aug 31 2013: (Start)
Written as irregular triangle with row lengths A000079:
  1;
  1, 2;
  1, 2, 3, 4;
  1, 2, 3, 4, 5, 6, 7, 8;
  1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16;
  ...
Row sums give A007582.
(End)

Links

Crossrefs

Cf. A053644, A053645.
Cf. A092754.

Programs

  • Haskell
    a062050 n = if n < 2 then n else 2 * a062050 n' + m - 1
            where (n',m) = divMod n 2
-- Reinhard Zumkeller, May 07 2012
  • Maple
    A062050 := proc(n) option remember; if n < 4 then return [1, 1, 2][n] fi;
2*A062050(floor(n/2)) + irem(n,2) - 1 end:
seq(A062050(n), n=1..89); # Peter Luschny, Apr 27 2020
  • Mathematica
    Flatten[Table[Range[2^n],{n,0,6}]] (* Harvey P. Dale, Oct 12 2015 *)
  • PARI
    a(n)=floor(n+1-2^logint(n,2))
  • PARI
    a(n)= n - 1<Ruud H.G. van Tol, Dec 13 2024
  • Python
    def A062050(n): return n-(1<Chai Wah Wu, Jan 22 2023

Formula

a(n) = A053645(n) + 1.
a(n) = n - msb(n) + 1 (where msb(n) = A053644(n)).
a(n) = 1 + n - 2^floor(log(n)/log(2)). - Benoit Cloitre, Feb 06 2003; corrected by Joseph Biberstine (jrbibers(AT)indiana.edu), Nov 25 2008
G.f.: 1/(1-x) * ((1-x+x^2)/(1-x) - Sum_{k>=1} 2^(k-1)*x^(2^k)). - Ralf Stephan, Apr 18 2003
a(1) = 1, a(2*n) = 2*a(n) - 1, a(2*n+1) = 2*a(n). - Ralf Stephan, Oct 06 2003
A005836(a(n+1)) = A107681(n). - Reinhard Zumkeller, May 20 2005
a(n) = if n < 2 then n else 2*a(floor(n/2)) - 1 + n mod 2. - Reinhard Zumkeller, May 07 2012
Without the constant 1, Ralf Stephan's g.f. becomes A(x) = x/(1-x)^2 - (1/(1-x)) * Sum_{k>=1} 2^(k-1)*x^(2^k) and satisfies the functional equation A(x) - 2*(1+x)*A(x^2) = x*(1 - x - x^2)/(1 - x^2). - Petros Hadjicostas, Apr 27 2020
For n > 0: a(n) = (A006257(n) + 1) / 2. - Frank Hollstein, Oct 25 2021

A151666 Number of partitions of n into distinct powers of 4.

Original entry on oeis.org

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

Views

Author

N. J. A. Sloane, May 30 2009

Keywords

Links

Crossrefs

For generating functions Prod_{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. A039966, A151667, A000695, A269707.

Programs

  • Haskell
    a151666 n = fromEnum (n < 2 || m < 2 && a151666 n' == 1)
   where (n', m) = divMod n 4
-- Reinhard Zumkeller, Dec 03 2011
  • Mathematica
    terms = 105;
kmax = Log[4, terms] // Ceiling;
CoefficientList[Product[1+x^(4^k), {k, 0, kmax}] + O[x]^(kmax terms), x][[1 ;; terms]] (* Jean-François Alcover, Jul 31 2018 *)

Formula

G.f.: Prod_{k >= 0 } (1+x^(4^k)). Exponents give A000695.
G.f. A(x) satisfies: A(x) = (1 + x) * A(x^4). - Ilya Gutkovskiy, Aug 12 2019
Previous Showing 51-60 of 238 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.