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 11-20 of 44 results. Next

A249681 A084937(3n+2).

Original entry on oeis.org

2, 4, 8, 6, 10, 16, 14, 22, 12, 18, 20, 26, 28, 32, 38, 34, 44, 40, 24, 30, 36, 46, 52, 56, 50, 58, 62, 64, 68, 74, 48, 42, 54, 60, 70, 76, 82, 86, 80, 88, 92, 94, 98, 100, 104, 66, 72, 78, 84, 96, 90, 106, 112, 110, 116, 118, 122, 124, 128, 134, 136, 142, 130
Offset: 0

Views

Author

N. J. A. Sloane, Nov 03 2014

Keywords

Crossrefs

Programs

A249575 Indices of prime powers in A084937.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 16, 17, 19, 21, 22, 24, 27, 28, 30, 31, 34, 37, 41, 43, 45, 51, 54, 55, 58, 60, 61, 67, 69, 73, 76, 79, 82, 83, 85, 91, 93, 94, 96, 97, 100, 103, 106, 115, 121, 124, 127, 130, 133, 136, 139, 142, 144, 145, 151, 160, 163
Offset: 1

Views

Author

Reinhard Zumkeller, Nov 04 2014

Keywords

Crossrefs

Programs

  • Haskell
    a249575 n = a249575_list !! (n-1)
    a249575_list = filter ((== 1) . a010055 . a084937) [1..]

A249602 Indices of primes in A084937.

Original entry on oeis.org

2, 3, 4, 6, 9, 10, 12, 13, 16, 19, 24, 27, 28, 30, 31, 37, 43, 45, 51, 54, 55, 58, 60, 61, 67, 73, 76, 79, 82, 85, 94, 96, 97, 100, 103, 106, 115, 121, 124, 130, 133, 136, 139, 142, 144, 145, 151, 160, 163, 166, 169, 175, 178, 187, 192, 198, 201, 202, 205
Offset: 1

Views

Author

N. J. A. Sloane, Nov 03 2014

Keywords

Crossrefs

Programs

  • Haskell
    a249602 n = a249602_list !! (n-1)
    a249602_list = filter ((== 1) . a010051' . a084937) [1..]
    -- Reinhard Zumkeller, Nov 04 2014

A249682 A084937(3n).

Original entry on oeis.org

3, 7, 11, 17, 21, 15, 25, 31, 37, 43, 33, 45, 39, 51, 61, 65, 67, 71, 77, 83, 85, 75, 81, 87, 93, 99, 105, 111, 115, 119, 125, 131, 133, 143, 123, 129, 141, 147, 153, 135, 159, 165, 171, 177, 175, 185, 187, 197, 205, 203, 217, 183, 201, 189, 195, 207, 213, 219
Offset: 1

Views

Author

N. J. A. Sloane, Nov 03 2014

Keywords

Crossrefs

Programs

A249603 A084937 mod 3.

Original entry on oeis.org

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

Views

Author

N. J. A. Sloane, Nov 03 2014

Keywords

Comments

By definition of A084937 this cannot contain two consecutive 0's. So the 0's divide this sequence into runs of 1's and 2's. In the first 100000 terms the following runs of 1's and 2's appear: 11, 12, 21, 22, 111, 112, 121, 122, 211, 212, 221, 222, 1112, 1121, 1122, 1212, 1221, 2111, 2112, 2121, 2211, 2212, 2221. There is no obvious structure. - N. J. A. Sloane, Nov 07 2014

Crossrefs

Cf. A084937.

Programs

A249689 Numbers k such that A084937(3k) > A084937(3k+1).

Original entry on oeis.org

16, 65, 91, 1514, 1993, 1994, 2452, 2722, 3047, 3214, 3931, 3957, 4356, 4366, 5191, 5581, 5805, 5806, 7519, 8871, 9228, 9752, 10036, 10037, 10039, 10040, 10963, 10964, 11278, 11279, 12015, 12281, 12595, 12665, 13262, 13618, 13648, 15102, 15103, 18529, 18991
Offset: 1

Views

Author

N. J. A. Sloane, Nov 12 2014

Keywords

Crossrefs

Programs

  • AWK
    # Using a-file for A084937.
    awk ' BEGIN {s = 0}
          NR%3 == 0 { s = $2 }
          NR%3 == 1 { t = $2
              if (s > t)
                 print( (NR-1)/3 )
            } ' a084937.txt | awk '{ print NR, $1}' >b249689.txt
    
  • Python
    from math import gcd
    A249689_list, l1, l2, s, b = [], 2, 1, 3, set()
    for n in range(3,10**4):
        i = s
        while True:
            if not i in b and gcd(i,l1) == 1 and gcd(i,l2) == 1:
                l2, l1 = l1, i
                b.add(i)
                if l2 > l1 and n % 3 == 1:
                    A249689_list.append((n-1)//3)
                while s in b:
                    b.remove(s)
                    s += 1
                break
            i += 1 # Chai Wah Wu, Dec 12 2014

A249694 a(n) = gcd(A084937(n), A084937(n+3)).

Original entry on oeis.org

1, 2, 1, 1, 4, 1, 1, 2, 1, 1, 2, 1, 1, 2, 3, 1, 2, 5, 1, 2, 1, 1, 2, 1, 1, 6, 1, 1, 2, 1, 1, 2, 3, 1, 2, 3, 1, 4, 3, 1, 2, 1, 1, 2, 1, 3, 2, 1, 3, 4, 1, 1, 8, 1, 1, 6, 1, 1, 6, 1, 1, 2, 5, 1, 2, 3, 1, 4, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 1, 4, 1, 1, 2
Offset: 1

Views

Author

Reinhard Zumkeller, Nov 04 2014

Keywords

Comments

gcd(A084937(n), A084937(n+1)) = gcd(A084937(n), A084937(n+2)) = 1 by definition of A084937.

Crossrefs

Cf. A084937.

Programs

  • Haskell
    a249694 n = a249694_list !! (n-1)
    a249694_list = zipWith gcd (drop 3 a084937_list) a084937_list

A011655 Period 3: repeat [0, 1, 1].

Original entry on oeis.org

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

Views

Author

Keywords

Comments

A binary m-sequence: expansion of reciprocal of x^2+x+1 (mod 2).
A Chebyshev transform of the Jacobsthal numbers A001045: if A(x) is the g.f. of a sequence, map it to ((1-x^2)/(1+x^2))*A(x/(1+x^2)). - Paul Barry, Feb 16 2004
This is the r = 1 member of the r-family of sequences S_r(n) defined in A092184 where more information can be found.
This is the Fibonacci sequence (A000045) modulo 2. - Stephen Jordan (sjordan(AT)mit.edu), Sep 10 2007
For n > 0: a(n) = A084937(n-1) mod 2. - Reinhard Zumkeller, Dec 16 2007
This is also the Lucas numbers (A000032) mod 2. In general, this is the parity of any Lucas sequence associated with any pair (P,Q) when P and Q are odd; i.e., a(n) = U_n(P,Q) mod 2 = V_n(P,Q) mod 2. See Ribenboim. - Rick L. Shepherd, Feb 07 2009
Starting with offset 1: (1, 1, 0, 1, 1, 0, ...) = INVERTi transform of the tribonacci sequence A001590 starting (1, 2, 3, 6, 11, 20, 37, ...). - Gary W. Adamson, May 04 2009
From Reinhard Zumkeller, Nov 30 2009: (Start)
Characteristic function of numbers coprime to 3.
a(n) = 1 - A079978(n); a(A001651(n)) = 1; a(A008585(n)) = 0;
A000212(n) = Sum_{k=0..n} a(k)*(n-k). (End)
Sum_{k>0} a(k)/k/2^k = log(7)/3. - Jaume Oliver Lafont, Jun 01 2010
The sequence is the principal Dirichlet character of the reduced residue system mod 3 (the other is A102283). Associated Dirichlet L-functions are L(2,chi) = Sum_{n>=1} a(n)/n^2 = 4*Pi^2/27 = A214549, and L(3,chi) = Sum_{n>=1} a(n)/n^3 = 1.157536... = -(psi''(1/3) + psi''(2/3))/54 where psi'' is the tetragamma function. [Jolley eq 309 and arXiv:1008.2547, L(m = 3, r = 1, s)]. - R. J. Mathar, Jul 15 2010
a(n+1), n >= 0, is the sequence of the row sums of the Riordan triangle A158454. - Wolfdieter Lang, Dec 18 2010
Removing the first two elements and keeping the offset at 0, this is a periodic sequence (1, 0, 1, 1, 0, 1, ...). Its INVERTi transform is (1, -1, 2, -2, 2, -2, ...) with period (2,-2). - Gary W. Adamson, Jan 21 2011
Column k = 1 of triangle in A198295. - Philippe Deléham, Jan 31 2012
The set of natural numbers, A000027: (1, 2, 3, ...); is the INVERT transform of the signed periodic sequence (1, 1, 0, -1, -1, 0, 1, 1, 0, ...). - Gary W. Adamson, Apr 28 2013
Any integer sequence s(n) = |s(n-1) - s(n-2)| (equivalently, max(s(n-1), s(n-2)) - min(s(n-1), s(n-2))) for n > i + 1 with s(i) = j and s(i+1) = k, where j and k are not both 0, is or eventually becomes a multiple of this sequence, namely, the sequence repeat gcd(j, k), gcd(j, k), 0 (at some offset). In particular, if j and k are coprime, then s(n) is or eventually becomes this sequence (see, e.g., A110044). - Rick L. Shepherd, Jan 21 2014
For n >= 1, a(n) is also the characteristic function for rational g-adic integers (+n/3)A001651).%20See%20the%20definition%20in%20the%20Mahler%20reference,%20p.%207%20and%20also%20p.%2010.%20-%20_Wolfdieter%20Lang">g and also (-n/3)_g for all integers g >= 2 without a factor 3 (A001651). See the definition in the Mahler reference, p. 7 and also p. 10. - _Wolfdieter Lang, Jul 11 2014
Characteristic function for A007908(n+1) being divisible by 3. a(n) = bit flipped A007908(n+1) (mod 3) = bit flipped A079978(n). - Wolfdieter Lang, Jun 12 2017
Also Jacobi or Kronecker symbol (n/9) (or (n/9^e) for all e >= 1). - Jianing Song, Jul 09 2018
The binomial trans. is 0, 1, 3, 6, 11, 21, 42, 85, 171, 342,.. (see A024495). - R. J. Mathar, Feb 25 2023

Examples

			G.f. = x + x^2 + x^4 + x^5 + x^7 + x^8 + x^10 + x^11 + x^13 + x^14 + x^16 + x^17 + ...
		

References

  • S. W. Golomb, Shift-Register Sequences, Holden-Day, San Francisco, 1967.
  • H. D. Lueke, Korrelationssignale, Springer 1992, pp. 43-48.
  • F. J. MacWilliams and N. J. A. Sloane, The Theory of Error-Correcting Codes, Elsevier/North Holland, 1978, p. 408.
  • K. Mahler, p-adic numbers and their functions, 2nd ed., Cambridge University press, 1981.
  • Paulo Ribenboim, The Little Book of Big Primes. Springer-Verlag, NY, 1991, p. 46. [Rick L. Shepherd, Feb 07 2009]

Crossrefs

Partial sums of A057078 give A011655(n+1).
Cf. A035191 (Mobius transform), A001590, A002487, A049347.
Cf. A000027, A000045, A004523 (partial sums), A057078 (first differences).
Cf. A007908, A079978 (bit flipped).
Cf. A011656 - A011751 for other binary m-sequences.
Cf. A002264.

Programs

Formula

G.f.: (x + x^2) / (1 - x^3) = Sum_{k>0} (x^k - x^(3*k)).
G.f.: x / (1 - x / (1 + x / (1 + x / (1 - 2*x / (1 + x))))). - Michael Somos, Apr 02 2012
a(n) = a(n+3) = a(-n), a(3*n) = 0, a(3*n + 1) = a(3*n + 2) = 1 for all n in Z.
a(n) = (1/2)*( (-1)^(floor((2n + 4)/3)) + 1 ). - Mario Catalani (mario.catalani(AT)unito.it), Oct 22 2003
a(n) = Fibonacci(n) mod 2. - Paul Barry, Nov 12 2003
a(n) = (2/3)*(1 - cos(2*Pi*n/3)). - Ralf Stephan, Jan 06 2004
a(n) = 1 - a(n-1)*a(n-2), a(n) = n for n < 2. - Reinhard Zumkeller, Feb 28 2004
a(n) = 2*(1 - T(n, -1/2))/3 with Chebyshev's polynomials T(n, x) of the first kind; see A053120. - Wolfdieter Lang, Oct 18 2004
a(n) = n*Sum_{k=0..floor(n/2)} (-1)^k*binomial(n-k, k)*A001045(n-2k)/(n-k). - Paul Barry, Oct 31 2004
a(n) = A002487(n) mod 2. - Paul Barry, Jan 14 2005
From Bruce Corrigan (scentman(AT)myfamily.com), Aug 08 2005: (Start)
a(n) = n^2 mod 3.
a(n) = (1/3)*(2 - (r^n + r^(2*n))) where r = (-1 + sqrt(-3))/2. (End)
From Michael Somos, Sep 23 2005: (Start)
Euler transform of length 3 sequence [ 1, -1, 1].
Moebius transform is length 3 sequence [ 1, 0, -1].
Multiplicative with a(3^e) = 0^e, a(p^e) = 1 otherwise. (End)
From Hieronymus Fischer, Jun 27 2007: (Start)
a(n) = (4/3)*(|sin(Pi*(n-2)/3)| + |sin(Pi*(n-1)/3)|)*|sin(Pi*n/3)|.
a(n) = ((n+1) mod 3 + 1) mod 2 = (1 - (-1)^(n - 3*floor((n+1)/3)))/2. (End)
a(n) = 2 - a(n-1) - a(n-2) for n > 1. - Reinhard Zumkeller, Apr 13 2008
a(2*n+1) = a(n+1) XOR a(n), a(2*n) = a(n), a(1) = 1, a(0) = 0. - Reinhard Zumkeller, Dec 27 2008
Sum_{n>=1} a(n)/n^s = (1-1/3^s)*Riemann_zeta(s), s > 1. - R. J. Mathar, Jul 31 2010
a(n) = floor((4*n-5)/3) mod 2. - Gary Detlefs, May 15 2011
a(n) = (a(n-1) - a(n-2))^2 with a(0) = 0, a(1) = 1. - Francesco Daddi, Aug 02 2011
Convolution of A040000 with A049347. - R. J. Mathar, Jul 21 2012
G.f.: Sum_{k>0} x^A001651(k). - L. Edson Jeffery, Dec 05 2012
G.f.: x/(G(0) - x^2) where G(k) = 1 - x/(x + 1/(1 - x/G(k+1))); (recursively defined continued fraction). - Sergei N. Gladkovskii, Feb 15 2013
For the general case: The characteristic function of numbers that are not multiples of m is a(n) = floor((n-1)/m) - floor(n/m) + 1, with m,n > 0. - Boris Putievskiy, May 08 2013
a(n) = sign(n mod 3). - Wesley Ivan Hurt, Jun 22 2013
a(n) = A000035(A000032(n)) = A000035(A000045(n)). - Omar E. Pol, Oct 28 2013
a(n) = (-n mod 3)^((n-1) mod 3). - Wesley Ivan Hurt, Apr 16 2015
a(n) = (2/3) * (1 - sin((Pi/6) * (4*n + 3))) for n >= 0. - Werner Schulte, Jul 20 2017
a(n) = a(n-1) XOR a(n-2) with a(0) = 0, a(1) = 1. - Chunqing Liu, Dec 18 2022
a(n) = floor((n+2)/3) - floor(n/3) = A002264(n+2) - A002264(n). - Aaron J Grech, Jul 30 2024
E.g.f.: 2*(exp(x) - exp(-x/2)*cos(sqrt(3)*x/2))/3. - Stefano Spezia, Mar 30 2025
Dirichlet g.f.: zeta(s)*(1-1/3^s). - R. J. Mathar, Aug 10 2025

Extensions

Better name from Omar E. Pol, Oct 28 2013

A329333 There is exactly one odd prime among the pairwise sums of any three consecutive terms: Lexicographically earliest sequence of distinct nonnegative integers with this property.

Original entry on oeis.org

0, 1, 2, 7, 3, 6, 4, 5, 8, 10, 11, 9, 12, 14, 15, 13, 18, 17, 19, 20, 21, 24, 16, 23, 25, 22, 26, 27, 28, 31, 29, 32, 33, 34, 30, 39, 37, 36, 38, 41, 40, 42, 43, 46, 35, 44, 47, 45, 50, 51, 48, 49, 56, 52, 53, 54, 57, 55, 58, 59, 68, 60, 63, 64, 61, 66, 62, 69, 67, 72, 71, 65, 74, 70, 75, 76, 77
Offset: 0

Views

Author

Keywords

Comments

This is conjectured and designed to be a permutation of the nonnegative integers, therefore the offset is taken to be zero.
Restricted to positive indices, this is a sequence of positive integers having the same property, then conjectured to be a permutation of the positive integers. (The word "odd" can be omitted in this case.)
If the word "odd" is dropped from the original definition, the sequence starts (0, 1, 3, 6, 2, 7), and then continues from a(6) = 4 onward as the present sequence. This is again conjectured to be a permutation of the nonnegative integers, and a permutation of the positive integers when restricted to the domain [1..oo). The latter however no longer has the property of lexicographic minimality.
See the OEIS wiki page for further considerations about existence, surjectivity and variants. - M. F. Hasler, Nov 24 2019

Examples

			For the first two terms there is no restriction regarding primality, so a(0) = 0, a(1) = 1. (If only positive values and indices are considered, then a(1) = 1 and a(2) = 2.)
Then a(2) must be such that among { 0+1, 0+a(2), 1+a(2) } there is exactly one odd prime, and 2 works.
Then a(3) must be such that among { 1+2, 1+a(3), 2+a(3) } there is only one (odd) prime. Since 1+2 = 3, the other two sums must both yield a composite. This excludes 3, 4, 5 and 6 and the smallest possibility is a(3) = 7.
And so on.
		

Crossrefs

For the primes that arise, or are missing, see A328997, A328998.
See A329450 for the variant having 0 primes among a(n+i) + a(n+j), 0 <= i < j < 3.
See A329452 for the variant having 2 primes among a(n+i) + a(n+j), 0 <= i < j < 4.
A084937, A305369 have comparable conditions on three consecutive terms.

Programs

  • Mathematica
    a[0]=0;a[1]=1;a[2]=2;a[n_]:=a[n]=(k=1;While[Length@Select[Plus@@@Subsets[{a[n-1],a[n-2],++k},{2}],PrimeQ]!=1||MemberQ[Array[a,n-1,0],k]];k);Array[a,100,0] (* Giorgos Kalogeropoulos, May 07 2021 *)
  • PARI
    A329333(n,show=0,o=0,p=0,U=[])={for(n=o,n-1, show&&print1(o","); U=setunion(U,[o]); while(#U>1&&U[1]==U[2]-1,U=U[^1]); for(k=U[1]+1,oo, setsearch(U,k)|| if(isprime(o+p), isprime(o+k)|| isprime(p+k), isprime(o+k)==isprime(p+k)&&p)||[o&&p=o, o=k, break]));o} \\ Optional args: show = 1: print all values up to a(n); o = 1: start with a(1) = 1; p = 1: compute the variant with a(2) = 3. See the wiki page for more general code which returns the whole vector: Use S(n_max,1,3,1) or S(n_max,1,3,2,[0,1]); S(n_max,1,3,0) gives the variant (0, 1, 3, ...)

Extensions

Entry revised by N. J. A. Sloane, Nov 14 2019 and M. F. Hasler, Nov 15 2019

A090252 The Two-Up sequence: a(n) is the least positive number not already used that is coprime to the previous floor(n/2) terms.

Original entry on oeis.org

1, 2, 3, 5, 4, 7, 9, 11, 13, 17, 8, 19, 23, 25, 21, 29, 31, 37, 41, 43, 47, 53, 16, 59, 61, 67, 71, 73, 55, 79, 27, 49, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 26, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 85, 121, 223, 227, 57, 229
Offset: 1

Views

Author

Amarnath Murthy, Nov 27 2003

Keywords

Comments

a(n) is coprime to the next n terms. - David Wasserman, Oct 24 2005
All values up to a(1000000) are either prime powers or semiprimes; this suggests the sequence is unlikely to be a permutation of the integers.
It appears that a(n) is even iff n = 3*2^k-1 for some k (A083356). - N. J. A. Sloane, Nov 01 2014
The even terms in the present sequence are listed in A354255.
We have a(1) = 1 and a(2) = 2. At step k >= 2, the sequence is extended by adding two terms: a(2*k-1) = smallest unused number which is relatively prime to a(k), a(k+1), ..., a(2*k-2), and a(2*k) = smallest unused number which is relatively prime to a(k), a(k+1), ..., a(2*k-1). So at step k=2 we add a(3)=3, a(4)=5; at step k=3 we add a(5)=4, a(6)=7; and so on. - N. J. A. Sloane, May 21 2022
Comments from N. J. A. Sloane, May 23 2022: (Start)
Conjecture 1. A090252 is a subsequence of A354144 (prime powers and semiprimes).
Conjecture 2. The terms of A354144 that are missing from A090252 are 6, 10, 14, 15, 22, 33, 34, 35, 38, 39, 46, 51, 58, 62, 65, 69, 74, 77, 82, 86, 87, 91, 93, 94, 95, 106, 111, 115, 118, 119, 122, 123, 129, 133, 134, 141, 142, 143, 145, 146, 155, 158, 166, 177, 178, 183, 185, 187, 194, 201, 202, 203, 209, 213, 214, 215, 218, 219, 221, ...
But since there is no proof that any one of these numbers is really missing, this list cannot yet have an entry in the OEIS.
Let S_p = list of indices of terms in A090252 that are divisible by the prime p.
Conjecture 3. For a prime p, there are constants v_1, v_2, ..., v_K and c such that
S_p = { v_1, v_2, ..., v_k, lambda*2^i - 1, i >= c}.
For example, from Michael S. Branicky's 10000-term b-file, it appears that:
S_2 = { 3*2^k-1, k >= 0 } cf. A083329
S_3 = { 2^k-1, k >= 2 } cf. A000225
S_5 = { 4 then 15*2^k-1 k >= 0 } cf. A196305
S_7 = { 6, 15, then 33*2^k-1, k >= 0 }
S_11 = { 8, 29, then 61*2^k-1, k >= 0 }
S_13 = { 9, 47, 97*2^n-1, n >= 0 }
S_17 = { 10, 59, 121*2^n-1, n >= 0 }
S_19 = { 12, 63, 129*2^n-1, n >= 0 }
S_23 = { 13, 65, 133*2^n-1, n >= 0 }
S_29 = { 16, 121, 245*2^n-1, n >= 0 }
S_31 = { 17, 131, 265*2^n-1, n >= 0 }
The initial primes p and the corresponding values of lambda are:
p: 2 3 5 7 11 13 17 19 23 29 31
lambda:..3...1..15..33...61...97..121..129..133..245..265
(This sequence of lambdas does not seem to have any simpler explanation, is not in the OEIS, and cannot be since the terms shown are all conjectural.)
Conjecture 2 is a consequence of Conjecture 3. For example, 6 does not appear in A090252, since the sets S_2 and S_3 are disjoint.
Also 10 does not appear, since S_2 and S_5 are disjoint.
In fact 2*p for 3 <= p <= 11 does not appear, but 26 = 2*13 does appear since S_2 and S_13 have 47 in common.
Assuming the numbers that appear to be missing (see Conjecture 2) really are missing, the numbers that take a record number of steps to appear are 1, 2, 3, 4, 7, 8, 16, 26, 32, 64, 128, 206, 256, 478, 512, 933, ..., and the indices where they appear are 1, 2, 3, 5, 6, 11, 23, 47, 95, 191, 383, 767, 1535, 3071, 6143, 8191, .... These two sequences are not yet in the OEIS, and cannot be added since the terms are all conjectural.
(End)
From N. J. A. Sloane, Jun 06 2022 (Start)
Theorem: (a) a(n) <= prime(n-1) for all n >= 2 (cf. A354154).
(b) A stronger upper bound is the following. Let c(n) = A354166(n) denote the number of nonprime terms among a(1) .. a(n). Note c(1)=1. Then a(n) <= prime(n-c(n)) for n <> 7 and 14.
It appears that a(n) = prime(n-c(n)) for almost all n. That is, this is the equation to the line in the graph that contains most of the terms.
For example, a(34886) = 408710 (see the b-file) = prime(34886 - A354166(34886)) = prime(34886 - 374) = prime(34512) = 408710.
Another example: Consider Russ Cox's table of the first N = 5764982 terms. We see that a(5764982) = 99999989 = prime(5761455) = prime(N - 3527) which agrees with c(N) = 3527 (from the first Russ Cox link).
(End)
If we consider the May 23 2022 comment, note the conjectured indices show near complete overlap with terms of A081026: 1, 2, 3, 5, 6, 11, 23, 47, 95, 191, 383, 767, 1535, 3071, 6143, 8191. - Bill McEachen, Aug 09 2024

Crossrefs

See A247665 for the case when the numbers are required to be at least 2. A353730 is another version.
For a squarefree analog, see A354790, A354791, A354792.

Programs

Extensions

More terms from David Wasserman, Oct 24 2005
Previous Showing 11-20 of 44 results. Next