%I #25 Sep 08 2022 08:45:55
%S 44,46,80,88,102,104,108,226,234,238,246,272,290,308,310,328,334,358,
%T 370,426,456,480,514,526,530,586,588,614,720,766,790,842,846,848,872,
%U 880,884,896,898,900,934,940,974,980,1040,1076,1078,1088,1110,1160,1208
%N Numbers n such that n^2 + {1,3,7} are semiprimes.
%C This is to A182238 as A001358 semiprimes are to A000040 primes.
%H Alois P. Heinz, <a href="/A182261/b182261.txt">Table of n, a(n) for n = 1..1000</a>
%F { n : {n^2+1, n^2+3, n^2+7} in A001358 }.
%e 44 is in the sequence because (44^2) + 1 = 1937 = 13 * 149, (44^2) + 3 = 1939 = 7 * 277, and (442) + 7 = 1943 = 29 * 67.
%p a:= proc(n) option remember; local k;
%p for k from 1+a(n-1) while map(x-> not isprime(k^2+x) and
%p add(i[2], i=ifactors(k^2+x)[2])=2, [1, 3, 7])<>[true$3]
%p do od; k
%p end: a(0):=0:
%p seq(a(n), n=1..50); # _Alois P. Heinz_, Apr 22 2012
%t okQ[n_] := AllTrue[n^2 + {1, 3, 7}, PrimeOmega[#] == 2&];
%t Select[Range[2000], okQ] (* _Jean-François Alcover_, Jun 01 2022 *)
%o (Magma) IsSemiprime:=func<n | &+[m[2]: m in Factorization(n)] eq 2>; [n: n in [2..1225] | forall{n^2+i: i in [1,3,7] | IsSemiprime(n^2+i)}]; // _Bruno Berselli_, Apr 22 2012
%Y Cf. A001358, A182238.
%K nonn,easy
%O 1,1
%A _Jonathan Vos Post_, Apr 21 2012