close
login
A387122
Greatest prime divisor of 2*prime(n) + 3.
0
7, 3, 13, 17, 5, 29, 37, 41, 7, 61, 13, 11, 17, 89, 97, 109, 11, 5, 137, 29, 149, 23, 13, 181, 197, 41, 19, 31, 17, 229, 257, 53, 277, 281, 43, 61, 317, 47, 337, 349, 19, 73, 11, 389, 397, 401, 17, 449, 457, 461, 67, 37, 97, 101, 47, 23, 541, 109, 557, 113, 569, 31
OFFSET
1,1
FORMULA
a(n) = A006530(A129069(n)).
EXAMPLE
a(1) = Gpf(2*2+3) = Gpf(7) = 7.
a(2) = Gpf(2*3+3) = Gpf(9) = 3.
MAPLE
f:= n -> max(numtheory:-factorset(2*ithprime(n)+3)):
map(f, [$1..100]); # Robert Israel, Nov 12 2025
MATHEMATICA
a[n_] := FactorInteger[2*Prime[n] + 3][[-1, 1]]; Array[a, 65] (* Amiram Eldar, Nov 12 2025 *)
PROG
(PARI) a(n) = vecmax(factor(2*prime(n) + 3)[, 1]);
CROSSREFS
Sequence in context: A361416 A076569 A061194 * A248280 A255772 A195159
KEYWORD
nonn
AUTHOR
Alain Rocchelli, Nov 12 2025
STATUS
approved