cp's OEIS Frontend

A108278 Numbers k such that k^2-1 and k^2+1 are semiprimes.

%I A108278 #18 Sep 08 2022 08:45:19
%S A108278 12,30,42,60,102,108,198,312,462,522,600,810,828,1020,1050,1062,1278,
%T A108278 1452,1488,1872,1950,2028,2130,2142,2340,2790,2802,2970,3000,3120,
%U A108278 3252,3300,3330,3672,3930,4020,4092,4230,4548,4800,5280,5640,5652,5658,6198
%N A108278 Numbers k such that k^2-1 and k^2+1 are semiprimes.
%C A108278 Subsequence of A069062. - _Michel Marcus_, Jan 22 2016
%C A108278 Subsequence of A014574. - _Robert Israel_, Jan 24 2016
%H A108278 Robert Israel, <a href="/A108278/b108278.txt">Table of n, a(n) for n = 1..10000</a>
%e A108278 a(1)=12 because 12^2-1=143=11*13 and 12^2+1=145=5*29 are both semiprimes.
%p A108278 filter:= n -> isprime(n+1) and isprime(n-1) and numtheory:-bigomega(n^2+1)=2:
%p A108278 select(filter, [seq(i,i=2..1000, 2)]); # _Robert Israel_, Jan 24 2016
%t A108278 Select[Range[7000], PrimeOmega[#^2 - 1] == PrimeOmega[#^2 + 1]== 2 &] (* _Vincenzo Librandi_, Jan 22 2016 *)
%o A108278 (Magma) IsSemiprime:=func< n | &+[k[2]: k in Factorization(n)] eq 2 >; [ n: n in [4..7000] | IsSemiprime(n^2+1) and IsSemiprime(n^2-1) ]; // _Vincenzo Librandi_, Jan 22 2016
%o A108278 (PARI) isok(n) = (bigomega(n^2-1) == 2) && (bigomega(n^2+1) == 2); \\ _Michel Marcus_, Jan 22 2016
%Y A108278 Cf. A001358 (semiprimes), A069062 (k^2-1 and k^2+1 have the same number of divisors), A014574 (average of twin prime pairs).
%K A108278 nonn
%O A108278 1,1
%A A108278 _Hugo Pfoertner_, May 30 2005