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!)
A306601 Let b(1) = 3 and let b(n+1) be the least prime expressible as k*(b(n)-1)*b(n)-1; this sequence gives the values of k in order. 1

%I #92 Dec 10 2023 18:09:20

%S 1,1,2,4,8,16,5,360,142,104,34,1904,3127,253,1219,8755,16222,7672,

%T 22515

%N Let b(1) = 3 and let b(n+1) be the least prime expressible as k*(b(n)-1)*b(n)-1; this sequence gives the values of k in order.

%C The corresponding primes in order are 5, 19, 683, 1863223, P14, P29, P57, P117, P235, P472, P945, P1893, P3789, P7581, P15164, P30332, P60668, P121339, P242682.

%C After each iteration the number of decimal digits is roughly twice that of the previous iteration. These primes can generally be easily certified using the N+1 method since all the prime factors for N+1 are known.

%F Nested f(k) = k*(p-1)*p-1 for p=3. After each iteration the last obtained f(k) is substituted for p. The primes can be certified using OpenPFGW by adding each previous iteration to the helper file.

%e For p = 3, the smallest k for which f(k) = k*(p-1)*p-1 is prime is 1:

%e f(1) = k*(p-1)*p-1 = 1*(3-1)*3-1 = 5.

%e This sets p = 5 for the next iteration for which the smallest k for which f(k) is prime is 1:

%e f(1) = k*(p-1)*p-1 = 1*(5-1)*5-1 = 19.

%e This sets p = 19 for the next iteration for which the smallest k for which f(k) is prime is 2:

%e f(2) = k*(p-1)*p-1 = 2*(19-1)*19-1 = 683.

%e This sets p = 683 for the next iteration for which the smallest k for which f(k) is prime is 4:

%e f(4) = k*(p-1)*p-1 = 4*(683-1)*683-1 = 1863223.

%e This sets p = 1863223 for the next iteration for which the smallest k for which f(k) is prime is 8:

%e f(8) = k*(p-1)*p-1 = 8*(1863223-1)*1863223-1 = P14.

%o (PARI) p=3; k=1; while(1, runningP=k*(p-1)*p-1; if(ispseudoprime(runningP), print1(k,", "); k=1; p=runningP;, k=k+1))

%o The largest prime (P242682) can be generated by using the code:

%o (PARI) k=[1, 1, 2, 4, 8, 16, 5, 360, 142, 104, 34, 1904, 3127, 253, 1219, 8755, 16222, 7672, 22515]; p=3; for(i=1, #k, p=k[i]*(p-1)*p-1); print("\n", p, "\n")

%Y Cf. A000058.

%K nonn,hard,more

%O 1,3

%A _Rashid Naimi_, Apr 10 2019

%E Definition clarified by _Charlie Neder_, Jun 03 2019

%E a(17) from _Rashid Naimi_, Aug 23 2019

%E a(18) from _Rashid Naimi_, Oct 22 2019

%E a(19) from _Rashid Naimi_, Aug 01 2020

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 March 28 12:26 EDT 2024. Contains 371254 sequences. (Running on oeis4.)