cp's OEIS Frontend

A241975 Numbers n such that n^4 - n^3 - n - 1 is a semiprime.

%I A241975 #23 Sep 08 2022 08:46:07
%S A241975 4,6,10,14,16,20,36,40,54,56,66,84,90,94,116,126,146,150,156,160,170,
%T A241975 204,210,260,264,306,340,350,386,396,406,420,464,474,496,570,634,674,
%U A241975 696,700,716,740,764,780,816,826,864,890,966,1054,1070,1094,1106,1144
%N A241975 Numbers n such that n^4 - n^3 - n - 1 is a semiprime.
%C A241975 Since n^4 - n^3 - n - 1 = (n^2 + 1)*(n^2 - n - 1), these are also numbers n such that n^2 + 1 and n^2 - n - 1 are both prime. Numbers in the intersection of A005574 and A002328. - _Derek Orr_, Aug 10 2014 [Sequence numbers corrected by _Jens Kruse Andersen_, Aug 11 2014]
%H A241975 Vincenzo Librandi, <a href="/A241975/b241975.txt">Table of n, a(n) for n = 1..400</a>
%t A241975 Select[Range[2000], PrimeOmega[#^4 - #^3 - # - 1]==2 &]
%o A241975 (Magma) IsSemiprime:=func<n | &+[ k[2]: k in Factorization(n) ] eq 2 >; [ n: n in [2..1500] | IsSemiprime(n^4 - n^3 - n - 1)];
%o A241975 (PARI) for(n=1,10^4,if(isprime(n^2+1)&&isprime(n^2-n-1),print1(n,", "))) \\ _Derek Orr_, Aug 10 2014
%Y A241975 Cf. A001358, A002327, A002496, A085722, A241670, A243436.
%K A241975 nonn,less
%O A241975 1,1
%A A241975 _Vincenzo Librandi_, Aug 10 2014