|
COMMENTS
|
A122150[n] = Numerator[ Sum[ (-1)^(k+1) * 1/2^Prime[k], {k,1,n} ] ] begins {1,1,5,19,305,1219,19505,78019,1248305,79891519,319566077,20452228927,327235662833,...}. Indices of primes in A122150[n] are listed in A122151[n] = {3,4,13,60,66,75,175,...}.
|
|
MATHEMATICA
|
Do[f=Numerator[Sum[(-1)^(k+1)*1/2^Prime[k], {k, 1, n}]]; If[PrimeQ[f], Print[{n, f}]], {n, 1, 1000}]
|