OFFSET
1,1
MATHEMATICA
a=Table[If[PrimeQ[4*n+1]==True&&PrimeQ[4*n+3]==False, 1, 0], {n, 0, 200}]
b=Table[If[PrimeQ[4*n+1]==False&&PrimeQ[4*n+3]==True, 1, 0], {n, 0, 200}]
d=Table[If[PrimeQ[a[[n]]*(4*n+1)+b[[n]]*(4*n+3)]==True, a[[n]]*(4*n+1)+b[[n]]*(4*n+3), 0], {n, 1, 200}]
PROG
(PARI) av(terms) = my(v=[]); for(n=0, oo, if(ispseudoprime(4*n+1) && !ispseudoprime(4*n+3), v=concat(v, [1]), v=concat(v, [0])); if(#v==terms, break)); v
bv(terms) = my(v=[]); for(n=0, oo, if(!ispseudoprime(4*n+1) && ispseudoprime(4*n+3), v=concat(v, [1]), v=concat(v, [0])); if(#v==terms, break)); v
a(n) = my(avv=av(n), bvv=bv(n), x=avv[#avv], y=bvv[#bvv], val=x*(4*n+1)+y*(4*n+3)); if(ispseudoprime(val), return(val), return(0)) \\ Felix Fröhlich, Jul 03 2021; Program adapted from Mathematica code
CROSSREFS
KEYWORD
nonn,uned,obsc,less
AUTHOR
Roger L. Bagula, Jun 14 2004
STATUS
approved