OFFSET
1,2
MATHEMATICA
Select[Range[200], PrimeQ[Numerator[Abs[(1 - Sum[Prime[k] + 1, {k, 1, #}])/Product[Prime[k] + 1, {k, 1, #}] ]]]&]
PROG
(PARI) isok(n) = isprime(abs(numerator(matdet(matrix(n, n, i, j, if(i==j, prime(i)/(1+prime(i)), 1)))))); \\ Michel Marcus, May 10 2020
CROSSREFS
Cf. A120292 = absolute value of the numerator of the determinant of n X n matrix M with M[i, j] = prime[i]/(1 + prime[i]) if i = j, and 1 otherwise.
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Jul 04 2008
STATUS
approved