|
PROG
|
(PARI) limit=L; v=vector(limit); for(n=1, limit, forprime(p=2, n/2, for(k=1, logint(n\2, 2), for(j=k, logint(n, 2), forprime(q=p+1, n/2, for(t=1, logint(n\2, 3), for(u=t, logint(n, 3), if(n-p^k==p^j&&n-q^t==q^u, v[n]=1)))))))); for(n=1, limit, if(v[n], print1(n", ")))
|