OFFSET
1,1
COMMENTS
The prime factors of x^3 - x^2 - 2x + 1 come exclusively from this sequence. - Charles R Greathouse IV, Mar 18 2022
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
FORMULA
A045472 UNION {7}. - R. J. Mathar, Oct 18 2008
a(n) ~ 3n log n. - Charles R Greathouse IV, Mar 17 2022
MATHEMATICA
Select[Prime[Range[200]], MemberQ[{0, 1, 6}, Mod[#, 7]]&] (* Vincenzo Librandi, Aug 13 2012 *)
PROG
(Magma) [ p: p in PrimesUpTo(900) | p mod 7 in {0, 1, 6} ]; // Vincenzo Librandi, Aug 13 2012
(PARI) select(p->abs(centerlift(Mod(p, 7)))<2, primes(100)) \\ Charles R Greathouse IV, Mar 17 2022
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved