OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
285 is in the sequence because 285 = 3*5*19 and (3+19) - 5 = 17 is prime, where 5 ==1 mod 4 and 3, 19 ==3 mod 4.
MAPLE
with(numtheory):for n from 2 to 1000 do:x:=factorset(n):n1:=nops(x):s1:=0:s3:=0:for m from 1 to n1 do: if irem(x[m], 4)=1 then s1:=s1+x[m]:else if irem(x[m], 4)=3 then s3:=s3+x[m]:else fi:fi:od:x:=abs(s1-s3):if s1>0 and s1>0 and s3>0 and type (x, prime)=true then printf(`%d, `, n):else fi:od:
MATHEMATICA
aQ[n_] := Module[{p = FactorInteger[n][[;; , 1]]}, (t1 = Total[Select[p, Mod[#, 4] == 1 &]]) > 0 && (t2 = Total[Select[p, Mod[#, 4] == 3 &]]) > 0 && PrimeQ@Abs[t1 - t2]]; Select[Range[630], aQ] (* Amiram Eldar, Sep 09 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Aug 28 2012
STATUS
approved