下半素数[n_,k_:1]:=块[{c=0,sgn=Sign[k]},sp=n+sgn;当[c<Abs[k]时,当[PrimeOmega[sp]!=2,如果[sgn<0,sp--,sp++]];如果[sgn<0,sp--,sp++];c++];sp+如果[sgn<0,1,-1]];f[n_]:=Times@@NestList[NextSemiPrime,2^2,n-1];k=1;lst={};当[k<3501时,如果[PrimeQ[f[k]-1],打印[k];附录[lst,k]];k++];第一次