(PARI) a=[]; for(n=1, 10, forprime(p=if(a, a[#a]+2, 3), , setsearch(a, p)&&next; for(i=1, 2^#a-1, isprime(normlp(vecextract(a, i), 1)+p)&&next(2)); a=concat(a, p); print1(p", "))) \\ Very simplistic, should at least avoid an odd number of terms in the partial sum of earlier terms. \\ M. F. Hasler, Jan 29 2015