OFFSET
1,2
MATHEMATICA
f=1; Do[ p=Prime[n]; f=f - (-1)^(n+1)*1/p^2; g=Numerator[f] ; If[ PrimeQ[g], Print[ {n, g} ] ], {n, 1, 100} ]
PROG
(Python) # uses A136370gen() and imports from A136370
from sympy import isprime
def agen(): yield from (k for k, ak in enumerate(A136370gen(), 1) if isprime(ak))
print(list(islice(agen(), 5))) # Michael S. Branicky, Jun 26 2022
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Alexander Adamchuk, Dec 27 2007
EXTENSIONS
More terms added and edited by Alexander Adamchuk, Sep 15 2010
a(10) from Robert Price, Aug 29 2019
STATUS
approved