|
|
A288212
|
|
Start with k=2*n, and until k+1 is prime, apply the map k -> k*(least prime factor of (k+1)); then a(n) = k+1, or 0 if k+1 never reaches a prime.
|
|
1
|
|
|
3, 5, 7, 1321, 11, 13, 43, 17, 19, 61, 23, 1321, 79, 29, 31, 97, 3571, 37, 571, 41, 43, 4621, 47, 337, 151, 53, 271, 21217561, 59, 61, 47059, 29761, 67, 1021, 71, 73, 223, 6917, 79, 241, 83, 421, 6221671, 89, 631, 277, 23971, 97, 1471, 101, 103, 313, 107, 109, 331
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
a(126) is unknown; it's known to either be 0 or have more than 152 digits.
a(126) has 277 digits. - J. Lowell, Jan 01 2023
|
|
LINKS
|
|
|
EXAMPLE
|
a(17) = 3571 because:
34 + 1 = 35 is composite with smallest prime factor 5, and 34*5 = 170;
170 + 1 = 171 is composite with smallest prime factor 3, and 170*3 = 510;
510 + 1 = 511 is composite with smallest prime factor 7, and 510*7 = 3570;
3570 + 1 = 3571 is prime.
If A359444 has a term k that is one less than a prime, then a(161) = k + 1; otherwise, a(161) = 0.
|
|
PROG
|
(PARI) a(n) = {my(x = 2*n); while (! isprime(x+1), x = x*vecmin(factor(x+1)[, 1]); ); x+1; } \\ Michel Marcus, Jun 07 2017
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|