OFFSET
1,1
FORMULA
PROG
(PARI) a(n) = my(p = prime(n)+3, v = p/(2^valuation(p, 2))) ; if (v == 1, 1, factor(v)[1, 1]); \\ Michel Marcus, Aug 05 2021
(Python)
from sympy import factorint, prime
def A023578(n): return min((p for p in factorint(prime(n)+3) if p > 2), default=1) # Chai Wah Wu, Feb 03 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved