%I #16 Dec 07 2021 10:59:43
%S 1,4,1,10,7,4,5,2,1,8,9,4,3,2,1,6,5,4,3,2,3,2,5,6,3,6,1,2,13,8,7,6,3,
%T 6,3,4,9,2,1,6,7,4,11,8,13,6,3,2,3,2,3,8,3,2,11,8,1,6,3,6,9,2,3,2,5,4,
%U 3,8,1,2,3,6,3,8,3,2,5,2,7,2,3,2,5,6,5,6
%N a(n) = n^2 - A099611(n).
%H Hugo Pfoertner, <a href="/A349809/b349809.txt">Table of n, a(n) for n = 4..10000</a>
%o (PARI) for(n=4,89,my(n2=n*n);forstep(k=n2,1,-1,if(k%2==1&&omega(k)==2&&bigomega(k)==2,print1(n2-k,", ");break))) \\ _Hugo Pfoertner_, Dec 06 2021
%o (Python)
%o from itertools import count
%o from sympy import factorint
%o def A349809(n):
%o for i in count(n**2-(n%2)-1,-2):
%o fs = factorint(i)
%o if len(fs) == 2 == sum(fs.values()):
%o return n**2-i # _Chai Wah Wu_, Dec 06 2021
%Y Cf. A099611.
%K nonn
%O 4,2
%A _N. J. A. Sloane_, Dec 06 2021