OFFSET
1,1
COMMENTS
Complement of A067377 in the primes. For the primes less than 10^6, these primes make up about 56%.
LINKS
Robert Israel, Table of n, a(n) for n = 1..5000
FORMULA
Prime(n) such that A307610(n) = 1. - Ray Chandler, Sep 21 2023
MAPLE
P:= [seq(ithprime(i), i=1..10^3)]:
S:= ListTools:-PartialSums([0, op(P)]):
sort(convert(convert(P, set) minus {seq(seq(S[i]-S[j], j=1..i-2), i=1..10^3+1)}, list)); # Robert Israel, May 09 2021
MATHEMATICA
lim = 1000; pFound = {}; ps = Prime[Range[PrimePi[lim]]]; sm = ps; i = 0; While[i++; j = 1; While[sm[[j]] = sm[[j]] + ps[[i + j]]; sm[[j]] <= lim, If[PrimeQ[sm[[j]]], AppendTo[pFound, sm[[j]]]]; j++]; j > 1]; Complement[ps, pFound]
CROSSREFS
KEYWORD
nonn
AUTHOR
T. D. Noe, Nov 03 2011
EXTENSIONS
Definition clarified by Jonathan Sondow, May 18 2013
STATUS
approved