OFFSET
1,2
COMMENTS
4*a(n) + 1 or (4*a(n) + 1)/3 is a prime number.
Compare with A380550, numbers not of the form i + 3*j + 4*i*j. See also A006093, numbers not of the form i + j + i*j and A005097, numbers not of the form i + j + 2*i*j. - Peter Bala, Jan 30 2025
LINKS
MAPLE
S := {}:
for n from 1 to 150 do
if isprime(4*n+1) then S := `union`(S, {n}) fi;
if type((4*n+1)*(1/3), integer) then if isprime((4*n+1)*(1/3)) then S := `union`(S, {n}) fi; fi;
end do:
S; # Peter Bala, Jan 30 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Davide Rotondo, Jan 27 2025
STATUS
approved