OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
EXAMPLE
2951=5^5-5^3-2*5^2+1 is in the sequence because 2951=13*227 and 13, 227 are primes;
58159=9^5-9^3-2*9^2+1 is in the sequence because 58159=19*3061 and 19, 3061 are primes, etc.
MATHEMATICA
f[n_] := Last/@FactorInteger[n] == {1, 1}; Select[Array[#^5 - #^3 - 2*#^2 + 1 &, 100], f[#] &] (* Vincenzo Librandi, Sep 21 2012 *)
Select[Table[n^5 - n^3 - 2 n^2 + 1, {n, 60}], PrimeOmega[#] == 2 &] (* Bruno Berselli, Sep 21 2012 *)
PROG
(Magma) IsSemiprime:=func<n | &+[d[2]: d in Factorization(n)] eq 2>; [s: n in [2..60] | IsSemiprime(s) where s is n^5-n^3-2*n^2+1]; // Bruno Berselli, Sep 21 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Jun 05 2010
EXTENSIONS
Checked by Jud McCranie, Jun 16 2010
STATUS
approved