OFFSET
1,1
COMMENTS
PROG
(PARI) lista(nn) = my(c, f, m=1, q=1, r=0, t, u=List([]), w=vector(nn)); for(n=2, nn, t=0; forprime(p=2, max(m, sqrtint(n)), c=n-1-w[p]; if(c>1&&!w[c], listput(u, c))); listsort(u, 1); for(i=1, #u, c=u[i]; f=factor(c)[, 1]; t=1; for(j=1, #f, if(n-c<=w[f[j]], t=0; break)); if(t, u=u[i+1..#u]; if(r, print1(n-1, ", ")); r=1; w[c]=1; for(j=1, #f, w[f[j]]=n); m=max(m, f[#f]); break)); if(!t, r=0; u=List([]); if(nn>q=nextprime(q+1), w[q]=n))); \\ Jinyuan Wang, Feb 16 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Feb 15 2025
EXTENSIONS
a(7)-a(8) from Michael De Vlieger, Feb 15 2025
a(9)-a(34) from Jinyuan Wang, Feb 16 2025
STATUS
approved
