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.

Showing 1-10 of 15 results. Next

A127252 Sequence composed of 1 and -1 with the -1's occurring at odious indexed positions given by A091855.

Original entry on oeis.org

1, -1, 1, 1, -1, 1, 1, -1, 1, 1, 1, -1, 1, -1, 1, 1, -1, 1, 1, -1, 1, -1, 1, 1, 1, -1, 1, 1, -1, 1, 1, -1, 1, 1, 1, -1, 1, -1, 1, 1, 1, -1, 1, 1, -1, 1, 1, -1, 1, -1, 1, 1, -1, 1, 1, -1, 1, 1, 1, -1, 1, -1, 1, 1, -1, 1, 1, -1, 1, -1, 1, 1, 1, -1, 1, 1, -1, 1, 1
Offset: 0

Views

Author

Paul Barry, Jan 10 2007

Keywords

Crossrefs

Row sums of A127251.
Partial sums are A127255.

Programs

  • Mathematica
    a[n_] := If[EvenQ[IntegerExponent[n, 2]] && OddQ[DigitCount[n, 2, 1]], -1, 1]; Array[a, 100, 0] (* Amiram Eldar, Aug 04 2023 *)

Extensions

More terms from Amiram Eldar, Aug 04 2023

A127254 (0,1) sequence whose zero positions are indexed by twice the odious numbers given by A091855.

Original entry on oeis.org

1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1
Offset: 0

Views

Author

Paul Barry, Jan 10 2007

Keywords

Comments

0's appear at positions indexed by 2*A091855(n).

Crossrefs

Row sums of A127253.

Programs

  • Mathematica
    a[n_] := If[EvenQ[n] && EvenQ[IntegerExponent[n/2, 2]] && OddQ[DigitCount[n/2, 2, 1]], 0, 1]; Array[a, 100, 0] (* Amiram Eldar, Aug 06 2023 *)
  • PARI
    A127254(n) = if(n%2,1,!(hammingweight(n)%2 && !(valuation(n/2, 2)%2))); \\ After PARI-code in A091855. - Antti Karttunen, Nov 02 2018

Extensions

More terms from Antti Karttunen, Nov 02 2018

A127250 Sequence consisting of 1,3 or 5 with 3's occurring at the odious indices given by A091855 and 5's occurring at twice these odious indices.

Original entry on oeis.org

1, 3, 5, 1, 3, 1, 1, 3, 5, 1, 1, 3, 1, 3, 5, 1, 3, 1, 1, 3, 1, 3, 5, 1, 1, 3, 5, 1, 3, 1, 1, 3, 5, 1, 1, 3, 1, 3, 5, 1, 1, 3, 5, 1, 3, 1, 1, 3, 1, 3, 5, 1, 3, 1, 1, 3, 5, 1, 1, 3, 1, 3, 5, 1, 3, 1, 1, 3, 1, 3, 5, 1, 1, 3, 5, 1, 3, 1, 1, 3, 1, 3, 5, 1, 3, 1, 1
Offset: 0

Views

Author

Paul Barry, Jan 10 2007

Keywords

Crossrefs

Row sums of triangle A127249.

Programs

  • Mathematica
    a[n_] := If[OddQ[DigitCount[n, 2, 1]], If[EvenQ[IntegerExponent[n, 2]], 3, 5], 1]; Array[a, 100, 0] (* Amiram Eldar, Aug 04 2023 *)

Extensions

More terms from Amiram Eldar, Aug 04 2023

A356133 Complement of A026430.

Original entry on oeis.org

2, 4, 7, 11, 13, 17, 20, 22, 25, 29, 32, 34, 38, 40, 43, 47, 49, 53, 56, 58, 62, 64, 67, 71, 74, 76, 79, 83, 85, 89, 92, 94, 97, 101, 104, 106, 110, 112, 115, 119, 122, 124, 127, 131, 133, 137, 140, 142, 146, 148, 151, 155, 157, 161, 164, 166, 169, 173, 176
Offset: 1

Views

Author

Clark Kimberling, Aug 04 2022

Keywords

Examples

			The partial sums of the Thue-Morse sequence A001285 = (1,2,2,1,2,1,1,...) are A026430 = (0,1,3,5,6,8,9,10,...), from which the missing positive integers are (2,4,7,11,...).
		

Crossrefs

Programs

  • Mathematica
    u = Accumulate[1 + ThueMorse /@ Range[0, 2^7]]; (* A026430 *)
    Complement[Range[Max[u]], u]  (* A356133 *)
  • PARI
    a(n) = 3*n - 1 - hammingweight(n-1)%2; \\ Kevin Ryde, Aug 04 2022
    
  • Python
    def A356133(n): return 3*n-(2 if (n-1).bit_count()&1 else 1) # Chai Wah Wu, Mar 01 2023

Formula

a(n) = 3n - A001285(n-1) for n >= 1.

A050292 a(2n) = 2n - a(n), a(2n+1) = 2n + 1 - a(n) (for n >= 0).

Original entry on oeis.org

0, 1, 1, 2, 3, 4, 4, 5, 5, 6, 6, 7, 8, 9, 9, 10, 11, 12, 12, 13, 14, 15, 15, 16, 16, 17, 17, 18, 19, 20, 20, 21, 21, 22, 22, 23, 24, 25, 25, 26, 26, 27, 27, 28, 29, 30, 30, 31, 32, 33, 33, 34, 35, 36, 36, 37, 37, 38, 38, 39, 40, 41, 41, 42, 43, 44, 44, 45, 46, 47, 47, 48, 48, 49, 49, 50, 51, 52, 52, 53, 54
Offset: 0

Views

Author

Keywords

Comments

Note that the first equation implies a(0)=0, so there is no need to specify an initial value.
Maximal cardinality of a double-free subset of {1, 2, ..., n}, or in other words, maximal size of a subset S of {1, 2, ..., n} with the property that if x is in S then 2x is not. a(0)=0 by convention.
Least k such that a(k)=n is equal to A003159(n).
To construct the sequence: let [a, b, c, a, a, a, b, c, a, b, c, ...] be the fixed point of the morphism a -> abc, b ->a, c -> a, starting from a(1) = a, then write the indices of a, b, c, that of a being written twice; see A092606. - Philippe Deléham, Apr 13 2004
Number of integers from {1,...,n} for which the subtraction of 1 changes the parity of the number of 1's in their binary expansion. - Vladimir Shevelev, Apr 15 2010
Number of integers from {1,...,n} the factorization of which over different terms of A050376 does not contain 2. - Vladimir Shevelev, Apr 16 2010
a(n) modulo 2 is the Prouhet-Thue-Morse sequence A010060. Each number n appears A026465(n+1) times. - Philippe Deléham, Oct 19 2011
Another way of stating the last two comments from Philippe Deléham: the sequence can be obtained by replacing each term of the Thue-Morse sequence A010060 by the run number that term is in. - N. J. A. Sloane, Dec 31 2013

Examples

			Examples for n = 1 through 8: {1}, {1}, {1,3}, {1,3,4}, {1,3,4,5}, {1,3,4,5}, {1,3,4,5,7}, {1,3,4,5,7}.
Binary expansion of 5 is 101, so Sum{i>=0} b_i*(-1)^i = 2. Therefore a(5) = 10/3 + 2/3 = 4. - _Vladimir Shevelev_, Apr 15 2010
		

References

  • S. R. Finch, Mathematical Constants, Cambridge, 2003, Section 2.26.
  • Wang, E. T. H. "On Double-Free Sets of Integers." Ars Combin. 28, 97-100, 1989.

Crossrefs

Programs

  • Haskell
    a050292 n = a050292_list !! (n-1)
    a050292_list = scanl (+) 0 a035263_list
    -- Reinhard Zumkeller, Jan 21 2013
    
  • Maple
    A050292:=n->add((-1)^k*floor(n/2^k), k=0..n); seq(A050292(n), n=0..100); # Wesley Ivan Hurt, Feb 14 2014
  • Mathematica
    a[n_] := a[n] = If[n < 2, 1, n - a[Floor[n/2]]]; Table[ a[n], {n, 1, 75}]
    Join[{0},Accumulate[Nest[Flatten[#/.{0->{1,1},1->{1,0}}]&,{0},7]]] (* Harvey P. Dale, Apr 29 2018 *)
  • PARI
    a(n)=if(n<2,1,n-a(floor(n/2)))
    
  • Python
    from sympy.ntheory import digits
    def A050292(n): return ((n<<1)+sum((0,1,-1,0)[i] for i in digits(n,4)[1:]))//3 # Chai Wah Wu, Jan 30 2025

Formula

Partial sums of A035263. Close to (2/3)*n.
a(n) = A123087(2*n) = n - A123087(n). - Max Alekseyev, Mar 05 2023
From Benoit Cloitre, Nov 24 2002: (Start)
a(1)=1, a(n) = n - a(floor(n/2));
a(n) = (2/3)*n + (1/3)*A065359(n);
more generally, for m>=0, a(2^m*n) - 2^m*a(n) = A001045(m)*A065359(n) where A001045(m) = (2^m - (-1)^m)/3 is the Jacobsthal sequence;
a(A039004(n)) = (2/3)*A039004(n);
a(2*A039004(n)) = 2*a(A039004(n));
a(A003159(n)) = n;
a(A003159(n)-1) = n-1;
a(n) mod 2 = A010060(n) the Thue-Morse sequence;
a(n+1) - a(n) = A035263(n+1);
a(n+2) - a(n) = abs(A029884(n)).
(End)
G.f.: (1/(x-1)) * Sum_{i>=0} (-1)^i*x^(2^i)/(x^(2^i)-1). - Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Feb 17 2003
a(n) = Sum_{k>=0} (-1)^k*floor(n/2^k). - Benoit Cloitre, Jun 03 2003
a(A091785(n)) = 2n; a(A091855(n)) = 2n-1. - Philippe Deléham, Mar 26 2004
a(2^n) = (2^(n+1) + (-1)^n)/3. - Vladimir Shevelev, Apr 15 2010
If n = Sum_{i>=0} b_i*2^i is the binary expansion of n, then a(n) = 2n/3 + (1/3)Sum_{i>=0} b_i*(-1)^i. Thus a(n) = 2n/3 + O(log(n)). - Vladimir Shevelev, Apr 15 2010
Moreover, the equation a(3m)=2m has infinitely many solutions, e.g., a(3*2^k)=2*2^k; on the other hand, a((4^k-1)/3)=(2*(4^k-1))/9+k/3, i.e., limsup |a(n)-2n/3| = infinity. - Vladimir Shevelev, Feb 23 2011
a(n) = Sum_{k>=0} A030308(n,k)*A001045(k+1). - Philippe Deléham, Oct 19 2011
From Peter Bala, Feb 02 2013: (Start)
Product_{n >= 1} (1 + x^((2^n - (-1)^n)/3 )) = (1 + x)^2(1 + x^3)(1 + x^5)(1 + x^11)(1 + x^21)... = 1 + sum {n >= 1} x^a(n) = 1 + 2x + x^2 + x^3 + 2x^4 + 2x^5 + .... Hence this sequence lists the numbers representable as a sum of distinct Jacobsthal numbers A001045 = [1, 1', 3, 5, 11, 21, ...], where we distinguish between the two occurrences of 1 by writing them as 1 and 1'. For example, 9 occurs twice in the present sequence because 9 = 5 + 3 + 1 and 9 = 5 + 3 + 1'. Cf. A197911 and A080277. See also A120385.
(End)

Extensions

Extended with formula by Christian G. Bower, Sep 15 1999
Corrected and extended by Reinhard Zumkeller, Aug 16 2006
Extended with formula by Philippe Deléham, Oct 19 2011
Entry revised to give a simpler definition by N. J. A. Sloane, Jan 03 2014

A036577 Ternary Thue-Morse sequence: closed under a->abc, b->ac, c->b.

Original entry on oeis.org

2, 1, 0, 2, 0, 1, 2, 1, 0, 1, 2, 0, 2, 1, 0, 2, 0, 1, 2, 0, 2, 1, 0, 1, 2, 1, 0, 2, 0, 1, 2, 1, 0, 1, 2, 0, 2, 1, 0, 1, 2, 1, 0, 2, 0, 1, 2, 0, 2, 1, 0, 2, 0, 1, 2, 1, 0, 1, 2, 0, 2, 1, 0, 2, 0, 1, 2, 0, 2, 1, 0, 1, 2, 1, 0, 2, 0, 1, 2, 0, 2, 1, 0, 2, 0, 1, 2, 1, 0, 1, 2, 0, 2, 1, 0, 1, 2, 1, 0, 2, 0, 1, 2, 1, 0
Offset: 1

Views

Author

Keywords

Comments

Number of 1's between successive 0's in A010060.
The infinite sequence is abcacbabcbac... which is encoded with a->2, b->1, c->0 to produce this integer sequence.
From Jeffrey Shallit, Dec 07 2019: (Start)
This word is sometimes called 'vtm'; see, for example, see the Blanchet-Sadri et al. reference.
It is a squarefree word containing no instances of the factor 010 or 212 (or cbc, aba in the encoding).
Berstel proved many different definitions (e.g., Braunholtz, Istrail) of the word are equivalent. (End)

Examples

			2*x + x^2 + 2*x^4 + x^6 + 2*x^7 + x^8 + x^10 + 2*x^11 + 2*x^13 + x^14 + ...
		

References

  • M. Lothaire, Combinatorics on Words. Addison-Wesley, Reading, MA, 1983, p. 26.

Crossrefs

See A007413, A036580 for other versions.

Programs

  • Mathematica
    (* ThueMorse is built-in since version 10.2, for lower versions it needs to be defined manually *) ThueMorse[n_] := Mod[DigitCount[n, 2, 1], 2]; Table[1 + ThueMorse[n] - ThueMorse[n-1], {n, 1, 100}]  (* Vladimir Reshetnikov, May 17 2016 *)
    Nest[Flatten[# /. {2 -> {2, 1, 0}, 1 -> {2, 0}, 0 -> {1}}] &, {2, 1, 0}, 7] (* Robert G. Wilson v, Jul 30 2018 *)
    Differences[ThueMorse[Range[0, 100]]] + 1 (* Paolo Xausa, Jul 17 2025 *)
  • PARI
    {a(n) = if( n<1, 0, if( valuation( n, 2)%2, 1, 1 - (-1)^subst( Pol( binary(n)), x, 1)))} /* Michael Somos, Aug 03 2011 */
    
  • Python
    def A036577(n): return (n.bit_count()&1)+((n-1).bit_count()&1^1) # Chai Wah Wu, Mar 03 2023

Formula

a(n) = A036585(n) - 1 = A029883(n) + 1.
a(n) = 3 - A007413(n). a(A036554(n)) = 1; a(A091785(n)) = 0; a(A091855(n)) = 2. - Philippe Deléham, Mar 20 2004
a(4*n + 2) = 1. a(2*n + 1) = 2 * A010059(n). a(4*n + 3) = 2 * A010060(n). - Michael Somos, Aug 03 2011
a(n) = A010060(2*n - 1) + A010060(2*n) = A115384(2*n) - A115384(2*n - 2). - Zhuorui He, Jul 11 2025

A157971 Odious twin locations: first members of pairs of consecutive odious numbers.

Original entry on oeis.org

1, 7, 13, 21, 25, 31, 37, 41, 49, 55, 61, 69, 73, 81, 87, 93, 97, 103, 109, 117, 121, 127, 133, 137, 145, 151, 157, 161, 167, 173, 181, 185, 193, 199, 205, 213, 217, 223, 229, 233, 241, 247, 253, 261, 265, 273, 279, 285, 289, 295
Offset: 1

Views

Author

John W. Layman, Mar 10 2009

Keywords

Comments

An odious number (A000069) is a nonnegative integer with an odd number of ones in its binary expansion.
In the reference it is shown that these odious twins alternate with the evil twins (see A157970), which are pairs of consecutive evil numbers (A001969) having even numbers of ones in their binary expansions.

Examples

			The sequence of odious numbers (A000069) begins 1,2,4,7,8,11,13,14,16,19,21,..., so the first few odious twins are at 1,7,13, ... .
		

Crossrefs

Programs

  • Mathematica
    SequencePosition[Table[If[OddQ[DigitCount[n,2,1]],1,0],{n,300}],{1,1}][[All,1]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Oct 05 2016 *)
  • PARI
    lista(nn) = select(n->((hammingweight(n) % 2) && (hammingweight(n+1) % 2)), vector(nn, i, i)); \\ Michel Marcus, Jul 10 2014

Formula

a(n) = A248057(n) - 1. - Amiram Eldar, Jun 16 2025

Extensions

Comment corrected by Jeff Aronson. - N. J. A. Sloane, Oct 04 2020

A029883 First differences of Thue-Morse sequence A001285.

Original entry on oeis.org

1, 0, -1, 1, -1, 0, 1, 0, -1, 0, 1, -1, 1, 0, -1, 1, -1, 0, 1, -1, 1, 0, -1, 0, 1, 0, -1, 1, -1, 0, 1, 0, -1, 0, 1, -1, 1, 0, -1, 0, 1, 0, -1, 1, -1, 0, 1, -1, 1, 0, -1, 1, -1, 0, 1, 0, -1, 0, 1, -1, 1, 0, -1, 1, -1, 0, 1, -1, 1, 0, -1, 0, 1, 0, -1, 1, -1, 0, 1, -1, 1, 0, -1, 1, -1, 0, 1, 0, -1, 0, 1, -1, 1, 0, -1, 0, 1, 0, -1, 1, -1, 0, 1, 0, -1
Offset: 1

Views

Author

N. J. A. Sloane, Dec 11 1999

Keywords

Comments

Also first differences of {0,1} Thue-Morse sequence A010060.- N. J. A. Sloane, Jan 05 2021
Fixed point of the morphism a->abc, b->ac, c->b, with a = 1, b = 0, c = -1, starting with a(1) = 1. - Philippe Deléham
From Thomas Anton, Sep 22 2020: (Start)
This sequence, interpreted as an infinite word, is squarefree.
Let & represent concatenation. For a word w of integers, let -w be the same word with each symbol negated. Then, starting with the empty word, this sequence can be obtained by iteratively applying the transformation T(w) = w & 1 & -w & 0 & -w & -1 & w. (End)

Crossrefs

Apart from signs, same as A035263. Cf. A001285, A010060, A036554, A091785, A091855.
a(n+1) = A036577(n) - 1 = A036585(n) - 2.

Programs

  • Mathematica
    Nest[ Function[ l, {Flatten[(l /. {0 -> {1, -1}, 1 -> {1, 0, -1}, -1 -> {0}})]}], {1}, 7] (* Robert G. Wilson v, Feb 26 2005 *)
    ThueMorse /@ Range[0, 105] // Differences (* Jean-François Alcover, Oct 15 2019 *)
  • PARI
    a(n)=if(n<1||valuation(n,2)%2,0,-(-1)^subst(Pol(binary(n)),x,1)) /* Michael Somos, Jul 08 2004 */
    
  • PARI
    a(n)=hammingweight(n)%2-hammingweight(n-1)%2 \\ Charles R Greathouse IV, Mar 26 2013
    
  • Python
    def A029883(n): return (bin(n).count('1')&1)-(bin(n-1).count('1')&1) # Chai Wah Wu, Mar 03 2023

Formula

Recurrence: a(4*n) = a(n), a(4*n+1) = a(2*n+1), a(4*n+2) = 0, a(4*n+3) = -a(2*n+1), starting a(1) = 1.
a(n) = 2 - A007413(n). a(A036554(n)) = 0; a(A091785(n)) = -1; a(A091855(n)) = 1. - Philippe Deléham, Mar 20 2004
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w) = -v+w+u^2-v^2+2*w^2-2*u*w. - Michael Somos, Jul 08 2004

Extensions

Edited by Ralf Stephan, Dec 09 2004

A248057 Positions of 1,1 in the Thue-Morse sequence (A010060).

Original entry on oeis.org

2, 8, 14, 22, 26, 32, 38, 42, 50, 56, 62, 70, 74, 82, 88, 94, 98, 104, 110, 118, 122, 128, 134, 138, 146, 152, 158, 162, 168, 174, 182, 186, 194, 200, 206, 214, 218, 224, 230, 234, 242, 248, 254, 262, 266, 274, 280, 286, 290, 296, 302, 310, 314, 322, 328
Offset: 1

Views

Author

Clark Kimberling, Sep 30 2014

Keywords

Comments

Every positive integer lies in exactly one of these four sequences: A248056, A091855, A091855, A248057.

Examples

			Thue-Morse sequence:  0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,..., so that a(1) = 2 and a(2) = 8.
		

Crossrefs

Programs

  • Mathematica
    z = 400; u = Nest[Flatten[# /. {0 -> {0, 1}, 1 -> {1, 0}}] &, {0}, 9] (* A010060 *)
    v = Rest[u]
    t1 = Table[If[u[[n]] == 0 && v[[n]] == 0, 1, 0], {n, 1, z}];
    t2 = Table[If[u[[n]] == 0 && v[[n]] == 1, 1, 0], {n, 1, z}];
    t3 = Table[If[u[[n]] == 1 && v[[n]] == 0, 1, 0], {n, 1, z}];
    t4 = Table[If[u[[n]] == 1 && v[[n]] == 1, 1, 0], {n, 1, z}];
    Flatten[Position[t1, 1]]  (* A248056 *)
    Flatten[Position[t2, 1]]  (* A091855 *)
    Flatten[Position[t3, 1]]  (* A091785 *)
    Flatten[Position[t4, 1]]  (* A248057 *)
    SequencePosition[ThueMorse[Range[400]],{1,1}][[All,2]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, May 16 2017 *)
  • PARI
    t(n)=hammingweight(n)%2;
    for(n=1,500,if(t(n)==1&&t(n-1)==1,print1(n,", "))); \\ Joerg Arndt, Mar 12 2022

Formula

a(n) = 2*A091855(n) for n >= 1.
a(n) = A157971(n) + 1. - Amiram Eldar, Jun 16 2025

A248056 Positions of 0,0 in the Thue-Morse sequence (A010060).

Original entry on oeis.org

6, 10, 18, 24, 30, 34, 40, 46, 54, 58, 66, 72, 78, 86, 90, 96, 102, 106, 114, 120, 126, 130, 136, 142, 150, 154, 160, 166, 170, 178, 184, 190, 198, 202, 210, 216, 222, 226, 232, 238, 246, 250, 258, 264, 270, 278, 282, 288, 294, 298, 306, 312, 318, 326, 330
Offset: 1

Views

Author

Clark Kimberling, Sep 30 2014

Keywords

Comments

Every positive integer lies in exactly one of these four sequences: A248056, A091855, A091855, A248057.

Examples

			Thue-Morse sequence:  0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,..., so that a(1) = 6 and a(2) = 10.
		

Crossrefs

Programs

  • Mathematica
    z = 400; u = Nest[Flatten[# /. {0 -> {0, 1}, 1 -> {1, 0}}] &, {0}, 9] (* A010060 *)
    v = Rest[u]
    t1 = Table[If[u[[n]] == 0 && v[[n]] == 0, 1, 0], {n, 1, z}];
    t2 = Table[If[u[[n]] == 0 && v[[n]] == 1, 1, 0], {n, 1, z}];
    t3 = Table[If[u[[n]] == 1 && v[[n]] == 0, 1, 0], {n, 1, z}];
    t4 = Table[If[u[[n]] == 1 && v[[n]] == 1, 1, 0], {n, 1, z}];
    Flatten[Position[t1, 1]]  (* A248056 *)
    Flatten[Position[t2, 1]]  (* A091855 *)
    Flatten[Position[t3, 1]]  (* A091785 *)
    Flatten[Position[t4, 1]]  (* A248057 *)
    SequencePosition[ThueMorse[Range[400]],{0,0}][[All,2]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Mar 02 2020 *)

Formula

a(n) = 2*A091785(n) for n >= 1.
a(n) = A157970(n) + 1. - Amiram Eldar, Jun 16 2025
Showing 1-10 of 15 results. Next