OFFSET
1,1
LINKS
Vincenzo Manto, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(n) isprime(n) and isprime(convert(NumberTheory:-PrimeFactors(n-1), `+`)) end proc:
select(filter, [seq(i, i=3..2000)]); # Robert Israel, May 29 2026
MATHEMATICA
Select[Prime[Range[300]], PrimeQ[Total[FactorInteger[# - 1][[All, 1]]]] &]
PROG
(Python)
from sympy import isprime, primerange, primefactors
def sequence(limit):
res = []
for p in primerange(2, limit):
if isprime(sum(primefactors(p - 1))):
res.append(p)
return res
print(sequence(2000))
CROSSREFS
KEYWORD
nonn
AUTHOR
Vincenzo Manto, May 28 2026
STATUS
approved
