The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A355176 a(n) is the smallest index k such that prime(n) divides both A090252(k) and A090252(2*k+1). 2
2, 3, 14, 32, 60, 96, 120, 128, 132, 244, 264, 388, 480, 484, 488, 2064, 1056, 571, 776, 960, 968, 976, 980, 2112, 2128, 1143, 1536, 1552, 1556, 1920, 3872, 1937, 3904, 3920, 1961, 4128, 4256, 3104, 6224, 3113, 3844, 3848, 7808, 7824, 7840, 8256, 8448, 8452 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
For n > 2, a(n) is not the smallest k such that prime(n) divides A090252(k), but it is the smallest k such that prime(n) divides both A090252(k) and A090252(2*k+1). If k_(0) = a(n) we may find either an infinite or finite range of indices where prime(n) divides A090252 using the recurrence k_(n) = 2*k_(n-1)+1, but there is a caveat: in very rare cases, some k values of this recurrence may be wrong by +-1, and the next iteration will then fit again. This uncertainty is caused by the fact that two terms of A090252 will be governed by the same floor(n/2) history. For yet unknown reasons, there may be an upper limit where such a recurrence may break.
This works because in A090252 the number of primes which do not divide the last floor(n/2) terms is growing faster than they are used up by this sequence. For each prime p then there exists an index k into A090252 where the supply of unused factors is so large that, when p becomes coprime to the last floor(n/2) terms, we can always immediately find a matching second prime to build a yet-unused semiprime or use p as a yet-unused power of itself.
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..476
FORMULA
A090252(a(n)) mod A000040(n) = 0 and a(n) is either even or A090252((a(n)-1)/2) mod A000040(n) > 0 is valid too.
A090252(2*a(n)+1) mod A000040(n) = 0.
A090252(f^m(a(n))) mod A000040(n) = 0, with f(x) = 2*x+1. The range of m is yet unknown.
EXAMPLE
prime(1) = 2 divides A090252(2) = 2, A090252(5) = 4, A090252(11) = 8, A090252(23) = 16, A090252(47) = 26, ... .
2*2+1 = 5; 2*5+1 = 11; 2*11+1 = 23; 2*23+1 = 47.
CROSSREFS
Sequence in context: A306267 A042551 A296149 * A346606 A059188 A080768
KEYWORD
nonn
AUTHOR
Thomas Scheuerle, Jun 22 2022
EXTENSIONS
a(41) and beyond (using Russ Cox's gzipped b-file at A090252) from Michael S. Branicky, Jun 23 2022
STATUS
approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified June 5 13:53 EDT 2024. Contains 373105 sequences. (Running on oeis4.)