OFFSET
1,1
COMMENTS
If p is a cuban prime (A002407) and p == 3 (mod 4) (A002145), then m = 3*p is a term. Indeed, there is k for which p = 1 + 3*k*(k + 1) and m = 3*p = 3 + 9*k*(k + 1) = (3*k + 2)^2 - (3*k + 2) + 1, so m is a term.
The sequence also includes terms that do not have this form: 133 = 12^2 - 12 + 1 = 7*19, 553 = 24^2 - 24 + 1 = 7*79, 1057 = 33^2 - 33 + 1 = 7*151, 1333 = 37^2 - 37 + 1= 31*43 and others.
EXAMPLE
MATHEMATICA
TR=40000; R1=Ceiling[(1+Sqrt[1-4(1-TR)])/2]; R2=TR/4; Intersection[Table[n^2-n+1, {n, 0, R1}], Select[4Range[5, R2]+1, PrimeNu[#]==2&&MoebiusMu[#]==1&&Mod[FactorInteger[#][[1, 1]], 4]!=1&]](* James C. McMahon, Feb 27 2024 *)
PROG
(Magma) pd:=PrimeDivisors; blum:=func<n|#Divisors(n) eq 4 and #pd(n) eq 2 and pd(n)[1] mod 4 eq 3 and pd(n)[2] mod 4 eq 3>; [n:n in [s^2-s+1:s in [2..2000]]|blum(n)];
CROSSREFS
KEYWORD
nonn
AUTHOR
Marius A. Burtea, Feb 27 2024
STATUS
approved