OFFSET
1,1
COMMENTS
From Jon E. Schoenfield, Jan 28 2024: (Start)
If it exists, a(13) > 10^2000.
Conjecture: a(12) = 1594323 is the final term of the sequence.
(End)
EXAMPLE
(prime = greatest odd divisor of a(n)-1; a(n); prime = greatest odd divisor of a(n)+1): (3; 4; 5), (5; 11; 3), (3; 13; 7), (11; 23; 3), (3; 25; 13), (13; 27; 7), (23; 47; 3), (5; 81; 41), (3; 193; 97), (191; 383; 3), (1093; 2187; 547), (797161; 1594323; 398581).
MATHEMATICA
q[n_] := PrimeQ[n/2^IntegerExponent[n, 2]]; Select[Range[2*10^6], PrimePowerQ[#] && And @@ q /@ {# - 1, # + 1} &] (* Amiram Eldar, Jan 28 2024 *)
PROG
(Magma) [k: k in [2..1600000] | #Divisors(2*k-2)-#Divisors(k-1) eq 2 and
#PrimeDivisors(k) eq 1 and #Divisors(2*k+2)-#Divisors(k+1) eq 2];
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Juri-Stepan Gerasimov, Jan 27 2024
STATUS
approved